特征值错误(Sigma,symmetric=TRUE):在';x';用zelig逻辑回归

特征值错误(Sigma,symmetric=TRUE):在';x';用zelig逻辑回归,r,logistic-regression,r-zelig,R,Logistic Regression,R Zelig,因此,我尝试进行倾向评分匹配,然后利用Zelig对匹配的数据进行逻辑回归 以下是我的数据的示例: 因此,我试图根据痴呆症、中风和年龄对患者进行匹配。然后,我试图看看,一旦患者配对,药物X是否会增加跌倒的风险。然后,我想对药物Y和药物Z进行同样的过程,看看每种药物是否都有关联 以下是我迄今为止用于执行匹配和逻辑回归分析的代码: data <- read_csv("Desktop/data.csv") View(data) attach(data) data[1:10,] m.out =

因此,我尝试进行倾向评分匹配,然后利用Zelig对匹配的数据进行逻辑回归

以下是我的数据的示例:

因此,我试图根据痴呆症、中风和年龄对患者进行匹配。然后,我试图看看,一旦患者配对,药物X是否会增加跌倒的风险。然后,我想对药物Y和药物Z进行同样的过程,看看每种药物是否都有关联

以下是我迄今为止用于执行匹配和逻辑回归分析的代码:

data <- read_csv("Desktop/data.csv")
View(data)
attach(data)
data[1:10,]

m.out = matchit(Fall ~ Dementia + Stroke + Age, method = "exact", data = data)

m.data = match.data(m.out)
library(Zelig)
z.out = zelig(Fall ~ Drug X + Dementia + Stroke + Age, model = "logit", data = m.data)

x.out0 <- setx(z.out, Drug X = 0)
x1.out0 <- setx(z.out, Drug X = 1)

s.out0 <- sim(z.out, x = x.out0, x1 = x1.out0)
知道我做错了什么吗?我已经在这里呆了一段时间,决定寻求帮助


提前谢谢

目前,Zelig已经有几年没有更新了,似乎有一些bug没有解决。一个问题是sim命令似乎对变量名和数据类型非常敏感。对于出现此错误的任何人,请尝试以下操作:

  • 变量
    druge X
    有一个Zelig可能无法处理的空间。尝试将其重命名为类似于
    drugh\u X

  • 如果有一个变量由于共线性而从回归模型中删除,则在转到
    sim()

  • 在运行模型之前,还要检查是否将任何字符变量转换为因子,否则可能会出现以下错误:

    对比误差
    Error in eigen(Sigma, symmetric = TRUE) : 
      infinite or missing values in 'x'