我想知道如果p值小于要在Sweave中使用的<0.001
时如何放置0.001
符号。>>
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) group <- gl(2, 10, 20, labels = c("Ctl","Trt")) weight <- c(ctl, trt) lm.D9 <- lm(weight ~ group) summary(lm.D9)$coef Estimate Std. Error t value Pr(>|t|) (Intercept) 4.8465 0.1557174 31.12368 4.185248e-17 group1 -0.1855 0.1557174 -1.19126 2.490232e-01
所需的输出
Estimate Std. Error t value Pr(>|t|) (Intercept) 4.8465 0.1557174 31.12368 <0.001 group1 -0.1855 0.1557174 -1.19126 0.249
我想知道如果p值小于0.001以便在Sweave中使用,如何放置<0.001符号。 ctl
[我使用两个主要功能,format.pval
,我从gforge撕下并调整了这个功能。
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm.D9 <- lm(weight ~ group)
tmp <- data.frame(summary(lm.D9)$coef)
tmp <- setNames(tmp, colnames(summary(lm.D9)$coef))
tmp[ , 4] <- format.pval(tmp[ , 4], eps = .001, digits = 2)
tmp
# Estimate Std. Error t value Pr(>|t|)
# (Intercept) 5.032 0.2202177 22.85012 <0.001
# groupTrt -0.371 0.3114349 -1.19126 0.25