To whom it may concern,
I happened to run the following R code just to check the layout of the
output, but found that the code doesn't work the way I thought it should
work.
''
> lm(rnorm(100) ~ rnorm(100))
Call:
lm(formula = rnorm(100) ~ rnorm(100))
Coefficients:
(Intercept)
-0.07966
Warning messages:
1: In model.matrix.default(mt, mf, contrasts) :
the response appeared on the right-hand side and was dropped
2: In model.matrix.default(mt, mf, contrasts) :
problem with term 1 in model.matrix: no columns are assigned
"
It appears that rnorm(100) produces the same array of numbers on both sides
of the ~ sign.
It can be further verified by having the same error message if we do x <-
rnorm(100) and lm(x ~ x).
I would expect the two rnorm(100) functions in the lm function return two
different arrays of numbers, but am open to hear reasons from the other
side.
Thanks,
--
*Xu Tian*
[[alternative HTML version deleted]]
______________________________________________
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
