面に色を付けてわかりやすくする

k<-3
t1<-acos(0.9)
t2<-acos(0.2)
CategoryVectorFull<-function(k){
	tmp<-CategoryVector(k)
	tmp2<-NULL
	n<-length(tmp[,1])
	for(i in 1:(n-1)){
		for(j in (i+1):n){
			x<-(tmp[i,]-tmp[j,])
			tmp2<-rbind(tmp2,c(x),c(-x))
		}
	}
	tmp2<-tmp2/sqrt(sum(tmp2[1,]*tmp2[1,]))
	tmp2
}

SimplexAroundTipFull<-function(k,t1){
	R2<-sin(t1)
	X<-CategoryVectorFull(k+1)*R2
	R3<-cos(t1)
	Y<-cbind(rep(R3,k+1),X)
	Y
}
t1<-pi/3

Y<-SimplexAroundTipFull(k-1,t1)
#Y<-SimplexAtTip(k-1,t1)
Y<-rbind(c(1,rep(0,k-1)),Y)
Y<-rbind(Y,-Y)

#Y<-SimplexAtTip(k,t1)
ee<-EularAngleGrid(3,200)
plot3d(ee,type="l")
Tiis<-RandomSphere(3,1,4)
Tiis<-rbind(Tiis,-Tiis)

cc<-ee%*%t(Y)
Mcc<-apply(cc,1,max)
Mid<-rep(0,length(cc[,1]))
for(i in 1:length(Mid)){
	Mid[i]<-which(cc[i,]==max(cc[i,]))
}
ee2<-ee/Mcc
xlim<-range(ee2)

plot3d(ee2,type="l",xlim=xlim,ylim=ylim,zlim=zlim,col=Mid)