# NOT RUN {
xy <- MSDecision
xy$sensoryOnset1 <- factor(xy$OnsetSymptom1=="SENSORY")
xy$brainstemOnset2 <- factor(xy$OnsetSymptom2=="BRAINSTEM")
xy$sensoryTrigger1 <- factor(xy$TriggerSymptom1=="SENSORY")
xy$brainstemTrigger2 <- factor(xy$TriggerSymptom2=="BRAINSTEM")
xy <- xy[, -c(5, 6, 8, 9, 15)]
xy[,1]<-log(xy[,1])
names(xy)[1] <- "y"
out <- ices(y ~ ., data = xy)
summary(out)
plot(out)
plot(out, normqq=TRUE)
plot(out, scaleloc=TRUE)
# }
Run the code above in your browser using DataLab