Rで「lm」関数を実行するにはどうすればよいのですか?


8

のような数式を実行してみましたy ~ 1。これは切片だけを与えますが、表示される値は切片推定(共変量なし)+ 1に等しいか、切片推定(共変量なし)だけですか?

任意の助けいただければ幸いです。


1
切片は単なる卑劣なものではない(y)
アーロン

唯一式参照するのでy、一定の1どのようにそれは可能性も任意の共変量について知っているだろうか?
whuber

@whuberわからない。Rでは、-1を指定すると、lmは定数なしで回帰を行います。1だけを入力することで、私には知られていない別のオプションをアクティブにすることができます。それが私が求めている理由です。
海の老人。

1
オプションかどうか、プロシージャに共変量を指定しない場合、それを使用する方法はありませんか?
whuber

@whuberあなたはおそらく正しいと思いますが、lm(式、データ)と言うことはできますが、データではlm関数は共変量を見つけることができます...
海の老人。

回答:


14

定数のみで線形回帰を行うと、切片値のみが取得されます。これは、実際には結果の平均にすぎません。

Rでは、

y <- rnorm(1000)
lm(y ~ 1)
# intercept = 0.00965
mean(y)
# Equal to 0.00965

単に平均値を計算するのではなく、回帰方法でそれを行う理由は、簡単な標準誤差を得るためです。

「インターセプト+1」の意味がわからないが、得られるのは結果の平均、つまり何も制御しない場合のyの期待値です。


1
Repmatに感謝します。また、これを行うもう1つの理由は、一部のパッケージではユーザーがベースモデルを定義する必要があることです(例: 'BayesVarSel')。
海の老人。
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.