#tree-E,P-E,st source("MST.R") par(ask=F) Np<-100 #MST何回か HN<-200 df<-3500 N<-100 #バリアントの数 i ka<-600 kr<-c(0,10^(2:6)) lkr<-length(kr) syu<-matrix(1,Np,lkr) syul<-matrix(0,Np,lkr) colnames(syu)<-kr #allele frequency afb<-0 for (i…