中介 R 包 p 值：获得更多有效数字的解决方法？

我正在运行多项中介分析，需要纠正多重比较。然而，中介包提供的 p 值似乎已四舍五入存储，正如 @Roland 关于此问题的非常有用的答案中所观察到的：

R中介包：逗号后面的数字 https://stackoverflow.com/questions/53850958/r-mediation-package-digit-behind-comma

有没有人找到一种解决方法来实际提取 p 值的更多有效数字？

这些都是用 2 位数字打印的，因为这就是该方法计算的全部内容。这是一项随机分析，因此 p 值只是满足某些条件的模拟的 50 个分数（在参考答案中）。如果你想要更多的数字，你需要运行sims设置为比 50 大得多的值。执行此操作后，您可以使用与该帖子中相同的技术或更简单地更改数字设置，如下所示：

1. 执行这个：

   print.summary.mediate

2. 使用fix(print.summary.mediate)编辑源，并更改行

   printCoefmat(smat, 数字 = 3)

到您想要的任何位数。我选择了6。

然后运行代码：

data(jobs)
b <- lm(job_seek ~ treat + econ_hard + sex + age, data=jobs)
c <- lm(depress2 ~ treat + job_seek + econ_hard + sex + age, data=jobs)
 
contcont <- mediate(b, c, sims=1000, treat="treat", mediator="job_seek")
summary(contcont)

我得到这个输出：

Causal Mediation Analysis 

Quasi-Bayesian Confidence Intervals

                 Estimate 95% CI Lower 95% CI Upper p-value
ACME           -0.0166933   -0.0404983      0.00735   0.168
ADE            -0.0412011   -0.1278011      0.04465   0.348
Total Effect   -0.0578944   -0.1449898      0.02716   0.196
Prop. Mediated  0.2313315   -1.7557486      1.89099   0.288

Sample Size Used: 899 


Simulations: 1000 

您只能在 p 值上看到 3 位数字，因为这就是您通过 1000 次模拟得到的全部结果：n/1000始终有 3 位数字。如果你选择sims不是整数， 你会看到更多的数字：但是后面的数字将毫无价值，它们除了像 123/456 这样的分数没有很好的小数扩展这一事实之外，并没有真正表示任何东西。

你的会有所不同，因为这是随机分析。