#
# Elston, D. A., R. Moss, T. Boulinier, C. Arrowsmith, and X.
# Lambin. 2001. "Analysis of Aggregation, a Worked Example: Numbers
# of Ticks on Red Grouse Chicks." Parasitology 122 (05): 563-569.
# doi:10.1017/S0031182001007740. <URL:
# http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=82701>.
#
set loc TICKS    qua number of ticks sampled
set loc HEIGHT   qua height above sea level (meters)
set loc YEAR     cat year 1995-7
set loc LOCATION cat geographic location code
set loc cHEIGHT  qua centred height
read pedigree inline
# BROOD INDEX             TICKS HEIGHT YEAR LOCATION cHEIGHT           
501   1     s1   d1   x   0     465    95   32       2.75930521091811  
501   2     s2   d2   x   0     465    95   32       2.75930521091811  
502   3     s3   d3   x   0     472    95   36       9.75930521091811  
503   4     s4   d4   x   0     475    95   37       12.7593052109181  
503   5     s5   d5   x   0     475    95   37       12.7593052109181  
503   6     s6   d6   x   3     475    95   37       12.7593052109181  
503   7     s7   d7   x   2     475    95   37       12.7593052109181  
504   8     s8   d8   x   0     488    95   44       25.7593052109181  
504   9     s9   d9   x   0     488    95   44       25.7593052109181  
504   10    s10  d10  x   2     488    95   44       25.7593052109181  
505   11    s11  d11  x   0     492    95   47       29.7593052109181  
505   12    s12  d12  x   0     492    95   47       29.7593052109181  
505   13    s13  d13  x   0     492    95   47       29.7593052109181  
506   14    s14  d14  x   0     490    95   45       27.7593052109181  
506   15    s15  d15  x   0     490    95   45       27.7593052109181  
506   16    s16  d16  x   0     490    95   45       27.7593052109181  
507   17    s17  d17  x   0     464    95   31       1.75930521091811  
507   18    s18  d18  x   0     464    95   31       1.75930521091811  
507   19    s19  d19  x   0     464    95   31       1.75930521091811  
507   20    s20  d20  x   1     464    95   31       1.75930521091811  
507   21    s21  d21  x   2     464    95   31       1.75930521091811  
509   22    s22  d22  x   1     457    95   28       -5.24069478908189 
510   23    s23  d23  x   0     457    95   28       -5.24069478908189 
511   24    s24  d24  x   0     457    95   28       -5.24069478908189 
511   25    s25  d25  x   5     457    95   28       -5.24069478908189 
512   26    s26  d26  x   8     451    95   26       -11.2406947890819 
512   27    s27  d27  x   3     451    95   26       -11.2406947890819 
512   28    s28  d28  x   4     451    95   26       -11.2406947890819 
513   29    s29  d29  x   7     437    95   17       -25.2406947890819 
513   30    s30  d30  x   0     437    95   17       -25.2406947890819 
513   31    s31  d31  x   4     437    95   17       -25.2406947890819 
514   32    s32  d32  x   4     430    95   14       -32.2406947890819 
514   33    s33  d33  x   1     430    95   14       -32.2406947890819 
514   34    s34  d34  x   0     430    95   14       -32.2406947890819 
514   35    s35  d35  x   0     430    95   14       -32.2406947890819 
514   36    s36  d36  x   3     430    95   14       -32.2406947890819 
514   37    s37  d37  x   1     430    95   14       -32.2406947890819 
515   38    s38  d38  x   6     427    95   13       -35.2406947890819 
515   39    s39  d39  x   0     427    95   13       -35.2406947890819 
515   40    s40  d40  x   1     427    95   13       -35.2406947890819 
515   41    s41  d41  x   0     427    95   13       -35.2406947890819 
516   42    s42  d42  x   2     419    95   7        -43.2406947890819 
516   43    s43  d43  x   7     419    95   7        -43.2406947890819 
516   44    s44  d44  x   31    419    95   7        -43.2406947890819 
517   45    s45  d45  x   34    411    95   4        -51.2406947890819 
517   46    s46  d46  x   17    411    95   4        -51.2406947890819 
517   47    s47  d47  x   16    411    95   4        -51.2406947890819 
518   48    s48  d48  x   66    419    95   7        -43.2406947890819 
518   49    s49  d49  x   49    419    95   7        -43.2406947890819 
518   50    s50  d50  x   82    419    95   7        -43.2406947890819 
518   51    s51  d51  x   85    419    95   7        -43.2406947890819 
518   52    s52  d52  x   64    419    95   7        -43.2406947890819 
519   53    s53  d53  x   11    424    95   11       -38.2406947890819 
519   54    s54  d54  x   14    424    95   11       -38.2406947890819 
519   55    s55  d55  x   4     424    95   11       -38.2406947890819 
519   56    s56  d56  x   10    424    95   11       -38.2406947890819 
520   57    s57  d57  x   3     427    95   13       -35.2406947890819 
520   58    s58  d58  x   15    427    95   13       -35.2406947890819 
520   59    s59  d59  x   8     427    95   13       -35.2406947890819 
521   60    s60  d60  x   9     422    95   9        -40.2406947890819 
521   61    s61  d61  x   11    422    95   9        -40.2406947890819 
521   62    s62  d62  x   7     422    95   9        -40.2406947890819 
521   63    s63  d63  x   13    422    95   9        -40.2406947890819 
522   64    s64  d64  x   3     503    95   53       40.7593052109181  
523   65    s65  d65  x   0     496    95   51       33.7593052109181  
523   66    s66  d66  x   0     496    95   51       33.7593052109181  
523   67    s67  d67  x   1     496    95   51       33.7593052109181  
523   68    s68  d68  x   1     496    95   51       33.7593052109181  
525   69    s69  d69  x   0     507    95   54       44.7593052109181  
525   70    s70  d70  x   0     507    95   54       44.7593052109181  
525   71    s71  d71  x   0     507    95   54       44.7593052109181  
525   72    s72  d72  x   0     507    95   54       44.7593052109181  
526   73    s73  d73  x   0     496    95   51       33.7593052109181  
526   74    s74  d74  x   0     496    95   51       33.7593052109181  
526   75    s75  d75  x   0     496    95   51       33.7593052109181  
526   76    s76  d76  x   1     496    95   51       33.7593052109181  
526   77    s77  d77  x   1     496    95   51       33.7593052109181  
526   78    s78  d78  x   0     496    95   51       33.7593052109181  
528   79    s79  d79  x   2     466    95   33       3.75930521091811  
528   80    s80  d80  x   0     466    95   33       3.75930521091811  
528   81    s81  d81  x   3     466    95   33       3.75930521091811  
531   82    s82  d82  x   1     488    95   44       25.7593052109181  
533   83    s83  d83  x   7     442    95   19       -20.2406947890819 
533   84    s84  d84  x   2     442    95   19       -20.2406947890819 
533   85    s85  d85  x   16    442    95   19       -20.2406947890819 
533   86    s86  d86  x   12    442    95   19       -20.2406947890819 
533   87    s87  d87  x   0     442    95   19       -20.2406947890819 
535   88    s88  d88  x   1     442    95   19       -20.2406947890819 
537   89    s89  d89  x   0     533    95   63       70.7593052109181  
537   90    s90  d90  x   0     533    95   63       70.7593052109181  
537   91    s91  d91  x   1     533    95   63       70.7593052109181  
537   92    s92  d92  x   0     533    95   63       70.7593052109181  
537   93    s93  d93  x   0     533    95   63       70.7593052109181  
537   94    s94  d94  x   1     533    95   63       70.7593052109181  
537   95    s95  d95  x   0     533    95   63       70.7593052109181  
538   96    s96  d96  x   0     533    95   63       70.7593052109181  
539   97    s97  d97  x   0     515    95   59       52.7593052109181  
539   98    s98  d98  x   0     515    95   59       52.7593052109181  
539   99    s99  d99  x   0     515    95   59       52.7593052109181  
539   100   s100 d100 x   0     515    95   59       52.7593052109181  
540   101   s101 d101 x   5     518    95   60       55.7593052109181  
540   102   s102 d102 x   2     518    95   60       55.7593052109181  
542   103   s103 d103 x   2     493    95   48       30.7593052109181  
542   104   s104 d104 x   1     493    95   48       30.7593052109181  
542   105   s105 d105 x   1     493    95   48       30.7593052109181  
548   106   s106 d106 x   0     468    95   34       5.75930521091811  
548   107   s107 d107 x   0     468    95   34       5.75930521091811  
548   108   s108 d108 x   1     468    95   34       5.75930521091811  
549   109   s109 d109 x   1     476    95   38       13.7593052109181  
549   110   s110 d110 x   1     476    95   38       13.7593052109181  
549   111   s111 d111 x   0     476    95   38       13.7593052109181  
549   112   s112 d112 x   0     476    95   38       13.7593052109181  
550   113   s113 d113 x   5     446    95   22       -16.2406947890819 
550   114   s114 d114 x   3     446    95   22       -16.2406947890819 
553   115   s115 d115 x   2     460    95   30       -2.24069478908189 
559   116   s116 d116 x   2     525    95   62       62.7593052109181  
559   117   s117 d117 x   1     525    95   62       62.7593052109181  
601   118   s118 d118 x   1     410    96   3        -52.2406947890819 
601   119   s119 d119 x   0     410    96   3        -52.2406947890819 
601   120   s120 d120 x   2     410    96   3        -52.2406947890819 
601   121   s121 d121 x   5     410    96   3        -52.2406947890819 
601   122   s122 d122 x   1     410    96   3        -52.2406947890819 
601   123   s123 d123 x   2     410    96   3        -52.2406947890819 
601   124   s124 d124 x   2     410    96   3        -52.2406947890819 
602   125   s125 d125 x   3     417    96   6        -45.2406947890819 
602   126   s126 d126 x   14    417    96   6        -45.2406947890819 
602   127   s127 d127 x   11    417    96   6        -45.2406947890819 
602   128   s128 d128 x   9     417    96   6        -45.2406947890819 
602   129   s129 d129 x   4     417    96   6        -45.2406947890819 
602   130   s130 d130 x   10    417    96   6        -45.2406947890819 
602   131   s131 d131 x   33    417    96   6        -45.2406947890819 
602   132   s132 d132 x   19    417    96   6        -45.2406947890819 
602   133   s133 d133 x   16    417    96   6        -45.2406947890819 
603   134   s134 d134 x   16    430    96   14       -32.2406947890819 
603   135   s135 d135 x   13    430    96   14       -32.2406947890819 
603   136   s136 d136 x   11    430    96   14       -32.2406947890819 
603   137   s137 d137 x   7     430    96   14       -32.2406947890819 
603   138   s138 d138 x   4     430    96   14       -32.2406947890819 
603   139   s139 d139 x   11    430    96   14       -32.2406947890819 
604   140   s140 d140 x   1     456    96   27       -6.24069478908189 
604   141   s141 d141 x   1     456    96   27       -6.24069478908189 
604   142   s142 d142 x   4     456    96   27       -6.24069478908189 
605   143   s143 d143 x   6     457    96   28       -5.24069478908189 
605   144   s144 d144 x   2     457    96   28       -5.24069478908189 
605   145   s145 d145 x   7     457    96   28       -5.24069478908189 
605   146   s146 d146 x   8     457    96   28       -5.24069478908189 
605   147   s147 d147 x   14    457    96   28       -5.24069478908189 
606   148   s148 d148 x   6     430    96   14       -32.2406947890819 
606   149   s149 d149 x   13    430    96   14       -32.2406947890819 
606   150   s150 d150 x   5     430    96   14       -32.2406947890819 
606   151   s151 d151 x   8     430    96   14       -32.2406947890819 
606   152   s152 d152 x   13    430    96   14       -32.2406947890819 
606   153   s153 d153 x   17    430    96   14       -32.2406947890819 
606   154   s154 d154 x   5     430    96   14       -32.2406947890819 
606   155   s155 d155 x   1     430    96   14       -32.2406947890819 
606   156   s156 d156 x   1     430    96   14       -32.2406947890819 
606   157   s157 d157 x   2     430    96   14       -32.2406947890819 
607   158   s158 d158 x   7     423    96   10       -39.2406947890819 
608   159   s159 d159 x   7     421    96   8        -41.2406947890819 
608   160   s160 d160 x   11    421    96   8        -41.2406947890819 
609   161   s161 d161 x   1     525    96   62       62.7593052109181  
609   162   s162 d162 x   0     525    96   62       62.7593052109181  
610   163   s163 d163 x   5     509    96   55       46.7593052109181  
610   164   s164 d164 x   4     509    96   55       46.7593052109181  
611   165   s165 d165 x   0     499    96   52       36.7593052109181  
611   166   s166 d166 x   0     499    96   52       36.7593052109181  
611   167   s167 d167 x   0     499    96   52       36.7593052109181  
612   168   s168 d168 x   7     503    96   53       40.7593052109181  
612   169   s169 d169 x   5     503    96   53       40.7593052109181  
612   170   s170 d170 x   3     503    96   53       40.7593052109181  
612   171   s171 d171 x   1     503    96   53       40.7593052109181  
612   172   s172 d172 x   6     503    96   53       40.7593052109181  
614   173   s173 d173 x   1     492    96   47       29.7593052109181  
614   174   s174 d174 x   2     492    96   47       29.7593052109181  
615   175   s175 d175 x   14    491    96   46       28.7593052109181  
615   176   s176 d176 x   5     491    96   46       28.7593052109181  
615   177   s177 d177 x   27    491    96   46       28.7593052109181  
616   178   s178 d178 x   1     475    96   37       12.7593052109181  
616   179   s179 d179 x   2     475    96   37       12.7593052109181  
616   180   s180 d180 x   3     475    96   37       12.7593052109181  
616   181   s181 d181 x   0     475    96   37       12.7593052109181  
617   182   s182 d182 x   1     479    96   40       16.7593052109181  
617   183   s183 d183 x   0     479    96   40       16.7593052109181  
617   184   s184 d184 x   5     479    96   40       16.7593052109181  
617   185   s185 d185 x   5     479    96   40       16.7593052109181  
617   186   s186 d186 x   8     479    96   40       16.7593052109181  
617   187   s187 d187 x   21    479    96   40       16.7593052109181  
618   188   s188 d188 x   15    472    96   36       9.75930521091811  
618   189   s189 d189 x   15    472    96   36       9.75930521091811  
618   190   s190 d190 x   6     472    96   36       9.75930521091811  
618   191   s191 d191 x   19    472    96   36       9.75930521091811  
618   192   s192 d192 x   14    472    96   36       9.75930521091811  
621   193   s193 d193 x   1     485    96   42       22.7593052109181  
621   194   s194 d194 x   1     485    96   42       22.7593052109181  
621   195   s195 d195 x   3     485    96   42       22.7593052109181  
621   196   s196 d196 x   2     485    96   42       22.7593052109181  
621   197   s197 d197 x   3     485    96   42       22.7593052109181  
623   198   s198 d198 x   2     495    96   50       32.7593052109181  
623   199   s199 d199 x   5     495    96   50       32.7593052109181  
624   200   s200 d200 x   0     472    96   36       9.75930521091811  
624   201   s201 d201 x   6     472    96   36       9.75930521091811  
624   202   s202 d202 x   3     472    96   36       9.75930521091811  
625   203   s203 d203 x   1     458    96   29       -4.24069478908189 
625   204   s204 d204 x   0     458    96   29       -4.24069478908189 
625   205   s205 d205 x   1     458    96   29       -4.24069478908189 
625   206   s206 d206 x   6     458    96   29       -4.24069478908189 
625   207   s207 d207 x   1     458    96   29       -4.24069478908189 
626   208   s208 d208 x   85    449    96   24       -13.2406947890819 
626   209   s209 d209 x   45    449    96   24       -13.2406947890819 
626   210   s210 d210 x   68    449    96   24       -13.2406947890819 
626   211   s211 d211 x   84    449    96   24       -13.2406947890819 
626   212   s212 d212 x   50    449    96   24       -13.2406947890819 
628   213   s213 d213 x   13    442    96   19       -20.2406947890819 
628   214   s214 d214 x   1     442    96   19       -20.2406947890819 
629   215   s215 d215 x   19    448    96   23       -14.2406947890819 
629   216   s216 d216 x   26    448    96   23       -14.2406947890819 
629   217   s217 d217 x   9     448    96   23       -14.2406947890819 
629   218   s218 d218 x   2     448    96   23       -14.2406947890819 
629   219   s219 d219 x   4     448    96   23       -14.2406947890819 
629   220   s220 d220 x   3     448    96   23       -14.2406947890819 
630   221   s221 d221 x   22    448    96   23       -14.2406947890819 
630   222   s222 d222 x   32    448    96   23       -14.2406947890819 
631   223   s223 d223 x   5     403    96   1        -59.2406947890819 
631   224   s224 d224 x   21    403    96   1        -59.2406947890819 
631   225   s225 d225 x   26    403    96   1        -59.2406947890819 
631   226   s226 d226 x   13    403    96   1        -59.2406947890819 
631   227   s227 d227 x   23    403    96   1        -59.2406947890819 
632   228   s228 d228 x   42    411    96   4        -51.2406947890819 
632   229   s229 d229 x   38    411    96   4        -51.2406947890819 
632   230   s230 d230 x   61    411    96   4        -51.2406947890819 
632   231   s231 d231 x   79    411    96   4        -51.2406947890819 
632   232   s232 d232 x   39    411    96   4        -51.2406947890819 
632   233   s233 d233 x   41    411    96   4        -51.2406947890819 
634   234   s234 d234 x   15    415    96   5        -47.2406947890819 
634   235   s235 d235 x   23    415    96   5        -47.2406947890819 
634   236   s236 d236 x   14    415    96   5        -47.2406947890819 
635   237   s237 d237 x   7     427    96   13       -35.2406947890819 
636   238   s238 d238 x   24    424    96   11       -38.2406947890819 
638   239   s239 d239 x   3     525    96   62       62.7593052109181  
638   240   s240 d240 x   1     525    96   62       62.7593052109181  
640   241   s241 d241 x   2     521    96   61       58.7593052109181  
640   242   s242 d242 x   1     521    96   61       58.7593052109181  
640   243   s243 d243 x   0     521    96   61       58.7593052109181  
641   244   s244 d244 x   3     518    96   60       55.7593052109181  
641   245   s245 d245 x   8     518    96   60       55.7593052109181  
642   246   s246 d246 x   1     495    96   50       32.7593052109181  
642   247   s247 d247 x   2     495    96   50       32.7593052109181  
642   248   s248 d248 x   0     495    96   50       32.7593052109181  
643   249   s249 d249 x   8     495    96   50       32.7593052109181  
643   250   s250 d250 x   3     495    96   50       32.7593052109181  
643   251   s251 d251 x   14    495    96   50       32.7593052109181  
643   252   s252 d252 x   16    495    96   50       32.7593052109181  
643   253   s253 d253 x   18    495    96   50       32.7593052109181  
643   254   s254 d254 x   11    495    96   50       32.7593052109181  
643   255   s255 d255 x   13    495    96   50       32.7593052109181  
645   256   s256 d256 x   6     460    96   30       -2.24069478908189 
645   257   s257 d257 x   7     460    96   30       -2.24069478908189 
645   258   s258 d258 x   10    460    96   30       -2.24069478908189 
647   259   s259 d259 x   5     442    96   19       -20.2406947890819 
647   260   s260 d260 x   7     442    96   19       -20.2406947890819 
648   261   s261 d261 x   25    443    96   20       -19.2406947890819 
648   262   s262 d262 x   11    443    96   20       -19.2406947890819 
648   263   s263 d263 x   6     443    96   20       -19.2406947890819 
648   264   s264 d264 x   4     443    96   20       -19.2406947890819 
648   265   s265 d265 x   7     443    96   20       -19.2406947890819 
650   266   s266 d266 x   4     425    96   12       -37.2406947890819 
650   267   s267 d267 x   6     425    96   12       -37.2406947890819 
650   268   s268 d268 x   2     425    96   12       -37.2406947890819 
651   269   s269 d269 x   5     439    96   18       -23.2406947890819 
651   270   s270 d270 x   3     439    96   18       -23.2406947890819 
651   271   s271 d271 x   7     439    96   18       -23.2406947890819 
652   272   s272 d272 x   3     444    96   21       -18.2406947890819 
701   273   s273 d273 x   1     450    97   25       -12.2406947890819 
701   274   s274 d274 x   4     450    97   25       -12.2406947890819 
701   275   s275 d275 x   4     450    97   25       -12.2406947890819 
701   276   s276 d276 x   2     450    97   25       -12.2406947890819 
701   277   s277 d277 x   5     450    97   25       -12.2406947890819 
702   278   s278 d278 x   3     446    97   22       -16.2406947890819 
702   279   s279 d279 x   0     446    97   22       -16.2406947890819 
702   280   s280 d280 x   3     446    97   22       -16.2406947890819 
702   281   s281 d281 x   1     446    97   22       -16.2406947890819 
702   282   s282 d282 x   2     446    97   22       -16.2406947890819 
702   283   s283 d283 x   3     446    97   22       -16.2406947890819 
704   284   s284 d284 x   1     472    97   36       9.75930521091811  
704   285   s285 d285 x   0     472    97   36       9.75930521091811  
704   286   s286 d286 x   4     472    97   36       9.75930521091811  
704   287   s287 d287 x   0     472    97   36       9.75930521091811  
704   288   s288 d288 x   0     472    97   36       9.75930521091811  
705   289   s289 d289 x   0     472    97   36       9.75930521091811  
706   290   s290 d290 x   0     460    97   30       -2.24069478908189 
706   291   s291 d291 x   0     460    97   30       -2.24069478908189 
706   292   s292 d292 x   0     460    97   30       -2.24069478908189 
708   293   s293 d293 x   1     442    97   19       -20.2406947890819 
708   294   s294 d294 x   3     442    97   19       -20.2406947890819 
708   295   s295 d295 x   4     442    97   19       -20.2406947890819 
708   296   s296 d296 x   0     442    97   19       -20.2406947890819 
708   297   s297 d297 x   4     442    97   19       -20.2406947890819 
708   298   s298 d298 x   2     442    97   19       -20.2406947890819 
709   299   s299 d299 x   0     525    97   62       62.7593052109181  
709   300   s300 d300 x   0     525    97   62       62.7593052109181  
709   301   s301 d301 x   1     525    97   62       62.7593052109181  
710   302   s302 d302 x   0     533    97   63       70.7593052109181  
710   303   s303 d303 x   1     533    97   63       70.7593052109181  
710   304   s304 d304 x   2     533    97   63       70.7593052109181  
710   305   s305 d305 x   0     533    97   63       70.7593052109181  
710   306   s306 d306 x   0     533    97   63       70.7593052109181  
710   307   s307 d307 x   0     533    97   63       70.7593052109181  
711   308   s308 d308 x   1     513    97   57       50.7593052109181  
711   309   s309 d309 x   0     513    97   57       50.7593052109181  
711   310   s310 d310 x   0     513    97   57       50.7593052109181  
711   311   s311 d311 x   1     513    97   57       50.7593052109181  
711   312   s312 d312 x   1     513    97   57       50.7593052109181  
711   313   s313 d313 x   0     513    97   57       50.7593052109181  
711   314   s314 d314 x   0     513    97   57       50.7593052109181  
712   315   s315 d315 x   0     514    97   58       51.7593052109181  
712   316   s316 d316 x   0     514    97   58       51.7593052109181  
713   317   s317 d317 x   0     511    97   56       48.7593052109181  
713   318   s318 d318 x   0     511    97   56       48.7593052109181  
713   319   s319 d319 x   1     511    97   56       48.7593052109181  
713   320   s320 d320 x   0     511    97   56       48.7593052109181  
713   321   s321 d321 x   0     511    97   56       48.7593052109181  
714   322   s322 d322 x   0     511    97   56       48.7593052109181  
714   323   s323 d323 x   0     511    97   56       48.7593052109181  
714   324   s324 d324 x   1     511    97   56       48.7593052109181  
714   325   s325 d325 x   0     511    97   56       48.7593052109181  
714   326   s326 d326 x   0     511    97   56       48.7593052109181  
714   327   s327 d327 x   0     511    97   56       48.7593052109181  
715   328   s328 d328 x   0     496    97   51       33.7593052109181  
715   329   s329 d329 x   0     496    97   51       33.7593052109181  
715   330   s330 d330 x   0     496    97   51       33.7593052109181  
715   331   s331 d331 x   0     496    97   51       33.7593052109181  
716   332   s332 d332 x   1     494    97   49       31.7593052109181  
716   333   s333 d333 x   0     494    97   49       31.7593052109181  
717   334   s334 d334 x   0     494    97   49       31.7593052109181  
717   335   s335 d335 x   0     494    97   49       31.7593052109181  
717   336   s336 d336 x   0     494    97   49       31.7593052109181  
718   337   s337 d337 x   2     411    97   4        -51.2406947890819 
718   338   s338 d338 x   4     411    97   4        -51.2406947890819 
718   339   s339 d339 x   4     411    97   4        -51.2406947890819 
718   340   s340 d340 x   1     411    97   4        -51.2406947890819 
718   341   s341 d341 x   2     411    97   4        -51.2406947890819 
719   342   s342 d342 x   2     411    97   4        -51.2406947890819 
719   343   s343 d343 x   1     411    97   4        -51.2406947890819 
719   344   s344 d344 x   4     411    97   4        -51.2406947890819 
719   345   s345 d345 x   1     411    97   4        -51.2406947890819 
719   346   s346 d346 x   3     411    97   4        -51.2406947890819 
719   347   s347 d347 x   3     411    97   4        -51.2406947890819 
720   348   s348 d348 x   2     423    97   10       -39.2406947890819 
720   349   s349 d349 x   0     423    97   10       -39.2406947890819 
721   350   s350 d350 x   3     424    97   11       -38.2406947890819 
722   351   s351 d351 x   2     403    97   1        -59.2406947890819 
722   352   s352 d352 x   0     403    97   1        -59.2406947890819 
722   353   s353 d353 x   1     403    97   1        -59.2406947890819 
723   354   s354 d354 x   0     409    97   2        -53.2406947890819 
723   355   s355 d355 x   3     409    97   2        -53.2406947890819 
723   356   s356 d356 x   1     409    97   2        -53.2406947890819 
724   357   s357 d357 x   10    434    97   16       -28.2406947890819 
724   358   s358 d358 x   4     434    97   16       -28.2406947890819 
724   359   s359 d359 x   1     434    97   16       -28.2406947890819 
725   360   s360 d360 x   2     477    97   39       14.7593052109181  
725   361   s361 d361 x   0     477    97   39       14.7593052109181  
725   362   s362 d362 x   1     477    97   39       14.7593052109181  
727   363   s363 d363 x   0     472    97   36       9.75930521091811  
728   364   s364 d364 x   0     468    97   34       5.75930521091811  
728   365   s365 d365 x   0     468    97   34       5.75930521091811  
729   366   s366 d366 x   0     470    97   35       7.75930521091811  
730   367   s367 d367 x   0     486    97   43       23.7593052109181  
731   368   s368 d368 x   0     495    97   50       32.7593052109181  
731   369   s369 d369 x   0     495    97   50       32.7593052109181  
731   370   s370 d370 x   0     495    97   50       32.7593052109181  
731   371   s371 d371 x   0     495    97   50       32.7593052109181  
731   372   s372 d372 x   0     495    97   50       32.7593052109181  
732   373   s373 d373 x   0     483    97   41       20.7593052109181  
732   374   s374 d374 x   0     483    97   41       20.7593052109181  
732   375   s375 d375 x   2     483    97   41       20.7593052109181  
733   376   s376 d376 x   2     442    97   19       -20.2406947890819 
733   377   s377 d377 x   2     442    97   19       -20.2406947890819 
734   378   s378 d378 x   0     457    97   28       -5.24069478908189 
736   379   s379 d379 x   4     457    97   28       -5.24069478908189 
736   380   s380 d380 x   5     457    97   28       -5.24069478908189 
737   381   s381 d381 x   2     457    97   28       -5.24069478908189 
737   382   s382 d382 x   3     457    97   28       -5.24069478908189 
737   383   s383 d383 x   2     457    97   28       -5.24069478908189 
737   384   s384 d384 x   2     457    97   28       -5.24069478908189 
737   385   s385 d385 x   1     457    97   28       -5.24069478908189 
737   386   s386 d386 x   2     457    97   28       -5.24069478908189 
737   387   s387 d387 x   0     457    97   28       -5.24069478908189 
738   388   s388 d388 x   2     464    97   31       1.75930521091811  
738   389   s389 d389 x   0     464    97   31       1.75930521091811  
738   390   s390 d390 x   1     464    97   31       1.75930521091811  
739   391   s391 d391 x   0     433    97   15       -29.2406947890819 
739   392   s392 d392 x   3     433    97   15       -29.2406947890819 
739   393   s393 d393 x   1     433    97   15       -29.2406947890819 
739   394   s394 d394 x   0     433    97   15       -29.2406947890819 
740   395   s395 d395 x   0     442    97   19       -20.2406947890819 
740   396   s396 d396 x   0     442    97   19       -20.2406947890819 
741   397   s397 d397 x   0     433    97   15       -29.2406947890819 
741   398   s398 d398 x   1     433    97   15       -29.2406947890819 
741   399   s399 d399 x   0     433    97   15       -29.2406947890819 
742   400   s400 d400 x   0     430    97   14       -32.2406947890819 
742   401   s401 d401 x   0     430    97   14       -32.2406947890819 
743   402   s402 d402 x   2     450    97   25       -12.2406947890819 
743   403   s403 d403 x   0     450    97   25       -12.2406947890819 
;;;;
run
set locus logTICKS qua
let logTICKS = log10(TICKS+1)
set iter 1000
set chains 4
fpm TICKS c nqtl 0 lik poisson cov YEAR + HEIGHT
fpm TICKS c s nqtl 0 lik poisson cov YEAR + HEIGHT

echo ~~~~~~~~~~~~~~~~~~~~~~
echo Same model using glmer
echo ~~~~~~~~~~~~~~~~~~~~~~
echo Random effects:
echo  Groups Name        Variance Std.Dev.
echo  INDEX  (Intercept) 0.2954   0.5435  
echo  BROOD  (Intercept) 0.8254   0.9085  
echo Number of obs: 403, groups:  INDEX, 403; BROOD, 118
echo 
echo Fixed effects:
echo              Estimate Std. Error z value Pr(>|z|)    
echo (Intercept)  0.409987   0.190356   2.154 0.031256 *  
echo YEAR96       1.148917   0.246303   4.665 3.09e-06 ***
echo YEAR97      -0.990614   0.273460  -3.623 0.000292 ***
echo cHEIGHT     -0.024049   0.003045  -7.898 2.83e-15 ***


