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)