data(Snacks)
Snacks2 <- Snacks[order(Snacks$snack,Snacks$week,Snacks$texture),]
attach(Snacks2)
snack <- factor(snack)
week2 <- week^2
week3 <- week^3
fit <- ssym.l(log(texture),formula.mu=~snack + week2 + week3,formula.phi=~snack,
ncs=week,family='Normal',local.influence=TRUE)
summary(fit)
Run the code above in your browser using DataLab