フェノタイプの単位とジェノタイプの単位は違って、フェノタイプ=2とジェノタイプ=2の意味は違うのだけど、MSTで処理するときに、同じものとしてあつかわれているので、ちょっと微妙になっているかもしれない。困った。リスクの与え方を工夫しないといけない…

source("MST.R") par(ask=F) Np<-100 HN<-500 df<-2000 N<-280 #バリアントの数 i ka<-1000 coa<-c(0,1000,2000,3000,4000,5000) lcoa<-length(coa) syu<-matrix(200,Np,lcoa) syu2<-matrix(0,Np,lcoa) colnames(syu)<-coa nrvcv<-matrix(200,Np,lcoa) afb<-…