1,没有任何影响,因为实际上这个代码里没有用到这个G矩阵,这个是在GBLUP里的,不影响
2,这个是我们的疏忽,感谢提醒,使用这两行代码替换
sets <- as.numeric(cut(1:length(y), breaks = 10, labels = FALSE))
sets <- sets[order(runif(nrow(X)))]
3,不绝对,和几个因素有关, 其中最重要的是h2和Ne
要是都是平均的,一般为了精度最低差不多是300个样本,
高h2低Ne的看情况可能可以缩减到200个