r - 仅当 >0.001 时才打印实际 P 值

Question

考虑这两个 t 检验的 P 值

set.seed(1)
x <- c(rnorm(50,1), rnorm(50, 2))
y <- (c(rep("a", 50), rep("b", 50)))

t.test(x ~ y)$p.value

[1] 1.776808e-07

set.seed(2)
x <- c(rnorm(50), rnorm(50))
y <- (c(rep("a", 50), rep("b", 50)))

t.test(x ~ y)$p.value

[1] 0.3922354

第一个 P 值 <0.001，第二个 P 值 >0.001。如果 P 值恰好 <0.001，我可以让 R 将 P 值打印为 <0.001。如果 P 值恰好 >0.001，我可以让 R 打印出实际的 P 值吗？所以第一个 t 检验的结果应该打印为 "<0.001"并且第二次 t 检验的结果应打印为 0.3922354 .

我正在使用 knitr为我的论文将 R 代码转换为 latex 。我的规则是只有 P 值 >0.001 应该作为实际 P 值打印。

Answer 1

只需使用 if 语句:

p_val <- t.test(x,y)$p.value;
if(p_val>=0.001) {
  print(p_val)
} else {
  print("<0.001")
}