w<-t(matrix(b[1:6],2,3)) rho<-b[7] X<-cbind(1,seq(1, 20, 1),min(birthY)) U<-cbind(X%*%w,0) ival<-log(exp(U[,1]/rho)+exp(U[,2]/rho)) #クラスター選択確率 CL1<-exp(rho*ival)/(exp(rho*ival)+exp(U[,3])) CL2<-exp(U[,3])/(exp(rho*ival)+exp(U[,3])…

とある飲み会で人名ビンゴをやることになった。飲み会の人数は12人なので、参加者は12人の中から9人選び、3×3のマスの中に名前を記入していく。12人のうち2人は人名ビンゴの企画者なので、ビンゴには参加しない。 （つまり参加者は10人）何回か抽選を行い、…