1

背景 在我的论文中,我进行了一个成对的自助餐厅实验,将成对的水果呈现给鸟类,并记录他们选择的水果。从研究中,我发现分析这些数据的合适模型是 Bradley Terry 模型。我不是统计学家,对 R 相对较新,所以任何帮助都将不胜感激

使用包 BradleyTerry2 的问题

我的数据框看起来像这样,其中fruit1 和 Fruit2 代表玩家 1 和 2,以及 player1 vsplayer 2 的相应胜负频率

关注https://cran.r-project.org/web/packages/BradleyTerry2/vignettes/BradleyTerry.pdf

我正在使用代码

R>prefModel<-BTm(cbind(win1,win2),fruit1,fruit2,~fruit,id="fruit",data=df)

这会产生一个错误

Error in Diff(player1, player2, formula, id, data, separate.ability, refcat,  : 
'player1$fruit' and 'player2$fruit' must be factors with the same levels

尝试其他方法也会产生相同的错误

BTm(outcome=cbind(win1,win2),player1=fruit1,player2=fruit2,
formula=~fruit,id="fruit",)

(!is.factor(df$fruit1) || !is.factor(df$fruit2) ||
!identical(levels(df$fruit1), levels(df$fruit2)))

返回 [真]

4

0 回答 0