rm(list=ls()) setwd("C:/Users/shuzhang/Desktop/Research/Zsl/PS/spatial/paper/Scientific_Report/program/") source("subfunction.R") pc0=read.table(file="pc2a.txt",header=FALSE, skip=1) pc0=as.matrix(pc0) sample=800 K=20 M=0.01 L=100 rep=100 rep1=5 mm=4 N=L*rep pc=0 xx=matrix(0,sample,rep*L) pvper=matrix(0,N,mm) power=matrix(0,6,mm+3) nn0=0 for(model in c(1,2)) { for(RR in c(3)) { nn0=nn0+1 for(i in 1:rep) { print(paste("model=",model,"RR=",RR,"rep=",i)) n1=(i-1)*L+1 n2=i*L xx[,n1:n2]=generate_genotype(sample, M, K, L) } percent=0.3 y=trait_nulla(sample, K, model, RR,2,4) y=quan_to_qual(y,percent) pvper[,1]=uncorrect_single(y,xx) pvper[,2]=GC(pvper[,1]) pvper[,3]=pc_linear_single(y, xx, pc0[,1:10]) h1=new_choose_h(y, xx, pc0, 10) pvper[,4]=pc_nonp_single(y, xx, pc0, 10, h1) for(jj in 1:3) { nn=jj+(nn0-1)*3 if(jj==1) alpha0=0.05 else if(jj==2) alpha0=0.01 else alpha0=0.001 power[nn,1]=model power[nn,2]=RR power[nn,3]=alpha0 for(k in 1:mm) { power[nn,(k+3)]=sum(pvper[,k]