powered by
rZ_nrm_fc(Z, EZ, rho,s2, Y)
## The function is currently defined as rZ_nrm_fc<-function(Z,EZ,rho,s2,Y) { ZS<-simY_nrm(EZ,rho,s2) diag(ZS)<-rnorm(nrow(Y),diag(EZ),sqrt(s2)) Z[is.na(Y)]<-ZS[is.na(Y)] Z }
Run the code above in your browser using DataLab