我正在尝试从数据集中模拟r中样本比例的采样分布。这是我到目前为止的代码:
library(MASS)
phat=replicate(100, sample (Aids2$status, 50, replace = FALSE, prob = NULL), simplify="matrix")
table(phat)["A"]/50
我认为这会将数组每列中“A”的观察次数除以50,所以我可以得到100%的“A”,我可以将它们放在直方图中,但它只取总数在整个阵列中观察“A”并将其除以50得到像39.8这样的东西并不是很有用。我认为我通过将phat数组转换为表格然后尝试在结果表上使用prop.table来获得最接近的结果,但这不起作用。这是代码:
newtable = ftable(phat)
prop.table(newtable, margin=NULL)
再次,我只需要在每个50个观察样本中获得一个百分比的“A”就可以制作直方图,所以如果有人对此有任何想法,那将非常感激。