我的数据是一组32个样本的32个3D地标坐标的数组（p x k x n）。我有几个因素来定义这些标本。我现在想用这些因素中的一个来划分这组标本，然后再应用mshape。
我可以使用split(F.gpa$coords, F.gpa$status)将坐标分成相应的组，但是这样我只得到所述组中的地标列表，但是失去了结构（我再次需要一个（pxkxn）矩阵，宁愿保留dimnames和其他因素）。
如果有任何帮助，那就是我数据的结构（我对R来说比较新，所以我可能不太擅长解释）：
$ coords: num [1:32, 1:3, 1:36] 0.1604 0.1371 0.1212 0.1058 0.0888 ...
..- attr(*, "dimnames")=List of 3
.. ..$ : NULL
.. ..$ : NULL
.. ..$ : chr [1:36] "106L" "106R" "124L" "124R" ...
$ Csize : Named num [1:36] 26.3 25.9 19.6 21 23.9 ...
..- attr(*, "names")= chr [1:36] "106L" "106R" "124L" "124R" ...
$ ind : Factor w/ 18 levels "22","44","55",..: 7 7 8 8 9 9 10 10 11 11 ...
$ side : Factor w/ 2 levels "L","R": 1 2 1 2 1 2 1 2 1 2 ...
$ sex : Factor w/ 2 levels "1","2": 1 1 2 2 2 2 1 1 1 1 ...
$ status: Factor w/ 3 levels "1","2","3": 1 1 1 1 2 2 1 1 1 1 ...
数据太大而无法在此处发布，但可以找到here
我很感激任何线索，谢谢你，
麦克拉