分からない問題はここに書いてね457

分からない問題はここ ..

48:１３２人目の素数さん
19/12/29 14:24:28.40 36JBVOYE.net
>>39
ありがとうございます。やはりCで書かれておりましたか。
以前にも教えていただいた方かと存じます。
これをRに移植するのは私の手には負えないので、Rで最初から作ってみました。
整数表示はできておりません。
ra2p <- function(x){ # 部屋割り配列からその確率を返す room allocation to probability
n=length(x)
y=numeric(3) # ルーレット後にA,B,Cに埋まった人数の数列
z=0 # ルーレット後に満室の数
p=numeric(n) # 収容可能な部屋から特定の部屋が選ばれる確率の配列
p[1]=1/3
for(i in 1:(n-1)){
xi=x[1:i] # i番目までの部分数列
y=c(sum(xi==1),sum(xi==2),sum(xi==3)) # A,B,cの収容人数
z=sum(y==n/3) # 満室の数　0 ～３
p[i+1]=c(1/3,1/2,1,1)[z+1] ;p # 1/(残り部屋種類数)　1/3, 1/2, 1のいずれか
}
prod(p)
}
# sum(apply(dat3,1,ra2p)) # 総和１を確認
srp <- function(x){ # a,bが同室の部屋割り配列の確率を総和する
a=x[1]
b=x[2]
sum(apply(dat3[dat3[,a]==dat3[,b],],1,ra2p))
}
p=combn(n,2,srp) # 組合わせ�

次ページ

続きを表示

1を表示