包Optmatch错误消息

我试图使用对功能在ř包optmatch上一些示例数据：包Optmatch错误消息

> tmp 
    id dose  Pred 
1 1 2 0.08494142 
2 1 2 -0.02784262 
3 1 2 -0.05284574 
4 1 2 -0.07428686 
5 1 3 0.04335151

可使用此代码重新创建：

tmp <- data.frame(id = rep(1,5), dose = c(rep(2,4),3), 
        Pred = c(0.08494142, -0.02784262, 
          -.05284574, -0.07428686, 
          0.04335151))

当我运行代码：

pair(x = tmp$Pred, z = tmp$dose, data=tmp, remove.unmatchables = TRUE)

我得到以下错误：

Error in toZ(z) : Treatment indicator must have exactly 2 levels not 1

然而，如可以从数据帧中可以看出，处理指示变量剂量仅具有两个等级：2或3

使用功能getAnywhere我能够打印出功能TOZ（其由一对称为），其推测生成错误的部分：

> getAnywhere(toZ) 
[Lines omitted] 
function (x) 
{ 
[Lines of code omitted] 
if (length(unique(x)) != 2) { 
     stop(paste("Treatment indicator must have exactly 2 levels not", 
      length(unique(x)))) 
}

为了证实我的怀疑，我运行以下命令：

> length(unique(tmp$dose)) 
[1] 2

确认处理指标确实具有完全相同2级。我很困惑。任何帮助将不胜感激。

来源

2017-02-21 Colin

怀疑错误来自'remove.unmatchables = TRUE'。你的单身人士，治疗观察正在被放弃。 – lmo

这是optmatch中的bug。我们有一个修复，将在下一个发布的版本。目前，解决方法是从数字处理指标中构造一个逻辑变量。例如。

tmp$dose2 <- tmp$dose == 3

来源

2017-02-21 18:28:08

包Optmatch错误消息

回答

相关问题