二元变量的条件概率

Question

我有这套：

    1 2 3 4
r01 1 0 1 1
r02 1 1 0 1
r03 0 0 1 0
r04 0 0 1 0

代码：

initial <- data.frame(c(1,1,0,0), c(0,1,0,0), c(1,0,1,1),c(1,1,0,0), row.names = c("r01","r02","r03","r04"))
colnames(initial) <- c(1:4)

我需要获取条件概率，例如如果我选择r03 作为选择对象概率条件表，r01 将是：

          |r03=0 |r03=1
    r01=0 | 1/3  | 0/1
    r01=1 | 2/3  | 1/1

其他示例，如果我选择r03 并评估r02：

          |r03=0 |r03=1
    r02=0 | 0/3  | 1/1
    r02=1 | 3/3  | 0/1

我不知道如何实现这个表或者只是获取这个矩阵值。

Answer 1

我认为您正在寻找这个：

initial <- data.frame(c(1,1,0,0), c(0,1,0,0), c(1,0,1,1),c(1,1,0,0), row.names = c("r01","r02","r03","r04"))
colnames(initial) <- c(1:4)

tin <- data.frame(matrix(t(initial), ncol = 4, byrow = FALSE))
colnames(tin) <- rownames(initial)

tab <- table(tin$r01, tin$r03) #here specify values and get count
prop.table(tab,2)

你也可以加一些rownames我colnames

编辑： 固定概率并按列创建矩阵