我有6个大小为100 * 50的矩阵(例如m1-m6)和6个大小为100 * 100(m7-m12)的矩阵。这些矩阵嵌套有两个因子,如F和G.我想创建一个反映这种设计的数据结构:
F1 G1 m1
F1 G2 m2
F1 G3 m3
F1 G1 m4
F1 G2 m5
F1 G3 m6
F2 G1 m7
F2 G2 m8
F2 G3 m9
F2 G1 m10
F2 G2 m11
F2 G3 m12
我想将这些结构用于ANOVA并绘制结果。每个矩阵由超过100次重复的值组成。 m1-m12矩阵存储在csv
个文件中。因此,我需要导入矩阵并创建数据结构。我尝试array
或list
但我找不到有效的并且正确的方法。
有什么想法吗?
答案 0 :(得分:1)
您可以使用array
,首先我创建数据结构:
A1 = array(0,dim=c(100,50,6))
A2 = array(0,dim=c(100,100,6))
我假设您有2个文件名列表,list_files1
和list_files2
:
lapply(seq_along(list_files1),function(x){
A1[,,x] <- read.csv(list_files1[[x]])
})
lapply(seq_along(list_files2),function(x){
A2[,,x] <- read.csv(list_files2[[x]])
})