# NOT RUN {
data(luzdat)
y <- luzdat$y
x1 <- luzdat$x1 ; x1 <- factor(x1) ; x1 <- C(x1, treatment)
x2 <- luzdat$x2
x3 <- (luzdat$x2)^2
luz <- data.frame(y, x1, x2, x3)
elliptical.fitt <- elliptical(y ~ x1+x2+x3, family = Student(df=5),
data = luz)
elliptical.diag(elliptical.fitt)
# }
Run the code above in your browser using DataLab