facedata<-read.table('facedata.txt') names(facedata)<-c('B10%', 'B20%','B40%','C10%','C20%','C40%','P10%','P20%','P40%') stars(facedata,ncol=6, lwd=2, key.loc=c(4,4)) stars(facedata,ncol=6, lwd=2, key.loc=c(8,-1), axes=T, labels=labels) lines(c(2,6), c(9,9)) text(4.1,9.2,'Without Auxiliary variables') lines(c(8.5,13), c(9,9)) text(10.6,9.2,'With Auxiliary variables') text(2,7.5,'MI') labels<-c('MCAR/MI','MCAR/FIML','MCAR/EM','MCAR/MI','MCAR/FIML','MCAR/EM','MAR/MI','MAR/FIML','MAR/EM','MAR/MI','MAR/FIML','MAR/EM','MNAR/MI','MNAR/FIML','MNAR/EM','MNAR/MI','MNAR/FIML','MNAR/EM') stars(facedata,ncol=6, lwd=2, key.loc=c(6.3,-1), xlim=c(-2,15)) stars(facedata,ncol=6, lwd=2, key.loc=c(4.3,-1), xlim=c(-2,15),axes=T, ylim=c(-3, 10))->mnar lines(c(2.1,6.3), c(8,8)) text(4.2,8.2,'Without Auxiliary variables') lines(c(8.4,12.5), c(8,8)) text(10.5,8.2,'With Auxiliary variables') text(2.1, 7.6, 'MI') text(4.2, 7.6, 'FIML') text(6.3, 7.6, 'EM') text(8.4, 7.6, 'MI') text(10.5, 7.6, 'FIML') text(12.6, 7.6, 'EM') ## add row names text(0,6.3,'MCAR') text(0,4.2,'MAR') text(0,2.1,'MNAR') ## add legends text(7, 0.2, 'B: Bias', adj=0) text(7, -.3, 'C: Coverage probability', adj=0) text(7, -.7, 'P: Power', adj=0)