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

[過去ﾛｸﾞ] 分からない問題はここに書いてね457 (1002ﾚｽ)
上下前次1-新
通常表示 512ﾊﾞｲﾄ分割ﾚｽ栞
抽出解除ﾚｽ栞

このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索歴削→次ｽﾚ栞削→次ｽﾚ過去ﾛｸﾞﾒﾆｭｰ

46: １３２人目の素数さん [sage] 2019/12/29(日) 14:24:28 ID:36JBVOYE

>>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) # 組合わせごとに実行
person=combn(6,2)
data.frame(person[1,],person[2,],同室確率=p)

結果は一致したようです。
person.1... person.2...     同室確率
1            1           2 0.3333333333
2            1           3 0.2222222222
3            1           4 0.1851851852
4            1           5 0.1296296296
5            1           6 0.1296296296
6            2           3 0.2222222222
7            2           4 0.1851851852
8            2           5 0.1296296296
9            2           6 0.1296296296
10           3           4 0.2407407407
11           3           5 0.1574074074
12           3           6 0.1574074074
13           4           5 0.1944444444
14           4           6 0.1944444444
15           5           6 0.3888888889

http://rio2016.5ch.net/test/read.cgi/math/1577457155/46

上下前次1-新書関写板覧索設栞歴

ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ AAｻﾑﾈｲﾙ

ぬこの手ぬこTOP 0.030s