用R语言运行多重插补的代码,运行过程非常慢,几个小时过去了, 本以为出来1-5个数据集就会结束,刚刚发现又出现iter2 imp1新的插补数据集,懵圈中 1.请问老师多重插补这个代码是否有误呢?(我是网上搜索的),我的样本量是2046人,其实我主要想插补的缺失值只有年龄(缺3个),hcy(缺74个) 2.运行这么久时间合理嘛?
r.squared <- NULLif(!is.null(test.preds) && !is.null(test.labels)) {# store predicted points:test.df <- data.frame("Prediction"=test.preds,"Outcome"=test.labels,"DataSet"="test")# store residuals for predictions on the test datatest.residuals <-test.labels -test.predstest.res.df <...
等人考虑的一种机制,在Y和X2中缺少一些值: 根据Sullivan等人的说法,#make缺少一些数据 z1 - x1 / 0.2 ^ 0.5 r_y - 1 *(runif(n)expit(2.5 + 2 * z1)) r_x2 - 1 *(runif(n)expit(2.5 + 2 * z1)) obsData - fullData obsData $ y [r_y == 0] - NA obsData $ x2 [r_x2 ...
内容提示: R 语言用多重插补法估算相对风险数据分析报告 原文链接:http://tecdat.cn/?p=6379 在这里,我将用 R 中的一个小模拟示例进行说明。首先,我们使用 X1 和 X2 双变量法线和 Y 模拟大型数据集,其中 Y 遵循给定 X1 和 X2 的逻辑模型。 首先,我们模拟一个非常大的完整数据集: #simulate 完整数据 ...
R语言用线性模型进行臭氧预测: 加权泊松回归,普通最小二乘,加权负二项式模型,多重插补缺失值|附代码数据, 最近我们被客户要求撰写关于线性模型的研究报告,包括一些图形和统计输出。 在这篇文章中,我将从一个基本的线性模型开始,然后从那里尝试找到一个更