我需要对我的小样本量(n = 46
)进行精确logistic回归。我有两个因变量(as
和fof
),我将为每个因变量建立两个单独的模型。我有17个自变量,例如:连续(minonground
、mental
、...)、分类(pdd
分为三组)和一些二进制(injury
、multiple
、...)。当我为数据集编写elrm
时,遇到了两个问题:
1.我收到一个无法识别n
变量的错误。我不知道互联网上所有示例中存在的二项式试验的数量是多少。我只有一行每个患者。
1.当结局变量为fof
且所有其他变量均为预测变量时,"兴趣"意味着什么?
我的代码:
library(elrm)
fit1 <- elrm(fof/n~ minonground+injury+pdd, interest = ~minonground, r=4, iter = 2000, burnIn = 100,dataset = com, alpha = 0.05)
我尝试用数字代替n
,但是失败了。事实上,我不知道它是什么,也不知道应该如何为我的数据集定义它。
1条答案
按热度按时间brccelvz1#
精确逻辑回归的任务是了解哪个预测因子对响应有影响,尤其是对于小样本。因此,在您的情况下,您希望了解
minonground
是否显著,是否考虑了其他变量的影响。因此,您将interest = ~minonground
指定为感兴趣的变量,但公式fof / n ~ minonground + injury + pdd
中的所有其他变量构成模型,您对此不感兴趣。它们是有害参数。在公式中,
x/y
表示success/total
,因此在您的数据集中,您应该有fof
列(成功数量)和n
列(病例总数)(n
应该是列,而不是数字)。如果您有一行用于每名患者并观察到一些病例,则在fof
中应该是二进制,n = 1
用于每行。因此,我根据您提到的结构生成了一些随机样本+添加了
n
列,您可以看到minoground
不显著(p > 0.8
,置信区间包括零,这不足为奇-它只是一个随机样本)。见下文:
输出:
电话: