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