1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
|
R version 2.4.0 Under development (unstable) (2006-08-11 r38834)
Copyright (C) 2006 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(foreign)
> pc5 <- read.dta("pc5.dta")
Warning message:
cannot read factor labels from Stata 5 files in: read.dta("pc5.dta")
> summary(pc5)
age alb alkphos ascites
Min. :26.28 Min. :1.960 Min. : -9 Min. :-9.000
1st Qu.:42.71 1st Qu.:3.330 1st Qu.: 423 1st Qu.: 0.000
Median :50.47 Median :3.570 Median : 1009 Median : 0.000
Mean :50.45 Mean :3.543 Mean : 1525 Mean :-2.074
3rd Qu.:56.99 3rd Qu.:3.820 3rd Qu.: 1677 3rd Qu.: 0.000
Max. :78.44 Max. :4.640 Max. :13862 Max. : 1.000
bili chol edema edemarx
Min. : 0.300 Min. : -9.0 Min. :0.00000 Min. :0.00000
1st Qu.: 0.700 1st Qu.: -9.0 1st Qu.:0.00000 1st Qu.:0.00000
Median : 1.200 Median : 256.0 Median :0.00000 Median :0.00000
Mean : 2.688 Mean : 245.8 Mean :0.09256 Mean :0.07815
3rd Qu.: 3.000 3rd Qu.: 346.0 3rd Qu.:0.00000 3rd Qu.:0.00000
Max. :28.000 Max. :1775.0 Max. :1.00000 Max. :1.00000
hepmeg time plate protime
Min. :-9.000 Min. : 41 Min. : -9.0 Min. : 9.00
1st Qu.: 0.000 1st Qu.:1434 1st Qu.:181.0 1st Qu.:10.00
Median : 0.000 Median :1487 Median :249.0 Median :10.60
Mean :-1.750 Mean :1760 Mean :252.0 Mean :10.67
3rd Qu.: 1.000 3rd Qu.:2153 3rd Qu.:321.0 3rd Qu.:11.00
Max. : 1.000 Max. :4795 Max. :721.0 Max. :18.00
sex sgot spiders stage
Min. :-9.000 Min. : -9.00 Min. :-9.000 Min. :-9.0000
1st Qu.: 0.000 1st Qu.: 44.20 1st Qu.: 0.000 1st Qu.: 1.0000
Median : 1.000 Median : 88.35 Median : 0.000 Median : 3.0000
Mean :-1.434 Mean : 88.22 Mean :-1.914 Mean : 0.1351
3rd Qu.: 1.000 3rd Qu.:130.20 3rd Qu.: 0.000 3rd Qu.: 3.0000
Max. : 1.000 Max. :457.25 Max. : 1.000 Max. : 4.0000
cens rx trig copper
Min. :0.0000 Min. :-9.0000 Min. : -9.00 Min. : -9.00
1st Qu.:0.0000 1st Qu.: 1.0000 1st Qu.: -9.00 1st Qu.: 10.00
Median :0.0000 Median : 1.0000 Median : 85.00 Median : 48.00
Mean :0.2443 Mean :-0.9788 Mean : 81.44 Mean : 65.26
3rd Qu.:0.0000 3rd Qu.: 2.0000 3rd Qu.:126.00 3rd Qu.: 88.00
Max. :1.0000 Max. : 2.0000 Max. :598.00 Max. :588.00
id first late t0
Min. : 1.0 Min. :0.0000 Min. :0.0000 Min. : 0.0
1st Qu.: 99.5 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.: 0.0
Median :195.0 Median :0.0000 Median :0.0000 Median : 0.0
Mean :199.6 Mean :0.3657 Mean :0.3657 Mean : 543.8
3rd Qu.:302.5 3rd Qu.:1.0000 3rd Qu.:1.0000 3rd Qu.:1487.0
Max. :418.0 Max. :1.0000 Max. :1.0000 Max. :1487.0
> str(pc5)
'data.frame': 659 obs. of 24 variables:
$ age : num 58.8 56.4 56.4 70.1 54.7 ...
$ alb : num 2.60 4.14 4.14 3.48 2.54 ...
$ alkphos: num 1718 7395 7395 516 6122 ...
$ ascites: num 1 0 0 0 0 0 0 0 0 0 ...
$ bili : num 14.5 1.1 1.1 1.4 1.8 ...
$ chol : num 261 302 302 176 244 244 279 279 248 248 ...
$ edema : num 1 0 0 1 1 1 0 0 0 0 ...
$ edemarx: num 1 0 0 0.5 0.5 0.5 0 0 0 0 ...
$ hepmeg : num 1 1 1 0 1 1 1 1 1 1 ...
$ time : num 400 1487 4500 1012 1487 ...
$ plate : num 190 221 221 151 183 183 136 136 -9 -9 ...
$ protime: num 12.2 10.6 10.6 12.0 10.3 ...
$ sex : num 1 1 1 0 1 1 1 1 1 1 ...
$ sgot : num 137.9 113.5 113.5 96.1 60.6 ...
$ spiders: num 1 1 1 0 1 1 1 1 0 0 ...
$ stage : num 4 3 3 4 4 4 3 3 3 3 ...
$ cens : num 1 0 0 1 0 1 0 0 0 1 ...
$ rx : num 1 1 1 1 1 1 2 2 2 2 ...
$ trig : num 172 88 88 55 92 92 72 72 63 63 ...
$ copper : num 156 54 54 210 64 64 143 143 50 50 ...
$ id : num 1 2 2 3 4 4 5 5 6 6 ...
$ first : num 0 1 0 0 1 0 1 0 1 0 ...
$ late : num 0 0 1 0 0 1 0 1 0 1 ...
$ t0 : num 0 0 1487 0 0 ...
- attr(*, "datalabel")= chr ""
- attr(*, "time.stamp")= chr "14 Feb 1997 14:22"
- attr(*, "formats")= chr "%9.0g" "%9.0g" "%9.0g" "%9.0g" ...
- attr(*, "types")= int 102 102 102 102 102 102 102 102 102 102 ...
- attr(*, "val.labels")= chr "" "" "" "" ...
- attr(*, "var.labels")= chr "" "" "" "" ...
- attr(*, "version")= int 5
> compressed<-read.dta("compressed.dta")
> summary(compressed)
age alb alkphos ascites
Min. :26.28 Min. :1.960 Min. : -9 Min. :-9.000
1st Qu.:42.71 1st Qu.:3.330 1st Qu.: 423 1st Qu.: 0.000
Median :50.47 Median :3.570 Median : 1009 Median : 0.000
Mean :50.45 Mean :3.543 Mean : 1525 Mean :-2.074
3rd Qu.:56.99 3rd Qu.:3.820 3rd Qu.: 1677 3rd Qu.: 0.000
Max. :78.44 Max. :4.640 Max. :13862 Max. : 1.000
bili chol edema edemarx
Min. : 0.300 Min. : -9.0 Min. :0.00000 Min. :0.00000
1st Qu.: 0.700 1st Qu.: -9.0 1st Qu.:0.00000 1st Qu.:0.00000
Median : 1.200 Median : 256.0 Median :0.00000 Median :0.00000
Mean : 2.688 Mean : 245.8 Mean :0.09256 Mean :0.07815
3rd Qu.: 3.000 3rd Qu.: 346.0 3rd Qu.:0.00000 3rd Qu.:0.00000
Max. :28.000 Max. :1775.0 Max. :1.00000 Max. :1.00000
hepmeg time plate protime
Min. :-9.000 Min. : 41 Min. : -9.0 Min. : 9.00
1st Qu.: 0.000 1st Qu.:1434 1st Qu.:181.0 1st Qu.:10.00
Median : 0.000 Median :1487 Median :249.0 Median :10.60
Mean :-1.750 Mean :1760 Mean :252.0 Mean :10.67
3rd Qu.: 1.000 3rd Qu.:2153 3rd Qu.:321.0 3rd Qu.:11.00
Max. : 1.000 Max. :4795 Max. :721.0 Max. :18.00
sex sgot spiders stage
Min. :-9.000 Min. : -9.00 Min. :-9.000 Min. :-9.0000
1st Qu.: 0.000 1st Qu.: 44.20 1st Qu.: 0.000 1st Qu.: 1.0000
Median : 1.000 Median : 88.35 Median : 0.000 Median : 3.0000
Mean :-1.434 Mean : 88.22 Mean :-1.914 Mean : 0.1351
3rd Qu.: 1.000 3rd Qu.:130.20 3rd Qu.: 0.000 3rd Qu.: 3.0000
Max. : 1.000 Max. :457.25 Max. : 1.000 Max. : 4.0000
cens rx trig copper
Min. :0.0000 Min. :-9.0000 Min. : -9.00 Min. : -9.00
1st Qu.:0.0000 1st Qu.: 1.0000 1st Qu.: -9.00 1st Qu.: 10.00
Median :0.0000 Median : 1.0000 Median : 85.00 Median : 48.00
Mean :0.2443 Mean :-0.9788 Mean : 81.44 Mean : 65.26
3rd Qu.:0.0000 3rd Qu.: 2.0000 3rd Qu.:126.00 3rd Qu.: 88.00
Max. :1.0000 Max. : 2.0000 Max. :598.00 Max. :588.00
id first late t0
Min. : 1.0 Min. :0.0000 Min. :0.0000 Min. : 0.0
1st Qu.: 99.5 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.: 0.0
Median :195.0 Median :0.0000 Median :0.0000 Median : 0.0
Mean :199.6 Mean :0.3657 Mean :0.3657 Mean : 543.8
3rd Qu.:302.5 3rd Qu.:1.0000 3rd Qu.:1.0000 3rd Qu.:1487.0
Max. :418.0 Max. :1.0000 Max. :1.0000 Max. :1487.0
> all.equal(summary(pc5), summary(compressed))
[1] TRUE
> sun6 <- read.dta("sun6.dta")
> summary(sun6)
age alb alkphos ascites
Min. :26.28 Min. :1.960 Min. : -9 Min. :-9.000
1st Qu.:42.71 1st Qu.:3.330 1st Qu.: 423 1st Qu.: 0.000
Median :50.47 Median :3.570 Median : 1009 Median : 0.000
Mean :50.45 Mean :3.543 Mean : 1525 Mean :-2.074
3rd Qu.:56.99 3rd Qu.:3.820 3rd Qu.: 1677 3rd Qu.: 0.000
Max. :78.44 Max. :4.640 Max. :13862 Max. : 1.000
bili chol edema edemarx
Min. : 0.300 Min. : -9.0 Min. :0.00000 Min. :0.00000
1st Qu.: 0.700 1st Qu.: -9.0 1st Qu.:0.00000 1st Qu.:0.00000
Median : 1.200 Median : 256.0 Median :0.00000 Median :0.00000
Mean : 2.688 Mean : 245.8 Mean :0.09256 Mean :0.07815
3rd Qu.: 3.000 3rd Qu.: 346.0 3rd Qu.:0.00000 3rd Qu.:0.00000
Max. :28.000 Max. :1775.0 Max. :1.00000 Max. :1.00000
hepmeg time plate protime
Min. :-9.000 Min. : 41 Min. : -9.0 Min. : 9.00
1st Qu.: 0.000 1st Qu.:1434 1st Qu.:181.0 1st Qu.:10.00
Median : 0.000 Median :1487 Median :249.0 Median :10.60
Mean :-1.750 Mean :1760 Mean :252.0 Mean :10.67
3rd Qu.: 1.000 3rd Qu.:2153 3rd Qu.:321.0 3rd Qu.:11.00
Max. : 1.000 Max. :4795 Max. :721.0 Max. :18.00
sex sgot spiders stage
Min. :-9.000 Min. : -9.00 Min. :-9.000 Min. :-9.0000
1st Qu.: 0.000 1st Qu.: 44.20 1st Qu.: 0.000 1st Qu.: 1.0000
Median : 1.000 Median : 88.35 Median : 0.000 Median : 3.0000
Mean :-1.434 Mean : 88.22 Mean :-1.914 Mean : 0.1351
3rd Qu.: 1.000 3rd Qu.:130.20 3rd Qu.: 0.000 3rd Qu.: 3.0000
Max. : 1.000 Max. :457.25 Max. : 1.000 Max. : 4.0000
cens rx trig copper
Min. :0.0000 Min. :-9.0000 Min. : -9.00 Min. : -9.00
1st Qu.:0.0000 1st Qu.: 1.0000 1st Qu.: -9.00 1st Qu.: 10.00
Median :0.0000 Median : 1.0000 Median : 85.00 Median : 48.00
Mean :0.2443 Mean :-0.9788 Mean : 81.44 Mean : 65.26
3rd Qu.:0.0000 3rd Qu.: 2.0000 3rd Qu.:126.00 3rd Qu.: 88.00
Max. :1.0000 Max. : 2.0000 Max. :598.00 Max. :588.00
id first late t0
Min. : 1.0 Min. :0.0000 Min. :0.0000 Min. : 0.0
1st Qu.: 99.5 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.: 0.0
Median :195.0 Median :0.0000 Median :0.0000 Median : 0.0
Mean :199.6 Mean :0.3657 Mean :0.3657 Mean : 543.8
3rd Qu.:302.5 3rd Qu.:1.0000 3rd Qu.:1.0000 3rd Qu.:1487.0
Max. :418.0 Max. :1.0000 Max. :1.0000 Max. :1487.0
> str(sun6)
'data.frame': 659 obs. of 24 variables:
$ age : num 58.8 56.4 56.4 70.1 54.7 ...
$ alb : num 2.60 4.14 4.14 3.48 2.54 ...
$ alkphos: num 1718 7395 7395 516 6122 ...
$ ascites: int 1 0 0 0 0 0 0 0 0 0 ...
$ bili : num 14.5 1.1 1.1 1.4 1.8 ...
$ chol : int 261 302 302 176 244 244 279 279 248 248 ...
$ edema : int 1 0 0 1 1 1 0 0 0 0 ...
$ edemarx: num 1 0 0 0.5 0.5 0.5 0 0 0 0 ...
$ hepmeg : int 1 1 1 0 1 1 1 1 1 1 ...
$ time : int 400 1487 4500 1012 1487 1925 1487 1504 1487 2503 ...
$ plate : int 190 221 221 151 183 183 136 136 -9 -9 ...
$ protime: num 12.2 10.6 10.6 12.0 10.3 ...
$ sex : int 1 1 1 0 1 1 1 1 1 1 ...
$ sgot : num 137.9 113.5 113.5 96.1 60.6 ...
$ spiders: int 1 1 1 0 1 1 1 1 0 0 ...
$ stage : int 4 3 3 4 4 4 3 3 3 3 ...
$ cens : int 1 0 0 1 0 1 0 0 0 1 ...
$ rx : int 1 1 1 1 1 1 2 2 2 2 ...
$ trig : int 172 88 88 55 92 92 72 72 63 63 ...
$ copper : int 156 54 54 210 64 64 143 143 50 50 ...
$ id : int 1 2 2 3 4 4 5 5 6 6 ...
$ first : int 0 1 0 0 1 0 1 0 1 0 ...
$ late : int 0 0 1 0 0 1 0 1 0 1 ...
$ t0 : int 0 0 1487 0 0 1487 0 1487 0 1487 ...
- attr(*, "datalabel")= chr ""
- attr(*, "time.stamp")= chr "14 Sep 2000 10:40"
- attr(*, "formats")= chr "%9.0g" "%9.0g" "%9.0g" "%9.0g" ...
- attr(*, "types")= int 102 102 102 98 102 105 98 102 98 105 ...
- attr(*, "val.labels")= chr "" "" "" "" ...
- attr(*, "var.labels")= chr "" "" "" "" ...
- attr(*, "version")= int 6
- attr(*, "label.table")=List of 24
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
> all.equal(summary(sun6),summary(pc5))
[1] TRUE
> df<-read.dta("datefactor.dta")
> summary(df)
hlth159 id adate
EXCELLENT: 0 Min. : 1.00 Min. :1960-01-02
VERY GOOD: 6 1st Qu.: 8.25 1st Qu.:1960-01-10
GOOD : 8 Median :15.50 Median :1960-01-17
FAIR : 4 Mean :15.50 Mean :1960-01-17
POOR : 1 3rd Qu.:22.75 3rd Qu.:1960-01-24
NA's :11 Max. :30.00 Max. :1960-01-31
> data(esoph)
> write.dta(esoph,esophile<-tempfile())
Warning message:
abbreviating variable names in: write.dta(esoph, esophile <- tempfile())
> esoph2<-read.dta(esophile)
> all.equal(ordered(esoph2$alcgp),esoph$alcgp)
[1] TRUE
> write.dta(esoph,esophile,convert.factors="string")
Warning message:
abbreviating variable names in: write.dta(esoph, esophile, convert.factors = "string")
> esoph2<-read.dta(esophile)
> all.equal(as.character(esoph$alcgp),esoph2$alcgp)
[1] TRUE
> write.dta(esoph,esophile,convert.factors="code")
Warning message:
abbreviating variable names in: write.dta(esoph, esophile, convert.factors = "code")
> esoph2<-read.dta(esophile)
> all.equal(as.numeric(esoph$alcgp),as.numeric(esoph2$alcgp))
[1] TRUE
>
> se<-read.dta("stata7se.dta")
> print(se)
race number
1 white 2
2 asian 5
3 hispanic 5
4 white 4
5 black 5
6 white 6
7 black 7
> v8<-read.dta("stata8mac.dta")
> print(v8)
race number
1 white 2
2 asian 5
3 hispanic 5
4 white 4
5 black 5
6 white 6
7 black 7
>
> stata8<-read.dta("auto8.dta",missing=TRUE,convert.underscore=FALSE)
> str(stata8)
'data.frame': 74 obs. of 13 variables:
$ make : chr "AMC Concord" "AMC Pacer" "AMC Spirit" "Buick Century" ...
$ price : int 4099 4749 3799 4816 7827 5788 4453 5189 10372 4082 ...
$ mpg : int 22 17 22 20 15 18 26 20 16 19 ...
$ rep78 : int 3 3 NA 3 4 3 NA 3 3 3 ...
$ headroom : num 2.5 3 3 4.5 4 4 3 2 3.5 3.5 ...
$ trunk : int 11 11 12 16 20 21 10 16 17 13 ...
$ weight : int 2930 3350 2640 3250 4080 3670 2230 3280 3880 3400 ...
$ length : int 186 173 168 196 222 218 170 200 207 200 ...
$ turn : int 40 40 35 40 43 43 34 42 43 42 ...
$ displacement: int 121 258 121 196 350 231 304 196 231 231 ...
$ gear_ratio : num 3.58 2.53 3.08 2.93 2.41 ...
$ foreign : Factor w/ 2 levels "Domestic","Foreign": 1 1 1 1 1 1 1 1 1 1 ...
$ testmiss : num NA NA NA NA NA NA NA NA NA NA ...
- attr(*, "datalabel")= chr "1978 Automobile Data"
- attr(*, "time.stamp")= chr "20 May 2003 14:39"
- attr(*, "formats")= chr "%-18s" "%8.0gc" "%8.0g" "%8.0g" ...
- attr(*, "types")= int 18 252 252 252 254 252 252 252 252 252 ...
- attr(*, "val.labels")= chr "" "" "" "" ...
- attr(*, "var.labels")= chr "Make and Model" "Price" "Mileage (mpg)" "Repair Record 1978" ...
- attr(*, "version")= int -8
- attr(*, "label.table")=List of 13
..$ origin: Named num 0 1
.. ..- attr(*, "names")= chr "Domestic" "Foreign"
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
..$ : NULL
- attr(*, "missing")=List of 13
..$ make : NULL
..$ price : num NA NA NA NA NA NA NA NA NA NA ...
..$ mpg : num NA NA NA NA NA NA NA NA NA NA ...
..$ rep78 : num NA NA 0 NA NA NA 0 NA NA NA ...
..$ headroom : num NA NA NA NA NA NA NA NA NA NA ...
..$ trunk : num NA NA NA NA NA NA NA NA NA NA ...
..$ weight : num NA NA NA NA NA NA NA NA NA NA ...
..$ length : num NA NA NA NA NA NA NA NA NA NA ...
..$ turn : num NA NA NA NA NA NA NA NA NA NA ...
..$ displacement: num NA NA NA NA NA NA NA NA NA NA ...
..$ gear_ratio : num NA NA NA NA NA NA NA NA NA NA ...
..$ foreign : num NA NA NA NA NA NA NA NA NA NA ...
..$ testmiss : num 18 18 18 18 18 18 18 18 18 18 ...
>
> q()
|
