-1 
set.seed(98234)
y <- rnorm(100)
x <- rnorm(100)
lm0 <- lm(y ~ x)
lm1 <- lm(I(y) ~ I(x))

すべてが完全に正常に動作し、「lm0」が発生すると予想されることに同意できると思います。lm1に等しいlm0(係数による判断)。そうです

set.seed(98234)
lm3 <- lm(I(rnorm(100)) ~ rnorm(100))
set.seed(98234)
lm4 <- lm(rnorm(100) ~ I(rnorm(100)))

しかしI()、式の両側または両側にある場合、上記の結果は得られません。

set.seed(98234)
lm2 <- lm(I(rnorm(100)) ~ I(rnorm(100)))
set.seed(98234)
lm5 <- lm(rnorm(100) ~ rnorm(100))

理由はありますか？

4

1 に答える 1

1

引数の左側はコンポーネントに分解されないため、Iグループ化記号は不要です。

lm(I(a+b+c+d/e)~f)

と同じだろう

lm(a+b+c+d/e~f)

一方

lm(a~b+c+d/e)

と同じではありません

lm(a~I(b+c+d/e))

あるいは

lm(a~b+c+I(d/e))
于 2012-12-14T22:29:15.300 に答える