Last data update: 2014.03.03

R: Classic gamma regression. Identity link for the mean
gammahetero2R Documentation

Classic gamma regression. Identity link for the mean

Description

Performs the Classic Gamma Regression for joint modeling of mean and shape parameters.

Usage

gammahetero2(formula1, formula2)

Arguments

formula1

object of class formula. It describes yi and xi for the mean equation of the gamma regression.

formula2

object of class formula. It describes zi for the shape equation of the gamma regression.

Details

The classic gamma regression allow the joint modeling of mean and shape parameters of a gamma distributed variable, as is proposed in Cepeda (2001), using the Fisher Scoring algorithm, with log link for the mean and log link for the shape.

Value

object of class Gammareg with the following:

X

object of class matrix, with the variables for modelling the mean.

Z

object of class matrix, with the variables for modelling the shape.

beta

object of class matrix with the estimated coefficients of beta.

gamma

object of class matrix with the estimated coefficients of gamma.

ICB

object of class matrix with the estimated confidence intervals of beta.

ICG

object of class matrix with the estimated confidence intervals of gamma.

CovarianceMatrixbeta

object of class matrix with the estimated covariances of beta.

CovarianceMatrixgamma

object of class matrix with the estimated covariances of gamma.

AIC

the AIC criteria

iteration

numbers of iterations to convergence

convergence

value of convergence obtained

Author(s)

Martha Corrales martha.corrales@usa.edu.co Edilberto Cepeda-Cuervo ecepedac@unal.edu.co,

References

1. Cepeda-Cuervo, E. (2001). Modelagem da variabilidade em modelos lineares generalizados. Unpublished Ph.D. tesis. Instituto de Matem<c3><a1>ticas. Universidade Federal do R<c3><ad>o do Janeiro. //http://www.docentes.unal.edu.co/ecepedac/docs/MODELAGEM20DA20VARIABILIDADE.pdf. http://www.bdigital.unal.edu.co/9394/. 2. McCullagh, P. and Nelder, N.A. (1989). Generalized Linear Models. Second Edition. Chapman and Hall.

Examples


# Simulation Example

X1 <- rep(1,500)
X2 <- runif(500,0,30)
X3 <- runif(500,0,15)
X4 <- runif(500,10,20)
mui <- 15 + 2*X2 + 3*X3
alphai <- exp(0.2 + 0.1*X2 + 0.3*X4)
Y <- rgamma(500,shape=alphai,scale=mui/alphai)
X <- cbind(X1,X2,X3)
Z <- cbind(X1,X2,X4)
formula.mean= Y~X2+X3
formula.shape= ~X2+X4
a=gammahetero2(formula.mean,formula.shape)
a

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

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(Gammareg)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Gammareg/gammahetero2.Rd_%03d_medium.png", width=480, height=480)
> ### Name: gammahetero2
> ### Title: Classic gamma regression. Identity link for the mean
> ### Aliases: gammahetero2
> ### Keywords: Gamma regression Fisher Scoring Joint modeling Mean and shape
> ###   parameters
> 
> ### ** Examples
> 
> 
> # Simulation Example
> 
> X1 <- rep(1,500)
> X2 <- runif(500,0,30)
> X3 <- runif(500,0,15)
> X4 <- runif(500,10,20)
> mui <- 15 + 2*X2 + 3*X3
> alphai <- exp(0.2 + 0.1*X2 + 0.3*X4)
> Y <- rgamma(500,shape=alphai,scale=mui/alphai)
> X <- cbind(X1,X2,X3)
> Z <- cbind(X1,X2,X4)
> formula.mean= Y~X2+X3
> formula.shape= ~X2+X4
> a=gammahetero2(formula.mean,formula.shape)
> a
$X
    (Intercept)          X2          X3
1             1 24.72253110  3.17756913
2             1  7.07405104  4.94910154
3             1 10.00790352  8.19931691
4             1 10.12412773  4.96111449
5             1  6.89208160 14.56202674
6             1 25.60260931 11.08233996
7             1 21.99736667 10.11866940
8             1 17.88023686  0.50860617
9             1 23.04712733  9.90414722
10            1 10.12713857  9.05841584
11            1 16.13933552  2.16983055
12            1 13.27626317  4.54434514
13            1 13.18377552  5.24641080
14            1 24.81757191  3.48331838
15            1  4.39448632  5.85993372
16            1 27.01520659  6.74568447
17            1  2.83342048  5.56282125
18            1 26.90295449  8.02228663
19            1 13.86241719 13.36493096
20            1 16.61213602 14.25865321
21            1  2.54739666  5.74126596
22            1  5.88260275  0.52135271
23            1 25.67845593 10.09646401
24            1 24.21961809  9.45740145
25            1 10.04279952 12.94263834
26            1  9.76154584  8.96705287
27            1 11.79945016 11.61041088
28            1 10.93303695 13.69649905
29            1  3.64479174  4.12583416
30            1  0.25151455  7.58590367
31            1 27.13902405  9.43099465
32            1  0.52369982  3.85917518
33            1 28.92391022 10.14455920
34            1 15.62674368 13.77887961
35            1 21.55948666  2.46444030
36            1 19.86926021  8.88371680
37            1  3.97427840  1.07032557
38            1  2.63540732  7.26909313
39            1 13.73258693  5.80199512
40            1 14.98342639 11.47877364
41            1 14.29279030 13.75053190
42            1 20.88979448 12.21630076
43            1  4.35563548 13.90801317
44            1  3.90902523  8.60517203
45            1  4.61900427  4.28138781
46            1 24.93680947  8.47613175
47            1 28.88759418  8.23591551
48            1  2.02154807  6.63556091
49            1 24.83720761  5.88153364
50            1  0.16396840 14.67120895
51            1 13.83883427 12.12652579
52            1 27.57210989 11.06608160
53            1  1.49722709 10.47351608
54            1  2.40985300  4.61915852
55            1 21.49295525 13.36754497
56            1 22.69841801  8.98673603
57            1  1.94895925  4.39678649
58            1  2.23141362  9.69273811
59            1 22.66253099 13.15617039
60            1 23.37023410  2.99646092
61            1  6.35717608 13.73950246
62            1 19.59995912  5.64931051
63            1 17.24968369 13.67550411
64            1 15.38818212  1.69369614
65            1 21.22613921  8.83442353
66            1  6.77979921  5.30556647
67            1  7.92420222  7.41452066
68            1 10.03525541 10.03362786
69            1  8.27169803  0.61895201
70            1  3.61965260  5.26730978
71            1 12.92176917 10.94960087
72            1 23.71847446  5.65432297
73            1 25.39190743 13.02445371
74            1  0.62090801 12.21240921
75            1  6.31442562  8.87919034
76            1 26.53439480  6.58539460
77            1 27.48904197  2.64648968
78            1 17.80098016 13.12192080
79            1  4.29562906 11.67173546
80            1 16.54324720  9.66006150
81            1 10.73209742  7.99956690
82            1  7.13659343  4.45237717
83            1 12.31524947  7.84185008
84            1  6.70797131  1.44426655
85            1 27.43820701  8.11948307
86            1 13.39230055  2.72098058
87            1 21.57409098 14.93269603
88            1  3.94527385  1.29563467
89            1  2.20709723  4.95552912
90            1  8.12248744 10.44547238
91            1 19.29064296  8.23809615
92            1 20.90471363  2.71779816
93            1 11.91917336  0.31193964
94            1 14.46105857 12.75200849
95            1  1.38115130  2.64204312
96            1 28.33883988 14.53022478
97            1 20.21450771  9.20533000
98            1  6.70517870  6.87401759
99            1 27.08232331 14.00824778
100           1 19.21269606  4.07866304
101           1 26.70268477  1.31359926
102           1  1.16627432  9.07832254
103           1 27.96514377  9.24446513
104           1  0.21500076 12.91213485
105           1 10.75766448  6.58085706
106           1 13.12050846  2.89449959
107           1  1.57355832  2.87163618
108           1  9.29538527  0.08099745
109           1 26.26807471  2.00364447
110           1  9.74017067  5.85966463
111           1 16.70221499  2.91836270
112           1 13.11648876  8.90903003
113           1 15.82924896  6.18227976
114           1 18.59954804 12.89855523
115           1  6.45501269 12.51171266
116           1 27.84218557  0.27991190
117           1 17.10159975  3.84501457
118           1 24.30868549  9.63824051
119           1 17.01677156  4.72231891
120           1  3.40334528  2.54029597
121           1  5.53478618  4.44744935
122           1 18.77708562 12.23892968
123           1 22.67428813  0.89466205
124           1 16.72377637  9.17452747
125           1  9.77678758  7.77388646
126           1  5.67670837  5.55530851
127           1  0.07211418  2.20111944
128           1 16.05733008 14.48350125
129           1  1.43336098 11.94014373
130           1  4.60081269  6.99719315
131           1  9.73013312  6.72932751
132           1 10.47236488  8.57846342
133           1  9.96437280  8.53608201
134           1 23.35217732  1.03566693
135           1 18.87619826  4.96637337
136           1 12.23189945  3.63654893
137           1 14.86156901  4.35220496
138           1 23.63052885  9.73273101
139           1 16.71979667  0.60836318
140           1 14.93658144 12.31537338
141           1 28.56726173 12.16439913
142           1 24.94012931  1.06723211
143           1 15.49415249  5.84961341
144           1 16.97706906 10.65166940
145           1 28.96205420 12.17898349
146           1 15.71221253  7.71161008
147           1 25.29777259  0.33712323
148           1 23.61929812  2.29550365
149           1  9.33856281  0.17663234
150           1 29.19288259  1.23073316
151           1 28.01400196  8.28669582
152           1 11.81503752  2.98427460
153           1  5.60530061  2.98944565
154           1 28.48947077 10.90896973
155           1  4.85282763  4.68274883
156           1 22.34695001 14.54715399
157           1 24.14828588 11.42307769
158           1 17.64847732 10.15308689
159           1 27.08401588  0.87261470
160           1 28.51078386  4.33107962
161           1  5.92444717  3.03471336
162           1 12.72590685 11.95983257
163           1  3.93794781 14.57516124
164           1 15.49459361  2.23772660
165           1  8.89701352 10.18968401
166           1  3.72436841  5.93302440
167           1  2.41288047  0.72074082
168           1  4.31144123  4.06240629
169           1 21.46084698 10.35579315
170           1  3.34014048  5.76298767
171           1 22.28174626  4.68270538
172           1  0.07971427  0.40371464
173           1 13.66240869  6.63654114
174           1  3.37509154  7.51783147
175           1 12.09906204 12.30483700
176           1 13.09438448  1.66164186
177           1 18.38770032  2.68269597
178           1  0.90526420  2.34427300
179           1 24.53641264  3.15504559
180           1 11.07248885  2.15695786
181           1 29.75123278  9.50873039
182           1  3.51813120  4.92506926
183           1 11.31896202  7.59320197
184           1 24.99027881 10.21371104
185           1  9.23333831 14.37266533
186           1  4.31053132 11.44940239
187           1  9.88480875  6.85257759
188           1  9.28301391  4.73624435
189           1  6.72554405  0.29056148
190           1 27.48191994 10.78560056
191           1 20.11630363 11.75296232
192           1 26.19192549  2.32025230
193           1  7.17281753 10.21699657
194           1 29.06524045  3.97117311
195           1 11.41594983  3.71145479
196           1 21.79984578 12.62826322
197           1 23.32687551 14.23318555
198           1 12.91897079  4.50529402
199           1 12.76246575  1.47058612
200           1 25.11776769  6.52134905
201           1 11.88029100 10.61272869
202           1 17.88826080  0.99177650
203           1  1.45431655 11.13678549
204           1 21.62140878  4.87419768
205           1 28.00591232 10.24168771
206           1 28.44356806  9.28473471
207           1 19.38723746  9.33924228
208           1 19.61220003  5.56913084
209           1 20.05806120 12.27606201
210           1  0.99825994  2.92284539
211           1 29.38576008  3.45816928
212           1  3.19759092  4.03503390
213           1  3.74973696  3.45189958
214           1 24.36117251 14.92968958
215           1  4.64522211  2.49755448
216           1 12.89546262  8.45050221
217           1 28.65181496  8.68883238
218           1  5.58822893  3.69651087
219           1  3.00535957 11.01761193
220           1 23.48511413 12.49502304
221           1  2.25061891  2.71109115
222           1  4.83245644  9.73293687
223           1  5.30237749 10.44391424
224           1 13.08378044  8.35534685
225           1  1.79869444  6.38626599
226           1 26.93214454  9.48034368
227           1 24.80540332  2.84677512
228           1 17.16293783  6.74693326
229           1 28.74116243  4.00902387
230           1 20.19843590  1.70742153
231           1  9.55075190  5.82777785
232           1  4.87268739  9.46636986
233           1  0.49033069 11.71773266
234           1  5.01403851  5.97963751
235           1  5.40694085 14.45666718
236           1 23.54595847  1.32773008
237           1  3.72915963  1.77115489
238           1 10.45373394  1.74792602
239           1 16.54957034 10.80645100
240           1 26.72166242  0.29593645
241           1  0.88368159 14.48641400
242           1 16.73458605  6.15067296
243           1 11.30916021 14.35379131
244           1  3.23337893  4.47202599
245           1 15.03776601 10.61764326
246           1 10.07554060  0.28843453
247           1  8.80227552  4.44765071
248           1 15.55532217  8.26161488
249           1  4.66681670 10.52738857
250           1  6.52839452  8.40378586
251           1 29.38683630  6.57358926
252           1 21.12527108 11.11895079
253           1  0.67814983 14.13817482
254           1 11.55289834  3.70659754
255           1  0.60880536  6.60979498
256           1 25.85585239  9.30498582
257           1 29.56353621  7.71236590
258           1  8.37897586  9.35238370
259           1  3.04939166  7.34689210
260           1 16.27215765 14.16034338
261           1  5.63052282 14.27705494
262           1 13.88824330  4.64389004
263           1 25.50527879  9.02132935
264           1 17.69220178  0.68361167
265           1  4.35289280 12.00012198
266           1 15.68959943  8.77847162
267           1 12.47080107 12.59590900
268           1 24.02576832  4.26339576
269           1 20.70828730 14.98221126
270           1 11.07329531  0.88323591
271           1  9.27819083 13.00552077
272           1  7.51417221  2.71435300
273           1 16.49200589  9.31680765
274           1 17.75953191  2.18687940
275           1 27.99685839  5.98548116
276           1 17.45847064 12.30914582
277           1 21.59845255  7.71864010
278           1 27.14801322 14.32558087
279           1 22.25875234  6.52020349
280           1 26.68913641  1.31345824
281           1 27.22712232  8.41727299
282           1 22.02145531 12.65289778
283           1  1.78460702 14.32971816
284           1 19.31797493  6.40640841
285           1 10.88016083 13.47664873
286           1 12.59904403  7.19203773
287           1 18.77755239  2.33098168
288           1 26.11253844 10.81969284
289           1  7.85159575 14.85792159
290           1 26.79609704 12.57109349
291           1 24.07107668  8.02307130
292           1  1.70475374 13.98404160
293           1 28.67198238  3.91298946
294           1  9.39763114  5.86147252
295           1 20.88695939  3.11279334
296           1 27.76709385  5.64156905
297           1 24.66185275 11.40024816
298           1 10.31414258 13.04854255
299           1  8.71189392  0.44590361
300           1 26.24108496  7.43228538
301           1 14.15107340  9.59112350
302           1 10.09628222  2.13013814
303           1  0.15107819  5.55878078
304           1 19.99346992 14.24584041
305           1 23.91534724  8.01224280
306           1 18.14463506 14.10827067
307           1 18.62603526 12.68939494
308           1 18.52428061  5.03223120
309           1 15.99289888  4.66958368
310           1 11.67586131 13.27891125
311           1  4.55759491  0.72701366
312           1 18.37493862  5.72008486
313           1  7.41999096  7.13332071
314           1  6.25341196  3.45715163
315           1 13.57254742 10.38290986
316           1 18.11996279  0.03025505
317           1 14.70339363 11.78330565
318           1  5.92862839  5.48925999
319           1 28.91499694  9.39482240
320           1 16.80856196 10.46889102
321           1 21.93594280 11.34724045
322           1 15.36431154  9.87077114
323           1 11.84119534 10.42811748
324           1 19.52206433  7.10867310
325           1 20.47725979  9.56632072
326           1 12.35378132  3.76903325
327           1 14.66797197 13.42878386
328           1 25.97003023 12.07425491
329           1 22.33673230 13.36949973
330           1 21.60996387 12.86034709
331           1 16.48966789  6.84034275
332           1 21.23278299  4.81869968
333           1 23.10720964  8.56001595
334           1 20.95316447 11.76963020
335           1  4.43269350 10.84755115
336           1 15.06014995 14.93322180
337           1  7.15298616  8.59213375
338           1 13.97470719  1.74504979
339           1  8.44645983 11.73836993
340           1  5.15987466 11.81352250
341           1 14.43019241  3.35383386
342           1 11.21752386  8.91908712
343           1 14.38105728 14.59014616
344           1 28.42061440 11.35938733
345           1 19.50704519  7.65565850
346           1 27.71453034  2.66760315
347           1 13.89465578 12.08540430
348           1  0.38070746  8.77947082
349           1  5.22947607 11.64295658
350           1 24.83918875 13.28929354
351           1 26.43600933  4.68466518
352           1 19.22990807  2.79865589
353           1 22.29846352 14.98370244
354           1  9.05562675 10.82459138
355           1 21.65174157 13.21567907
356           1 18.67198503  7.16163390
357           1 10.26918695  3.81359876
358           1 22.57927887  0.95725822
359           1 13.67168896 10.78917263
360           1  4.50336697  8.26094576
361           1 23.97654619  0.37970431
362           1  1.22466082  6.00578571
363           1 24.45052168  9.66202742
364           1 23.54317752 11.15029749
365           1 16.94828406  7.74655102
366           1 14.65426841 10.13885858
367           1 14.61753216 11.96208493
368           1  5.46446643 12.65795462
369           1 14.91135502 11.62732365
370           1 27.63798337 13.12877602
371           1  6.26957035  8.18210706
372           1 28.01163680 11.41976730
373           1 29.53891478  5.76101677
374           1  6.82935687 14.46819808
375           1 26.19689218  6.50537717
376           1 26.94746240  3.72147965
377           1 22.83205192 13.01368579
378           1 10.15505623 10.00042914
379           1 19.22787069 14.61646023
380           1  4.87280410  8.31492064
381           1 24.15230116  3.20422879
382           1 16.33196251 13.56532889
383           1 12.32337110  6.78906753
384           1  0.50099012 11.98541906
385           1 13.73262028  2.09504980
386           1 23.84278004  1.54251339
387           1 27.24253460 13.64825137
388           1  9.60615029 11.90121837
389           1  5.61683530 12.73069681
390           1  4.87567028 11.21175991
391           1 15.81467470  4.49690170
392           1  4.43847307 12.36121430
393           1  4.45156549 11.56038563
394           1 14.44912300  8.99867677
395           1 16.41580652  0.55750398
396           1 26.45224677  8.33331870
397           1 12.01957168  0.87585998
398           1 27.73371991  1.11151017
399           1 19.55383019 11.20261962
400           1  2.84149921  6.71448084
401           1  4.57131299  5.67445595
402           1 12.17223248  6.73773403
403           1 15.57346312  6.87685686
404           1 13.82520655  8.16437110
405           1 10.12103918 14.96589966
406           1  9.53955410  5.31456511
407           1 25.60160339 11.87825520
408           1 11.07221988 13.48140130
409           1  2.31305166  7.15585699
410           1  6.13828668  1.85168831
411           1  9.31020336 12.80583341
412           1 12.80527603  5.61381118
413           1 10.65806824  6.92626990
414           1 24.80908586  5.59281519
415           1  3.98190353  8.37612707
416           1  8.32006980 12.97262608
417           1  9.69285802  7.79907558
418           1 28.96168556 12.40023894
419           1  1.39309075  9.92596089
420           1  3.99216511  3.14111403
421           1 25.65617071 14.54858344
422           1 29.48714362 11.41436542
423           1  8.83724497  1.12145677
424           1 24.25193456  3.56878950
425           1 10.35582889 12.44622927
426           1 25.59441816  5.71182627
427           1 13.98030250  9.56751245
428           1  7.54372286 13.10038536
429           1  8.08380258 10.98667977
430           1 21.21541270 13.85511902
431           1 10.30630505  2.55104290
432           1 24.69932938  5.54495459
433           1 24.55247132  7.33624593
434           1 23.39587597  5.52771749
435           1 12.58292401 13.94711231
436           1 20.40068790  5.57562720
437           1  8.74694078  2.70569006
438           1 18.28908734  4.91113054
439           1  4.00198726  2.71955424
440           1  1.05375289 11.91755556
441           1  2.87580305  4.54080996
442           1 19.35738263 14.28297825
443           1 16.30702590  4.18066339
444           1 25.39569848 11.50130469
445           1 29.73890583  0.18902658
446           1 10.90074735  6.68435660
447           1 21.16904110  4.25306113
448           1 21.59592701  2.13074893
449           1  3.29415409 14.03188485
450           1 21.19532043  5.64125938
451           1  5.35415160 13.96893765
452           1 26.08873240  5.51464496
453           1  6.02810107  1.54902114
454           1 12.44526061  3.99849458
455           1  3.20169472  5.34536859
456           1 22.09032357  5.18461678
457           1 28.97981264  8.29068351
458           1 11.58388639 13.42926796
459           1  2.60665853 14.88982124
460           1 28.83005932  8.57835823
461           1 28.78123133  6.48516816
462           1 10.50209803 11.62877546
463           1 20.30417501  1.42502082
464           1 17.36412058  2.23657974
465           1  7.55637879 11.96641325
466           1 25.37770704 12.11422487
467           1 14.04302486  3.81716435
468           1  6.18051535  8.74866773
469           1 12.46610570  2.25319544
470           1 24.15661255  7.07480285
471           1 13.68010196 10.97855800
472           1 27.52172696  7.77756838
473           1 22.41514356  0.81859051
474           1  7.29576522  4.77842357
475           1 22.67068102  2.63326235
476           1  6.98522558 12.86073055
477           1 18.75794888  3.24543731
478           1 15.84404469  1.86385120
479           1 13.24406296 11.95902743
480           1 23.79571911  4.55730224
481           1 21.47306963  1.89963567
482           1 10.33275858  6.23884265
483           1 21.40716831  4.90055819
484           1 14.87419252  0.05640801
485           1  4.72184242 14.13570285
486           1  5.70793379  9.33000354
487           1  5.36796084 14.65925941
488           1  4.08048956  9.24160928
489           1 21.57885778  4.50060224
490           1  4.93713136  1.33061654
491           1 23.08441229 12.68901398
492           1  3.57821164 12.88585716
493           1 17.59118190  6.98036631
494           1 15.71197667  1.77413311
495           1 13.13536860  4.51592736
496           1 15.44981388 10.17295322
497           1 12.95604119 13.94139470
498           1 11.19623518 13.35003465
499           1 15.04995548  7.69060813
500           1 14.16441970  6.99815379
attr(,"assign")
[1] 0 1 2

$Z
    (Intercept)          X2       X4
1             1 24.72253110 16.19027
2             1  7.07405104 18.42628
3             1 10.00790352 19.36589
4             1 10.12412773 16.32179
5             1  6.89208160 15.68601
6             1 25.60260931 10.76560
7             1 21.99736667 15.15421
8             1 17.88023686 19.44353
9             1 23.04712733 10.04723
10            1 10.12713857 18.89696
11            1 16.13933552 13.35606
12            1 13.27626317 16.64784
13            1 13.18377552 14.29273
14            1 24.81757191 16.53257
15            1  4.39448632 18.07197
16            1 27.01520659 14.10879
17            1  2.83342048 13.10117
18            1 26.90295449 19.60441
19            1 13.86241719 10.98553
20            1 16.61213602 12.32830
21            1  2.54739666 17.22927
22            1  5.88260275 15.99635
23            1 25.67845593 15.49630
24            1 24.21961809 18.18733
25            1 10.04279952 19.55958
26            1  9.76154584 17.99382
27            1 11.79945016 15.20294
28            1 10.93303695 18.10999
29            1  3.64479174 11.45313
30            1  0.25151455 14.38223
31            1 27.13902405 14.59115
32            1  0.52369982 15.25214
33            1 28.92391022 13.47154
34            1 15.62674368 14.20057
35            1 21.55948666 14.38443
36            1 19.86926021 12.55566
37            1  3.97427840 15.16774
38            1  2.63540732 12.42377
39            1 13.73258693 14.42726
40            1 14.98342639 14.57043
41            1 14.29279030 14.35044
42            1 20.88979448 19.50229
43            1  4.35563548 18.19820
44            1  3.90902523 14.20005
45            1  4.61900427 10.21901
46            1 24.93680947 16.59840
47            1 28.88759418 14.26173
48            1  2.02154807 11.18099
49            1 24.83720761 15.97045
50            1  0.16396840 15.27805
51            1 13.83883427 14.90798
52            1 27.57210989 13.94518
53            1  1.49722709 17.25840
54            1  2.40985300 10.55961
55            1 21.49295525 19.30409
56            1 22.69841801 16.45718
57            1  1.94895925 10.37129
58            1  2.23141362 17.32626
59            1 22.66253099 18.12682
60            1 23.37023410 17.64295
61            1  6.35717608 16.93980
62            1 19.59995912 10.61766
63            1 17.24968369 11.12465
64            1 15.38818212 17.39332
65            1 21.22613921 11.34797
66            1  6.77979921 16.77022
67            1  7.92420222 16.19426
68            1 10.03525541 12.99085
69            1  8.27169803 11.27939
70            1  3.61965260 19.26367
71            1 12.92176917 11.93635
72            1 23.71847446 11.75331
73            1 25.39190743 15.96168
74            1  0.62090801 12.02004
75            1  6.31442562 18.60314
76            1 26.53439480 13.55693
77            1 27.48904197 10.00585
78            1 17.80098016 15.78815
79            1  4.29562906 14.33619
80            1 16.54324720 15.42256
81            1 10.73209742 13.03951
82            1  7.13659343 15.92314
83            1 12.31524947 16.98066
84            1  6.70797131 19.79251
85            1 27.43820701 18.75661
86            1 13.39230055 10.97334
87            1 21.57409098 10.91009
88            1  3.94527385 18.35099
89            1  2.20709723 13.15072
90            1  8.12248744 15.07096
91            1 19.29064296 11.82010
92            1 20.90471363 17.70462
93            1 11.91917336 11.99011
94            1 14.46105857 18.42356
95            1  1.38115130 14.92371
96            1 28.33883988 16.60897
97            1 20.21450771 14.69829
98            1  6.70517870 11.39707
99            1 27.08232331 10.39271
100           1 19.21269606 15.07913
101           1 26.70268477 10.77975
102           1  1.16627432 19.92901
103           1 27.96514377 13.84657
104           1  0.21500076 13.81663
105           1 10.75766448 12.40766
106           1 13.12050846 19.29226
107           1  1.57355832 18.10184
108           1  9.29538527 17.67297
109           1 26.26807471 15.39002
110           1  9.74017067 18.41570
111           1 16.70221499 16.07380
112           1 13.11648876 15.27738
113           1 15.82924896 19.54311
114           1 18.59954804 18.75743
115           1  6.45501269 18.61968
116           1 27.84218557 18.33588
117           1 17.10159975 13.22954
118           1 24.30868549 15.85464
119           1 17.01677156 19.02580
120           1  3.40334528 11.97287
121           1  5.53478618 18.46267
122           1 18.77708562 17.71224
123           1 22.67428813 13.83164
124           1 16.72377637 12.35172
125           1  9.77678758 11.11058
126           1  5.67670837 15.03214
127           1  0.07211418 15.37998
128           1 16.05733008 14.14373
129           1  1.43336098 14.87125
130           1  4.60081269 14.75958
131           1  9.73013312 17.00414
132           1 10.47236488 18.16461
133           1  9.96437280 13.06820
134           1 23.35217732 13.11177
135           1 18.87619826 10.80808
136           1 12.23189945 14.04539
137           1 14.86156901 14.08419
138           1 23.63052885 17.06683
139           1 16.71979667 13.70694
140           1 14.93658144 13.76466
141           1 28.56726173 14.06356
142           1 24.94012931 17.27515
143           1 15.49415249 14.19119
144           1 16.97706906 11.01732
145           1 28.96205420 14.53108
146           1 15.71221253 17.81142
147           1 25.29777259 15.30815
148           1 23.61929812 13.69507
149           1  9.33856281 18.72297
150           1 29.19288259 11.57232
151           1 28.01400196 11.38900
152           1 11.81503752 17.49778
153           1  5.60530061 10.75768
154           1 28.48947077 15.18413
155           1  4.85282763 12.87153
156           1 22.34695001 19.51592
157           1 24.14828588 11.55500
158           1 17.64847732 14.76040
159           1 27.08401588 12.44893
160           1 28.51078386 17.46516
161           1  5.92444717 16.63261
162           1 12.72590685 13.67538
163           1  3.93794781 18.58969
164           1 15.49459361 16.61230
165           1  8.89701352 13.30316
166           1  3.72436841 17.31881
167           1  2.41288047 10.34141
168           1  4.31144123 18.06772
169           1 21.46084698 12.25055
170           1  3.34014048 11.89053
171           1 22.28174626 16.43798
172           1  0.07971427 18.25177
173           1 13.66240869 14.55288
174           1  3.37509154 17.97972
175           1 12.09906204 19.71738
176           1 13.09438448 11.12564
177           1 18.38770032 10.86406
178           1  0.90526420 12.28241
179           1 24.53641264 19.47463
180           1 11.07248885 17.37872
181           1 29.75123278 19.12459
182           1  3.51813120 19.95068
183           1 11.31896202 11.94532
184           1 24.99027881 15.36263
185           1  9.23333831 18.39574
186           1  4.31053132 18.91188
187           1  9.88480875 18.88327
188           1  9.28301391 19.27461
189           1  6.72554405 18.64191
190           1 27.48191994 12.12843
191           1 20.11630363 10.19680
192           1 26.19192549 13.24297
193           1  7.17281753 14.80643
194           1 29.06524045 17.55250
195           1 11.41594983 17.40299
196           1 21.79984578 17.15397
197           1 23.32687551 18.35467
198           1 12.91897079 16.37776
199           1 12.76246575 19.53340
200           1 25.11776769 16.75712
201           1 11.88029100 11.31211
202           1 17.88826080 18.84513
203           1  1.45431655 16.62768
204           1 21.62140878 17.17863
205           1 28.00591232 10.59971
206           1 28.44356806 18.61046
207           1 19.38723746 11.54922
208           1 19.61220003 14.84114
209           1 20.05806120 17.31062
210           1  0.99825994 13.60334
211           1 29.38576008 18.41106
212           1  3.19759092 15.76667
213           1  3.74973696 12.07206
214           1 24.36117251 16.78401
215           1  4.64522211 18.92882
216           1 12.89546262 19.79099
217           1 28.65181496 13.82700
218           1  5.58822893 19.06452
219           1  3.00535957 15.63875
220           1 23.48511413 17.46122
221           1  2.25061891 13.04117
222           1  4.83245644 15.74218
223           1  5.30237749 11.85141
224           1 13.08378044 19.66038
225           1  1.79869444 10.44291
226           1 26.93214454 14.14182
227           1 24.80540332 10.34122
228           1 17.16293783 17.38747
229           1 28.74116243 13.51417
230           1 20.19843590 1