<0.001

Question

के साथ प्रिंटिंग पी-वैल्यू मुझे आश्चर्य है कि <0.001 प्रतीक कैसे डालें यदि पी-मान से Sweave में उपयोग किया जा सकता है। किसी भी तरह की सहायता की हम सराहना करेंगे। धन्यवाद

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 

Answer 1

दो मुख्य कार्य है कि मैं का उपयोग करें, format.pval और यह एक है कि मैं एक ईकाई और फेरबदल से फट रहे हैं।

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 

मैं यह पसंद है क्योंकि यह pvalues ​​से परिशुद्धता को हटा> .1 (या अगर आप कुछ अलग चाहते हैं जो कुछ भी आप की तरह सीमा; वह यह है कि digits की परवाह किए बिना, यह केवल दो दशमलव स्थानों होता रहता है तो मान है> .1), पीछे शून्यों को रखता है (नीचे उदाहरण देखें), और < में जोड़ता है जैसे आप कुछ स्तर की सटीकता (यहां 0.001) चाहते हैं।

pvalr <- function(pvals, sig.limit = .001, digits = 3, html = FALSE) { 

    roundr <- function(x, digits = 1) { 
    res <- sprintf(paste0('%.', digits, 'f'), x) 
    zzz <- paste0('0.', paste(rep('0', digits), collapse = '')) 
    res[res == paste0('-', zzz)] <- zzz 
    res 
    } 

    sapply(pvals, function(x, sig.limit) { 
    if (x < sig.limit) 
     if (html) 
     return(sprintf('&lt; %s', format(sig.limit))) else 
      return(sprintf('< %s', format(sig.limit))) 
    if (x > .1) 
     return(roundr(x, digits = 2)) else 
     return(roundr(x, digits = digits)) 
    }, sig.limit = sig.limit) 
} 

और उदाहरण:

pvals <- c(.133213, .06023, .004233, .000000134234) 
pvalr(pvals, digits = 3) 
# [1] "0.13" "0.060" "0.004" "< 0.001"