原始数据,作图目的缺失值的可视化
set.seed(123456)
age<-round(abs(rnorm(200, mean = 67, sd = 19)))
set.seed(12345)
sex<-rbinom(200, 1, 0.45)
set.seed(12356)
sex.miss.tag<-rbinom(200, 1, 0.3) #MCAR
sex[sex.miss.tag==1]<-NA
sex[sex==1]<-"male"
sex[sex==0]<-"female"
set.seed(12456)
lac<-round(abs(rnorm(200, mean = 3, sd = 4)),1)
set.seed(13456)
lac.miss.tag<-rbinom(200, 1, 0.3)
lac[lac<=3&lac.miss.tag==1]<-NA # NMAR
set.seed(23456)
wbc<-round(abs(rnorm(200, mean = 12, sd = 4)),1)
set.seed(123)
wbc.miss.tag<-rbinom(200, 1, 0.3)
wbc[wbc.miss.tag==1]<-NA
set.seed(1234)
crp<-round(abs(rnorm(200, mean = 50, sd = 100)),1)
set.seed(3456)
crp.miss.tag<-rbinom(200, 1, 0.4)
crp[wbc<=12&crp.miss.tag==1]<-NA # MAR
data<-data.frame(age,sex,lac,wbc,crp)
画图代码
library(VIM)
matrixplot(data,interactive=T)
R gui中可以根据每列的数据大小排序互动作画,Rstudio只能看到最终的作图,是我的设置问题吗