2017-08-02 24 views
1

加权比例我已经建立使用3种状态(A,B,C)如下的转移概率矩阵:获取中的R

transition <-prop.table(with(data, table(data$old, 
       data$new)), 2) 

例如,如果你想获得对于A的概率 - > B,你会计算你看到B遵循A的次数,并将它除以你看到任何状态遵循A的次数。现在假设每一行数据都有一定的权重/重要性。我将如何修改上述以获得加权概率转换矩阵?

回答

0

你可以这样做......

transition <- prop.table(tapply(data$weight, list(data$old, data$new), sum), 2) 

其中data$weight为权数,以每行数据的列。

tapplylength将重现您拥有的。将其更改为sum将为每个组合添加权重,而不仅仅是对它们进行计数。