logistf

    1热度

    1回答

    我正在尝试data.frame逻辑回归(11359行,137列)。 data.frame包含Y(一个因变量)和预测变量(136个独立变量)。所有变量都是二进制的。 我创建基于 “my_data” data.frame公式是f = as.formula(paste('y ~', paste(colnames(my_data)[c(3:52, 54:133, 138:143)], collapse =

    3热度

    1回答

    我遇到了一些与我的数据有关的问题,需要一些帮助。 我试图运行有/无变量作为响应变量和几个解释变量(时间,地点,存在/缺勤数据,丰度数据)的glm分析。 首先,我试图使用GLM()函数,但是我有关于glm.fit()2个警告: 1:glm.fit:算法没有收敛 2:glm.fit:嵌合概率数值0或发生了1次 经过一番调查后,我发现问题很可能是准完全分离,因此决定使用brglm和/或logistf。

    0热度

    2回答

    我一直对下面的代码有问题; “权重=重量”显示为未使用的参数。我应该如何解决这个问题? x_0 <- rbinom(1,100, 0.01) x_1 <- rbinom(1,100, 0.1) x <- c(0,0,1,1) y <- c(0,1,0,1) weight <- c(100-x_0, x_0, 100-x_1, x_1) result <- logistf(y ~ x,

    1热度

    1回答

    我无法弄清楚如何从logistf()回归模型中提取标准错误“sd(coef)”信息。这些模型类logistf和数据可以提取这样的手册指出的: 以下通用方法可用于logistf`s输出对象:打印,汇总,COEF,vcov,confint,方差分析,extractAIC ,add1,drop1,个人资料,条款,nobs。 但是,标准错误并不存在。se(coef)的str(summary(fit))中没