Dear all,
We are trying the code from the paper
(http://www.jstatsoft.org/v40/i01/) to build stars for each point on
the map (plyr package). It uses ggplot2 package and works fine for
ozone data but we can't get it work for our data (probably due to
misunderstanding of some issues).
Our data is SO2 time series for each point on the map with 0.25
degrees step. It has both NA, negative and positive values. Its length
is about 7 years with daily time step.
Here is a small sample for point (48, 28):
time_unix date value
1 1.096589e+12 2004-10-01 NA
2 1.096675e+12 2004-10-02 NA
3 1.096762e+12 2004-10-03 NA
4 1.096848e+12 2004-10-04 0.0061420719
5 1.096934e+12 2004-10-05 0.1289926320
6 1.097021e+12 2004-10-06 NA
7 1.097107e+12 2004-10-07 0.8269771338
8 1.097194e+12 2004-10-08 -0.2335753143
9 1.097280e+12 2004-10-09 -0.0868944377
10 1.097366e+12 2004-10-10 NA
11 1.097453e+12 2004-10-11 1.8495836258
12 1.097539e+12 2004-10-12 1.5986032486
13 1.097626e+12 2004-10-13 0.8445309997
14 1.097712e+12 2004-10-14 0.6054053903
15 1.097798e+12 2004-10-15 1.6727224588
16 1.097885e+12 2004-10-16 0.3150097430
17 1.097971e+12 2004-10-17 NA
18 1.098058e+12 2004-10-18 0.9932233691
19 1.098144e+12 2004-10-19 -0.7627518177
20 1.098230e+12 2004-10-20 -1.1208484173
21 1.098317e+12 2004-10-21 NA
22 1.098403e+12 2004-10-22 NA
23 1.098490e+12 2004-10-23 NA
24 1.098576e+12 2004-10-24 0.3474421799
25 1.098662e+12 2004-10-25 -0.1100406647
26 1.098749e+12 2004-10-26 -0.6024221778
27 1.098835e+12 2004-10-27 NA
28 1.098922e+12 2004-10-28 NA
29 1.099008e+12 2004-10-29 NA
30 1.099094e+12 2004-10-30 NA
31 1.099181e+12 2004-10-31 -0.7240343094
32 1.099267e+12 2004-11-01 NA
33 1.099354e+12 2004-11-02 NA
34 1.099440e+12 2004-11-03 NA
35 1.099526e+12 2004-11-04 NA
36 1.099613e+12 2004-11-05 NA
37 1.099699e+12 2004-11-06 NA
38 1.099786e+12 2004-11-07 NA
39 1.099872e+12 2004-11-08 NA
40 1.099958e+12 2004-11-09 NA
41 1.100045e+12 2004-11-10 NA
42 1.100131e+12 2004-11-11 NA
43 1.100218e+12 2004-11-12 NA
44 1.100304e+12 2004-11-13 NA
45 1.100390e+12 2004-11-14 NA
46 1.100477e+12 2004-11-15 NA
47 1.100563e+12 2004-11-16 NA
48 1.100650e+12 2004-11-17 NA
49 1.100736e+12 2004-11-18 NA
50 1.100822e+12 2004-11-19 NA
51 1.100909e+12 2004-11-20 NA
52 1.100995e+12 2004-11-21 NA
53 1.101082e+12 2004-11-22 NA
54 1.101168e+12 2004-11-23 NA
55 1.101254e+12 2004-11-24 NA
56 1.101341e+12 2004-11-25 NA
57 1.101427e+12 2004-11-26 NA
58 1.101514e+12 2004-11-27 NA
59 1.101600e+12 2004-11-28 NA
60 1.101686e+12 2004-11-29 NA
61 1.101773e+12 2004-11-30 NA
62 1.101859e+12 2004-12-01 NA
63 1.101946e+12 2004-12-02 NA
64 1.102032e+12 2004-12-03 NA
65 1.102118e+12 2004-12-04 NA
66 1.102205e+12 2004-12-05 NA
67 1.102291e+12 2004-12-06 NA
68 1.102378e+12 2004-12-07 NA
69 1.102464e+12 2004-12-08 NA
70 1.102550e+12 2004-12-09 NA
71 1.102637e+12 2004-12-10 NA
72 1.102723e+12 2004-12-11 NA
73 1.102810e+12 2004-12-12 NA
74 1.102896e+12 2004-12-13 NA
75 1.102982e+12 2004-12-14 NA
76 1.103069e+12 2004-12-15 NA
77 1.103155e+12 2004-12-16 NA
78 1.103242e+12 2004-12-17 NA
79 1.103328e+12 2004-12-18 NA
80 1.103414e+12 2004-12-19 NA
81 1.103501e+12 2004-12-20 NA
82 1.103587e+12 2004-12-21 NA
83 1.103674e+12 2004-12-22 NA
84 1.103760e+12 2004-12-23 NA
85 1.103846e+12 2004-12-24 NA
86 1.103933e+12 2004-12-25 NA
87 1.104019e+12 2004-12-26 NA
88 1.104106e+12 2004-12-27 NA
89 1.104192e+12 2004-12-28 NA
90 1.104278e+12 2004-12-29 NA
91 1.104365e+12 2004-12-30 NA
92 1.104451e+12 2004-12-31 NA
93 1.104538e+12 2005-01-01 NA
94 1.104624e+12 2005-01-02 NA
95 1.104710e+12 2005-01-03 NA
96 1.104797e+12 2005-01-04 NA
97 1.104883e+12 2005-01-05 NA
98 1.104970e+12 2005-01-06 NA
99 1.105056e+12 2005-01-07 NA
100 1.105142e+12 2005-01-08 NA
101 1.105229e+12 2005-01-09 NA
102 1.105315e+12 2005-01-10 NA
103 1.105402e+12 2005-01-11 NA
104 1.105488e+12 2005-01-12 NA
105 1.105574e+12 2005-01-13 NA
106 1.105661e+12 2005-01-14 NA
107 1.105747e+12 2005-01-15 NA
108 1.105834e+12 2005-01-16 NA
109 1.105920e+12 2005-01-17 NA
110 1.106006e+12 2005-01-18 NA
111 1.106093e+12 2005-01-19 NA
112 1.106179e+12 2005-01-20 NA
113 1.106266e+12 2005-01-21 NA
114 1.106352e+12 2005-01-22 NA
115 1.106438e+12 2005-01-23 NA
116 1.106525e+12 2005-01-24 NA
117 1.106611e+12 2005-01-25 NA
118 1.106698e+12 2005-01-26 NA
119 1.106784e+12 2005-01-27 NA
120 1.106870e+12 2005-01-28 NA
121 1.106957e+12 2005-01-29 NA
122 1.107043e+12 2005-01-30 NA
123 1.107130e+12 2005-01-31 NA
124 1.107216e+12 2005-02-01 NA
125 1.107302e+12 2005-02-02 NA
126 1.107389e+12 2005-02-03 NA
127 1.107475e+12 2005-02-04 NA
128 1.107562e+12 2005-02-05 NA
129 1.107648e+12 2005-02-06 NA
130 1.107734e+12 2005-02-07 NA
131 1.107821e+12 2005-02-08 NA
132 1.107907e+12 2005-02-09 NA
133 1.107994e+12 2005-02-10 NA
134 1.108080e+12 2005-02-11 NA
135 1.108166e+12 2005-02-12 NA
136 1.108253e+12 2005-02-13 NA
137 1.108339e+12 2005-02-14 NA
138 1.108426e+12 2005-02-15 NA
139 1.108512e+12 2005-02-16 NA
140 1.108598e+12 2005-02-17 NA
141 1.108685e+12 2005-02-18 NA
142 1.108771e+12 2005-02-19 NA
143 1.108858e+12 2005-02-20 NA
144 1.108944e+12 2005-02-21 NA
145 1.109030e+12 2005-02-22 NA
146 1.109117e+12 2005-02-23 NA
147 1.109203e+12 2005-02-24 NA
148 1.109290e+12 2005-02-25 NA
149 1.109376e+12 2005-02-26 NA
150 1.109462e+12 2005-02-27 NA
151 1.109549e+12 2005-02-28 NA
152 1.109635e+12 2005-03-01 NA
153 1.109722e+12 2005-03-02 NA
154 1.109808e+12 2005-03-03 NA
155 1.109894e+12 2005-03-04 NA
156 1.109981e+12 2005-03-05 NA
157 1.110067e+12 2005-03-06 NA
158 1.110154e+12 2005-03-07 NA
159 1.110240e+12 2005-03-08 NA
160 1.110326e+12 2005-03-09 NA
161 1.110413e+12 2005-03-10 NA
162 1.110499e+12 2005-03-11 NA
163 1.110586e+12 2005-03-12 NA
164 1.110672e+12 2005-03-13 0.7098249197
165 1.110758e+12 2005-03-14 -0.5577268004
166 1.110845e+12 2005-03-15 -0.1150340438
167 1.110931e+12 2005-03-16 -2.6746203899
168 1.111018e+12 2005-03-17 1.9455249310
169 1.111104e+12 2005-03-18 NA
170 1.111190e+12 2005-03-19 NA
171 1.111277e+12 2005-03-20 NA
172 1.111363e+12 2005-03-21 NA
173 1.111450e+12 2005-03-22 NA
174 1.111536e+12 2005-03-23 NA
175 1.111622e+12 2005-03-24 -1.4572575092
176 1.111709e+12 2005-03-25 NA
177 1.111795e+12 2005-03-26 NA
178 1.111882e+12 2005-03-27 -0.4886645079
179 1.111968e+12 2005-03-28 NA
180 1.112054e+12 2005-03-29 -0.2307986468
181 1.112141e+12 2005-03-30 NA
182 1.112227e+12 2005-03-31 -0.2024482042
183 1.112314e+12 2005-04-01 NA
184 1.112400e+12 2005-04-02 NA
185 1.112486e+12 2005-04-03 -0.9777982235
186 1.112573e+12 2005-04-04 -1.9230638742
187 1.112659e+12 2005-04-05 -0.3645920753
188 1.112746e+12 2005-04-06 1.0298686028
189 1.112832e+12 2005-04-07 0.5949411392
190 1.112918e+12 2005-04-08 -1.1224043369
191 1.113005e+12 2005-04-09 -1.4254966974
192 1.113091e+12 2005-04-10 NA
193 1.113178e+12 2005-04-11 NA
194 1.113264e+12 2005-04-12 NA
195 1.113350e+12 2005-04-13 NA
196 1.113437e+12 2005-04-14 NA
197 1.113523e+12 2005-04-15 NA
198 1.113610e+12 2005-04-16 NA
199 1.113696e+12 2005-04-17 -0.6095501184
200 1.113782e+12 2005-04-18 -0.4926427007
201 1.113869e+12 2005-04-19 -0.7748999000
202 1.113955e+12 2005-04-20 NA
203 1.114042e+12 2005-04-21 0.5438590050
204 1.114128e+12 2005-04-22 NA
205 1.114214e+12 2005-04-23 0.6878462434
206 1.114301e+12 2005-04-24 NA
207 1.114387e+12 2005-04-25 0.4751869142
208 1.114474e+12 2005-04-26 NA
209 1.114560e+12 2005-04-27 NA
210 1.114646e+12 2005-04-28 NA
211 1.114733e+12 2005-04-29 NA
212 1.114819e+12 2005-04-30 NA
213 1.114906e+12 2005-05-01 -0.9722039104
214 1.114992e+12 2005-05-02 NA
215 1.115078e+12 2005-05-03 -0.3842672706
216 1.115165e+12 2005-05-04 -0.1891176254
217 1.115251e+12 2005-05-05 NA
218 1.115338e+12 2005-05-06 NA
219 1.115424e+12 2005-05-07 NA
220 1.115510e+12 2005-05-08 -0.4297648370
221 1.115597e+12 2005-05-09 0.7431994677
222 1.115683e+12 2005-05-10 NA
223 1.115770e+12 2005-05-11 NA
224 1.115856e+12 2005-05-12 -0.3754873574
225 1.115942e+12 2005-05-13 1.0822590590
226 1.116029e+12 2005-05-14 0.1991832405
227 1.116115e+12 2005-05-15 -0.0029508672
228 1.116202e+12 2005-05-16 0.1759385169
229 1.116288e+12 2005-05-17 -0.7524182200
230 1.116374e+12 2005-05-18 0.0271098670
Here is the code wich forms input matrix to make_stars function:
transpose_stars = function(lat_b, long_b,lat_t,long_t,step){
matrix <- readSO2Data(lat_b, long_b) // read an arbitrary point for
calculating matrix size
nlat <- (lat_t-lat_b)/step+1
nlong <- (long_t-long_b)/step+1
ntime <- length(matrix$value)
list_lat<-list()
list_long<-list()
list_time<-list()
for( m in 1:nlat){
list_lat[m]<-toString(lat_b + step*(m-1))
}
for( m in 1:nlong){
list_long[m]<-toString(long_b + step*(m-1))
}
for( m in 1:ntime){
list_time[m]<-toString(m)
}
test<-array(data=NA,dim=c(nlat,nlong,ntime),list(lat = list_lat,long
list_long,time =list_time))
lat_loop<-lat_b
long_loop<-long_b
for(i in 1:nlat){
for(j in 1:nlong){
matrix <- readSO2Data(lat_loop+step*(i-1), long_loop+step*(j-1))
test[i,j,]<-matrix$value
}
}
return(test)
}
Slightly changed make_stars code (with na.rm = TRUE):
make_stars <- function(data, time, value) {
data[, c(time, value)] <- lapply(data[, c(time, value)], function(x) {
x <- as.numeric(x)
(x - min(x,na.rm = TRUE)) / diff(range(x,na.rm = TRUE))
})
ddply(data, .(lat, long), function(df) {
df$x <- df[, value] * cos(df[, time] * 2 * pi + pi / 2)
df$y <- df[, value] * sin(df[, time] * 2 * pi + pi / 2)
df[order(df[, time]), ]
})
}
The data after transpose_stars(sample of 7 dates)
, , time = 1
long
lat 28 28.25 28.5
48 NA NA NA
48.25 NA NA NA
48.5 NA NA NA
, , time = 2
long
lat 28 28.25 28.5
48 NA NA NA
48.25 NA NA NA
48.5 NA NA NA
, , time = 3
long
lat 28 28.25 28.5
48 NA NA NA
48.25 NA NA NA
48.5 NA NA NA
, , time = 4
long
lat 28 28.25 28.5
48 0.006142072 NA NA
48.25 0.986131787 0.9736934 NA
48.5 0.669782579 0.6609899 -0.8407977
, , time = 5
long
lat 28 28.25 28.5
48 0.1289926 0.1274767 -1.0643847
48.25 0.2148655 0.2121553 -0.9771701
48.5 -0.1101724 -1.2263418 -0.6803514
, , time = 6
long
lat 28 28.25 28.5
48 NA NA NA
48.25 NA NA NA
48.5 NA NA NA
, , time = 7
long
lat 28 28.25 28.5
48 0.8269771 0.90059173 0.65524501
48.25 1.1812921 -0.49642536 -0.48989806
48.5 -0.3733957 0.03384778 0.03337998
The data after applying melt()
lat long time value
1 48.00 28.00 1 NA
2 48.25 28.00 1 NA
3 48.50 28.00 1 NA
4 48.00 28.25 1 NA
5 48.25 28.25 1 NA
6 48.50 28.25 1 NA
7 48.00 28.50 1 NA
8 48.25 28.50 1 NA
9 48.50 28.50 1 NA
10 48.00 28.00 2 NA
11 48.25 28.00 2 NA
12 48.50 28.00 2 NA
13 48.00 28.25 2 NA
14 48.25 28.25 2 NA
15 48.50 28.25 2 NA
16 48.00 28.50 2 NA
17 48.25 28.50 2 NA
18 48.50 28.50 2 NA
19 48.00 28.00 3 NA
20 48.25 28.00 3 NA
21 48.50 28.00 3 NA
22 48.00 28.25 3 NA
23 48.25 28.25 3 NA
24 48.50 28.25 3 NA
25 48.00 28.50 3 NA
26 48.25 28.50 3 NA
27 48.50 28.50 3 NA
28 48.00 28.00 4 0.0061420719
29 48.25 28.00 4 0.9861317873
30 48.50 28.00 4 0.6697825789
31 48.00 28.25 4 NA
32 48.25 28.25 4 0.9736934304
33 48.50 28.25 4 0.6609899402
34 48.00 28.50 4 NA
35 48.25 28.50 4 NA
36 48.50 28.50 4 -0.8407977223
37 48.00 28.00 5 0.1289926320
38 48.25 28.00 5 0.2148654610
39 48.50 28.00 5 -0.1101723686
40 48.00 28.25 5 0.1274766922
41 48.25 28.25 5 0.2121553123
42 48.50 28.25 5 -1.2263418436
43 48.00 28.50 5 -1.0643846989
44 48.25 28.50 5 -0.9771701097
45 48.50 28.50 5 -0.6803514361
46 48.00 28.00 6 NA
47 48.25 28.00 6 NA
48 48.50 28.00 6 NA
49 48.00 28.25 6 NA
50 48.25 28.25 6 NA
51 48.50 28.25 6 NA
52 48.00 28.50 6 NA
53 48.25 28.50 6 NA
54 48.50 28.50 6 NA
55 48.00 28.00 7 0.8269771338
56 48.25 28.00 7 1.1812920570
57 48.50 28.00 7 -0.3733956814
58 48.00 28.25 7 0.9005917311
59 48.25 28.25 7 -0.4964253604
60 48.50 28.25 7 0.0338477790
61 48.00 28.50 7 0.6552450061
62 48.25 28.50 7 -0.4898980558
63 48.50 28.50 7 0.0333799757
64 48.00 28.00 8 -0.2335753143
65 48.25 28.00 8 -0.6291310191
66 48.50 28.00 8 1.3873713017
67 48.00 28.25 8 0.1986568868
68 48.25 28.25 8 -0.6211956143
69 48.50 28.25 8 1.3691585064
70 48.00 28.50 8 0.3291694224
71 48.25 28.50 8 0.3399876952
72 48.50 28.50 8 1.8036993742
73 48.00 28.00 9 -0.0868944377
74 48.25 28.00 9 -1.5597105026
75 48.50 28.00 9 -0.6437677741
76 48.00 28.25 9 -0.0858732387
77 48.25 28.25 9 -1.5400375128
78 48.50 28.25 9 -0.6353166699
79 48.00 28.50 9 0.6432881355
80 48.25 28.50 9 2.5510394573
81 48.50 28.50 9 -0.6265360117
82 48.00 28.00 10 NA
83 48.25 28.00 10 NA
84 48.50 28.00 10 NA
85 48.00 28.25 10 NA
86 48.25 28.25 10 NA
87 48.50 28.25 10 NA
88 48.00 28.50 10 NA
89 48.25 28.50 10 NA
90 48.50 28.50 10 NA
91 48.00 28.00 11 1.8495836258
92 48.25 28.00 11 1.1024228334
93 48.50 28.00 11 0.3401820362
94 48.00 28.25 11 -0.3195888400
95 48.25 28.25 11 0.6344406009
96 48.50 28.25 11 0.4953756630
97 48.00 28.50 11 -0.3156543970
98 48.25 28.50 11 0.6260985732
99 48.50 28.50 11 0.4885291457
100 48.00 28.00 12 1.5986032486
101 48.25 28.00 12 -0.9401375055
102 48.50 28.00 12 -0.2272775322
103 48.00 28.25 12 0.0605854318
104 48.25 28.25 12 -0.0968705639
105 48.50 28.25 12 -1.7634104490
106 48.00 28.50 12 0.6993049383
107 48.25 28.50 12 -0.0955968499
108 48.50 28.50 12 -1.7390385866
109 48.00 28.00 13 0.8445309997
110 48.25 28.00 13 0.0624825843
111 48.50 28.00 13 1.3104938269
112 48.00 28.25 13 2.9179391861
113 48.25 28.25 13 1.6861394644
114 48.50 28.25 13 1.2932901382
115 48.00 28.50 13 0.9382953048
116 48.25 28.50 13 0.3322718143
117 48.50 28.50 13 1.7414759398
118 48.00 28.00 14 0.6054053903
119 48.25 28.00 14 0.4198371172
120 48.50 28.00 14 -0.4944095612
121 48.00 28.25 14 0.5982905626
122 48.25 28.25 14 0.4145416021
123 48.50 28.25 14 0.2337413281
124 48.00 28.50 14 0.5909250975
125 48.25 28.50 14 0.5290067792
126 48.50 28.50 14 0.2871938348
127 48.00 28.00 15 1.6727224588
128 48.25 28.00 15 1.8313188553
129 48.50 28.00 15 NA
130 48.00 28.25 15 -0.1610381603
131 48.25 28.25 15 1.8082200289
132 48.50 28.25 15 NA
133 48.00 28.50 15 -0.1590556353
134 48.25 28.50 15 NA
135 48.50 28.50 15 NA
136 48.00 28.00 16 0.3150097430
137 48.25 28.00 16 0.8319899440
138 48.50 28.00 16 NA
139 48.00 28.25 16 -0.4835641086
After make_stars
lat long time value x y
1 48.00 28.00 0.0000000000 NA NA NA
2 48.00 28.00 0.0003547357 NA NA NA
3 48.00 28.00 0.0007094714 NA NA NA
4 48.00 28.00 0.0010642072 0.241280057 -0.0016133338 0.2412746630
5 48.00 28.00 0.0014189429 0.248470868 -0.0022152080 0.2484609927
6 48.00 28.00 0.0017736786 NA NA NA
7 48.00 28.00 0.0021284143 0.289325988 -0.0038691052 0.2893001162
8 48.00 28.00 0.0024831501 0.227248681 -0.0035454110 0.2272210229
9 48.00 28.00 0.0028378858 0.235834352 -0.0042049306 0.2357968618
10 48.00 28.00 0.0031926215 NA NA NA
11 48.00 28.00 0.0035473572 0.349182204 -0.0077821740 0.3490954729
12 48.00 28.00 0.0039020929 0.334491572 -0.0082001000 0.3343910435
13 48.00 28.00 0.0042568287 0.290353468 -0.0077649967 0.2902496188
14 48.00 28.00 0.0046115644 0.276356731 -0.0080064025 0.2762407289
15 48.00 28.00 0.0049663001 0.338829991 -0.0105711976 0.3386650446
16 48.00 28.00 0.0053210358 0.259359005 -0.0086695484 0.2592140671
17 48.00 28.00 0.0056757716 NA NA NA
18 48.00 28.00 0.0060305073 0.299056877 -0.0113287915 0.2988422229
19 48.00 28.00 0.0063852430 0.196274398 -0.0078723508 0.1961164590
20 48.00 28.00 0.0067399787 0.175313933 -0.0074220694 0.1751567530
21 48.00 28.00 0.0070947144 NA NA NA
22 48.00 28.00 0.0074494502 NA NA NA
23 48.00 28.00 0.0078041859 NA NA NA
24 48.00 28.00 0.0081589216 0.261257373 -0.0133872368 0.2609141559
25 48.00 28.00 0.0085136573 0.234479534 -0.0125370060 0.2341441336
26 48.00 28.00 0.0088683930 0.205658971 -0.0114537499 0.2053397773
27 48.00 28.00 0.0092231288 NA NA NA
28 48.00 28.00 0.0095778645 NA NA NA
29 48.00 28.00 0.0099326002 NA NA NA
30 48.00 28.00 0.0102873359 NA NA NA
31 48.00 28.00 0.0106420717 0.198540650 -0.0132657501 0.1980969699
32 48.00 28.00 0.0109968074 NA NA NA
33 48.00 28.00 0.0113515431 NA NA NA
34 48.00 28.00 0.0117062788 NA NA NA
35 48.00 28.00 0.0120610145 NA NA NA
36 48.00 28.00 0.0124157503 NA NA NA
37 48.00 28.00 0.0127704860 NA NA NA
38 48.00 28.00 0.0131252217 NA NA NA
39 48.00 28.00 0.0134799574 NA NA NA
40 48.00 28.00 0.0138346932 NA NA NA
41 48.00 28.00 0.0141894289 NA NA NA
42 48.00 28.00 0.0145441646 NA NA NA
43 48.00 28.00 0.0148989003 NA NA NA
44 48.00 28.00 0.0152536360 NA NA NA
45 48.00 28.00 0.0156083718 NA NA NA
46 48.00 28.00 0.0159631075 NA NA NA
47 48.00 28.00 0.0163178432 NA NA NA
48 48.00 28.00 0.0166725789 NA NA NA
49 48.00 28.00 0.0170273147 NA NA NA
50 48.00 28.00 0.0173820504 NA NA NA
51 48.00 28.00 0.0177367861 NA NA NA
52 48.00 28.00 0.0180915218 NA NA NA
53 48.00 28.00 0.0184462575 NA NA NA
54 48.00 28.00 0.0188009933 NA NA NA
55 48.00 28.00 0.0191557290 NA NA NA
56 48.00 28.00 0.0195104647 NA NA NA
57 48.00 28.00 0.0198652004 NA NA NA
58 48.00 28.00 0.0202199361 NA NA NA
59 48.00 28.00 0.0205746719 NA NA NA
60 48.00 28.00 0.0209294076 NA NA NA
61 48.00 28.00 0.0212841433 NA NA NA
62 48.00 28.00 0.0216388790 NA NA NA
63 48.00 28.00 0.0219936148 NA NA NA
64 48.00 28.00 0.0223483505 NA NA NA
65 48.00 28.00 0.0227030862 NA NA NA
66 48.00 28.00 0.0230578219 NA NA NA
67 48.00 28.00 0.0234125576 NA NA NA
68 48.00 28.00 0.0237672934 NA NA NA
69 48.00 28.00 0.0241220291 NA NA NA
70 48.00 28.00 0.0244767648 NA NA NA
71 48.00 28.00 0.0248315005 NA NA NA
72 48.00 28.00 0.0251862363 NA NA NA
73 48.00 28.00 0.0255409720 NA NA NA
74 48.00 28.00 0.0258957077 NA NA NA
75 48.00 28.00 0.0262504434 NA NA NA
76 48.00 28.00 0.0266051791 NA NA NA
77 48.00 28.00 0.0269599149 NA NA NA
78 48.00 28.00 0.0273146506 NA NA NA
79 48.00 28.00 0.0276693863 NA NA NA
80 48.00 28.00 0.0280241220 NA NA NA
81 48.00 28.00 0.0283788578 NA NA NA
82 48.00 28.00 0.0287335935 NA NA NA
.........................................
10314 48.25 28.00 0.6573252927 NA NA NA
10315 48.25 28.00 0.6576800284 NA NA NA
10316 48.25 28.00 0.6580347641 NA NA NA
10317 48.25 28.00 0.6583894998 0.153376685 0.1286619830 -0.0834895293
10318 48.25 28.00 0.6587442355 NA NA NA
10319 48.25 28.00 0.6590989713 NA NA NA
10320 48.25 28.00 0.6594537070 NA NA NA
10321 48.25 28.00 0.6598084427 NA NA NA
10322 48.25 28.00 0.6601631784 NA NA NA
10323 48.25 28.00 0.6605179142 NA NA NA
10324 48.25 28.00 0.6608726499 NA NA NA
10325 48.25 28.00 0.6612273856 NA NA NA
10326 48.25 28.00 0.6615821213 NA NA NA
10327 48.25 28.00 0.6619368570 NA NA NA
10328 48.25 28.00 0.6622915928 NA NA NA
10329 48.25 28.00 0.6626463285 NA NA NA
10330 48.25 28.00 0.6630010642 NA NA NA
10331 48.25 28.00 0.6633557999 NA NA NA
10332 48.25 28.00 0.6637105357 NA NA NA
10333 48.25 28.00 0.6640652714 NA NA NA
10334 48.25 28.00 0.6644200071 NA NA NA
10335 48.25 28.00 0.6647747428 NA NA NA
10336 48.25 28.00 0.6651294785 NA NA NA
10337 48.25 28.00 0.6654842143 NA NA NA
10338 48.25 28.00 0.6658389500 NA NA NA
10339 48.25 28.00 0.6661936857 NA NA NA
10340 48.25 28.00 0.6665484214 NA NA NA
10341 48.25 28.00 0.6669031571 NA NA NA
10342 48.25 28.00 0.6672578929 NA NA NA
10343 48.25 28.00 0.6676126286 NA NA NA
10344 48.25 28.00 0.6679673643 NA NA NA
10345 48.25 28.00 0.6683221000 NA NA NA
10346 48.25 28.00 0.6686768358 NA NA NA
10347 48.25 28.00 0.6690315715 NA NA NA
10348 48.25 28.00 0.6693863072 NA NA NA
10349 48.25 28.00 0.6697410429 NA NA NA
10350 48.25 28.00 0.6700957786 NA NA NA
10351 48.25 28.00 0.6704505144 NA NA NA
10352 48.25 28.00 0.6708052501 NA NA NA
10353 48.25 28.00 0.6711599858 NA NA NA
10354 48.25 28.00 0.6715147215 NA NA NA
10355 48.25 28.00 0.6718694573 NA NA NA
10356 48.25 28.00 0.6722241930 NA NA NA
10357 48.25 28.00 0.6725789287 NA NA NA
10358 48.25 28.00 0.6729336644 NA NA NA
10359 48.25 28.00 0.6732884001 NA NA NA
10360 48.25 28.00 0.6736431359 NA NA NA
10361 48.25 28.00 0.6739978716 NA NA NA
10362 48.25 28.00 0.6743526073 NA NA NA
10363 48.25 28.00 0.6747073430 NA NA NA
10364 48.25 28.00 0.6750620788 NA NA NA
10365 48.25 28.00 0.6754168145 NA NA NA
10366 48.25 28.00 0.6757715502 NA NA NA
10367 48.25 28.00 0.6761262859 NA NA NA
10368 48.25 28.00 0.6764810216 NA NA NA
10369 48.25 28.00 0.6768357574 NA NA NA
10370 48.25 28.00 0.6771904931 NA NA NA
10371 48.25 28.00 0.6775452288 NA NA NA
10372 48.25 28.00 0.6778999645 NA NA NA
10373 48.25 28.00 0.6782547002 NA NA NA
10374 48.25 28.00 0.6786094360 NA NA NA
10375 48.25 28.00 0.6789641717 NA NA NA
10376 48.25 28.00 0.6793189074 NA NA NA
10377 48.25 28.00 0.6796736431 NA NA NA
10378 48.25 28.00 0.6800283789 NA NA NA
10379 48.25 28.00 0.6803831146 NA NA NA
10380 48.25 28.00 0.6807378503 NA NA NA
10381 48.25 28.00 0.6810925860 NA NA NA
10382 48.25 28.00 0.6814473217 NA NA NA
10383 48.25 28.00 0.6818020575 NA NA NA
10384 48.25 28.00 0.6821567932 NA NA NA
10385 48.25 28.00 0.6825115289 NA NA NA
10386 48.25 28.00 0.6828662646 NA NA NA
10387 48.25 28.00 0.6832210004 NA NA NA
10388 48.25 28.00 0.6835757361 NA NA NA
10389 48.25 28.00 0.6839304718 NA NA NA
10390 48.25 28.00 0.6842852075 NA NA NA
10391 48.25 28.00 0.6846399432 NA NA NA
10392 48.25 28.00 0.6849946790 NA NA NA
10393 48.25 28.00 0.6853494147 NA NA NA
10394 48.25 28.00 0.6857041504 NA NA NA
10395 48.25 28.00 0.6860588861 NA NA NA
10396 48.25 28.00 0.6864136219 NA NA NA
10397 48.25 28.00 0.6867683576 NA NA NA
10398 48.25 28.00 0.6871230933 NA NA NA
10399 48.25 28.00 0.6874778290 NA NA NA
10400 48.25 28.00 0.6878325647 NA NA NA
10401 48.25 28.00 0.6881873005 NA NA NA
10402 48.25 28.00 0.6885420362 NA NA NA
10403 48.25 28.00 0.6888967719 NA NA NA
10404 48.25 28.00 0.6892515076 NA NA NA
10405 48.25 28.00 0.6896062433 NA NA NA
10406 48.25 28.00 0.6899609791 NA NA NA
10407 48.25 28.00 0.6903157148 NA NA NA
10408 48.25 28.00 0.6906704505 NA NA NA
10409 48.25 28.00 0.6910251862 NA NA NA
10410 48.25 28.00 0.6913799220 NA NA NA
10411 48.25 28.00 0.6917346577 NA NA NA
10412 48.25 28.00 0.6920893934 NA NA NA
10413 48.25 28.00 0.6924441291 NA NA NA
10414 48.25 28.00 0.6927988648 NA NA NA
10415 48.25 28.00 0.6931536006 NA NA NA
10416 48.25 28.00 0.6935083363 NA NA NA
10417 48.25 28.00 0.6938630720 NA NA NA
10418 48.25 28.00 0.6942178077 NA NA NA
10419 48.25 28.00 0.6945725435 NA NA NA
10420 48.25 28.00 0.6949272792 NA NA NA
10421 48.25 28.00 0.6952820149 NA NA NA
10422 48.25 28.00 0.6956367506 NA NA NA
10423 48.25 28.00 0.6959914863 NA NA NA
10424 48.25 28.00 0.6963462221 NA NA NA
10425 48.25 28.00 0.6967009578 NA NA NA
10426 48.25 28.00 0.6970556935 NA NA NA
10427 48.25 28.00 0.6974104292 NA NA NA
10428 48.25 28.00 0.6977651650 NA NA NA
10429 48.25 28.00 0.6981199007 NA NA NA
10430 48.25 28.00 0.6984746364 NA NA NA
10431 48.25 28.00 0.6988293721 NA NA NA
10432 48.25 28.00 0.6991841078 NA NA NA
10433 48.25 28.00 0.6995388436 NA NA NA
10434 48.25 28.00 0.6998935793 NA NA NA
10435 48.25 28.00 0.7002483150 NA NA NA
10436 48.25 28.00 0.7006030507 NA NA NA
10437 48.25 28.00 0.7009577864 NA NA NA
10438 48.25 28.00 0.7013125222 NA NA NA
10439 48.25 28.00 0.7016672579 NA NA NA
10440 48.25 28.00 0.7020219936 NA NA NA
10441 48.25 28.00 0.7023767293 NA NA NA
10442 48.25 28.00 0.7027314651 NA NA NA
10443 48.25 28.00 0.7030862008 NA NA NA
10444 48.25 28.00 0.7034409365 NA NA NA
10445 48.25 28.00 0.7037956722 NA NA NA
10446 48.25 28.00 0.7041504079 NA NA NA
10447 48.25 28.00 0.7045051437 0.227387721 0.2181606880 -0.0641177820
10448 48.25 28.00 0.7048598794 NA NA NA
10449 48.25 28.00 0.7052146151 NA NA NA
10450 48.25 28.00 0.7055693508 NA NA NA
10451 48.25 28.00 0.7059240866 NA NA NA
10452 48.25 28.00 0.7062788223 NA NA NA
10453 48.25 28.00 0.7066335580 NA NA NA
10454 48.25 28.00 0.7069882937 NA NA NA
10455 48.25 28.00 0.7073430294 NA NA NA
10456 48.25 28.00 0.7076977652 0.217451042 0.2098151531 -0.0571187985
10457 48.25 28.00 0.7080525009 NA NA NA
10458 48.25 28.00 0.7084072366 0.242139028 0.2339174158 -0.0625615796
10459 48.25 28.00 0.7087619723 0.166997248 0.1614227761 -0.0427874772
10460 48.25 28.00 0.7091167081 0.253136578 0.2448306734 -0.0643122731
10461 48.25 28.00 0.7094714438 NA NA NA
10462 48.25 28.00 0.7098261795 NA NA NA
10463 48.25 28.00 0.7101809152 0.187984340 0.1821314900 -0.0465428011
10464 48.25 28.00 0.7105356509 NA NA NA
10465 48.25 28.00 0.7108903867 NA NA NA
10466 48.25 28.00 0.7112451224 NA NA NA
10467 48.25 28.00 0.7115998581 0.146712962 0.1424632838 -0.0350557550
10468 48.25 28.00 0.7119545938 0.192607598 0.1871306509 -0.0456048931
10469 48.25 28.00 0.7123093295 NA NA NA
10470 48.25 28.00 0.7126640653 NA NA NA
10471 48.25 28.00 0.7130188010 NA NA NA
10472 48.25 28.00 0.7133735367 0.076736324 0.0747132903 -0.0175039340
10473 48.25 28.00 0.7137282724 NA NA NA
10474 48.25 28.00 0.7140830082 NA NA NA
10475 48.25 28.00 0.7144377439 0.343848538 0.3353004722 -0.0761932435
10476 48.25 28.00 0.7147924796 NA NA NA
10477 48.25 28.00 0.7151472153 NA NA NA
10478 48.25 28.00 0.7155019510 NA NA NA
10479 48.25 28.00 0.7158566868 NA NA NA
10480 48.25 28.00 0.7162114225 NA NA NA
10481 48.25 28.00 0.7165661582 NA NA NA
10482 48.25 28.00 0.7169208939 NA NA NA
10483 48.25 28.00 0.7172756297 NA NA NA
10484 48.25 28.00 0.7176303654 NA NA NA
10485 48.25 28.00 0.7179851011 NA NA NA
10486 48.25 28.00 0.7183398368 0.180105287 0.1765534760 -0.0355919177
10487 48.25 28.00 0.7186945725 0.235603990 0.2310609038 -0.0460445305
10488 48.25 28.00 0.7190493083 NA NA NA
10489 48.25 28.00 0.7194040440 NA NA NA
10490 48.25 28.00 0.7197587797 0.165559861 0.1625801344 -0.0312692718
10491 48.25 28.00 0.7201135154 0.200346696 0.1968247306 -0.0374008537
10492 48.25 28.00 0.7204682512 NA NA NA
10493 48.25 28.00 0.7208229869 NA NA NA
10494 48.25 28.00 0.7211777226 NA NA NA
10495 48.25 28.00 0.7215324583 0.146792124 0.1444501918 -0.0261164628
10496 48.25 28.00 0.7218871940 NA NA NA
10497 48.25 28.00 0.7222419298 NA NA NA
10498 48.25 28.00 0.7225966655 NA NA NA
10499 48.25 28.00 0.7229514012 0.213383214 0.2103090006 -0.0360904436
10500 48.25 28.00 0.7233061369 0.222296905 0.2191775297 -0.0371096299
10501 48.25 28.00 0.7236608726 NA NA NA
10502 48.25 28.00 0.7240156084 0.168157496 0.1659213190 -0.0273323866
10503 48.25 28.00 0.7243703441 NA NA NA
10504 48.25 28.00 0.7247250798 0.156438489 0.1544699680 -0.0247392356
10505 48.25 28.00 0.7250798155 NA NA NA
10506 48.25 28.00 0.7254345513 NA NA NA
10507 48.25 28.00 0.7257892870 0.186324323 0.1841726470 -0.0282345410
10508 48.25 28.00 0.7261440227 NA NA NA
10509 48.25 28.00 0.7264987584 0.204793234 0.2025646052 -0.0301305379
10510 48.25 28.00 0.7268534941 NA NA NA
10511 48.25 28.00 0.7272082299 0.226793849 0.2244723167 -0.0323670999
10512 48.25 28.00 0.7275629656 NA NA NA
10513 48.25 28.00 0.7279177013 0.231039458 0.2288191789 -0.0319533203
10514 48.25 28.00 0.7282724370 NA NA NA
10515 48.25 28.00 0.7286271728 NA NA NA
10516 48.25 28.00 0.7289819085 NA NA NA
10517 48.25 28.00 0.7293366442 NA NA NA
10518 48.25 28.00 0.7296913799 NA NA NA
10519 48.25 28.00 0.7300461156 NA NA NA
10520 48.25 28.00 0.7304008514 0.242148703 0.2403149645 -0.0297441120
10521 48.25 28.00 0.7307555871 NA NA NA
10522 48.25 28.00 0.7311103228 NA NA NA
10523 48.25 28.00 0.7314650585 0.188514360 0.1872374345 -0.0219045023
10524 48.25 28.00 0.7318197943 0.219057301 0.2176296807 -0.0249684458
10525 48.25 28.00 0.7321745300 NA NA NA
10526 48.25 28.00 0.7325292657 NA NA NA
10527 48.25 28.00 0.7328840014 0.285932831 0.2842809467 -0.0306908328
.........................................................................................................
What can be wrong?
Thank you.
Antonio
-------------- next part --------------
A non-text attachment was scrubbed...
Name: screenshot.png
Type: image/png
Size: 6508 bytes
Desc: not available
URL:
<https://stat.ethz.ch/pipermail/r-help/attachments/20121003/9ea48757/attachment-0002.png>
