2012-06-19 16 views
20

私たちは、私は、関数の要約、名前と係数を知っているからLMオブジェクトを持って、標準エラーにR:標準エラー出力LMオブジェクトから

lm_aaa<- lm(aaa~x+y+z) 

を抽出したいです。 しかし、概要が標準エラーに手動でアクセスする唯一の方法だと思われます。 どのように出力するだけですか?

ありがとうございました!

答えて

16

summary関数からの出力は、単にR リストあります。したがって、すべての標準リスト操作を使用できます。例:

#some data (taken from Roland's example) 
x = c(1,2,3,4) 
y = c(2.1,3.9,6.3,7.8) 

#fitting a linear model 
fit = lm(y~x) 
m = summary(fit) 

mオブジェクトまたはリストには多数の属性があります。あなたは、ブラケットや名前のアプローチを使用してそれらにアクセスすることができます

m$sigma 
m[[6]] 

について知るための便利な機能は、strあります。この関数は、すなわち

str(m) 
+0

しかし、@ csgillespieが指すものは、モデルの**残余**標準偏差であり、個々の係数の標準偏差ではありません。 'm $ sigma'関数は' sigma(fit) 'に対応しています。[here](https://stat.ethz.ch/R-manual/R-devel/library/stats/html/sigma.html)を参照してください。私は質問が本当に**個**係数の標準偏差に関するものだと信じています。 –

8
#some data 
x<-c(1,2,3,4) 
y<-c(2.1,3.9,6.3,7.8) 

#fitting a linear model 
fit<-lm(y~x) 

#look at the statistics summary 
summary(fit) 

#get the standard error of the slope 
se_slope<-summary(fit)$coef[[4]] 
#the index depends on the model and which se you want to extract 

#get the residual standard error 
rse<-summary(fit)$sigma 
11

すべてのパラメータのための標準的なエラーのリストを取得するには、属性オブジェクトの要約を提供して他の人が指摘したように、あなたは

summary(lm_aaa)$coefficients[, 2] 

を使用することができ、str(lm_aaa)が教えてくれますあなたのモデルから抽出できるすべての情報を取得できます。

0

あなたはモデルの標準誤差/偏差を取得する必要はありませんが、その代わり個々の係数の標準誤差/偏差は、の詳細については

# some data (taken from Roland's example) 
x = c(1, 2, 3, 4) 
y = c(2.1, 3.9, 6.3, 7.8) 

# fitting a linear model 
fit = lm(y ~ x) 

# get vector of all standard errors of the coefficients 
coef(summary(fit))[, "Std. Error"] 

を使用している場合モデルの標準誤差/偏差については、hereを参照してください。係数の標準誤差/偏差の詳細については、hereを参照してください。

関連する問題