「統計学関連なんでもあり」の過去ログ--- 046

No.21130　lmオブジェクトから説明変数名を取り出すには？　　【葉隠】　2014/07/03(Thu) 12:27

いつもお世話になっております。

lmオブジェクトから「説明変数名」を取り出す方法をご教示頂けませんでしょうか？

回帰係数は下記のように取り出せたのですが，「説明変数名」では上手くいきません。

ans <- lm(Y ~ X1 + X2, data=df)
coef <- data.frame(coef(ans))
k <- coef[i,]

variablesというキーワードまではたどり着いたのですが。。。

どうか，ご教示いただけますと幸いです。

No.21132　Re: lmオブジェクトから説明変数名を取り出すには？　　【葉隠】　2014/07/03(Thu) 14:02

連続投稿失礼いたします。

何とか下記のように取り出せてはいるのですが，不要な部分も付いてきてしまいます。
（"6 Levels: (Intercept) X2..."の部分）

何か上手い方法は無いでしょうか。。。

vars <- data.frame(names(coef(ans)))
vars.X1 <- vars[2,]

> vars.X1
[1] X1
3 Levels: (Intercept) X2...

No.21133　Re: lmオブジェクトから説明変数名を取り出すには？　　【青木繁伸】　2014/07/03(Thu) 16:21

ans <- lm( )
names(coef(ans)[-1])
で得られるベクトルが説明変数の名前になるでしょう

No.21139　Re: lmオブジェクトから説明変数名を取り出すには？　　【葉隠】　2014/07/04(Fri) 13:28

青木先生

ご教示ありがとうございます。

説明変数を個別に取り出したかったのですが，下記のようにすることで，上手くいきました。

本当にありがとうございます。

今後ともご指導のほど，お願い致します。

ans <- lm( )
list <- names(coef(ans)[-1])

X1 <- list[1]
X2 <- list[2]