R中成对偏好数据的Bradley Terry模型

Question

背景 在我的论文中，我进行了成对自助餐厅实验，将成对的水果呈现给鸟类，并记录它们选择的鸟类。从研究中，我发现分析此数据的合适模型是Bradley Terry模型。我不是统计学家，而且对R还是比较陌生，因此将不胜感激

问题 使用软件包BradleyTerry2

my data frame looks like this in which fruit1 and Fruit2 represent Player 1 and 2 and the coresponding frequencys of wins and looses for 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)))

返回[TRUE]