全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
13200 24
2013-05-09
最近在做缺失数据填补的分析,一直在用matlab做,实现了EM填补和一些插值填补,但是写不出MI填补的代码
看了看相关论文,MI填补基本都用SAS做的,就想试一试
我想请问一下:
1、SAS对数据集本身有什么要求吗?比如格式和大小,我处理的数据集都算是比较大的;
2、我只知道SAS里有现成的MI过程,那应用这个过程麻烦吗?我毫无SAS基础
3、填补完的数据集可以导出成txt或者逗号分隔值文件吗?
谢谢各位大神!~
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-5-9 17:47:29
1. NO;
2. Proc MI, pls refer to sas help and documentation for help;
3. Yes.
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-9 19:38:13
boe 发表于 2013-5-9 17:47
1. NO;
2. Proc MI, pls refer to sas help and documentation for help;
3. Yes.
谢谢~那再请问我导入数据集的时候,在我要对这个数据集使用MI过程的情况下,需要对这个数据集做那些预处理呢?还是我什么都不用做像matlab一样直接导入就可以了呢?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-9 19:45:26
一夜可以走多远 发表于 2013-5-9 19:38
谢谢~那再请问我导入数据集的时候,在我要对这个数据集使用MI过程的情况下,需要对这个数据集做那些预处理 ...
不需要预处理,直接导入即可
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-9 20:37:46
曾经做过一个,程序供你参考。本论坛也有SAS MI过程。
data ttest1_1;
        input p0 m0 p1 m1 p2 m2 p3 m3@@;
datalines;
144        11.1        74        10.8        58        13.8        48        10.3
156        10.9        174        10.9        180        10        244        10.7
94        11.7        77        13        77        10.4        69        13
205        9.7        200        10        220        10.4        197        10.1
90        10        56        12.5        37        10.7        39        10.3
247        13        229        12.9        248        12.5        277        11.2
149        9.5        134        9.7        115        9.7        124        9.8
156        11.5        79        11.4        60        12.4        59        13.6
270        10.7        232        10.8        209        11        220        10.9
146        11.1        157        10.8        135        12.1        145        13.3
108        11.3        112        13.4        99        13.1        84        13.1
90        11.1        98        10.2        120        10        146        10.2
67        10.2        61        9.8        39        10.3        48        10.4
646        11.5        425        11.3        372        11.9        295        11.9
121        11.8        107        10.2        86        10.4        92        10.8
261        10        179        10.6        179        11.7        167        12
91        12.3        230        10.9        194        11.3        202        10.4
181        11        165        10.9        148        10.3        177        10.7
229        10        234        10.3        218        11        296        11.5
229        10        234        10.3        218        11        296        11.5
186        10.2        184        9.8        137        9.6        83        9.9
160        13.1        199        12.1        207        12.1        221        12.2
124        10.5        125        10.6        101        11.1        84        12.1
298        10.1        221        9.9        206        10.2        128        10.1
173        10.4        148        10.3        173        10.4        192        10.4
24        8.2        14        11.4        5        .        34        11.8
92        10.8        72        11        64        13        89        12.6
164        12.8        108        12.2        76        12.1        48        12.8
146        12.3        133        12.3        120        11.8        127        10.4
162        12.3        135        12        124        12.3        130        13.1
133        11.5        154        11        192        10.4        249        10.4
106        14.5        88        14        71        11.2        60        11.6
139        10.9        151        10.7        145        11.2        146        11.1
114        11.4        83        12        72        12.5        33        9.1
139        12.5        127        12.9        86        14        34        11.8
159        13.4        164        12.8        146        13.9        175        13.9
107        12.3        188        10.4        178        10.7        182        12.1
298        12.8        226        12.8        198        12.6        179        10.1
50        11.9        56        10.6        51        11.7        73        10.9
180        10.6        160        11.5        122        11.6        114        12.6
76        12.1        41        11.2        32        11.9        34        11.8
69        10.1        67        11.6        71        11.2        100        14.3
189        12.2        208        11.1        239        11        249        10.8
126        13.8        126        12.9        103        14.9        114        13.5
188        10.8        161        11.4        107        11.5        73        10.6
113        9.9        44        11.8        42        11        62        11.6
116        14.1        83        13.7        73        12.6        54        11.5
114        11.4        59        10.2        88        11.4        82        11
74        12.4        62        11.6        81        11.4        54        13.3
146        10.2        124        10.5        102        11.8        116        10.3
244        12.7        219        12.8        217        12.4        357        12.1
478        10.3        392        10        361        10        310        10
94        13.1        68        12.1        59        12.2        59        11.8
173        11        135        11.3        107        11.4        117        12.3
160        12.7        159        11.9        137        12.4        135        12.6
153        10.5        170        11.2        166        11.4        144        11.3
688        11        564        11        390        11.1        293        11.4
120        11.8        125        11.4        134        11.3        159        11.9
57        10.5        48        11        35        11.4        31        9.7
104        13.5        73        14.2        70        13.1        77        13.2
141        12.1        111        11.9        114        12.5        138        13.3
112        12.7        121        12.6        128        11.9        151        13.5
211        13.4        200        13.7        191        14.3        159        10
170        10.6        139        10.9        88        10.5        65        12.6
126        12.9        124        13.1        91        12.4        71        13
132        11.7        159        11.8        195        11.6        250        10.8
114        12.5        100        11.2        103        11.8        144        11.8
48        15.3        49        13        61        13.7        71        12.3
92        11.1        85        10.8        94        10.9        105        11.6
163        10.4        119        11.9        136        11.2        127        11.2
105        12.6        101        9.9        75        9.9        99        9.5
129        10.7        103        12.8        108        10.3        114        11.1
78        14.5        61        10.5        93        14        106        13
177        11.4        180        11.2        174        10.4        206        10
147        13.8        81        11        63        12.2        48        10.4
94        10.1        79        10.6        86        10.7        88        11.9
204        8.5        226        8.2        105        8.3        99        8.7
126        8.1        118        7.7        82        7.8        85        7.5
100        10.4        82        10.7        76        10.1        102        9.3
144        10.5        131        9.8        125        10.3        151        10.9
112        7.2        73        6.6        60        6.7        58        6.9
94        9.5        89        9.8        89        9.4        123        9.3
196        8.9        162        9.6        114        10        146        10.2
47        10.6        51        13.7        39        12.8        30        10
189        12.6        187        12.5        163        12.9        161        9.9
63        9.3        47        9.6        45        8.9        59        9.3
98        12.3        70        14.7        61        15.6        52        13.7
111        12.9        86        11.6        68        13.5        54        10.2
97        9.3        103        9.2        84        10.7        78        10.3
299        11.9        298        11.3        391        12        389        11.1
161        13.4        180        12.8        176        11.3        232        12.1
72        11.3        80        11.4        69        11.2        82        11.2
150        14.8        183        13.3        114        13.3        108        14.1
99        15.2        78        12.6        97        15.7        93        13.1
138        14        106        14.5        104        13        133        12.8
179        13.2        149        13.1        121        13.1        95        13.3
261        9.6        213        10.3        208        10.1        174        10.3
189        10.6        146        10.3        195        10.4        148        11.4
116        12        91        11.5        107        13.4        132        14.2
224        8.3        193        8        200        10        147        8.6
175        10.8        155        11.4        131        10.2        98        10.6
143        11.4        104        13.6        107        13.7        118        14.5
103        10.6        81        11.1        51        13.9        41        13.9
180        11.7        144        11.8        117        12.1        106        11.5
161        9.3        127        9.1        65        10        46        9.8
120        13.6        .        .        .        .        71        13.6
133        10.5        104        10.6        .        .        92        12
112        13.5        .        .        36        11        .        .
242        12.8        173        12.1        .        .        .        .
134        9.7        127        15.2        .        .        87        10.3
148        12.2        44        11.4        .        .        15        13.3
177        12.4        .        .        68        10.3        96        10.2
61        10.2        56        11.1        46        11.2        .        .
125        11.4        131        12.2        162        11.7        .        .
352        11.4        292        11.3        .        .        159        12.6
85        12        81        11.4        .        .        64        12.5
145        11.7        130        12.3        96        11.7        .        .
114        11.6        .        .        .        .        162        13
161        12.4        .        .        44        11.3        45        13.7
139        12.9        99        13.5        .        .        49        10.2
86        11.6        43        11.6        .        .        31        9.7
189        11.6        .        .        152        11.2        .        .
35        11.4        5        .        69        11.6        .        .
210        11        104        11.5        67        10.4        72        11.1
103        12.8        .        .        118        13        108        13.1
173        11        175        10.9        .        .        88        12.5
378        9.8        .        .        .        .        160        13.1
37        10.8        22        9.1        .        .        5        .
245        9.8        .        .        .        .        237        10.5
199        11.1        .        .        202        11.4        210        12.4
94        11.9        82        12.4        81        14        .        .
.        .        .        .        47        8.9        32        9.9
.        .        185        9.7        170        11.2        .        .
175        13.1        243        13.6        .        .        107        13.1
80        12.8        .        .        74        12.4        .        .
100        12        .        .        33        10.9        .        .
158        13.9        43        11.6        12        16.7        .        .
79        12.7        53        13.2        46        13        .        .
.        .        71        11.3        99        11.1        114        10.5
213        12.7        .        .        94        13        .        .
259        8.3        .        .        247        10.9        195        10.8
.        .        83        10.8        .        .        94        10.6
274        8.9        259        9.5        266        9.1        .        .
398        10.6        .        .        358        11.5        354        11.6
.        .        135        14.1        128        13.3        .        .
62        9.7        37        10        32        10        47        10.6
39        10.3        22        13.6        .        .        23        8.7
;
proc mi data=ttest1_1 seed=2012 simple out=outmpv;
mcmc chain=multiple initial=EM;
var p0 m0 p1 m1 p2 m2 p3 m3;
run;
proc print data=outmpv;
run;
data ttest1_2;
        set outmpv;
        p01=p0-p1;
        p12=p1-p2;
        p23=p2-p3;
        p03=p0-p3;
        m01=m0-m1;
        m12=m1-m2;
        m23=m2-m3;
        m03=m0-m3;
proc ttest data=ttest1_2;
        var p01;
        var p12;
        var p23;
        var p03;
        var m01;
        var m12;
        var m23;
        var m03;
run;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-9 21:35:38
boe 发表于 2013-5-9 19:45
不需要预处理,直接导入即可
谢谢~~
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群