その3の続き。 L<- 1:500 # 尤度比 a<- seq(from=0,to=1,by=0.01) # 事前確率 f<- function(L,a){ return((L*a)/(L*a+1-a)) } theta<- 40 phi<- 35 library(rgl) m<- expand.grid(L,a) M<- matrix(f(m$Var1,m$Var2),nr=length(L),nc=length(a)) persp(M,col=…

その1とかその2とか条件付き確率とかそれについてプロットしたものについて続き。 A：擬父が子の真の父であるという事象 B：親子鑑定で擬父が真の父だという結果 親子鑑定結果が出たとき、擬父が本当に真の父か、という確率が問題となる。 ベイズの定理から …