data("Snacks", package="ssym")
fit <- ssym.l(log(texture) ~ type + ncs(week) | type, data=Snacks,
family='Student', xi=15)
summary(fit)
Run the code above in your browser using DataLab