как использовать команду cbind для получения значений переменных для каждого значения rho

Question

Я работаю с r, и я хочу использовать команду cbind, чтобы получить результаты для каждого значения rho. Но получить одно значение для каждого значения rho вместо n = 3

muy=100; mux=50; sigy=10; sigx=7
n=3; sim=5    
rho=seq(-1,1,0.5)
ybark <- numeric(sim)
for(j in 1:length(rho)) {
  for(i2 in 1:sim) {
    z1 = rnorm(n)
    z2= rnorm(n)
    x = sqrt(1-rho[j]^2)*sigx*z1[j] + rho[j]*sigx*z2[j] + mux
    y = sigy*z2[j] + muy
    ybark[i2] = (mean(y)*mean(x))/mux
  }
}

Answer 1

muy=100 ;mux=50 ; sigy=10; sigx=7
n=3;sim=5;rho=seq(-1,1,0.5)
ybark<-list()
for(j in 1:length(rho)){
    for(i2 in 1:sim)
    {
        z1 = rnorm(n)
        z2= rnorm(n)
        x = sqrt(1-rho[j]^2)*sigx*z1[j] + rho[j]*sigx*z2[j] + mux
        y = sigy*z2[j] + muy
        ybark[[paste(j,i2,sep="_")]] = (mean(y)*mean(x))/mux
    }
}
n<-length(ybark)/length(rho);grp<-rep(1:n,each=sim)
tapply(ybark,grp,function(x)matrix(x,nrow=1))

это то, что вы хотите?

$`1`
     [,1]     [,2]     [,3]     [,4]     [,5]   
[1,] 99.88803 92.99583 101.0367 98.73194 99.4599

$`2`
     [,1]     [,2]     [,3]     [,4]     [,5]    
[1,] 111.9889 87.19593 111.1355 86.76358 91.04678

$`3`
     [,1]     [,2]    [,3]     [,4]     [,5]   
[1,] 135.2721 93.4944 71.71243 121.4449 84.8635

$`4`
     [,1] [,2] [,3] [,4] [,5]
[1,] NA   NA   NA   NA   NA  

$`5`
     [,1] [,2] [,3] [,4] [,5]
[1,] NA   NA   NA   NA   NA