Skip to content

Instantly share code, notes, and snippets.

@ericpgreen

ericpgreen/stata_mi_example.R

Created April 8, 2020 20:03
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save ericpgreen/a1a640184db13a070753288782f76d06 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save ericpgreen/a1a640184db13a070753288782f76d06 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
stata_mi_example.R
df_mi <- structure(list(v1 = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L,
24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L,
37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L,
50L, 51L, 52L, 53L, 54L, 55L, 56L, 57L, 58L, 59L, 60L, 61L, 62L,
63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L,
76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 88L,
89L, 90L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L,
41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 91L, 92L, 93L,
94L, 95L, 96L, 97L, 98L, 99L, 100L, 41L, 42L, 43L, 44L, 45L,
46L, 47L, 48L, 49L, 50L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L,
99L, 100L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L,
91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 41L, 42L,
43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 91L, 92L, 93L, 94L, 95L,
96L, 97L, 98L, 99L, 100L, 41L, 42L, 43L, 44L, 45L, 46L, 47L,
48L, 49L, 50L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L
), id = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L,
13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L,
26L, 27L, 28L, 29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L,
39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L,
52L, 53L, 54L, 55L, 56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L,
65L, 66L, 67L, 68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L,
78L, 79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 88L, 89L, 90L,
91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 41L, 42L,
43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 91L, 92L, 93L, 94L, 95L,
96L, 97L, 98L, 99L, 100L, 41L, 42L, 43L, 44L, 45L, 46L, 47L,
48L, 49L, 50L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L,
41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 91L, 92L, 93L,
94L, 95L, 96L, 97L, 98L, 99L, 100L, 41L, 42L, 43L, 44L, 45L,
46L, 47L, 48L, 49L, 50L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L,
99L, 100L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L,
91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L), trt = c(0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L), x1 = c(0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L,
0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L,
1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L,
1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L,
0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L,
1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L,
0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L,
1L), y = c(1, 3, 1, 2, 1, 3, 3, 2, 2, 3, 3, 1, 1, 1, 2, 2, 2,
2, 3, 1, 5, 3, 3, 3, 3, 4, 3, 3, 4, 4, 4, 3, 6, 5, 3, 6, 5, 4,
4, 4, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 7, 8, 7, 6, 8,
7, 6, 6, 8, 7, 7, 8, 8, 7, 7, 7, 7, 6, 7, 7, 9, 9, 8, 10, 10,
9, 10, 10, 9, 10, 10, 8, 8, 8, 8, 10, 9, 10, 8, 8, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, 2.30907, 6.40722, 5.57695, -1.40922,
2.20291, 5.12494, 9.07028, 10.4139, 4.82098, 8.97737, -1.8358,
-0.307942, -2.76266, 3.04212, 4.82593, 2.86925, 8.98986, 6.18658,
2.50114, 12.5518, 5.65282, -2.42181, 4.89983, 6.27595, 6.96756,
8.34183, 11.6484, 10.7812, 4.89868, 10.6572, 6.47331, 5.90388,
12.2386, 8.55759, 4.80869, 7.95166, 12.7987, 7.50973, 9.95841,
8.38024, 6.01908, 2.00816, 1.95061, 1.21553, 1.10895, 6.41019,
7.15468, 7.4428, 9.07871, 3.4927, 5.06634, 1.95723, 3.36867,
8.81967, 4.19371, 5.46829, 6.78036, 4.89139, 8.33048, 8.78568,
7.04015, 3.66143, 1.45926, 5.99856, 2.91786, 8.73688, 7.61774,
10.5881, 2.97864, 4.66389, 9.06656, 5.5244, 6.61624, 8.77095,
6.29498, 3.21623, 10.311, 6.86147, 0.993279, 10.5456, 5.84009,
8.93219, 2.33313, 3.76373, -0.608017, 10.0227, 2.9371, 8.07382,
6.45556, 7.73459, 4.5517, 4.75401, 5.23356, 3.93879, 7.21409,
5.46405, 5.8554, 3.15436, 5.49279, 8.38921), X_mi_m = c(0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L), X_mi_id = c(1L, 2L, 3L, 4L, 5L, 6L,
7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L,
20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L,
33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L,
46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 55L, 56L, 57L, 58L,
59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L,
72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L,
85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L, 93L, 94L, 95L, 96L, 97L,
98L, 99L, 100L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L,
50L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 41L,
42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 91L, 92L, 93L, 94L,
95L, 96L, 97L, 98L, 99L, 100L, 41L, 42L, 43L, 44L, 45L, 46L,
47L, 48L, 49L, 50L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L,
100L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 91L,
92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 41L, 42L, 43L,
44L, 45L, 46L, 47L, 48L, 49L, 50L, 91L, 92L, 93L, 94L, 95L, 96L,
97L, 98L, 99L, 100L), X_mi_miss = c(0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA)), class = "data.frame", row.names = c(NA, -200L))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment