如何访问数组的元素？

Question

我正在处理包含对象的数组。当我比较两个不同数组的元素时，它在cell_list [x4，3，1]中给出了一个错误：错误的数量。的代码段如下所示，

cnt2 <- 25 
cell_list <- array(0, c(cnt2,6,1)) 
nbr_list <- array(0, c(25,6,1)) 
insert_nbrlist = function(nbr_list, cell_list, cnt1) 
    { 
    for(x3 in 1:cnt1) 
     for(x4 in 1:cnt2) 
     { 
     if((nbr_list[x3,3,1] == cell_list[x4,3,1]) && (nbr_list[x3,4,1]== cell_list[x4,4,1])) 
     { 
      theta_diff <- cell_list[x4,1,1] - nbr_list[x3,1,1] 
      phi_diff <- cell_list[x4,2,1] - nbr_list[x3,2,1] 
     } 
     else 
     { 
      cnt2 <- cnt2+1 
      cell_list <- rbind(cell_list[,,1],nbr_list[x3,,1]) 
     } 
    } 
    return(cell_list) 
    } 

这里insert_nbrlist是由主程序调用的函数。 Inner for循环将被执行一次，但是对于第二次执行，即当x4 = 2时，它给出了错误cell_list [x4,3,1]中的错误：不正确的维数。请有人能帮助我吗？ 这里cnt1被传递给函数insert_nbrlist。 cnt1是nbr_list的table1中存在的行数并且是恒定的。 cnt2表示cell_list的table1中存在的行数，并且它不断增加。 cnt2我没有传递给函数，因为它在该函数中正在增加。调用该函数不会产生任何错误，但是该错误是在函数insert_nbrlist中产生的。

我想以另一种方式提出这个问题。我有一个数组，它初始化为cell_list < - array（0，c（25,6,1]）。在程序中，我再向cell_list的table1中添加一行。向数组添加一行后，我无法我得到以下错误，错误在cell_list [x4,3,1]：不正确的维数。如何在向该数组添加一行后访问数组中的元素

我重写代码与主程序

func <- function() 
{ 

    cnt2 <- 25 

    insert_nbrlist = function(nbr_list, cell_list) 
    { 
    for(x3 in 1:25) 
     for(x4 in 1:cnt2) 
     { 
     if((nbr_list[x3,3,1] == cell_list[x4,3,1]) && (nbr_list[x3,4,1]== cell_list[x4,4,1])) 
     { 
      theta_diff <- cell_list[x4,1,1] - nbr_list[x3,1,1] 
      phi_diff <- cell_list[x4,2,1] - nbr_list[x3,2,1] 
     } 
     else 
     { 
      cnt2 <- cnt2+1 
      cell_list <- rbind(cell_list[,,1],nbr_list[x3,,1]) 
     } 
    } 
    return(cell_list) 
    } 

    cnt <- 0 
    m <- 50 
    r1 <- 10 
    dtheta <- 0.5 
    dphi <- 0.5 

    cell_list <- array(0, c(cnt2,6,1)) 
    nbr_list <- array(0, c(25,6,1)) 

    repeat 
    { 
    cnt <- cnt+1 
    s <- runif(2,-pi/2 : pi/2) 

    theta <- s[1] 
    phi <- s[2] 
    cnt1 <- 0 
    for(x1 in -2:2) 
     for(x2 in -2:2) 
     { 
     cnt1 <- cnt1+1 
     theta1 <- theta+x1*dtheta 
     phi1 <- phi +x2*dphi 

     x <- round(r1*(sin(theta1) + sin(theta1+phi1))+m/2, digits=0) 
     y <- round(r1*(cos(theta1)+ cos(theta1+phi1))+m/2,digits=0) 

     if((theta1 == theta) && (phi1 == phi)) 
      nbr_list[cnt1,,1] <- c(theta1, phi1, x, y, 1,1) 
     else 
      nbr_list[cnt1,,1] <- c(theta1, phi1, x, y, 1,0) 
     } 

     if(cnt==1) 
     cell_list <- nbr_list 

     else 
     cell_list <- insert_nbrlist(nbr_list, cell_list) 

     if(cnt == 50) 
     break 

    } 

}  

func() 

我想补充简单的代码。

test<- function() 
{ 
    sink("G:/rvma/test/test",append=FALSE) 
    nbr_list <- array(0, c(5,5,1)) 
    nbr_list[1,,1] <- 15 
    nbr_list[2,,1] <- -15 
    print(nbr_list) 
    n <- c(1,2,3,4,5) 
    nbr_list <- rbind(nbr_list[,,1],n) 
    print(dim(nbr_list)) 
    print(nbr_list) 
    print(nbr_list[,,1]) 
    sink() 
}  

我能够在rbind之前访问数组中的元素。在rbind之后，数组的维度从5X5X1更改为6X5X1，如果指定了print（nbr_list），我可以打印数组，但如果指定print（nbr_list [,, 1]），则无法打印相同的数组，在这里，我得到nbr_list [，，1]中的错误错误：不正确的维数。即使我无法在rbind之后访问元素。

Answer 1

我的猜测是你的问题是这条线在else块：

cnt2 <- cnt2+1 

这会推高cnt2值超过25，而且因为它是一个25x6x1阵列，这将是失界限。

编辑

我怀疑rbind是不是一个3-d阵列上正确工作。有一个包，abind你可以使用，而不是rbind。

更多on this answer。