On Jul 11, 2005, at 11:47 AM, S.O. Nyangoma wrote: > Hi there, > If I do an lm, I get p-vlues as > > p-value: < 2.2e-16 > > Suppose am interested in exact value such as > > p-value = 1.6e-16 (note = and not <) > > How do I go about it? > > stephen >
I think you're seeing a very small p-value after it has been formatted by format.pval. Consider: options(digits = 6) format.pval(0) # [1] "< 2.2e-16" format.pval(0, eps = 1e-20) # [1] "< 1e-20" you could access the p-values you're interested following this example: x <- 1:100 y <- x + rnorm(100, 0, 1) f <- lm(y ~ x) ff <- summary(f) print(ff) ff$coefficients Pvals <- ff$coefficients[, 4] print(Pvals, digits = 20) Hope this helps, Stephen ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
