Last data update: 2014.03.03

R: Classic gamma regression. Log link for the mean
gammahetero1R Documentation

Classic gamma regression. Log link for the mean

Description

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

Usage

gammahetero1(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 <- log(runif(500,0,30))
X3 <- log(runif(500,0,15))
X4 <- log(runif(500,10,20))
mui <- exp(-5 + 0.2*X2 - 0.03*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=gammahetero1(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/gammahetero1.Rd_%03d_medium.png", width=480, height=480)
> ### Name: gammahetero1
> ### Title: Classic gamma regression. Log link for the mean
> ### Aliases: gammahetero1
> ### Keywords: Gamma regression Fisher Scoring Joint modeling Mean and shape
> ###   parameters
> 
> ### ** Examples
> 
> 
> # Simulation Example
> 
> X1 <- rep(1,500)
> X2 <- log(runif(500,0,30))
> X3 <- log(runif(500,0,15))
> X4 <- log(runif(500,10,20))
> mui <- exp(-5 + 0.2*X2 - 0.03*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=gammahetero1(formula.mean,formula.shape)
> a
$X
    (Intercept)           X2           X3
1             1  2.460480688  1.120404490
2             1  3.071072147  0.167278094
3             1  2.440858211  2.154166883
4             1  2.618693808  2.294014411
5             1  1.826164137  2.323293227
6             1  2.638393996  1.935791800
7             1  2.935940018  0.596960173
8             1  2.988958046  2.030301800
9             1  3.240366465  2.450119521
10            1  2.835281440  1.335811846
11            1  3.055818422  1.930035679
12            1  0.091028114  2.477203625
13            1  3.363909543 -0.130843244
14            1  2.852081529  2.527786854
15            1  3.244143770  1.114333810
16            1  0.166773924  2.373344632
17            1  2.149602261  2.369002714
18            1  2.828009072  2.069024994
19            1  2.375838650  2.666658187
20            1  2.265606722  2.235112825
21            1  1.554575588  2.625095117
22            1  3.278303240  1.889365082
23            1  3.362375467  0.774328474
24            1  2.052741725  2.292696461
25            1  2.946978080  1.918096878
26            1  2.929980089  2.594624297
27            1  2.222360581  1.678551669
28            1  3.177142527  2.418650724
29            1  2.887638693  2.413255024
30            1  3.154703314  1.179292390
31            1  2.975710780 -0.105505594
32            1  0.985530729  1.488776048
33            1  2.448865888  1.556434911
34            1  1.831423716  2.510615933
35            1  2.701796527 -3.198418066
36            1  3.118742204  2.263996252
37            1  2.928075803  2.113261485
38            1  3.126200949  2.410952056
39            1 -0.859185349  1.285501214
40            1  3.374052489  2.491659971
41            1  2.737445376  2.658761498
42            1  0.910741737  2.510389581
43            1  3.080265271  2.342514336
44            1  3.238539599  2.425318183
45            1  3.189848140  1.807019812
46            1  3.232614692  1.831736668
47            1  2.780758406 -0.612308324
48            1  1.241923419  0.068563860
49            1  1.765220934  1.671352394
50            1  3.381927135  1.835672450
51            1  1.117655816 -0.227966175
52            1  3.042464106  2.504055098
53            1  2.172354800  2.429977522
54            1  3.057388738  1.820470992
55            1  2.926433371 -0.216628440
56            1  2.967875342  2.680590416
57            1  2.037444571  2.080182241
58            1  2.065586255  1.868325115
59            1  3.304697837 -0.029889715
60            1  1.036689768  1.794501836
61            1 -0.318951234  2.673101571
62            1  0.375193800  2.340049922
63            1  2.847319432  2.511664153
64            1  2.740197891  2.319850656
65            1  2.226536462  2.042099250
66            1  2.646354506  2.684371086
67            1  2.266911991  1.490969523
68            1  2.937220415  0.476458586
69            1  2.408346128  1.828713392
70            1  1.706080125  1.502022371
71            1  3.056765862  2.449261448
72            1  2.286234464  1.088270817
73            1  2.923591233  0.575852819
74            1  1.896383715  1.205640769
75            1  3.123462155  1.801003964
76            1  2.117160349  2.671575073
77            1  3.007337828  2.671479235
78            1  3.154686214  1.884236192
79            1  1.731989718  2.707674801
80            1  3.091378653  2.336377621
81            1  2.783963382  2.273444703
82            1  2.967726296  2.222501062
83            1  3.082657791  0.595668389
84            1  0.532691688  2.559538467
85            1  2.569812935  2.508630113
86            1  2.310591331  2.603095861
87            1  2.936082347  2.173105450
88            1  1.878904773  0.834001009
89            1  0.320255619  2.567750880
90            1  2.956004820  2.511460668
91            1  3.383270997  2.495437419
92            1  2.230798887  2.135124457
93            1  2.433218333  2.452194515
94            1  3.216669901  1.420024484
95            1  2.480656789  1.389948693
96            1  2.693471380  1.886521772
97            1  3.164106769 -3.401434174
98            1  2.908439282 -0.343389370
99            1  3.361252241  0.777300863
100           1  3.147999131  2.504313661
101           1 -0.703820559  2.232839150
102           1  2.239839639  1.330740851
103           1  2.268575109  1.366000439
104           1  1.564032778  2.588393913
105           1  2.404692474  2.387857335
106           1  2.762801434  1.665790410
107           1  2.061666252  2.412437027
108           1  2.161150920  2.630852313
109           1  2.855975092 -0.070135411
110           1  3.207020299  2.453447576
111           1  2.384599617 -0.043812022
112           1  3.334348890  2.626846729
113           1  2.773900698  1.541380039
114           1  1.621917883  1.607260098
115           1  2.242866077  1.978151397
116           1  0.575079351  0.155986499
117           1  1.405602820  1.199239177
118           1  2.018847650  0.828286111
119           1  1.722424815  1.970745162
120           1  1.208933683  0.216658902
121           1  3.233818120  2.674226652
122           1  3.240030357  2.610220491
123           1  3.311632897  2.047514689
124           1  2.083146204  2.606142455
125           1  2.376692210  2.285372387
126           1  1.275008747  2.701335089
127           1  2.831202109  1.985525491
128           1 -0.073983614  0.208599230
129           1  3.357065681  2.421596712
130           1  2.355172637  1.885895613
131           1  2.407472960  2.330554197
132           1  2.834613402 -5.346217326
133           1  2.996062571  2.250262501
134           1  3.346360428  2.620764597
135           1  2.433666344  2.551772332
136           1  3.220576572  2.183631640
137           1  3.328759902  2.030219131
138           1  1.789766611  2.684088661
139           1  3.027936974  2.035471493
140           1  0.676941112  2.004452773
141           1  3.393158993  1.813472336
142           1  3.260853974  2.377777893
143           1  2.666717714 -3.289485144
144           1  3.305098774  2.327108674
145           1  0.706266458  2.414336755
146           1  2.785883198  2.334710323
147           1  2.938536921  2.524733758
148           1  2.503005169  2.198525666
149           1  2.695300765  1.536051300
150           1  3.362062213  1.639822407
151           1  0.862644173  2.313999571
152           1  2.707831350  2.154516896
153           1  2.046967096  1.295952790
154           1  0.497790469  2.629217550
155           1  2.559682288  1.423949200
156           1  1.758879304  2.476703282
157           1  2.855825126  1.164700206
158           1  3.221774146  0.093638119
159           1  3.094589591  2.675039268
160           1  2.700893568  0.861934089
161           1  2.719173967 -0.007289349
162           1  2.442341816  2.539642936
163           1  2.495398457  0.679228212
164           1  1.762346127  0.709263718
165           1  2.937460220  2.307625531
166           1  3.166895092  1.066479491
167           1  1.263953605  2.184678098
168           1  3.362259222  0.356053377
169           1  1.045504910  2.338568705
170           1  2.915906143  0.214422858
171           1  3.146758238  2.605007885
172           1  2.291718595  2.479324924
173           1  3.227775234  2.323300102
174           1 -0.783912640 -1.430090963
175           1  3.172622970  2.508217373
176           1  2.114203182  2.579222883
177           1  1.992377156  2.544268034
178           1  2.864148451  2.512169287
179           1  2.116114141  2.196800117
180           1  0.274527917  2.456019590
181           1 -2.569331122  2.178823564
182           1  3.346411579  2.417298464
183           1  3.323137450  0.714891902
184           1  2.665984478  2.034958992
185           1  2.177240922  1.622411691
186           1  3.147762399  2.697312608
187           1  2.927397467  2.678955960
188           1  1.975934845 -0.797617297
189           1  2.895050444  1.894673963
190           1  2.884653981  2.564646032
191           1  0.676399728  2.605877565
192           1  3.354803752  2.257388054
193           1 -0.007892813  0.990646549
194           1  3.147483620 -0.078889889
195           1 -0.396902830  2.541473040
196           1  3.309155832  1.668340672
197           1  1.788526780  1.587347810
198           1  3.372676314  2.607749501
199           1  1.757613284  2.247428984
200           1  3.191592272  1.745226502
201           1  1.910018084  2.472320919
202           1  2.923756646  2.044357689
203           1  2.253487727  1.260896060
204           1  2.147294735  2.231090663
205           1  1.198908425  2.226143212
206           1  1.011679639  1.320233079
207           1  1.757819918  2.232934344
208           1  2.127352918  1.784196901
209           1  3.104226282  2.679561830
210           1  2.501228731  1.539144780
211           1  2.069247766  2.090627893
212           1  3.381589018  0.610992030
213           1  3.045839779  2.324564439
214           1  2.940878267  1.332290264
215           1  2.588592460  1.559279467
216           1  2.486228058  0.225366780
217           1  3.056784337  1.873743398
218           1  3.355536094  2.465893016
219           1  3.064819442  0.221387787
220           1  3.341604187  2.498034789
221           1  2.540894755  0.613624326
222           1  0.659056725  2.554751671
223           1  1.343191553  2.373074032
224           1  3.266073408  1.749091341
225           1  2.666358525  1.977559404
226           1  2.900256571  2.334911344
227           1  2.614035994  0.015921715
228           1  1.647565766  1.517055632
229           1  2.401490360  1.131763104
230           1  1.294732181  2.117990634
231           1  3.157991301  2.159991961
232           1  3.315108701  1.905301420
233           1  2.719300797  2.294276390
234           1  3.192396873  2.622789253
235           1  2.590963819  1.356043100
236           1  2.447221905  1.681592417
237           1  2.986701793  1.774017514
238           1  2.748099854  1.119869922
239           1  1.735663492  1.391567861
240           1  3.336334997  2.010552799
241           1  0.567850026  1.434806601
242           1  3.016947569  0.519094623
243           1  2.719986885  1.406290079
244           1  2.647634512  0.774331475
245           1  0.143897864  1.612916763
246           1  2.470169986  2.022019693
247           1  3.300728367  1.322078316
248           1  3.398464329  2.224892208
249           1  2.605823074  0.102700760
250           1  3.310616580  1.376985561
251           1  2.447148086  2.211286167
252           1  3.019134233  1.688149575
253           1  2.719719959  2.466391043
254           1  2.475962501  1.954298396
255           1  2.884451599 -1.863261944
256           1  1.592069036  2.401903261
257           1  2.333158593  1.995107339
258           1  3.266657992  2.444721896
259           1  3.053073330  2.227506507
260           1  2.240666667  0.055880563
261           1  3.252397766 -0.069526213
262           1  2.104556342  0.791207331
263           1  2.449187211  2.097976322
264           1  2.609277898  2.675811674
265           1  3.269678124  2.163963087
266           1  1.100730208  2.380722632
267           1  3.277476531 -1.664089953
268           1  3.050692178  0.813422408
269           1  2.787228728  2.390683807
270           1  2.993001397  2.405418844
271           1  2.222253576  1.975402187
272           1  2.283143934 -1.140040059
273           1  2.429452191  1.851009019
274           1  3.125581828  1.854836582
275           1  2.766462315  2.585199902
276           1  1.768827416  1.263959168
277           1  3.143360947  1.446975323
278           1  2.421469399  1.988512775
279           1  3.335701463  2.078519762
280           1  2.411731884  2.550554079
281           1  2.979372982  1.557808340
282           1  2.101155579  1.261832079
283           1  2.793837348  0.338726117
284           1  3.283910768  0.384504968
285           1  2.998793648  1.673881325
286           1  2.802454334  2.556401146
287           1  1.485204945  2.639569319
288           1  3.337648910  2.225650253
289           1  2.973040647  0.795144679
290           1  2.262494861  2.315622351
291           1  3.274096299  1.310890969
292           1  2.918646383  1.950097425
293           1  2.559053739 -0.505569276
294           1  3.350390148  2.382137344
295           1  2.419304192  2.011353788
296           1  1.957217944  2.598115603
297           1  2.959116501  2.539921822
298           1  2.792229597  1.881847209
299           1  3.144366861  1.284292832
300           1  1.785664920  1.901824301
301           1  2.397135588  2.049999961
302           1  2.887239447  0.858389997
303           1  0.038795304  0.757126015
304           1  2.833804828  1.603355788
305           1  3.320247208  2.278551374
306           1  2.072566734  2.187526546
307           1  2.332281708  1.007569014
308           1  2.584918820  2.598601641
309           1  1.349915663  2.574509321
310           1  2.791874994  1.079545499
311           1 -0.331121568  2.607663508
312           1 -2.035748030  2.095037283
313           1  1.194757754  1.797635465
314           1  3.389586772  2.581845410
315           1  1.363914383  2.302028888
316           1  2.313463936  2.549501311
317           1  3.148170053  2.465496566
318           1  3.150853524  0.666558151
319           1  2.781531339  2.307202930
320           1  1.332544592  0.016496183
321           1  2.420817498 -0.541894554
322           1  3.293641035  2.391697237
323           1  3.255547065  1.499868124
324           1  2.763210216  1.377744567
325           1  1.410009976  1.663000723
326           1  2.592573437  2.056015170
327           1  3.103339530  2.068362090
328           1  3.226161732  1.502849303
329           1  3.030091408  2.325075419
330           1  3.224536858  1.810700148
331           1  2.775301164  0.389893815
332           1  1.299748534  1.727256667
333           1  3.365733987  2.252407720
334           1  3.053984567  2.416361483
335           1  2.962382854  2.648455492
336           1  3.072555662  2.636210625
337           1  1.933103158  1.486923435
338           1  2.985277328  0.394119482
339           1  3.363893975  1.625751898
340           1  2.271356153  2.503321010
341           1  0.893770840  2.690357247
342           1  3.257737314  0.798852808
343           1  1.725649314  2.642459994
344           1  3.256423957  1.888593145
345           1  3.001769957  1.071294085
346           1  3.200356806  0.263348436
347           1  3.206719920 -0.049612036
348           1  0.461449930  2.469225298
349           1  2.844137576  2.285099531
350           1  3.193124124 -0.672288512
351           1  1.488738114  2.444935847
352           1  2.182161907  2.464719156
353           1  3.370136429  2.483955821
354           1  1.001550976  1.754060809
355           1  2.610241305 -1.128045954
356           1  2.068681910  0.883873050
357           1  0.404072068  2.546562299
358           1  2.524458483  1.516434202
359           1  2.223835770  2.374101851
360           1  1.748102017  2.497042052
361           1  1.460324354  1.186644846
362           1  0.561896551  1.555553257
363           1  3.335020786  2.213675121
364           1  2.163845210  2.585626685
365           1  1.522595156 -0.524945518
366           1  3.199911723  1.710412612
367           1  2.945271201 -0.312621208
368           1  1.582826896  1.321613197
369           1  2.919416908  2.039404923
370           1  0.357593936 -0.033686338
371           1  3.041312564  1.510474951
372           1  2.338400442  0.457775495
373           1  3.070358841  1.870537093
374           1  3.391105332  2.513523007
375           1  1.471423031  2.433933332
376           1  1.988687911 -0.202612053
377           1  3.247397756  0.422361364
378           1  0.877834904  1.713613829
379           1  2.162680294  1.960524322
380           1 -2.649916428  2.135430519
381           1  2.859585732  2.151393047
382           1  3.213638694  1.436840515
383           1  0.984016020  2.649315211
384           1  2.273648613  2.111860766
385           1  2.907326727  2.178779346
386           1  1.512502706  2.400597371
387           1  1.958592890  2.437673492
388           1  2.868824440  1.260306784
389           1  2.490552651  2.181240356
390           1  3.017593078  2.569130966
391           1  2.636978272  1.686480132
392           1  2.492517297  2.241176118
393           1  3.295607813  1.258667994
394           1  3.075887746  1.173315405
395           1  2.499588740  2.355653533
396           1  2.540360888  1.463631885
397           1  2.668663723 -1.609298619
398           1  1.516225745  2.309683740
399           1  2.053939272  1.154277261
400           1  3.123710863  1.575748251
401           1  2.025954427  2.382127933
402           1  2.557242348  1.687575981
403           1  2.974072080  0.244488267
404           1 -0.699946190  2.177579539
405           1  2.181227564  2.574938188
406           1  2.700994229  0.414527107
407           1  2.693757952  2.536723177
408           1  3.065399296  2.566953409
409           1  2.187411211  2.228065838
410           1  3.295199997  2.320405931
411           1  3.110832683  1.712702314
412           1  2.403760357  1.613834154
413           1  2.748275818  1.258197778
414           1  1.565769615  1.890968511
415           1  2.741997114 -0.282306067
416           1  0.649910132  2.209922665
417           1  3.117808319  1.278159261
418           1  2.989393032 -0.446922850
419           1 -3.173005494  2.304632687
420           1  1.494358536  0.884797814
421           1  2.364028937  2.264624062
422           1  2.709289627  2.563992964
423           1  3.074090203  1.836978699
424           1  1.941166660  1.499679451
425           1  3.206267821  2.582279996
426           1  2.363173820  2.339736500
427           1  2.105023935  2.292542434
428           1  3.398022707  2.387039081
429           1  1.684824999  0.752870682
430           1  1.852303453  2.032968412
431           1  1.287162286  2.663951935
432           1  1.072025457  0.604090952
433           1  2.530068659  1.800241394
434           1 -1.333571802  2.636358931
435           1  3.123540010  1.657026864
436           1  3.138957535  2.335175056
437           1  2.234040745  1.639792981
438           1  3.207572724  2.330970745
439           1  3.049647605  2.496227791
440           1  3.276264703 -1.690186245
441           1  3.047572921  2.427728934
442           1  3.196693006  2.182477447
443           1  1.394150197  2.023929558
444           1  2.313363826  2.247763648
445           1  2.299394600  1.571194713
446           1  2.143083832  2.583103668
447           1  3.348067692  1.566614382
448           1  3.349740205  1.225505029
449           1  3.360500203  2.458157678
450           1  2.790346899  1.333462315
451           1  3.107005587  0.396187821
452           1  1.601673207  1.696913943
453           1  2.681032017  1.743841282
454           1  1.610606946  2.303390109
455           1  2.977389156  2.128371732
456           1  3.210224147  2.135381647
457           1  1.016460258  2.344981483
458           1  1.255553849  1.121598208
459           1  3.177046504  2.649674682
460           1  3.217192080  1.951129327
461           1  2.911192755  1.634274499
462           1  2.879244056  2.618805659
463           1  2.179225989  2.043727591
464           1  2.173901335  1.597724710
465           1  2.383076827  1.375017083
466           1  3.246193983  2.093651897
467           1  1.889075232  2.339904370
468           1  3.320559653  2.407799663
469           1  3.205279758  2.540947469
470           1  2.821683308  1.898644646
471           1  2.418659383  2.315682354
472           1  3.222463973  2.191343428
473           1  2.575822788  1.512907290
474           1  3.038305168  1.148629228
475           1  2.575726585  1.267769826
476           1  3.240169547  2.548373595
477           1  1.665022848  1.107764601
478           1  3.271516200  1.944457208
479           1  2.421184102  1.219873991
480           1  3.239299620  2.515358909
481           1  3.287769954  2.009348665
482           1  2.898186595  2.640159644
483           1  2.207611951  2.263051266
484           1  2.576924420  2.525082407
485           1  1.831211072  0.544520530
486           1  2.565928761  2.491760090
487           1  1.617068592  0.594944731
488           1  2.638257666  0.914483984
489           1  3.123883560  2.027030742
490           1  2.531115010  2.461721649
491           1  1.212824990  1.201666320
492           1  0.538463973  0.587009581
493           1  3.358942158  1.738412230
494           1  2.568258228  1.653537835
495           1  2.204635802  0.929262802
496           1  3.031996416  2.062046414
497           1  2.760103328  1.466372981
498           1  1.579998176  2.422334874
499           1  3.395522357  1.873074347
500           1  2.056175810 -0.247798663
attr(,"assign")
[1] 0 1 2

$Z
    (Intercept)           X2       X4
1             1  2.460480688 2.735550
2             1  3.071072147 2.747352
3             1  2.440858211 2.439790
4             1  2.618693808 2.876898
5             1  1.826164137 2.844138
6             1  2.638393996 2.625225
7             1  2.935940018 2.908613
8             1  2.988958046 2.572912
9             1  3.240366465 2.711681
10            1  2.835281440 2.510048
11            1  3.055818422 2.886118
12            1  0.091028114 2.696716
13            1  3.363909543 2.894725
14            1  2.852081529 2.428150
15            1  3.244143770 2.838094
16            1  0.166773924 2.625847
17            1  2.149602261 2.943121
18            1  2.828009072 2.718452
19            1  2.375838650 2.872509
20            1  2.265606722 2.990748
21            1  1.554575588 2.653091
22            1  3.278303240 2.476100
23            1  3.362375467 2.940425
24            1  2.052741725 2.861607
25            1  2.946978080 2.724284
26            1  2.929980089 2.896996
27            1  2.222360581 2.685013
28            1  3.177142527 2.906905
29            1  2.887638693 2.502420
30            1  3.154703314 2.981337
31            1  2.975710780 2.371129
32            1  0.985530729 2.683553
33            1  2.448865888 2.497984
34            1  1.831423716 2.986070
35            1  2.701796527 2.415438
36            1  3.118742204 2.716570
37            1  2.928075803 2.827599
38            1  3.126200949 2.522805
39            1 -0.859185349 2.936775
40            1  3.374052489 2.576069
41            1  2.737445376 2.531094
42            1  0.910741737 2.370615
43            1  3.080265271 2.704788
44            1  3.238539599 2.841428
45            1  3.189848140 2.372029
46            1  3.232614692 2.536877
47            1  2.780758406 2.780941
48            1  1.241923419 2.924120
49            1  1.765220934 2.503191
50            1  3.381927135 2.729047
51            1  1.117655816 2.386323
52            1  3.042464106 2.885083
53            1  2.172354800 2.831557
54            1  3.057388738 2.798840
55            1  2.926433371 2.919968
56            1  2.967875342 2.947760
57            1  2.037444571 2.380580
58            1  2.065586255 2.812716
59            1  3.304697837 2.662491
60            1  1.036689768 2.319312
61            1 -0.318951234 2.585073
62            1  0.375193800 2.634904
63            1  2.847319432 2.856227
64            1  2.740197891 2.746857
65            1  2.226536462 2.710501
66            1  2.646354506 2.363469
67            1  2.266911991 2.371539
68            1  2.937220415 2.756212
69            1  2.408346128 2.866101
70            1  1.706080125 2.614560
71            1  3.056765862 2.654678
72            1  2.286234464 2.402316
73            1  2.923591233 2.940560
74            1  1.896383715 2.844914
75            1  3.123462155 2.494109
76            1  2.117160349 2.592996
77            1  3.007337828 2.695101
78            1  3.154686214 2.556280
79            1  1.731989718 2.828712
80            1  3.091378653 2.809817
81            1  2.783963382 2.785531
82            1  2.967726296 2.536280
83            1  3.082657791 2.742773
84            1  0.532691688 2.574384
85            1  2.569812935 2.376426
86            1  2.310591331 2.561168
87            1  2.936082347 2.430295
88            1  1.878904773 2.967608
89            1  0.320255619 2.883064
90            1  2.956004820 2.891493
91            1  3.383270997 2.346693
92            1  2.230798887 2.630993
93            1  2.433218333 2.916717
94            1  3.216669901 2.423384
95            1  2.480656789 2.488199
96            1  2.693471380 2.889477
97            1  3.164106769 2.508236
98            1  2.908439282 2.348897
99            1  3.361252241 2.699288
100           1  3.147999131 2.320822
101           1 -0.703820559 2.857370
102           1  2.239839639 2.497852
103           1  2.268575109 2.540791
104           1  1.564032778 2.872060
105           1  2.404692474 2.747259
106           1  2.762801434 2.930373
107           1  2.061666252 2.758175
108           1  2.161150920 2.651403
109           1  2.855975092 2.873028
110           1  3.207020299 2.962778
111           1  2.384599617 2.711983
112           1  3.334348890 2.369505
113           1  2.773900698 2.458485
114           1  1.621917883 2.578025
115           1  2.242866077 2.615477
116           1  0.575079351 2.312383
117           1  1.405602820 2.902039
118           1  2.018847650 2.585774
119           1  1.722424815 2.779573
120           1  1.208933683 2.675746
121           1  3.233818120 2.568471
122           1  3.240030357 2.382846
123           1  3.311632897 2.728318
124           1  2.083146204 2.991124
125           1  2.376692210 2.961496
126           1  1.275008747 2.434378
127           1  2.831202109 2.483450
128           1 -0.073983614 2.868362
129           1  3.357065681 2.765098
130           1  2.355172637 2.859220
131           1  2.407472960 2.801704
132           1  2.834613402 2.760148
133           1  2.996062571 2.927712
134           1  3.346360428 2.416456
135           1  2.433666344 2.850020
136           1  3.220576572 2.915335
137           1  3.328759902 2.847072
138           1  1.789766611 2.594650
139           1  3.027936974 2.992236
140           1  0.676941112 2.312545
141           1  3.393158993 2.538321
142           1  3.260853974 2.315757
143           1  2.666717714 2.344006
144           1  3.305098774 2.847374
145           1  0.706266458 2.888703
146           1  2.785883198 2.632484
147           1  2.938536921 2.623419
148           1  2.503005169 2.895838
149           1  2.695300765 2.783061
150           1  3.362062213 2.578509
151           1  0.862644173 2.867781
152           1  2.707831350 2.727226
153           1  2.046967096 2.977961
154           1  0.497790469 2.788997
155           1  2.559682288 2.830984
156           1  1.758879304 2.954443
157           1  2.855825126 2.447229
158           1  3.221774146 2.426539
159           1  3.094589591 2.626264
160           1  2.700893568 2.946201
161           1  2.719173967 2.373834
162           1  2.442341816 2.880379
163           1  2.495398457 2.686756
164           1  1.762346127 2.759111
165           1  2.937460220 2.619481
166           1  3.166895092 2.886059
167           1  1.263953605 2.360536
168           1  3.362259222 2.930907
169           1  1.045504910 2.879252
170           1  2.915906143 2.970631
171           1  3.146758238 2.896649
172           1  2.291718595 2.315782
173           1  3.227775234 2.949960
174           1 -0.783912640 2.781046
175           1  3.172622970 2.862263
176           1  2.114203182 2.958445
177           1  1.992377156 2.462531
178           1  2.864148451 2.850404
179           1  2.116114141 2.804877
180           1  0.274527917 2.633488
181           1 -2.569331122 2.521749
182           1  3.346411579 2.681137
183           1  3.323137450 2.771395
184           1  2.665984478 2.817918
185           1  2.177240922 2.816556
186           1  3.147762399 2.590687
187           1  2.927397467 2.787244
188           1  1.975934845 2.522652
189           1  2.895050444 2.350877
190           1  2.884653981 2.622808
191           1  0.676399728 2.648467
192           1  3.354803752 2.787732
193           1 -0.007892813 2.750109
194           1  3.147483620 2.847032
195           1 -0.396902830 2.424965
196           1  3.309155832 2.796240
197           1  1.788526780 2.