general.min2<-function(n1=n1,n2=n2,n3=n3,n4=n4,beta=beta,x=x,N.null=N.null,N.pow=N.power,level=level){

 N=N.power
 if(beta==0) N=N.null

MIN=CATT=Pearson=NULL

for(i in 1:N){
genJ.data(n1=250,n2=250,n3=250,n4=250,beta=beta,x=x)->mydata
T1=CATK2(mydata,x0)
T2=mychisq(mydata)
CATT=c(CATT,T1)
Pearson=c(Pearson,T2)
P1=1-pchisq(T1,df=1)
P2=1-pchisq(T2,df=J-1)
MIN=c(MIN,min(P1,P2))
}

if(beta==0) {cmin2=quantile(MIN2,level);print(cmin2=cmin2)}

pow.Ps=mean(Pearson>qchisq(0.95,df=J-1))
pow.CAT=mean(CATT>qchisq(0.95,df=1))
pow.MIN=mean(MIN<cmin2)

cbind(Pearson=pow.Ps, CAT=pow.CAT,MIN2=pow.MIN)

}