rの多変量時系列データをシミュレートするrコードを知りたいです。正確には、Genenet Package にある時系列データのサブセットをシミュレートしたいと思います。
library(Genenet)
data(arth800)
subset = c(60, 141, 260, 333, 365, 424, 441, 512, 521, 578, 789, 799)
arth12 = arth.expr[, subset]
シミュレートされたデータを次のようにしたいと思います
> str(arth12)
num [1:22, 1:12] 6.59 6.33 7.12 6.9 7.22 ...
- attr(*, "dimnames")=List of 2
..$ : chr [1:22] "0-1" "0-2" "1-1" "1-2" ...
..$ : chr [1:12] "265768_at" "263426_at" "260676_at" "258736_at"
> arth12[1:8,1:8]
265768_at 263426_at 260676_at 258736_at 257710_at 255764_at 255070_at 253425_at
0-1 6.589464 8.708049 6.802193 2.560715 5.958843 8.695228 8.556889 6.137504
0-2 6.329124 9.084011 6.923625 4.145677 5.632268 8.421644 9.479982 6.228819
1-1 7.116864 8.865733 7.795066 4.832890 5.289097 8.994071 8.355792 6.705978
1-2 6.898450 8.792790 7.941341 5.061776 5.741467 8.657854 9.210915 6.982994
2-1 7.224002 8.910793 7.991522 5.053111 5.402586 8.757223 8.178914 6.556123
2-2 6.993221 8.978710 8.093286 5.026800 5.812498 8.868823 8.981567 6.852998
4-1 7.266787 9.261390 8.396177 5.251719 6.133399 9.241745 7.596935 7.074677
4-2 7.483816 9.300124 8.137504 5.551516 5.794416 9.110092 7.914684 6.904484