@ND - lab:projects:12simulation_design_cautions
https://nd.psychstat.org/
2026-05-19T11:47:54+00:00@ND
https://nd.psychstat.org/
https://nd.psychstat.org/_media/wiki/dokuwiki.svgtext/html2016-01-24T14:48:26+00:00Anonymous (anonymous@undisclosed.example.com)index
https://nd.psychstat.org/lab/projects/12simulation_design_cautions/index
* Simple regression analysis with missing datatext/html2016-01-24T14:48:26+00:00Anonymous (anonymous@undisclosed.example.com)simple_regression_analysis_with_missing_data
https://nd.psychstat.org/lab/projects/12simulation_design_cautions/simple_regression_analysis_with_missing_data
B<-100
N<-1000
R<-0.3
runsim<-function(B, N, R){
## (1) Simulation based on different data sets
b.mcar.all<-NULL
b.mar.all<-NULL
for (i in 1:B){
x<-rnorm(N)
y<-1+2*x +rnorm(N)
## estimation based on MLE
## MAR data
y1<-y
c.value<-x[rank(x)==N*R]
y1[x < c.value]<-NA
b.mar<-coef(lm(y1~x))
b.mar.all<-rbind(b.mar.all, b.mar)
##MCAR data
x<-rnorm(N)
y<-1+2*x +rnorm(N)
## estimation based on MLE
#MCAR data
…