我需要帮助创建一个for循环来使用R填充5X5表。每行将是一个没有替换的观察。数字范围是1:75,并且我分别具有这些数字的概率。那么我该如何创建一个随机数生成代码,考虑每个数字的具体概率?
答案 0 :(得分:1)
以下是一些示例数据:
A <- seq_len(75)
B <- rpois(75, 3)
B <- B / sum(B)
所以现在B
是A
中每个元素的概率向量。
要提取25个样本,只需使用sample(A, size = 25, replace = FALSE, prob = B)
即可。像往常一样填充矩阵MAT <- matrix(sample(A, size = 25, replace = FALSE, prob = B), nrow = 5)
。