r - 使用模拟从 t 表中创建一条线（在 R 中）

Question

在为分布运行模拟之后，我将如何从 R 中的 t 表创建一条线？本质上，我想仅使用从正态分布中的随机样本计算的值来执行 qt 函数，而不是使用置信水平作为输入。

我已经对 at(degrees of freedom "df") 分布进行了 5000 次模拟，我想获取这些值并使用它们在 t 表中重新创建与 df 对应的行。试图弄清楚它有困难。

任何帮助是极大的赞赏！

Answer 1

quantile 函数将为您提供与我使用的 t 表中的值等效的值：

quants <- 1-c(0.2, 0.1, 0.05, 0.025, 0.01, 0.005, 0.001)
n <- 15

sims <- replicate( 5000, {x <- rnorm(n); sqrt(n)*mean(x)/sd(x)} )

(t.table <- quantile(sims, quants) )

## compare

rbind( t.table, qt( quants, n-1 ) )

pt(t.table, n-1)

这些是左边的区域，常见的 t 表给出了右边的区域，我们可以更改标签以匹配它names(t.table) <- 1-quants。更改quants以匹配您的特定表格（如果您想要更接近的匹配，请增加模拟次数）。

在这种情况下，可以通过一步模拟整个法线矩阵然后使用rowMeans.