[R] error code trying to extract second column from coeftest output

Petr PIKAL petr.pikal at precheza.cz
Thu May 17 15:28:35 CEST 2012


Hi

> 
> I want to use the standard error values in the summary that is produced 
using
> coeftest, but I am getting an error code- any ideas?

See what is structure of coeftest object by

str(coeftest(lmodT_WBHO))

and from this you shall deduct how to select second column.

Regards
Petr


> 
> > library(lmtest)
> > coeftest(lmodT_WBHO)
> 
> t test of coefficients:
> 
>     Estimate Std. Error t value  Pr(>|t|) 
> t1W  5.94819    0.17072 34.8410 < 2.2e-16 ***
> t2W  6.56216    0.17438 37.6322 < 2.2e-16 ***
> t3W  6.08252    0.16525 36.8082 < 2.2e-16 ***
> t4W  6.18041    0.17028 36.2949 < 2.2e-16 ***
> t1B  5.50000    0.50566 10.8768 < 2.2e-16 ***
> t2B  5.65000    0.53034 10.6535 < 2.2e-16 ***
> t3B  4.52381    0.51756  8.7406 < 2.2e-16 ***
> t4B  4.38095    0.51756  8.4646 < 2.2e-16 ***
> t1H  5.05000    0.53034  9.5221 < 2.2e-16 ***
> t2H  4.77778    0.55903  8.5465 < 2.2e-16 ***
> t3H  5.52632    0.54412 10.1564 < 2.2e-16 ***
> t4H  4.71429    0.63388  7.4372 2.236e-13 ***
> t1O  5.17647    0.57524  8.9988 < 2.2e-16 ***
> t2O  5.81818    0.50566 11.5060 < 2.2e-16 ***
> t3O  6.50000    0.63388 10.2543 < 2.2e-16 ***
> t4O  5.71429    0.63388  9.0147 < 2.2e-16 ***
> ---
> Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
> 
> > se1 <- coeftest(lmodT_WBHO)$coef[,2]
> Error in coeftest(lmodT_WBHO)$coef : 
>   $ operator is invalid for atomic vectors
> > 
> 
> 
> --
> View this message in context: http://r.789695.n4.nabble.com/error-code-
> trying-to-extract-second-column-from-coeftest-output-tp4630298.html
> Sent from the R help mailing list archive at Nabble.com.
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide 
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list