将采样数据存储在数组中

时间:2016-03-03 13:48:59

标签: arrays r storing-data

 No.       Austria Germany Hungary Poland Spain Italy Portugal Slovenia Belgium France Netherlands Switzerland  UK Norway Denmark Sweden
    1       4.4     4.3     6.8    5.7   4.4   4.6      4.6      4.8     4.4    4.4         4.4         3.0 5.1    4.9     4.4    4.3
    2       4.3     4.2     6.7    5.7   4.4   4.6      4.5      4.7     4.4    4.4         4.3         3.1 5.0    4.8     4.4    4.2
    3       4.4     4.3     6.6    5.6   4.4   4.6      4.5      4.6     4.4    4.4         4.4         3.0 5.0    4.9     4.4    4.3
    4       4.2     4.1     6.7    5.7   4.2   4.4      4.4      4.4     4.3    4.2         4.2         2.9 4.7    4.8     4.2    4.2
    5       4.3     4.2     6.9    5.9   4.3   4.5      4.5      4.5     4.4    4.3         4.3         3.1 4.7    4.7     4.3    4.3
    6       4.2     4.0     7.1    5.8   4.2   4.4      4.3      4.4     4.3    4.2         4.1         2.9 4.5    4.5     4.1    4.1
    7       4.2     4.0     7.6    5.8   4.1   4.4      4.3      4.3     4.2    4.1         4.0         3.1 4.6    4.4     4.1    4.0
    8       4.1     3.8     8.4    6.0   4.1   4.4      4.4      4.3     4.1    4.0         4.0         3.0 4.4    4.3     4.0    3.9
    9       4.3     4.0     8.0    6.0   4.3   4.5      4.5      4.5     4.4    4.3         4.2         3.2 4.6    4.5     4.3    4.1
    10      4.5     4.2     8.1    6.1   4.4   4.7      4.6      4.6     4.5    4.4         4.4         3.3 4.9    4.6     4.4    4.2

这16个欧洲国家/地区的虚拟数据可以通过构建q=4区域以多变量形式撰写:Region-1 (Austria,Germany,Hungry,Poland)Region-2 (Spain, Italy, Portugal, Slovenia)Region-3 (Belgium,France,Netherlands,Switzerland)Region-4 (UK, Norway, Denmark, Sweden)

对于上述数据,我想绘制大小为N=100行的n=5个引导样本,并将采样数据存储在DAT的4维数组(n*p*q*N)中其中p=4是国家/地区,q=4是区域数量。

我正在使用以下代码:

DATA16 <- read.table("20160302.txt",header=TRUE)
n <- 5; p <- 4; q <- 4; N <- 100
DAT    <- array(data=NA,dim=c(n,p,q,N))

for(i in 1:N){
for(j in 1:q){
for(k in 1:p){
            DAT[,,,i] <- DATA16[sample(10,n,replace=TRUE),]
}
}
}

此代码出错

Error in DAT[, , , i] <- DATA16[sample(100, n, replace = TRUE), ] : 
  incorrect number of subscripts

如何更正此代码?请帮帮我。

0 个答案:

没有答案