Hi, I am trying to run a conditional logistic model on a nested case-control
study using cph() and then estimate survival based on the model. The data came
from Prof Bryan Langholz website where he also has the SAS code to this, so I
am
trying to replicate the SAS results.
The data attached. Basically, the variables are:
rstime: risk set age
rsentry: fake entry time, just before rstime
setno: risk set identifier
cc: lung cancer death (0/1)
cr500: independent variable 1 (0/1), cumulative radon > 500WLM
smoke25: independent variable 2 (0/1), 1/2+ pack/dat at age 25
logw: weight to be used as an offset in the model
rstime2: indicator variable to indicate whether age is below 50, 50-69 or >=70,
used as strata in the model
The SAS code is (if this helps to address my question/problem)
proc phreg data=uminers;
model rstime*cc(0)=cr500 smoke25/ entry=rsentry offset=logw rl;
baseline out=ch covariates=covs cumhaz=cumhaz stdcumhaz=secumhaz
lowercumhaz=lci
uppercumhaz=uci/ nomean;
strata rstime2;
run;
My R code is:
library(survival)
library(Design)
> uminers<-read.table("uminers.txt",sep='\t',header=T,row.names=NULL)
> fit <- cph(Surv(rsentry,rstime,cc)~cr500+smoke25+strat(rstime2)
> +offset(logw),
>uminers, x=T, y=T)
> fit
Cox Proportional Hazards Model
cph(formula = Surv(rsentry, rstime, cc) ~ cr500 + smoke25 + strat(rstime2) +
offset(logw), data = uminers, x = T, y = T)
Obs Events Model L.R. d.f. P
774 258 89.98 2 0
Score Score P R2
82.61 0 0.211
Status
Stratum No Event Event
rstime2=1 100 50
rstime2=2 348 174
rstime2=3 68 34
coef se(coef) z p
cr500 1.491 0.184 8.08 6.66e-16
smoke25 0.461 0.182 2.53 1.16e-02
These results are the same as what SAS reported!
Next, I want to estimate the survival for stratum age 50-69 (rstime2=2), so I
used the following, this is where I got the error message:
> survest(fit, expand.grid(cr500=c(0,1),smoke25=c(0,1),rstime2=2, logw=0),
>times=unique(uminers$rstime[uminers$rstime2=='2']), conf.int=.95)
Error in factor(rep(1:nstrat, scount), labels = names(fit$strata)) :
invalid labels; length 3 should be 1 or 2
Anyone has any suggestions what went wrong?
Thanks!
John
setno rstime cc rsentry cr500 smoke25 logw rstime2
1 34.5 1 34.4999 1 1 6.315358002 1
1 34.5 0 34.4999 1 0 6.315358002 1
1 34.5 0 34.4999 0 1 6.315358002 1
2 36.42 1 36.4199 1 1 6.404125922 1
2 36.42 0 36.4199 0 0 6.404125922 1
2 36.42 0 36.4199 1 1 6.404125922 1
3 38.17 1 38.1699 1 0 6.470283375 1
3 38.17 0 38.1699 0 1 6.470283375 1
3 38.17 0 38.1699 0 1 6.470283375 1
4 39.92 1 39.9199 1 1 6.532334292 1
4 39.92 0 39.9199 0 0 6.532334292 1
4 39.92 0 39.9199 1 0 6.532334292 1
5 40 1 39.9999 1 1 6.532334292 1
5 40 0 39.9999 0 1 6.532334292 1
5 40 0 39.9999 1 1 6.532334292 1
6 40.08 1 40.0799 0 1 6.537657315 1
6 40.08 0 40.0799 1 1 6.537657315 1
6 40.08 0 40.0799 1 1 6.537657315 1
7 41.92 1 41.9199 1 1 6.582486713 1
7 41.92 0 41.9199 1 1 6.582486713 1
7 41.92 0 41.9199 0 0 6.582486713 1
8 42.08 1 42.0799 1 0 6.592130875 1
8 42.08 0 42.0799 1 1 6.592130875 1
8 42.08 0 42.0799 0 1 6.592130875 1
9 42.081 1 42.0809 1 1 6.593044534 1
9 42.081 0 42.0809 0 0 6.593044534 1
9 42.081 0 42.0809 1 0 6.593044534 1
10 42.42 1 42.4199 1 1 6.599870499 1
10 42.42 0 42.4199 1 1 6.599870499 1
10 42.42 0 42.4199 0 0 6.599870499 1
11 42.5 1 42.4999 0 1 6.598509029 1
11 42.5 0 42.4999 1 1 6.598509029 1
11 42.5 0 42.4999 0 1 6.598509029 1
12 42.92 1 42.9199 1 1 6.598054793 1
12 42.92 0 42.9199 1 1 6.598054793 1
12 42.92 0 42.9199 1 1 6.598054793 1
13 43.75 1 43.7499 1 1 6.612041035 1
13 43.75 0 43.7499 0 1 6.612041035 1
13 43.75 0 43.7499 0 1 6.612041035 1
14 44.33 1 44.3299 1 0 6.605748879 1
14 44.33 0 44.3299 0 1 6.605748879 1
14 44.33 0 44.3299 1 0 6.605748879 1
15 44.331 1 44.3309 1 1 6.6147256 1
15 44.331 0 44.3309 1 0 6.6147256 1
15 44.331 0 44.3309 1 1 6.6147256 1
16 44.67 1 44.6699 0 1 6.613384218 1
16 44.67 0 44.6699 0 1 6.613384218 1
16 44.67 0 44.6699 0 0 6.613384218 1
17 44.92 1 44.9199 1 0 6.611592906 1
17 44.92 0 44.9199 0 1 6.611592906 1
17 44.92 0 44.9199 0 1 6.611592906 1
18 45.17 1 45.1699 1 1 6.611144576 1
18 45.17 0 45.1699 1 0 6.611144576 1
18 45.17 0 45.1699 0 1 6.611144576 1
19 45.25 1 45.2499 0 0 6.610696045 1
19 45.25 0 45.2499 0 0 6.610696045 1
19 45.25 0 45.2499 1 0 6.610696045 1
20 45.42 1 45.4199 1 1 6.616511315 1
20 45.42 0 45.4199 1 0 6.616511315 1
20 45.42 0 45.4199 1 1 6.616511315 1
21 45.5 1 45.4999 1 1 6.613831546 1
21 45.5 0 45.4999 1 1 6.613831546 1
21 45.5 0 45.4999 0 0 6.613831546 1
22 45.58 1 45.5799 1 1 6.610696045 1
22 45.58 0 45.5799 0 1 6.610696045 1
22 45.58 0 45.5799 0 1 6.610696045 1
23 45.67 1 45.6699 1 0 6.608000625 1
23 45.67 0 45.6699 1 0 6.608000625 1
23 45.67 0 45.6699 0 0 6.608000625 1
24 45.75 1 45.7499 1 1 6.608450367 1
24 45.75 0 45.7499 1 0 6.608450367 1
24 45.75 0 45.7499 1 1 6.608450367 1
25 46 1 45.9999 1 0 6.605748879 1
25 46 0 45.9999 1 1 6.605748879 1
25 46 0 45.9999 1 1 6.605748879 1
26 46.001 1 46.0009 1 1 6.608899906 1
26 46.001 0 46.0009 0 1 6.608899906 1
26 46.001 0 46.0009 0 1 6.608899906 1
27 46.42 1 46.4199 1 1 6.604846759 1
27 46.42 0 46.4199 0 1 6.604846759 1
27 46.42 0 46.4199 0 1 6.604846759 1
28 47 1 46.9999 1 1 6.590758819 1
28 47 0 46.9999 1 1 6.590758819 1
28 47 0 46.9999 1 1 6.590758819 1
29 47.001 1 47.0009 1 1 6.596690846 1
29 47.001 0 47.0009 1 1 6.596690846 1
29 47.001 0 47.0009 0 1 6.596690846 1
30 47.17 1 47.1699 1 1 6.59121638 1
30 47.17 0 47.1699 0 1 6.59121638 1
30 47.17 0 47.1699 0 0 6.59121638 1
31 47.25 1 47.2499 1 1 6.588467867 1
31 47.25 0 47.2499 0 1 6.588467867 1
31 47.25 0 47.2499 1 0 6.588467867 1
32 47.251 1 47.2509 1 1 6.588926478 1
32 47.251 0 47.2509 0 0 6.588926478 1
32 47.251 0 47.2509 0 0 6.588926478 1
33 47.33 1 47.3299 1 1 6.584330881 1
33 47.33 0 47.3299 0 0 6.584330881 1
33 47.33 0 47.3299 0 0 6.584330881 1
34 47.42 1 47.4199 0 0 6.583870158 1
34 47.42 0 47.4199 1 1 6.583870158 1
34 47.42 0 47.4199 0 1 6.583870158 1
35 47.5 1 47.4999 1 0 6.578788142 1
35 47.5 0 47.4999 0 0 6.578788142 1
35 47.5 0 47.4999 1 1 6.578788142 1
36 47.83 1 47.8299 1 1 6.572748635 1
36 47.83 0 47.8299 0 0 6.572748635 1
36 47.83 0 47.8299 1 1 6.572748635 1
37 48 1 47.9999 1 1 6.570416 1
37 48 0 47.9999 1 1 6.570416 1
37 48 0 47.9999 0 0 6.570416 1
38 48.08 1 48.0799 1 1 6.569013803 1
38 48.08 0 48.0799 0 1 6.569013803 1
38 48.08 0 48.0799 0 1 6.569013803 1
39 48.17 1 48.1699 1 1 6.566203497 1
39 48.17 0 48.1699 0 0 6.566203497 1
39 48.17 0 48.1699 0 0 6.566203497 1
40 48.33 1 48.3299 1 1 6.558670504 1
40 48.33 0 48.3299 0 1 6.558670504 1
40 48.33 0 48.3299 1 0 6.558670504 1
41 48.5 1 48.4999 1 1 6.552032263 1
41 48.5 0 48.4999 1 1 6.552032263 1
41 48.5 0 48.4999 0 0 6.552032263 1
42 48.58 1 48.5799 1 1 6.548219103 1
42 48.58 0 48.5799 0 1 6.548219103 1
42 48.58 0 48.5799 0 1 6.548219103 1
43 48.67 1 48.6699 1 1 6.547263537 1
43 48.67 0 48.6699 0 0 6.547263537 1
43 48.67 0 48.6699 0 0 6.547263537 1
44 49.08 1 49.0799 1 1 6.534273217 1
44 49.08 0 49.0799 1 0 6.534273217 1
44 49.08 0 49.0799 0 0 6.534273217 1
45 49.081 1 49.0809 0 1 6.535241271 1
45 49.081 0 49.0809 0 0 6.535241271 1
45 49.081 0 49.0809 1 0 6.535241271 1
46 49.17 1 49.1699 1 1 6.531848973 1
46 49.17 0 49.1699 1 1 6.531848973 1
46 49.17 0 49.1699 0 1 6.531848973 1
47 49.25 1 49.2499 1 1 6.531848973 1
47 49.25 0 49.2499 1 1 6.531848973 1
47 49.25 0 49.2499 0 1 6.531848973 1
48 49.33 1 49.3299 1 1 6.528445128 1
48 49.33 0 49.3299 1 1 6.528445128 1
48 49.33 0 49.3299 0 1 6.528445128 1
49 49.67 1 49.6699 1 1 6.519638809 1
49 49.67 0 49.6699 1 1 6.519638809 1
49 49.67 0 49.6699 1 1 6.519638809 1
50 49.92 1 49.9199 1 1 6.513230111 1
50 49.92 0 49.9199 1 1 6.513230111 1
50 49.92 0 49.9199 1 1 6.513230111 1
51 50.08 1 50.0799 0 1 6.511249912 2
51 50.08 0 50.0799 0 1 6.511249912 2
51 50.08 0 50.0799 1 1 6.511249912 2
52 50.17 1 50.1699 0 1 6.509265785 2
52 50.17 0 50.1699 0 0 6.509265785 2
52 50.17 0 50.1699 1 1 6.509265785 2
53 50.171 1 50.1709 1 1 6.507277712 2
53 50.171 0 50.1709 0 0 6.507277712 2
53 50.171 0 50.1709 1 0 6.507277712 2
54 50.5 1 50.4999 1 0 6.492744758 2
54 50.5 0 50.4999 1 0 6.492744758 2
54 50.5 0 50.4999 0 0 6.492744758 2
55 50.58 1 50.5799 1 1 6.494257999 2
55 50.58 0 50.5799 1 1 6.494257999 2
55 50.58 0 50.5799 0 0 6.494257999 2
56 50.92 1 50.9199 0 1 6.482597538 2
56 50.92 0 50.9199 0 1 6.482597538 2
56 50.92 0 50.9199 0 1 6.482597538 2
57 51 1 50.9999 0 1 6.479021544 2
57 51 0 50.9999 0 1 6.479021544 2
57 51 0 50.9999 0 1 6.479021544 2
58 51.001 1 51.0009 1 1 6.479021544 2
58 51.001 0 51.0009 1 0 6.479021544 2
58 51.001 0 51.0009 0 0 6.479021544 2
59 51.08 1 51.0799 1 1 6.476459411 2
59 51.08 0 51.0799 1 0 6.476459411 2
59 51.08 0 51.0799 0 0 6.476459411 2
60 51.17 1 51.1699 1 0 6.475432717 2
60 51.17 0 51.1699 0 1 6.475432717 2
60 51.17 0 51.1699 0 0 6.475432717 2
61 51.25 1 51.2499 0 1 6.473890696 2
61 51.25 0 51.2499 1 1 6.473890696 2
61 51.25 0 51.2499 1 1 6.473890696 2
62 51.251 1 51.2509 1 0 6.474404967 2
62 51.251 0 51.2509 1 1 6.474404967 2
62 51.251 0 51.2509 1 1 6.474404967 2
63 51.33 1 51.3299 1 0 6.471315367 2
63 51.33 0 51.3299 1 1 6.471315367 2
63 51.33 0 51.3299 1 1 6.471315367 2
64 51.58 1 51.5799 1 1 6.457815681 2
64 51.58 0 51.5799 1 1 6.457815681 2
64 51.58 0 51.5799 1 1 6.457815681 2
65 51.67 1 51.6699 1 0 6.455198563 2
65 51.67 0 51.6699 1 1 6.455198563 2
65 51.67 0 51.6699 0 1 6.455198563 2
66 51.75 1 51.7499 1 0 6.452048954 2
66 51.75 0 51.7499 0 1 6.452048954 2
66 51.75 0 51.7499 0 1 6.452048954 2
67 51.83 1 51.8299 1 1 6.447833985 2
67 51.83 0 51.8299 0 1 6.447833985 2
67 51.83 0 51.8299 0 0 6.447833985 2
68 51.92 1 51.9199 0 1 6.442540166 2
68 51.92 0 51.9199 1 0 6.442540166 2
68 51.92 0 51.9199 0 0 6.442540166 2
69 52 1 51.9999 1 0 6.43775165 2
69 52 0 51.9999 0 0 6.43775165 2
69 52 0 51.9999 0 0 6.43775165 2
70 52.25 1 52.2499 0 1 6.428643631 2
70 52.25 0 52.2499 1 1 6.428643631 2
70 52.25 0 52.2499 0 1 6.428643631 2
71 52.42 1 52.4199 1 1 6.422164126 2
71 52.42 0 52.4199 0 1 6.422164126 2
71 52.42 0 52.4199 1 0 6.422164126 2
72 52.67 1 52.6699 1 0 6.40742989 2
72 52.67 0 52.6699 1 0 6.40742989 2
72 52.67 0 52.6699 0 0 6.40742989 2
73 52.75 1 52.7499 0 1 6.404677342 2
73 52.75 0 52.7499 0 1 6.404677342 2
73 52.75 0 52.7499 0 1 6.404677342 2
74 52.751 1 52.7509 0 1 6.404677342 2
74 52.751 0 52.7509 0 1 6.404677342 2
74 52.751 0 52.7509 0 1 6.404677342 2
75 53 1 52.9999 1 1 6.39804015 2
75 53 0 52.9999 1 1 6.39804015 2
75 53 0 52.9999 1 1 6.39804015 2
76 53.08 1 53.0799 0 0 6.393590754 2
76 53.08 0 53.0799 1 1 6.393590754 2
76 53.08 0 53.0799 0 1 6.393590754 2
77 53.25 1 53.2499 1 1 6.389121473 2
77 53.25 0 53.2499 0 1 6.389121473 2
77 53.25 0 53.2499 0 0 6.389121473 2
78 53.251 1 53.2509 1 1 6.390240667 2
78 53.251 0 53.2509 0 1 6.390240667 2
78 53.251 0 53.2509 1 0 6.390240667 2
79 53.252 1 53.2519 0 1 6.389681226 2
79 53.252 0 53.2519 0 1 6.389681226 2
79 53.252 0 53.2519 0 0 6.389681226 2
80 53.42 1 53.4199 1 1 6.384632127 2
80 53.42 0 53.4199 0 0 6.384632127 2
80 53.42 0 53.4199 0 0 6.384632127 2
81 54 1 53.9999 1 1 6.352629396 2
81 54 0 53.9999 1 1 6.352629396 2
81 54 0 53.9999 0 1 6.352629396 2
82 54.08 1 54.0799 1 0 6.34738921 2
82 54.08 0 54.0799 0 1 6.34738921 2
82 54.08 0 54.0799 0 1 6.34738921 2
83 54.25 1 54.2499 0 1 6.337415528 2
83 54.25 0 54.2499 0 0 6.337415528 2
83 54.25 0 54.2499 0 0 6.337415528 2
84 54.5 1 54.4999 1 1 6.323761412 2
84 54.5 0 54.4999 0 1 6.323761412 2
84 54.5 0 54.4999 0 1 6.323761412 2
85 54.58 1 54.5799 1 0 6.318968114 2
85 54.58 0 54.5799 1 1 6.318968114 2
85 54.58 0 54.5799 1 1 6.318968114 2
86 54.581 1 54.5809 1 0 6.317766191 2
86 54.581 0 54.5809 1 0 6.317766191 2
86 54.581 0 54.5809 1 1 6.317766191 2
87 54.67 1 54.6699 0 0 6.310524155 2
87 54.67 0 54.6699 1 0 6.310524155 2
87 54.67 0 54.6699 1 0 6.310524155 2
88 54.83 1 54.8299 1 0 6.305058001 2
88 54.83 0 54.8299 0 1 6.305058001 2
88 54.83 0 54.8299 0 1 6.305058001 2
89 54.92 1 54.9199 1 1 6.302008289 2
89 54.92 0 54.9199 0 0 6.302008289 2
89 54.92 0 54.9199 0 1 6.302008289 2
90 54.921 1 54.9209 1 1 6.300785795 2
90 54.921 0 54.9209 1 1 6.300785795 2
90 54.921 0 54.9209 1 1 6.300785795 2
91 55 1 54.9999 0 0 6.295266001 2
91 55 0 54.9999 1 1 6.295266001 2
91 55 0 54.9999 0 0 6.295266001 2
92 55.25 1 55.2499 1 1 6.283512077 2
92 55.25 0 55.2499 0 1 6.283512077 2
92 55.25 0 55.2499 1 0 6.283512077 2
93 55.251 1 55.2509 1 1 6.282266747 2
93 55.251 0 55.2509 0 1 6.282266747 2
93 55.251 0 55.2509 0 1 6.282266747 2
94 55.33 1 55.3299 1 0 6.27914662 2
94 55.33 0 55.3299 0 0 6.27914662 2
94 55.33 0 55.3299 1 0 6.27914662 2
95 55.5 1 55.4999 1 1 6.270988432 2
95 55.5 0 55.4999 0 1 6.270988432 2
95 55.5 0 55.4999 0 1 6.270988432 2
96 55.58 1 55.5799 1 1 6.267200549 2
96 55.58 0 55.5799 1 1 6.267200549 2
96 55.58 0 55.5799 0 0 6.267200549 2
97 55.67 1 55.6699 1 1 6.263398263 2
97 55.67 0 55.6699 0 1 6.263398263 2
97 55.67 0 55.6699 1 1 6.263398263 2
98 55.671 1 55.6709 1 1 6.262127614 2
98 55.671 0 55.6709 0 1 6.262127614 2
98 55.671 0 55.6709 0 1 6.262127614 2
99 55.672 1 55.6719 1 1 6.261491684 2
99 55.672 0 55.6719 1 1 6.261491684 2
99 55.672 0 55.6719 0 0 6.261491684 2
100 55.92 1 55.9199 1 1 6.251261416 2
100 55.92 0 55.9199 0 1 6.251261416 2
100 55.92 0 55.9199 0 1 6.251261416 2
101 55.921 1 55.9209 1 0 6.251261416 2
101 55.921 0 55.9209 1 0 6.251261416 2
101 55.921 0 55.9209 1 1 6.251261416 2
102 56.08 1 56.0799 0 1 6.246106765 2
102 56.08 0 56.0799 0 0 6.246106765 2
102 56.08 0 56.0799 1 1 6.246106765 2
103 56.081 1 56.0809 1 0 6.24481394 2
103 56.081 0 56.0809 1 1 6.24481394 2
103 56.081 0 56.0809 0 0 6.24481394 2
104 56.17 1 56.1699 1 1 6.239625861 2
104 56.17 0 56.1699 1 1 6.239625861 2
104 56.17 0 56.1699 1 0 6.239625861 2
105 56.171 1 56.1709 1 0 6.238324625 2
105 56.171 0 56.1709 1 0 6.238324625 2
105 56.171 0 56.1709 0 1 6.238324625 2
106 56.25 1 56.2499 1 0 6.234410726 2
106 56.25 0 56.2499 0 1 6.234410726 2
106 56.25 0 56.2499 0 1 6.234410726 2
107 56.33 1 56.3299 1 0 6.230481448 2
107 56.33 0 56.3299 1 1 6.230481448 2
107 56.33 0 56.3299 0 0 6.230481448 2
108 56.331 1 56.3309 0 0 6.2311374 2
108 56.331 0 56.3309 0 1 6.2311374 2
108 56.331 0 56.3309 1 1 6.2311374 2
109 56.42 1 56.4199 1 1 6.227853325 2
109 56.42 0 56.4199 1 1 6.227853325 2
109 56.42 0 56.4199 0 0 6.227853325 2
110 56.5 1 56.4999 1 1 6.21860012 2
110 56.5 0 56.4999 1 1 6.21860012 2
110 56.5 0 56.4999 1 1 6.21860012 2
111 56.501 1 56.5009 1 1 6.21860012 2
111 56.501 0 56.5009 0 1 6.21860012 2
111 56.501 0 56.5009 1 1 6.21860012 2
112 56.83 1 56.8299 1 1 6.207247744 2
112 56.83 0 56.8299 0 1 6.207247744 2
112 56.83 0 56.8299 0 0 6.207247744 2
113 56.92 1 56.9199 0 1 6.203884135 2
113 56.92 0 56.9199 1 1 6.203884135 2
113 56.92 0 56.9199 0 1 6.203884135 2
114 56.921 1 56.9209 1 1 6.205230937 2
114 56.921 0 56.9209 1 1 6.205230937 2
114 56.921 0 56.9209 1 1 6.205230937 2
115 57 1 56.9999 0 0 6.201185078 2
115 57 0 56.9999 0 1 6.201185078 2
115 57 0 56.9999 0 1 6.201185078 2
116 57.001 1 57.0009 0 1 6.200509174 2
116 57.001 0 57.0009 0 1 6.200509174 2
116 57.001 0 57.0009 1 1 6.200509174 2
117 57.33 1 57.3299 1 1 6.186208624 2
117 57.33 0 57.3299 1 1 6.186208624 2
117 57.33 0 57.3299 0 1 6.186208624 2
118 57.42 1 57.4199 1 1 6.181395964 2
118 57.42 0 57.4199 0 1 6.181395964 2
118 57.42 0 57.4199 0 1 6.181395964 2
119 57.5 1 57.4999 0 1 6.169610732 2
119 57.5 0 57.4999 0 0 6.169610732 2
119 57.5 0 57.4999 0 0 6.169610732 2
120 57.75 1 57.7499 1 0 6.159799862 2
120 57.75 0 57.7499 1 1 6.159799862 2
120 57.75 0 57.7499 0 1 6.159799862 2
121 57.92 1 57.9199 1 1 6.153441663 2
121 57.92 0 57.9199 0 1 6.153441663 2
121 57.92 0 57.9199 1 1 6.153441663 2
122 57.921 1 57.9209 1 1 6.153441663 2
122 57.921 0 57.9209 0 0 6.153441663 2
122 57.921 0 57.9209 0 0 6.153441663 2
123 58 1 57.9999 1 1 6.149180293 2
123 58 0 57.9999 1 1 6.149180293 2
123 58 0 57.9999 0 1 6.149180293 2
124 58.001 1 58.0009 1 1 6.149891784 2
124 58.001 0 58.0009 1 1 6.149891784 2
124 58.001 0 58.0009 0 0 6.149891784 2
125 58.002 1 58.0019 1 0 6.149180293 2
125 58.002 0 58.0019 0 0 6.149180293 2
125 58.002 0 58.0019 0 0 6.149180293 2
126 58.003 1 58.0029 1 1 6.148468296 2
126 58.003 0 58.0029 0 1 6.148468296 2
126 58.003 0 58.0029 0 0 6.148468296 2
127 58.08 1 58.0799 1 1 6.144185634 2
127 58.08 0 58.0799 1 1 6.144185634 2
127 58.08 0 58.0799 0 0 6.144185634 2
128 58.081 1 58.0809 0 1 6.144185634 2
128 58.081 0 58.0809 0 1 6.144185634 2
128 58.081 0 58.0809 1 1 6.144185634 2
129 58.17 1 58.1699 1 1 6.138446737 2
129 58.17 0 58.1699 1 1 6.138446737 2
129 58.17 0 58.1699 0 1 6.138446737 2
130 58.25 1 58.2499 1 1 6.135564891 2
130 58.25 0 58.2499 1 0 6.135564891 2
130 58.25 0 58.2499 1 0 6.135564891 2
131 58.251 1 58.2509 1 1 6.136286132 2
131 58.251 0 58.2509 1 1 6.136286132 2
131 58.251 0 58.2509 0 1 6.136286132 2
132 58.252 1 58.2519 1 1 6.135564891 2
132 58.252 0 58.2519 1 0 6.135564891 2
132 58.252 0 58.2519 1 0 6.135564891 2
133 58.42 1 58.4199 1 1 6.131226489 2
133 58.42 0 58.4199 1 1 6.131226489 2
133 58.42 0 58.4199 0 0 6.131226489 2
134 58.5 1 58.4999 1 1 6.115892125 2
134 58.5 0 58.4999 0 1 6.115892125 2
134 58.5 0 58.4999 1 0 6.115892125 2
135 58.58 1 58.5799 1 1 6.112944445 2
135 58.58 0 58.5799 0 1 6.112944445 2
135 58.58 0 58.5799 0 1 6.112944445 2
136 58.67 1 58.6699 1 0 6.109988049 2
136 58.67 0 58.6699 0 1 6.109988049 2
136 58.67 0 58.6699 0 1 6.109988049 2
137 58.671 1 58.6709 1 1 6.109247583 2
137 58.671 0 58.6709 0 1 6.109247583 2
137 58.671 0 58.6709 1 0 6.109247583 2
138 58.75 1 58.7499 1 1 6.104793232 2
138 58.75 0 58.7499 0 0 6.104793232 2
138 58.75 0 58.7499 1 1 6.104793232 2
139 58.751 1 58.7509 1 1 6.101812604 2
139 58.751 0 58.7509 0 0 6.101812604 2
139 58.751 0 58.7509 0 0 6.101812604 2
140 58.92 1 58.9199 1 1 6.094321933 2
140 58.92 0 58.9199 0 1 6.094321933 2
140 58.92 0 58.9199 0 1 6.094321933 2
141 59.17 1 59.1699 1 1 6.08221891 2
141 59.17 0 59.1699 0 1 6.08221891 2
141 59.17 0 59.1699 1 1 6.08221891 2
142 59.25 1 59.2499 1 0 6.078406477 2
142 59.25 0 59.2499 1 1 6.078406477 2
142 59.25 0 59.2499 1 1 6.078406477 2
143 59.42 1 59.4199 1 1 6.069196896 2
143 59.42 0 59.4199 0 0 6.069196896 2
143 59.42 0 59.4199 0 0 6.069196896 2
144 59.5 1 59.4999 0 1 6.055221513 2
144 59.5 0 59.4999 0 0 6.055221513 2
144 59.5 0 59.4999 0 0 6.055221513 2
145 59.67 1 59.6699 0 1 6.045005314 2
145 59.67 0 59.6699 0 0 6.045005314 2
145 59.67 0 59.6699 0 0 6.045005314 2
146 59.75 1 59.7499 1 0 6.040254711 2
146 59.75 0 59.7499 1 0 6.040254711 2
146 59.75 0 59.7499 1 1 6.040254711 2
147 59.83 1 59.8299 1 1 6.035481433 2
147 59.83 0 59.8299 1 1 6.035481433 2
147 59.83 0 59.8299 1 0 6.035481433 2
148 59.831 1 59.8309 1 1 6.033086222 2
148 59.831 0 59.8309 0 0 6.033086222 2
148 59.831 0 59.8309 0 1 6.033086222 2
149 59.832 1 59.8319 0 0 6.032286542 2
149 59.832 0 59.8319 0 1 6.032286542 2
149 59.832 0 59.8319 0 1 6.032286542 2
150 59.92 1 59.9199 1 1 6.029883657 2
150 59.92 0 59.9199 1 1 6.029883657 2
150 59.92 0 59.9199 0 0 6.029883657 2
151 60.5 1 60.4999 1 0 5.987289176 2
151 60.5 0 60.4999 0 1 5.987289176 2
151 60.5 0 60.4999 0 1 5.987289176 2
152 60.58 1 60.5799 1 1 5.980572106 2
152 60.58 0 60.5799 1 1 5.980572106 2
152 60.58 0 60.5799 1 0 5.980572106 2
153 60.67 1 60.6699 1 0 5.977196575 2
153 60.67 0 60.6699 0 1 5.977196575 2
153 60.67 0 60.6699 1 1 5.977196575 2
154 60.75 1 60.7499 0 1 5.972961076 2
154 60.75 0 60.7499 0 1 5.972961076 2
154 60.75 0 60.7499 1 0 5.972961076 2
155 60.83 1 60.8299 0 0 5.96870756 2
155 60.83 0 60.8299 0 1 5.96870756 2
155 60.83 0 60.8299 0 0 5.96870756 2
156 61.08 1 61.0799 1 0 5.954108761 2
156 61.08 0 61.0799 0 0 5.954108761 2
156 61.08 0 61.0799 1 0 5.954108761 2
157 61.17 1 61.1699 1 1 5.948034989 2
157 61.17 0 61.1699 0 1 5.948034989 2
157 61.17 0 61.1699 0 0 5.948034989 2
158 61.25 1 61.2499 1 1 5.940171253 2
158 61.25 0 61.2499 1 1 5.940171253 2
158 61.25 0 61.2499 1 0 5.940171253 2
159 61.251 1 61.2509 1 1 5.936656311 2
159 61.251 0 61.2509 1 0 5.936656311 2
159 61.251 0 61.2509 0 0 5.936656311 2
160 61.42 1 61.4199 1 1 5.92781452 2
160 61.42 0 61.4199 1 1 5.92781452 2
160 61.42 0 61.4199 1 1 5.92781452 2
161 61.5 1 61.4999 1 0 5.914403501 2
161 61.5 0 61.4999 1 1 5.914403501 2
161 61.5 0 61.4999 0 0 5.914403501 2
162 61.67 1 61.6699 1 1 5.896237697 2
162 61.67 0 61.6699 0 0 5.896237697 2
162 61.67 0 61.6699 0 1 5.896237697 2
163 61.75 1 61.7499 1 1 5.887029529 2
163 61.75 0 61.7499 0 1 5.887029529 2
163 61.75 0 61.7499 0 1 5.887029529 2
164 61.83 1 61.8299 1 1 5.877735782 2
164 61.83 0 61.8299 0 1 5.877735782 2
164 61.83 0 61.8299 0 0 5.877735782 2
165 61.92 1 61.9199 1 1 5.872117789 2
165 61.92 0 61.9199 1 1 5.872117789 2
165 61.92 0 61.9199 0 1 5.872117789 2
166 62 1 61.9999 1 1 5.864577697 2
166 62 0 61.9999 0 1 5.864577697 2
166 62 0 61.9999 0 1 5.864577697 2
167 62.08 1 62.0799 1 1 5.860786223 2
167 62.08 0 62.0799 0 0 5.860786223 2
167 62.08 0 62.0799 1 1 5.860786223 2
168 62.081 1 62.0809 0 1 5.85698032 2
168 62.081 0 62.0809 0 1 5.85698032 2
168 62.081 0 62.0809 0 1 5.85698032 2
169 62.17 1 62.1699 1 1 5.85220248 2
169 62.17 0 62.1699 0 0 5.85220248 2
169 62.17 0 62.1699 0 1 5.85220248 2
170 62.25 1 62.2499 1 1 5.848363703 2
170 62.25 0 62.2499 1 1 5.848363703 2
170 62.25 0 62.2499 0 1 5.848363703 2
171 62.251 1 62.2509 1 0 5.84547492 2
171 62.251 0 62.2509 0 1 5.84547492 2
171 62.251 0 62.2509 0 1 5.84547492 2
172 62.33 1 62.3299 1 1 5.84161018 2
172 62.33 0 62.3299 0 1 5.84161018 2
172 62.33 0 62.3299 1 1 5.84161018 2
173 62.58 1 62.5799 0 0 5.819093321 2
173 62.58 0 62.5799 1 0 5.819093321 2
173 62.58 0 62.5799 1 0 5.819093321 2
174 62.83 1 62.8299 1 1 5.798082043 2
174 62.83 0 62.8299 0 0 5.798082043 2
174 62.83 0 62.8299 1 1 5.798082043 2
175 62.831 1 62.8309 1 1 5.795044066 2
175 62.831 0 62.8309 0 1 5.795044066 2
175 62.831 0 62.8309 0 1 5.795044066 2
176 62.92 1 62.9199 1 1 5.788940283 2
176 62.92 0 62.9199 1 0 5.788940283 2
176 62.92 0 62.9199 1 1 5.788940283 2
177 63 1 62.9999 0 1 5.785874363 2
177 63 0 62.9999 1 0 5.785874363 2
177 63 0 62.9999 0 1 5.785874363 2
178 63.25 1 63.2499 1 1 5.766235489 2
178 63.25 0 63.2499 1 1 5.766235489 2
178 63.25 0 63.2499 0 1 5.766235489 2
179 63.251 1 63.2509 1 1 5.763099052 2
179 63.251 0 63.2509 0 1 5.763099052 2
179 63.251 0 63.2509 0 1 5.763099052 2
180 63.33 1 63.3299 1 1 5.759952746 2
180 63.33 0 63.3299 0 1 5.759952746 2
180 63.33 0 63.3299 1 1 5.759952746 2
181 63.42 1 63.4199 1 1 5.749392986 2
181 63.42 0 63.4199 0 0 5.749392986 2
181 63.42 0 63.4199 0 1 5.749392986 2
182 63.5 1 63.4999 0 1 5.737646989 2
182 63.5 0 63.4999 1 0 5.737646989 2
182 63.5 0 63.4999 0 1 5.737646989 2
183 63.67 1 63.6699 1 0 5.714832311 2
183 63.67 0 63.6699 0 0 5.714832311 2
183 63.67 0 63.6699 0 1 5.714832311 2
184 63.671 1 63.6709 0 1 5.713732806 2
184 63.671 0 63.6709 0 0 5.713732806 2
184 63.671 0 63.6709 1 0 5.713732806 2
185 64 1 63.9999 1 1 5.706002231 2
185 64 0 63.9999 0 1 5.706002231 2
185 64 0 63.9999 0 0 5.706002231 2
186 64.001 1 64.0009 1 0 5.702670746 2
186 64.001 0 64.0009 1 1 5.702670746 2
186 64.001 0 64.0009 1 0 5.702670746 2
187 64.33 1 64.3299 1 1 5.692609174 2
187 64.33 0 64.3299 0 1 5.692609174 2
187 64.33 0 64.3299 0 0 5.692609174 2
188 64.58 1 64.5799 1 1 5.676753802 2
188 64.58 0 64.5799 0 0 5.676753802 2
188 64.58 0 64.5799 1 1 5.676753802 2
189 64.581 1 64.5809 0 1 5.675611598 2
189 64.581 0 64.5809 0 0 5.675611598 2
189 64.581 0 64.5809 0 0 5.675611598 2
190 64.582 1 64.5819 1 0 5.674468087 2
190 64.582 0 64.5819 0 1 5.674468087 2
190 64.582 0 64.5819 0 1 5.674468087 2
191 64.583 1 64.5829 0 1 5.673323267 2
191 64.583 0 64.5829 1 1 5.673323267 2
191 64.583 0 64.5829 1 1 5.673323267 2
192 64.584 1 64.5839 1 1 5.672177135 2
192 64.584 0 64.5839 1 0 5.672177135 2
192 64.584 0 64.5839 1 0 5.672177135 2
193 64.67 1 64.6699 1 0 5.661802402 2
193 64.67 0 64.6699 0 1 5.661802402 2
193 64.67 0 64.6699 0 1 5.661802402 2
194 65.17 1 65.1699 0 1 5.619192406 2
194 65.17 0 65.1699 0 1 5.619192406 2
194 65.17 0 65.1699 1 1 5.619192406 2
195 65.25 1 65.2499 1 1 5.613128106 2
195 65.25 0 65.2499 0 0 5.613128106 2
195 65.25 0 65.2499 1 0 5.613128106 2
196 65.251 1 65.2509 0 1 5.607026806 2
196 65.251 0 65.2509 1 1 5.607026806 2
196 65.251 0 65.2509 1 1 5.607026806 2
197 65.33 1 65.3299 0 0 5.600888051 2
197 65.33 0 65.3299 0 1 5.600888051 2
197 65.33 0 65.3299 1 0 5.600888051 2
198 65.331 1 65.3309 0 1 5.598421959 2
198 65.331 0 65.3309 0 0 5.598421959 2
198 65.331 0 65.3309 0 1 5.598421959 2
199 65.67 1 65.6699 1 0 5.568344504 2
199 65.67 0 65.6699 1 0 5.568344504 2
199 65.67 0 65.6699 0 0 5.568344504 2
200 65.75 1 65.7499 1 1 5.555540232 2
200 65.75 0 65.7499 1 1 5.555540232 2
200 65.75 0 65.7499 1 1 5.555540232 2
201 66.08 1 66.0799 1 0 5.533389489 2
201 66.08 0 66.0799 0 0 5.533389489 2
201 66.08 0 66.0799 1 1 5.533389489 2
202 66.42 1 66.4199 1 1 5.512083756 2
202 66.42 0 66.4199 0 1 5.512083756 2
202 66.42 0 66.4199 1 1 5.512083756 2
203 66.421 1 66.4209 0 0 5.509388337 2
203 66.421 0 66.4209 0 0 5.509388337 2
203 66.421 0 66.4209 1 1 5.509388337 2
204 66.58 1 66.5799 1 0 5.490314189 2
204 66.58 0 66.5799 0 1 5.490314189 2
204 66.58 0 66.5799 1 1 5.490314189 2
205 66.67 1 66.6699 0 0 5.482026849 2
205 66.67 0 66.6699 1 0 5.482026849 2
205 66.67 0 66.6699 0 0 5.482026849 2
206 66.75 1 66.7499 1 1 5.476463552 2
206 66.75 0 66.7499 0 0 5.476463552 2
206 66.75 0 66.7499 0 1 5.476463552 2
207 67.25 1 67.2499 1 1 5.422008839 2
207 67.25 0 67.2499 1 0 5.422008839 2
207 67.25 0 67.2499 0 1 5.422008839 2
208 67.251 1 67.2509 1 1 5.413133041 2
208 67.251 0 67.2509 1 0 5.413133041 2
208 67.251 0 67.2509 0 0 5.413133041 2
209 67.33 1 67.3299 1 0 5.407171771 2
209 67.33 0 67.3299 1 1 5.407171771 2
209 67.33 0 67.3299 0 0 5.407171771 2
210 67.5 1 67.4999 1 1 5.373734006 2
210 67.5 0 67.4999 0 1 5.373734006 2
210 67.5 0 67.4999 0 1 5.373734006 2
211 67.83 1 67.8299 1 0 5.331107189 2
211 67.83 0 67.8299 0 0 5.331107189 2
211 67.83 0 67.8299 0 1 5.331107189 2
212 67.92 1 67.9199 1 1 5.321382639 2
212 67.92 0 67.9199 1 0 5.321382639 2
212 67.92 0 67.9199 0 0 5.321382639 2
213 68.33 1 68.3299 1 1 5.276412531 2
213 68.33 0 68.3299 1 1 5.276412531 2
213 68.33 0 68.3299 1 0 5.276412531 2
214 68.67 1 68.6699 0 0 5.23998179 2
214 68.67 0 68.6699 1 0 5.23998179 2
214 68.67 0 68.6699 1 1 5.23998179 2
215 68.75 1 68.7499 1 1 5.229324495 2
215 68.75 0 68.7499 0 0 5.229324495 2
215 68.75 0 68.7499 1 1 5.229324495 2
216 68.83 1 68.8299 1 0 5.214935758 2
216 68.83 0 68.8299 1 0 5.214935758 2
216 68.83 0 68.8299 0 1 5.214935758 2
217 69.25 1 69.2499 1 1 5.166688924 2
217 69.25 0 69.2499 0 0 5.166688924 2
217 69.25 0 69.2499 0 0 5.166688924 2
218 69.33 1 69.3299 1 0 5.157137753 2
218 69.33 0 69.3299 0 1 5.157137753 2
218 69.33 0 69.3299 1 1 5.157137753 2
219 69.331 1 69.3309 0 0 5.155216523 2
219 69.331 0 69.3309 1 1 5.155216523 2
219 69.331 0 69.3309 1 1 5.155216523 2
220 69.42 1 69.4199 1 0 5.149430586 2
220 69.42 0 69.4199 0 0 5.149430586 2
220 69.42 0 69.4199 1 0 5.149430586 2
221 69.5 1 69.4999 1 1 5.137757302 2
221 69.5 0 69.4999 1 1 5.137757302 2
221 69.5 0 69.4999 0 1 5.137757302 2
222 69.501 1 69.5009 0 1 5.131869159 2
222 69.501 0 69.5009 0 1 5.131869159 2
222 69.501 0 69.5009 1 1 5.131869159 2
223 69.67 1 69.6699 0 1 5.11599581 2
223 69.67 0 69.6699 1 1 5.11599581 2
223 69.67 0 69.6699 1 1 5.11599581 2
224 69.671 1 69.6709 1 1 5.107963638 2
224 69.671 0 69.6709 0 1 5.107963638 2
224 69.671 0 69.6709 0 0 5.107963638 2
225 70 1 69.9999 1 0 5.045573345 3
225 70 0 69.9999 1 0 5.045573345 3
225 70 0 69.9999 0 1 5.045573345 3
226 70.25 1 70.2499 1 1 4.999461993 3
226 70.25 0 70.2499 0 0 4.999461993 3
226 70.25 0 70.2499 1 0 4.999461993 3
227 70.33 1 70.3299 1 1 4.990432587 3
227 70.33 0 70.3299 1 0 4.990432587 3
227 70.33 0 70.3299 1 1 4.990432587 3
228 70.331 1 70.3309 0 0 4.988162438 3
228 70.331 0 70.3309 0 0 4.988162438 3
228 70.331 0 70.3309 1 1 4.988162438 3
229 70.5 1 70.4999 0 1 4.96284463 3
229 70.5 0 70.4999 1 0 4.96284463 3
229 70.5 0 70.4999 0 0 4.96284463 3
230 70.92 1 70.9199 1 0 4.905274778 3
230 70.92 0 70.9199 1 0 4.905274778 3
230 70.92 0 70.9199 0 0 4.905274778 3
231 71 1 70.9999 1 0 4.895349139 3
231 71 0 70.9999 0 0 4.895349139 3
231 71 0 70.9999 1 1 4.895349139 3
232 71.08 1 71.0799 1 1 4.882801923 3
232 71.08 0 71.0799 1 0 4.882801923 3
232 71.08 0 71.0799 0 1 4.882801923 3
233 71.17 1 71.1699 1 1 4.870095271 3
233 71.17 0 71.1699 0 0 4.870095271 3
233 71.17 0 71.1699 1 0 4.870095271 3
234 71.42 1 71.4199 1 1 4.849422701 3
234 71.42 0 71.4199 1 0 4.849422701 3
234 71.42 0 71.4199 0 0 4.849422701 3
235 71.421 1 71.4209 1 0 4.836281907 3
235 71.421 0 71.4209 1 0 4.836281907 3
235 71.421 0 71.4209 0 1 4.836281907 3
236 71.67 1 71.6699 0 1 4.767855768 3
236 71.67 0 71.6699 1 0 4.767855768 3
236 71.67 0 71.6699 0 0 4.767855768 3
237 71.75 1 71.7499 0 1 4.756459634 3
237 71.75 0 71.7499 1 0 4.756459634 3
237 71.75 0 71.7499 0 0 4.756459634 3
238 72.08 1 72.0799 1 0 4.703506087 3
238 72.08 0 72.0799 0 1 4.703506087 3
238 72.08 0 72.0799 0 1 4.703506087 3
239 72.081 1 72.0809 1 0 4.69440132 3
239 72.081 0 72.0809 0 1 4.69440132 3
239 72.081 0 72.0809 0 1 4.69440132 3
240 72.17 1 72.1699 1 0 4.682131227 3
240 72.17 0 72.1699 0 0 4.682131227 3
240 72.17 0 72.1699 0 1 4.682131227 3
241 72.33 1 72.3299 1 1 4.672828834 3
241 72.33 0 72.3299 0 1 4.672828834 3
241 72.33 0 72.3299 1 1 4.672828834 3
242 72.92 1 72.9199 1 0 4.560869927 3
242 72.92 0 72.9199 0 1 4.560869927 3
242 72.92 0 72.9199 1 1 4.560869927 3
243 73.08 1 73.0799 1 1 4.536177315 3
243 73.08 0 73.0799 1 0 4.536177315 3
243 73.08 0 73.0799 1 0 4.536177315 3
244 73.17 1 73.1699 0 1 4.529008825 3
244 73.17 0 73.1699 1 0 4.529008825 3
244 73.17 0 73.1699 1 0 4.529008825 3
245 73.25 1 73.2499 1 0 4.503506532 3
245 73.25 0 73.2499 1 1 4.503506532 3
245 73.25 0 73.2499 0 1 4.503506532 3
246 73.5 1 73.4999 1 1 4.454347296 3
246 73.5 0 73.4999 1 1 4.454347296 3
246 73.5 0 73.4999 1 0 4.454347296 3
247 73.58 1 73.5799 1 0 4.446565156 3
247 73.58 0 73.5799 1 0 4.446565156 3
247 73.58 0 73.5799 0 0 4.446565156 3
248 74.25 1 74.2499 1 1 4.335109715 3
248 74.25 0 74.2499 0 0 4.335109715 3
248 74.25 0 74.2499 1 0 4.335109715 3
249 74.251 1 74.2509 1 1 4.313033763 3
249 74.251 0 74.2509 1 1 4.313033763 3
249 74.251 0 74.2509 0 0 4.313033763 3
250 74.67 1 74.6699 1 0 4.204692619 3
250 74.67 0 74.6699 1 0 4.204692619 3
250 74.67 0 74.6699 0 0 4.204692619 3
251 75 1 74.9999 0 0 4.127134385 3
251 75 0 74.9999 1 0 4.127134385 3
251 75 0 74.9999 1 1 4.127134385 3
252 75.17 1 75.1699 1 1 4.071871706 3
252 75.17 0 75.1699 1 0 4.071871706 3
252 75.17 0 75.1699 0 1 4.071871706 3
253 76.67 1 76.6699 0 1 3.835861644 3
253 76.67 0 76.6699 1 1 3.835861644 3
253 76.67 0 76.6699 0 1 3.835861644 3
254 77.08 1 77.0799 0 0 3.79173684 3
254 77.08 0 77.0799 0 0 3.79173684 3
254 77.08 0 77.0799 0 0 3.79173684 3
255 79.17 1 79.1699 1 0 3.355735008 3
255 79.17 0 79.1699 0 0 3.355735008 3
255 79.17 0 79.1699 1 1 3.355735008 3
256 79.75 1 79.7499 1 1 3.258096538 3
256 79.75 0 79.7499 0 0 3.258096538 3
256 79.75 0 79.7499 0 0 3.258096538 3
257 80.17 1 80.1699 1 0 3.205452805 3
257 80.17 0 80.1699 1 0 3.205452805 3
257 80.17 0 80.1699 0 0 3.205452805 3
258 86.42 1 86.4199 1 1 1.540445041 3
258 86.42 0 86.4199 1 0 1.540445041 3
258 86.42 0 86.4199 0 0 1.540445041 3
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
