10,000 列 (SNP 頻度) の df があります。非反復ベクトルを使用してシミュレーション (因子分析) を実行する必要があります。これを行うには、列のサブセットで因子分析を実行し、10 個のグループに分割する必要があります。たとえば、列 1:10、11:20; 21:30。これを手動で指定するには時間がかかるため、それを行う簡単なスクリプトが必要です。これを書いたのですが、うまくいかないようです。各反復をいつ開始および停止するかをRに伝える方法がわかりません。
ind=seq(1,(ncol(df)-10),by=10)
for (i in ind) { start=i;end=i+9; rez = factanal(df,factors=1, start:end) }