利用R求p(x|y,z)？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我有一个三维矩阵XYZ
可以求得：
p(x)   <- rowSums(XYZ, dims = 1L)
p(y)   <- apply(XYZ[,,], 2, sum)
p(z)   <- colSums(XYZ, dims = 2L)

p(xz)  <- apply(XYZ, c(1,3), sum)
p(yz)  <- colSums(XYZ, dims = 1L)


请问大侠，怎么求p(x|y,z)的值，谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

#p_x_yz代表p(x|yz), p_xyz 代表p(x,y,z), p_yz代表p(yz)
#px_yz = p_xyz / p_yz

fun_x_yz <- function(x, tmp){
  x/tmp
}
p_x_yz <- apply(p_xyz, 1, fun_x_yz, p_yz)

p_x_yz <- array(t(p_x_yz), dim = dim(p_xyz))
p_x_yz[is.na(p_x_yz)] = 0
我自己的答案，如有不对，请指教，谢谢！