如何在r的逻辑回归方法中将p值格式化为4位数字

debugcn 发表于 Dev

king_arrival27

我有一个逻辑回归输出和赔率，格式为赔率：

  variables   oddsratio    2.5%      97.5%      p-value
      height  0.9810833  0.9724681  0.9898     2.197e-05 ***
      weight  0.9889900  0.9794387  0.9986     0.025356 *  
        BMI   1.0097044  1.0029179  1.0165     0.005005 **

但是我想要以十进制格式输出p值

pvalue
0.001
0.025
0.005

杰伊

对于仅控制台输出，请阅读问题注释中的@EliBerkow链接。对于表输出或在函数中使用，可以使用formatC()。

例

(output <- summary(glm(I(Sepal.Length > 5.7) ~ Sepal.Width + Species, iris, 
                       family=binomial()))$coefficients)
#                     Estimate Std. Error   z value     Pr(>|z|)
# (Intercept)       -15.846295   3.893382 -4.070059 4.700128e-05
# Sepal.Width         3.274112   0.970872  3.372342 7.453193e-04
# Speciesversicolor   7.149531   1.533740  4.661501 3.139108e-06
# Speciesvirginica    9.283911   1.630188  5.694992 1.233776e-08

全部包含四位数：

as.data.frame(apply(output, 2, formatC, format="f", digits=4))
#                   Estimate Std. Error z value Pr(>|z|)
# (Intercept)       -15.8463     3.8934 -4.0701   0.0000
# Sepal.Width         3.2741     0.9709  3.3723   0.0007
# Speciesversicolor   7.1495     1.5337  4.6615   0.0000
# Speciesvirginica    9.2839     1.6302  5.6950   0.0000

不同位数：

setNames(data.frame(output[,-4], formatC(output[,4], format="f", digits=4)), colnames(output))
#                     Estimate Std. Error   z value Pr(>|z|)
# (Intercept)       -15.846295   3.893382 -4.070059   0.0000
# Sepal.Width         3.274112   0.970872  3.372342   0.0007
# Speciesversicolor   7.149531   1.533740  4.661501   0.0000
# Speciesvirginica    9.283911   1.630188  5.694992   0.0000

本文收集自互联网，转载请注明来源。

如有侵权，请联系[email protected] 删除。