パワーを上げる友達

  • ほどほどの遠近関係の友達がいて、友達の多様度が高いとパワーが大きい

http://www.genome.med.kyoto-u.ac.jp/StatGenet/testRY20110208/df50.jpeg

Ks<-3
CheckDists<-seq(from=3,to=3,length.out=1)
nperm<-100

ks<-c(2,5,10,50)

#t1s<-acos(seq(from=0,to=1,length.out=20))
t1s<-seq(from=0,to=1,length.out=20)*pi/2

t2s<-seq(from=0,to=1,length.out=20)*pi/2
ppout2<-array(0,c(length(ks),length(t1s),length(t2s)))
for(i in 1:length(ks)){
	k<-ks[i]
	for(j in 1:length(t1s)){
		t1<-t1s[j]
		for(l in 1:length(t2s)){
			t2<-t2s[l]
			Y<-SimplexAroundTip(k,t1,t2)
			Kv<-Y[1,]
			spout<-SpherePower(k+1,Kv,Ks,CheckDists,Y,nperm)
			ppout2[i,j,l]<-c(spout$p.out)

		}
	}
}
plot(c(ppout2))
zlim<-range(ppout2)
filled.contour(ppout2[1,,],zlim=zlim)
filled.contour(ppout2[2,,],zlim=zlim)
filled.contour(ppout2[3,,],zlim=zlim)
filled.contour(ppout2[4,,],zlim=zlim)
  • 次元を上げると、コントラストが上る。また、「最強な友人」はより疎遠になる。疎遠とは言っても、半分よりは仲良しがよい?無限に幅広い友人関係のときには、「半分仲良し(仲がいいともわるいとも言えない、中立な関係)」がもっとも強力???
library(sphere)
Ks<-3
CheckDists<-seq(from=3,to=3,length.out=1)
nperm<-1000

ks<-2:50

#t1s<-acos(seq(from=0,to=1,length.out=20))
t1s<-seq(from=0,to=1,length.out=20)*pi/2

#t2s<-seq(from=0,to=1,length.out=20)*pi/2
ppout4<-array(0,c(length(ks),length(t1s)))
for(i in 1:length(ks)){
	k<-ks[i]
	for(j in 1:length(t1s)){
		t1<-t1s[j]
		#for(l in 1:length(t2s)){
			#t2<-t2s[l]
			#Y<-SimplexAroundTip(k,t1,t2)
			Y<-SimplexAtTip(k,t1)
			Kv<-Y[1,]
			spout<-SpherePower(k+1,Kv,Ks,CheckDists,Y,nperm)
			ppout4[i,j]<-c(spout$p.out)

		#}
	}
}
plot(c(ppout4))
matplot(ppout4,type="l")
matplot(t(ppout4),type="l")
image(t1s,ks,t(ppout4),xlab="r",ylab="df")

http://www.genome.med.kyoto-u.ac.jp/StatGenet/testRY20110208/dfvariationImage.jpeg