背景 在我的论文中,我进行了一个成对的自助餐厅实验,将成对的水果呈现给鸟类,并记录他们选择的水果。从研究中,我发现分析这些数据的合适模型是 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)))
返回 [真]