rm(list=ls()) setwd("C:/Users/shuzhang/Desktop/Research/Zsl/PS/spatial/paper/Scientific_Report/program/") source("subfunction.R") pc0=read.table(file="pc2.txt",header=FALSE, skip=1) pc0=as.matrix(pc0) sample=800 K=20 M=0.01 L=100 rep=100 rep1=5 model=2 RR=3 mm=4 N=L*rep pc=2 pvper=matrix(0,N,mm) power=matrix(0,6,mm+4) nn0=0 for(model in 1:2) { for(kk in c(10)) { LL=L*kk xx=matrix(0,sample,rep*LL) nn0=nn0+1 for(i in 1:rep) { print(paste("model=",model,"kk=",kk,"pc=",pc,"rep=",i)) n1=(i-1)*LL+1 n2=i*LL xx[,n1:n2]=generate_genotype(sample, M, K, LL) } y=trait_nulla(sample, K, model, RR,3,4) pvper[,1]=uncorrect_region1(y,xx,kk) pvper[,2]=GC(pvper[,1]) pvper[,3]=pc_linear_region1(y, xx, pc0[,1:10], kk) aa=sample(c(1:(rep*LL)),N,replace=FALSE) h1=new_choose_h(y, xx[,aa], pc0, 10) pvper[,4]=pc_nonp_region1(y, xx, pc0, 10, h1, kk) 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]=kk power[nn,3]=alpha0 power[nn,4]=pc for(k in 1:mm) { power[nn,(k+4)]=sum(pvper[,k]