Try: x = matrix(runif(30), 10, 3) newx = matrix(runif(15), 5, 3) y = rnorm(10)
fit = lm(y ~ x) newyhat = predict(fit, newdata=data.frame(x=I(newx))) HTH, Andy > From: Claus Gwiggner > > Hello, > > can someone give me an example of how to fit a linear model > where the > response is a matrix, please? (other than a loop over lm calls) > The ?lm points me to model.matrix but it is not clear to me > how to use a > model matrix in an lm call. > > fit<-lm(headmatrix ~ tailmatrix) causes problems as soon as I > use predict(fit,newframe). The variables in tailmatrix are > not found in > newframe. > > Thanks, > > Claus ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html