ソースメモ

• 次元で丸くするかボンフェロニか

df1<-1
df2s<-1:5
pdf1<-10^(-9:-1)
pp<-matrix(0,length(df2s),length(pdf1))
for(i in 1:length(df2s)){
	df2<-df2s[i]

	probs<-pchisq(qchisq(pdf1,df1,lower.tail=FALSE),df2,lower.tail=FALSE)
	plot(probs)
	probs

	plot(probs/pdf1)
	pp[i,]<-probs/pdf1
}
matplot(t(pp),type="l")

• さらにただのメモ

Ns<-1000
ptype<-c(0,1)
p<-sample(ptype,Ns,replace=TRUE)

Nm<-5
gtype<-c(0,1,2)
g0<-sample(gtype,Ns,replace=TRUE)
g<-matrix(rep(g0,Nm),Ns,Nm)
rt<-0.05
for(i in 2:Nm){
	for(j in 1:Ns){
		if(runif(1)<rt){
			g[j,i]<-sample(gtype,1)
		}
	}
}
#g<-matrix(sample(gtype,Nm*Ns,replace=TRUE),nrow=Ns)

O<-matrix(c(p,(p-1)*(-1)),nrow=2,byrow=TRUE)
O

Sp<-RegularSpherize(O)

W<-matrix(0,Nm,length(O))
for(i in 1:Nm){
	W[i,]<-c(g[,i],rep(0,Ns))
}

W
df <- (length(O[1, ]) - 1) * (length(O[, 1]) - 1)
bothside<-TRUE
Tiis <- matrix(0, Nm, df)
for (i in 1:Nm) {
	Tiis[i, ] <- TestDf1(O, matrix(W[i, ], nrow = length(O[,1]), byrow = TRUE), Sp= Sp, simple = TRUE)$Tiist
}
if (bothside) {
	Tiis <- rbind(Tiis, -Tiis)
}
#Kv<-Tiis[1,]
Kv<-apply(Tiis[1:Nm,],2,sum)
Kv<-Kv/sqrt(sum(Kv^2))
Kv<-c(1,rep(0,df-1))
Ks<-1:6
CheckDists<-sqrt(qchisq(10^(-8),df=1,lower.tail=FALSE))
CheckDists<-c(5)
CheckDist<-sqrt(Ns)
Ks<-CheckDist
nperm<-10000

pout<-SpherePower(df=df,Kv=Kv,Ks=Ks,CheckDists=CheckDists,Tiis=Tiis,nperm=nperm)
plot(pout$Ks,pout$p.out)

image(1-cor(g))

gmt<-GeneralMatrixTest(O=O,Sp=Sp,W=W,power=TRUE)

df1<-1
df2s<-1:5
pdf1<-10^(-9:-1)
pp<-matrix(0,length(df2s),length(pdf1))
for(i in 1:length(df2s)){
	df2<-df2s[i]

	probs<-pchisq(qchisq(pdf1,df1,lower.tail=FALSE),df2,lower.tail=FALSE)
	plot(probs)
	probs

	plot(probs/pdf1)
	pp[i,]<-probs/pdf1
}
matplot(t(pp),type="l")