目次はこちら こちらやこちらも参考 独立な分割表の確率行列のランクは１ Nr<-sample(2:10,1) Nc<-sample(2:10,1) library(MCMCpack) Pr<-rdirichlet(1,rep(1,Nr)) Pc<-rdirichlet(1,rep(1,Nc)) Pmat<-t(Pr) %*% Pc Pmat qr(Pmat)$rank m次元分割表の正確確…

n次元分割表でn軸の個のカテゴリ別に識別する、長さのベクトルを並べた行列をつくる > A [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 1 0 1 0 1 0 1 0 1 0 1 0 [2,] 0 1 0 1 0 1 0 1 0 1 0 1 [3,] 1 1 0 0 0 0 1 1 0 0 0 0 [4,] 0…