2016-06-30 7 views
0

私は一般的な質問があります。とにかく、私がRの回帰で使用された観測を識別(またはタグ付け)できるのですか?rの回帰サンプルを特定する

lligator = data.frame(lnLength = c(3.87, 3.61, NA, 3.43, 3.81, 3.83, 3.46, 3.76, 
3.50, 3.58, 4.19, 3.78, 3.71, 3.73, 3.78),lnWeight = c(4.87, 3.93, 6.46, 3.33, 4.38, 4.70, 3.50, 4.50,NA, 3.64, 5.90, 4.43, 4.38, 4.42, 4.25)) 

t.test=lm(lnWeight ~ lnLength, data = alligator) 

どの観測が使用されているかを示す別の列を持つデータフレームを作成したいと考えています。私はどのように

na.omit() 

na.exclude() 

.completecases 

仕事を知っていると私は、回帰を行うためにそれらを使用することができます。しかし、私が探しているのは、どの観測を使用するかを示す指標を作成する方法です。 Eに似たStataのユーザー何か(サンプル)

答えて

0

のそれらのために私が正しく理解していれば、あなたは回帰時に除外されたインデックスのベクトルを取得し、標識変数を計算するためにそれを使用するna.action()を使用することができます。

alligator$used <- !seq_len(nrow(alligator))%in%na.action(t.test); 
alligator; 
## lnLength lnWeight used 
## 1  3.87  4.87 TRUE 
## 2  3.61  3.93 TRUE 
## 3  NA  6.46 FALSE 
## 4  3.43  3.33 TRUE 
## 5  3.81  4.38 TRUE 
## 6  3.83  4.70 TRUE 
## 7  3.46  3.50 TRUE 
## 8  3.76  4.50 TRUE 
## 9  3.50  NA FALSE 
## 10  3.58  3.64 TRUE 
## 11  4.19  5.90 TRUE 
## 12  3.78  4.43 TRUE 
## 13  3.71  4.38 TRUE 
## 14  3.73  4.42 TRUE 
## 15  3.78  4.25 TRUE 

同等おそらく速い方法:

alligator$used <- TRUE; 
alligator$used[na.action(t.test)] <- FALSE; 
+0

はあなたに感謝します。だからna.action()もlmの別の戻りオブジェクトですか?それは正しい? – Yashar

+0

これは正しいです。返されるオブジェクトの名前付きリストコンポーネントの1つで、 '$ na.action'として直接アクセスできます。 – bgoldst