# 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.fitC <- elliptical(y ~ x1+x2+x3, family = Cauchy()
,data=luz)
family(elliptical.fitC)
# }
Run the code above in your browser using DataLab