2012-06-19 20 views
20

से मानक त्रुटि उत्पादन हम से एक एल एम वस्तु मिला है और मानक त्रुटिआर: एल एम वस्तु

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

मैं समारोह सारांश, नाम और गुणांक पता निकालना चाहते हैं। हालांकि, सारांश मानक त्रुटि को मैन्युअल रूप से एक्सेस करने का एकमात्र तरीका प्रतीत होता है। क्या आपको कोई विचार है कि मैं केवल आउटपुट कैसे कर सकता हूं?

धन्यवाद!

उत्तर

16

से summary समारोह सिर्फ एक अनुसंधान सूची है के उत्पादन में। तो आप सभी मानक सूची संचालन का उपयोग कर सकते हैं। उदाहरण के लिए:

#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 संदर्भित करता है ** अवशिष्ट ** मॉडल का मानक विचलन, व्यक्तिगत गुणांक का मानक विचलन नहीं। समारोह 'एम $ सिग्मा' सिग्मा (फिट) 'से मेल खाता है, [यहां] देखें (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 देखें।

संबंधित मुद्दे