ただのメモ

• 分割表の部分分割表制約の部分を取り出すと、その制約の具合がパターンになる
• このパターンから基底を取り出すとマルコフ基底ができる
• 使用する関数はこちらの記事で

Nv<-5
Rs<-sample(2:4,Nv,replace=TRUE)
Vs<-1:Nv
ns<-Rs
# Facesをランダムに作る
Nf<-3

Faces<-list()
maxN<-3
Ns<-sample(1:maxN,Nf,replace=TRUE)
for(i in 1:Nf){
	tmpVs<-sample(Vs,Ns[i])
	Faces[[i]]<-as.set(tmpVs)
}
# Facetsにする
Facets<-MakeFacets(Faces)
# グラフにしてみる
lenFacets[ii]<-length(Facets)
g<-GraphFacets(Nv,Facets,plot=FALSE)
##############
#Facets<-set(1:Nv)
# 観測表をランダムに作る
Obs<-array(runif(prod(ns)),ns)
Obs<-Obs/sum(Obs)

# 期待値表を作る
#Etable<-MakeExpected2(Obs,ns,Facets,plot=TRUE)
#Etable<-MakeExpected3(Obs,ns,Facets,plot=TRUE)
m.e.out<-MakeExpected4(Obs,ns,Facets)
Etable<-m.e.out$Etable
X<-m.e.out$X
Z<-m.e.out$Z
Xzero<-X[which(Z==0),]

corXzero<-cor(Xzero)

Thetas<-acos(corXzero)/pi

Thetas

image(Thetas)