矩阵中NA值需要将其替换为它所在列的平均值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

写完以后的结果总是出问题，算出平均值了，但是它会出问题，比如第3行第2列有个NA，他计算好第2列的平均值后，将其写到第3行第2列的NA了，但是把第3行第4,5列正常的数值也修改为刚刚求出的平均值了。。。不知道为什么啊。函数哪里写错了嘛？
for(i in 1:5){
a=mean(na.omit(data[,i]))
data[is.na(data[,i])]=a
}

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：平均值 Data mean 不知道 MIT 平均值

1. data[,is.na(data[,i])]=a

复制代码

少个, ,你那是提取的行，当然会出错

jiangbeilu 发表于 2015-10-28 23:56
少个, ,你那是提取的行，当然会出错

照你这样改之后报错了。
Error in `[<-`(`*tmp*`, , is.na(data[, i]), value = 5.84125874125874) :
  (subscript) logical subscript too long

for(i in 1:5){
  a=mean(na.omit(data[i,]))
  data[i,is.na(data[i,])]=a
}

creaoplus 发表于 2015-10-29 04:29
for(i in 1:5){
  a=mean(na.omit(data))
  data)]=a

你这个就完全错了额。。我要的是把NA替换成为当列的平均值哦

1. for(i in 1:5) {
   
2. a <- mean(data[, i], na.rm = T)
   
3. data[is.na(data[, i]), i] <- a
   
4. }

复制代码

上面是顺着你的思路做的改进，在我的电脑上已经通过。
你原来代码中的问题存在于data[is.na(data[,i])]=a这一句，你可以调试一下，is.na(data[,i])返回的是一个（FALSE FALSE TRUE FALSE FALSE）之类的向量，被data[is.na(data[,i])]引用后，会自动循环补齐ncol(data)次，变成选定了第i列有空值的那一整行，并进行赋值，于是就出现了你反映的问题，变成data[is.na(data[, i]), i] 就能特指第i列，就能解决了。

这个错误再次反映了——矩阵本质就是向量。如果用向量化的思维去解决这个问题，可以不用循环就能达成：（假设原数据是矩阵，且每列有且只有一个NA值）

1. data[is.na(data)] <- colMeans(data, na.rm = T)

复制代码

cheetahfly 发表于 2015-10-29 08:12
上面是顺着你的思路做的改进，在我的电脑上已经通过。
你原来代码中的问题存在于data)]=a这一句，你可以调 ...

太有启发了，学到了！

xuxinzju 发表于 2015-10-29 05:15
你这个就完全错了额。。我要的是把NA替换成为当列的平均值哦

对不起 把列看成行了。

假设数据框是dat，那么以下方法可以避免循环，虽然结果和楼上的是一样的
a<- colMeans(dat,na.rm = TRUE)
dat[is.na(dat)]<- a[ceiling(which(is.na(dat))/nrow(dat))]