我想获得一个矩阵数组，如下所示

Question

我想生成一个数据矩阵数组，其中每个矩阵将引用不同的场景，出于示例的目的，我只包含 1 个

p13=0.493;p43=0.325;p25=0.335;p35=0.574;p12=0.868
std_e2=sqrt(1-p12^2);std_e2
std_e3=sqrt(1-(p13^2+p43^2));std_e3
std_e5=sqrt(1-(p25^2+p35^2+2*p25*p35*(p13*p12)));std_e5
scenario_1<-matrix(c(3,0,0,0,0),ncol = 5,nrow = 1);scenario_1

genereting_fuction<- function(n,scenario){ 
  sample <- vector("list") 
  for (i in scenario){
    x1=rnorm(n)+scenario[i,1]
    x4=rnorm(n)+scenario[i,4]
    x2=x1*p12+std_e2*rnorm(n)+scenario[i,2]
    x3=x1*p13+x4*p43+std_e3*rnorm(n)+scenario[i,3]
    x5=x2*p25+x3*p35+std_e5*rnorm(n)+scenario[i,5]
    sample[[i]]=cbind(x1,x2,x3,x4,x5)
    colnames(sample[[i]])<-c("x1","x2","x3","x4","x5")
  }
  sample
} 
array_scenari<-array(dim=c(5,5,2));array_scenari

for(j in 1:nrow(scenario_1)){
  set.seed(1234)
  dati_prova<- sapply(rep(1,5), function(x) genereting_fuction(x,scenario_1),simplify = 'array');dati_prova
  dati_prova<-do.call(rbind.data.frame, dati_prova);dati_prova<-as.matrix(dati_prova);dati_prova
  dati_prova<-as.matrix(dati_prova)
  array_scenari[,,j]<-dati_prova[]
}

我无法理解为什么它不起作用并给我一个错误

Answer 1

如果没有有效的可重现示例，很难给出具体答案，因为您的函数包含外部定义的变量。然而，错误的根源是明确的。当你用 array() 创建一个空数组时，它有一个固定的维度：

matrix_classification <- array()

dim(matrix_classification)
#> [1] 1

如果你尝试写入它的第三个维度，你会得到一个错误：

k <- 1

matrix_classification[, , k] <- "x"
#> Error in matrix_classification[, , k] <- "x": incorrect number of subscripts

如果你想写入一个数组，你应该定义首先是它的尺寸。例如，以下代码创建一个空的 5 x 5 x 5 数组：

matrix_classification <- array("", dim = c(5, 5, 5))

dim(matrix_classification)
#> [1] 5 5 5

如果我们想将矩阵写入第 k 个切片，我们可以这样做：

matrix_classification[, , k] <- matrix(sample(letters, 25), nrow = 5)

matrix_classification[,,1]
#>      [,1] [,2] [,3] [,4] [,5]
#> [1,] "a"  "d"  "k"  "t"  "c" 
#> [2,] "f"  "b"  "n"  "m"  "s" 
#> [3,] "u"  "q"  "y"  "o"  "j" 
#> [4,] "l"  "g"  "h"  "w"  "v" 
#> [5,] "r"  "i"  "e"  "p"  "z"

^{创建于 2022 年 3 月 6 日，由reprex 包 (v2.0.1)}