Hi, folks,
x=1:10
y=rep(2:6,2)
lin=lm(y~x)
x=3:12
new=predict(lin,se.fit=T)
#se.fit: the standard error of the predicted means, namely, the square root
of Var( E[y|x] | x)
# How can I generate the variances of the new observations? Namely the
square root of var(y|x),
## Which I think should be much larger than the values from se.fit=T.
The reason why I need to know the estimations of expectation and also
variance is that Y=ln(Z) and I need to know the expectation of Z.
Thanks
Yi
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.