问题标签 [logistf]

我的数据遇到了一些问题，需要一些帮助。我正在尝试使用存在/不存在变量作为响应变量和几个解释变量（时间、位置、存在/不存在数据、丰度数据）运行 glm 分析。

首先，我尝试使用 glm() 函数，但是我有 2 个关于 glm.fit () 的警告：1：glm.fit：算法没有收敛 2：glm.fit：在经过一些调查后，出现数字 0 或 1 的拟合概率我发现问题很可能是准完全分离，因此决定使用 brglm 和/或 logistf。

• logistf：分析不运行运行 logistf() 时，我收到一条错误消息： chol.default(x) 中的错误：前导未成年人 39 不是肯定的 我在 Internet 上的理论和技术论文中查看了 logistf 包手册Heinze 和 Ploner 的，无法找到使用此功能的位置以及是否可以通过某些设置修复错误。

• brglm：分析运行但是我收到一条警告消息：在 fit.proc(x = X, y = Y, weights = weights, start = start, etastart # = etastart, : 达到迭代限制之前我找不到位置和原因此功能在运行包时使用，如果可以通过调整一些设置来修复它。

以更一般的方式，我想知道这些包的根本区别是什么。

我希望这足够有意义，如果这是我不知道的统计证据，我很抱歉。

这是我第一次提出问题，所以如果不应该这样，我深表歉意，并请您不要犹豫，让我知道。

谢谢您的帮助

霍奇特尔 C.

这是我的表格的摘录（由于表格太宽，我不得不截断行的长度：20 列）和我运行的不同公式：

I kept getting a problem for the following code; "weights=weight" was shown as an unused argument. How should I solve the problem?

Also, is there a way to perform the whole process shown above 30, 60, or 100 times and generate time (or count), x_0, x_1, and result for each time? Any suggestion would be great. Thanks.

我已经运行glm（用于逻辑回归），结果我得到算法发散。我想试试这个功能logistf来检查 Firth 的方法是否可以改进我的数据分析。让dat表示我的数据框并y ~ f1 + f2成为我的模型。试

或者

我总是得到错误

由于glm运行顺利，我看不出我的数据结构如何成为错误的根源（这就是我没有dat明确描述的原因）。你看到任何潜在的冲突吗？我非常感谢你，阿维图斯

我无法弄清楚如何从logistf()回归模型中提取标准错误“sd(coef)”信息。这些模型属于 logistf 类，并且手册指出可以通过这种方式提取数据：

以下通用方法可用于logistf 的输出对象：print、summary、coef、vcov、confint、anova、extractAIC、add1、drop1、profile、terms、nobs。

但是，标准错误不存在。在 str(summary(fit)) 中没有 se(coef) 的对象，我查看了源代码没有运气。

任何帮助，将不胜感激！

logistf::logistf我正在使用and函数编写一个logistf::forward函数。我将在这里给出一个使用包中的最小工作sex2示例logistf。

我得到以下输出：

步骤 0：启动模型

is.data.frame(data) 中的错误：找不到对象“datamod”。

有人对此有解释吗？

我已经使用 logistf 包为 R 中的惩罚逻辑模型运行了一些模型。但是，我希望为数据绘制一些森林图。

sjPlot 包： http: //www.strengejacke.de/sjPlot/custplot/

为 glm 输出提供了出色的功能，但没有为 logistf 功能提供功能。

有什么帮助吗？

我正在研究一个具有匹配对（1 到 5/6 匹配）的数据库，并尝试使用条件逻辑回归获得我的变量的优势比（95% 置信区间）。

此后，我发现我的两个预测变量具有准完全分离，因此我决定尝试惩罚最大似然估计方法（logistf）。

然而，出现了两个问题。首先，优势比（从logistf命令OR = 0.9（strata()考虑到我的条件逻辑回归中的匹配对）到logistf命令，所以我怀疑OR = 0.9获得的结果不正确，因为它没有将我的案例与其各自的 5/6 对照进行比较。

是否有对 OR 获得的差异的解释，有没有办法让我将strata()函数包含到我的logistf命令中？非常感谢您！

PS：我是 R 新手，所以如果问题没有明确表达或者您需要进一步澄清，请告诉我。

我正在使用logistfR 中的包使用 Firth 的校正对逻辑回归进行建模。当我使用命令DROP1/ADD1时，R 帮助说它基于惩罚似然比。但它所说的惩罚实际上是 Firth 的更正，还是只是一种惩罚的可能性？

在此先感谢您的帮助。

我在 R 上运行一个logistf模型

由惩罚性 ML 置信区间拟合的模型和由轮廓似然性拟合的 p 值 轮廓似然性轮廓似然性轮廓似然性轮廓似然性轮廓似然性轮廓似然性

协方差矩阵：

我一直在尝试提取 Chisq 值，但我不知道该怎么做！我可以轻松提取 coef、术语、p.values、上下 ci 等，但我无法找到提取 Chisq 值的方法

有什么建议么？

提前致谢

我正在尝试在数据框（dim：11359x139）上应用logistf（）。所有变量都是二进制的。我收到以下消息：“logistf.fit 中的错误（x = x, y = y, weight = weight, offset = offset, firth, : no memory available”。即使我只考虑 20 行和 139 个预测变量我得到相同的数据框是硬件问题还是我的错？