2012-02-05 9 views
12

मैं confint मदद पृष्ठआर में confint का उपयोग कर 95% सीआई गणना कैसे की जाती है?

> fit <- lm(100/mpg ~ disp + hp + wt + am, data=mtcars) 
> summary(fit) 

Call: 
lm(formula = 100/mpg ~ disp + hp + wt + am, data = mtcars) 

Residuals: 
    Min  1Q Median  3Q  Max 
-1.6923 -0.3901 0.0579 0.3649 1.2608 

Coefficients: 
      Estimate Std. Error t value Pr(>|t|) 
(Intercept) 0.740648 0.738594 1.003 0.32487 
disp  0.002703 0.002715 0.996 0.32832 
hp   0.005275 0.003253 1.621 0.11657 
wt   1.001303 0.302761 3.307 0.00267 ** 
am   0.155815 0.375515 0.415 0.68147 
--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 0.6754 on 27 degrees of freedom 
Multiple R-squared: 0.8527,  Adjusted R-squared: 0.8309 
F-statistic: 39.08 on 4 and 27 DF, p-value: 7.369e-11 

> confint(fit) 
        2.5 %  97.5 % 
(Intercept) -0.774822875 2.256118188 
disp  -0.002867999 0.008273849 
hp   -0.001400580 0.011949674 
wt   0.380088737 1.622517536 
am   -0.614677730 0.926307310 
> confint(fit, "wt") 
     2.5 % 97.5 % 
> wt 0.3800887 1.622518 

>confint.default(fit,"wt") 
     2.5 % 97.5 % 
wt 0.4079023 1.594704 

> 1.001303 + 1.96*0.302761 
[1] 1.594715 
> 1.001303 - 1.96*0.302761 
[1] 0.4078914 

तो 95% confint.default से प्राप्त सीआई asymptotic सामान्य पर आधारित है में प्रदान की जाती उदाहरण का उपयोग करें। confint के बारे में क्या?

धन्यवाद

उत्तर

17

आप प्रत्येक विधि के लिए कोड देख सकते हैं।

# View code for 'default' 
confint.default 
# View code of lm objects 
getAnywhere(confint.lm) 

अंतर यह है कि डिफ़ॉल्ट सामान्य quantiles का उपयोग करता है और रैखिक मॉडलों के लिए विधि टी quantiles बजाय का उपयोग करता प्रतीत होता है।

+0

धन्यवाद। यह स्पष्ट करता है। स्वतंत्रता की बड़ी संख्या में डिग्री के साथ, सामान्य क्वांटाइल और टी-क्वांटाइल के बीच अंतर नगण्य होना चाहिए। – Tony

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