r - 用计算值替换缺失值
问题描述
我正在尝试学习如何用计算值替换一个变量中的缺失数据。
我的数据集 ( bk3
) 看起来像:
ign: 80, 96, 75, 66, 53
Mean: 26, 24, 27, 34, 41
sd: 6, 7, NA, 8, 4
lci: 24, 25, 20, 32, 38
uci: 29, 26, 29, 33, 43
输入:
bk3 <- structure(list(ign = c(80L, 96L, 75L, 66L, 53L), mean = c(26L, 24L,
27L, 34L, 41L), sd = c(6L, 7L, NA, 8L, 4L), lci = c(24L, 25L, 20L,
32L, 38L), uci = c(29L, 26L, 29L, 33L, 43L)), .Names = c("ign",
"mean", "sd", "lci", "uci"), class = "data.frame", row.names = c(NA, -5L))
基本上,我使用 95% 置信区间 ( uci
, lci
) 和样本 n 的 ( ign
) 来计算缺失的 SD ( sd
)。
我尝试使用的代码是:
bk3$sd[is.na(bk3$sd)] <- (bk3$uci - bk3$lci) * sqrt(bk3$ign)/3.92
但我收到以下警告消息:
“要替换的项目数不是替换长度的倍数”
更新:我正在尝试创建一个函数,如果提供了适当的变量,它将自动执行此操作。我尝试将其设置为以下格式:
fillsd <- function(x, n, u, l)
{
i1 <- is.na(x)
i2 <- n > 59
x[i1 & i2] <- with(df, (u[i1 & i2] - l[i1 & i2]) * (sqrt(n[i1 &
i2])/3.92)) }
虽然函数“fillsd”似乎已在我的全局环境中正确保存,但当我尝试将它与以下代码一起使用时它不起作用:
fillsd(x="bk3$sd", n="bk3$ign", u="bk3$uci", l="bk3$lci")
该代码没有产生错误消息,但该函数似乎也没有做任何事情。这是我处理的第一个函数,我无法找到可比较的示例来知道代码的哪一部分是不正确的。如果您对如何完成这项工作有任何想法,请告诉我。谢谢!
解决方案
如果我们replace
将'sd'的NA元素与其他列的计算值对应的元素,那么逻辑索引应该在两边赋值。根据计算的性质,它给出的长度等于数据集的行数,而 lhs 只有较小的长度,因为我们只对具有 NA 元素的行进行子集化,这会导致长度不等,从而导致错误
i1 <- is.na(bk3$sd)
bk3$sd[i1] <- with(bk3, (uci[i1] - lci[i1]) * sqrt(ign[i1])/3.92)
但是,如果我们决定基于获取某些列的mean
of来获取摘要sum
,则它是一个数字,并且在 rhs 上没有逻辑索引是有意义的,因为值会被回收
数据
bk3 <- structure(list(ign = c(80, 96, 75, 66, 53), Mean = c(26, 24,
27, 34, 41), sd = c(6, 7, NA, 8, 4), lci = c(24, 25, 20, 32,
38), uci = c(29, 26, 29, 33, 43)), .Names = c("ign", "Mean",
"sd", "lci", "uci"), row.names = c(NA, -5L), class = "data.frame")
推荐阅读
- ruby - 扫描/匹配不正确的输入错误消息
- matlab - 如何在 MATLAB 中使用 fzero 获得给定方程的零点?
- javascript - 数据表中的列到行
- android - 无法解决:support-core-ui
- angular - 带参数 Angular 5 的数据绑定
- python - (Python Tkinter) 点击一个按钮,然后依次运行不同的函数,并在我的 GUI 中实时显示结果
- redis - Redis 扫描中的“或”匹配模式
- swift - 在 Swift 中选择时如何将子视图置于前面?
- html - 使浮动 div 填充空白
- javascript - REACT 16.3 - 将 props 数组渲染为地图(可选)