x <- c(1:30,1:30)
y <- c(1:30+rnorm(30,0,3),seq(10,22,12/29)+rnorm(30,0,3))
fact <- gl(2,30,labels=LETTERS[1:2])
regression <- least.rect(y~x|fact)
predict(regression,newdata=list(x=c(-2,18),fact=c("B","A")))
tab <- data.frame(x=c(-2,18),fact=c("B","A"))
predict(regression,newdata=tab)
Run the code above in your browser using DataLab