r语言多重插补

r语言多重插补

inputfile<-read.csv(“F:/rstudiospace/2019.3.7.csv”, header = TRUE)
View(inputfile)
is.na(inputfile) # 判断是否存在缺失
n <- sum(is.na(inputfile)) # 输出缺失值个数
n
par(mfrow = c(1, 2))
dotchart(inputfile$ET)

boxplot(inputfile$ET, horizontal = TRUE)

inputfile$ET[38] = NA
fix(inputfile)
library(lattice)
library(MASS)
library(nnet)
library(mice)
imp <- mice(inputfile, m = 4)
fit <- with(imp,lm(ET ~ time, data = inputfile))
pooled <- pool(fit)
summary(pooled)
result4 <- complete(imp, action = 3)

查看结果：
（1）imp i m p imp