假设我有n个统一(0,1)随机变量并调用y。然后我有b个统一(0,sum(y))随机变量叫做x。现在我根据y的累积和创建了一个带有终点的区间,并希望将b随机x存储在矩阵的y对应的n个区间中。我希望代码在x落入区间时将1的值添加到矩阵中。所以,假设我有两个y统一的rv,两者恰好都是0.5。所以我有一个区间为(0,0.5),另一个区间为(0.5,1)。现在我带两个x,分别是0.75和0.65。所以第一个应该落入第二个区间,第二个区域也应该落入第二个区间。所以我想要一个显示这个的矩阵。这个简短示例中的2X2矩阵在左上角和右下角将为零,在左下角和右上角为1 + 1(因为我只是希望能够查看任何行或列并查看值是否为在相应的行/列中)。如果这有意义,请告诉我!到目前为止,这是我的代码:
b<-5
n<-10
set.seed(440)
mymat <- matrix(0, n, n)
for(i in 1 : n) {
for(j in 1 : n) {
y <- runif(n)
end_points <- cumsum(c(0, y))
x <- runif(b, 0 , sum(y))
intervals <- findInterval(x, end_points)
}
}
因此,对于第一个x,我们在第1行第7列中将得到0 + 1,在第7行第1列中得到0 + 1.此外,我想将其保留为for循环。提前谢谢!