library(MASS)
n=100;p=10;per=0.1
X0=data=matrix(mvrnorm(n*p,0,1),n,p)
m=round(per*n*p,digits=0)
mr=sample(1:(n*p),m,replace=FALSE)
X0[mr]=NA;data0=X0
SR(data=data,data0=data0,real=FALSE,example=FALSE)
Run the code above in your browser using DataLab