已有的气象缺测数据中,多采用相关分析法和简单的线性插补法对缺测的气象数据进行插补,可能误差较大,本文本文主要采样missForest的R包来对缺测的气象数据进行填充,该包也可用于其他缺测数据的填充,但前提是缺测数据不可过多。
具体如下
library(missForest)
data<-read.csv('d:/missexample.csv',header = TRUE)
data.imp <- missForest(data)
data.imp$OOBerror
a1<-data.imp$ximp
write.csv(a1,file ='d:/result.csv',row.names = F )
最后的输出结果在d盘下的result部分,后续会更新下批量的循环插补,以及该方法和相关法、线性插值法的对比分析
更多需求,请查看个人介绍