[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