Yes, this is it (as would say michael)! Thank you guys!
Last question about another function on this list: imagine this list is my
data after your function for the regression model:
mydf <- data.frame(x=c(1:5), y=c(21:25),z=rnorm(1:5))
mylist <- rep(list(mydf),5)
Don't care about this fake data
Hello,
Your example doesn't run, lmList needs a 'data' argument.
set.seed(8109)
d <- rep(1:10, each=10)
x <- rnorm(100)
e <- rnorm(100)
y <- 2*x + e
dat <- data.frame(x=x, y=y, d=d, e=e)
model <- lmList(y ~ x | d, data = dat)
predict(model)
To the op: if you already have the larger data.farme
The easiest way may be to use lmList in the nlme library:
#simulate data
d<-rep(1:10,each=10)
x<-rnorm(100)
e<-rnorm(100)
y<-2*x+e
require(nlme) #or install and load package
lmList(y~x|d)
#predicted values are obtained with:
predict(lmList(y~x|d)
HTH,
Daniel
jeff6868 wrote
>
> Hi everyb
3 matches
Mail list logo
