这可能是一个简单的问题,但我就是无法解决。我有一个生物化学测试结果的数据框。其中一些测试由于检测限制base_crp而返回值。<3在继续之前我需要估算这些数据。我想正确地做到这一点,所以不仅仅是替换。
我尝试了zCompositions包中的 multLN ,但似乎认为所有<3值都是负数(错误提示X contains negative values)。似乎也没有太多文档——这是一个不起眼的包吗?
我还查看了LODI,但它希望我指定插补模型的协变量 - 有没有正确的方法来选择这些变量?不管怎样,我选择了 3 个理论上关联良好的,并使用了以下代码:
clmi.out <- clmi(formula = log(base_crp) ~ base_wcc + base_neut + base_lymph, df = all, lod = crplim, seed = 12345, n.imps = 5)
Run Code Online (Sandbox Code Playgroud)
其中 base_crp 是我要修复的变量。我用 替换了所有 <3NA并插入了一个新列all$crplim <- "3"。然而,这才刚刚回归
Error in sprintf("%s must be numeric.") : too few arguments。
即使我能让 LODI 工作,我也不确定它是否是正确的工具。我只是一名本科生,几乎没有统计背景,所以我不太明白我在做什么——我只是想要一些用数字填充列的东西,这样我就可以继续使用皮尔逊相关性和线性回归等。我真的很感激一些帮助。提前致谢。