Question

我试图使用R中的BradleyTerry2软件包输出布拉德利特里模型的能力估计值。我不断得到一个非常神秘的错误。他们的文档工作中的一些示例和其他示例返回了我使用我的数据时出现的相同错误。此代码使用文档中的示例分析之一。因此，如果你加载图书馆＆＃34; chameleon＆＃34;数据应该已经存在

install.packages("BradleyTerry2")
library (BradleyTerry2)
summary(chameleon.model <- BTm(player1 = winner, player2 = loser,formula = ~ prev.wins.2 + ch.res[ID] + prop.main[ID] + (1|ID), id = "ID",data = chameleons))
BTabilities(chameleon.model)

我得到的错误是

Error in X[, est, drop = FALSE] : (subscript) logical subscript too long

任何人都知道怎么做？

Answer 1

我维持BradleyTerry2。如果您在能力公式中具有特定于比赛的预测变量，则会出现此错误。它应该像?BTabilities：

中记录的那样工作

......能力是根据拟合的条款计算的仅涉及玩家协变量的模型（由模型公式中的'model $ id'）。因此参数在任何其他术语假定为零。

我们不知道这不起作用，所以感谢错误报告。在修复之前，您可以直接计算出能力和标准错误：

## get names of relevant coefficients
> nm <- grep("[ID]", names(coef(chameleon.model)),
>            fixed = TRUE, value = TRUE)
> nm
[1] "ch.res[ID]"    "prop.main[ID]"

## get names of corresponding predictors
> IDvar <- gsub("[ID]", "", nm, fixed = TRUE)
> IDvar
[1] "ch.res"    "prop.main"

## create coefficient vector and design matrix
> cf <- coef(chameleon.model)[nm]
> X <- as.matrix(chameleons$predictors[, IDvar])

## compute abilities
> abilities <- X %*% cf
> colnames(abilities) <- "abilities"

## compute standard errors
> V <- vcov(chameleon.model)[nm, nm]
> res <- cbind(abilities = abilities,
>               se = sqrt(diag(X %*% V %*% t(X))))
> head(res)
    abilities       se
C01  3.757229 1.655205
C02  2.404778 1.017782
C03  2.319346 1.133959
C04  1.892671 1.399391
C05  2.253472 1.101628
C06  2.015840 1.075806

这将为预测变量中具有缺失值的个体提供NA，但是他们的能力将单独建模并在模型摘要中返回。上面的代码也假定了连续的协变量，如果你的某些预测因子是因素，你需要更正式地创建一个模型矩阵，例如。

X <- model.matrix(reformulate(IDvar, intercept = FALSE),
                  c(chameleons$predictors))

计算BradleyTerry能力排名

1 个答案: