如何生成具有50个自由度的10个chi ^ 2（即“ chi平方”）分布变量的50个观测值？

Question

我正在尝试使用10个chi ^ 2分布变量和50个来生成50个观测值 自由程度。然后，我必须创建一个新变量，该变量是这10个chi ^ 2分布变量中每个变量的平均值，并创建该变量的直方图（这将是10个观测值的样本）。有什么想法吗？

Answer 1

rchisq(10, 50)将根据卡方分布生成10个具有50个自由度的随机数。我们可以做50个重复，然后取每个平均值。或者，我们可以只创建50 * 10个观测值，然后将它们分成10个一组。之后，只需运行hist来获取直方图即可。

obs <- matrix(rchisq(50 * 10, 50), nrow = 10)
hist(colMeans(obs))