\[ %% % Add your macros here; they'll be included in pdf and html output. %% \newcommand{\R}{\mathbb{R}} % reals \newcommand{\E}{\mathbb{E}} % expectation \renewcommand{\P}{\mathbb{P}} % probability \DeclareMathOperator{\logit}{logit} \DeclareMathOperator{\logistic}{logistic} \DeclareMathOperator{\SE}{SE} \DeclareMathOperator{\sd}{sd} \DeclareMathOperator{\var}{var} \DeclareMathOperator{\cov}{cov} \DeclareMathOperator{\cor}{cor} \DeclareMathOperator{\Normal}{Normal} \DeclareMathOperator{\MVN}{MVN} \DeclareMathOperator{\LogNormal}{logNormal} \DeclareMathOperator{\Poisson}{Poisson} \DeclareMathOperator{\Beta}{Beta} \DeclareMathOperator{\Binom}{Binomial} \DeclareMathOperator{\Gam}{Gamma} \DeclareMathOperator{\Exp}{Exponential} \DeclareMathOperator{\Cauchy}{Cauchy} \DeclareMathOperator{\Unif}{Unif} \DeclareMathOperator{\Dirichlet}{Dirichlet} \DeclareMathOperator{\Wishart}{Wishart} \DeclareMathOperator{\StudentsT}{StudentsT} \DeclareMathOperator{\Weibull}{Weibull} \newcommand{\given}{\;\vert\;} \]

Multivariate responses

Peter Ralph

1 February 2021 – Advanced Biological Statistics

When you’ve got more than one response variable

How often do you really just measure one thing?

Even if you’re interested in just one thing, how often is there one way to measure it?

Examples

  • how oral and axillary temperature depend on age and treatment

  • how gene expression in four tissue types varies with time of day, sex, and metabolic state

  • how jaw length and eye size vary with age and habitat

Things we might want to know

  • How do the responses depend on the predictors?

  • How do the responses covary after accounting for the predictors?

\(\ast\) “predictor” = “explanatory variable”

Multivariate model

Suppose the responses, \((X_1, \ldots, X_k)\) depend on the predictors \((Y_1, \ldots, Y_m)\).

\[\begin{aligned} \text{(response)}_1 &= \text{(intercept)}_1 + \text{(predictors)}_1 + \text{(noise)}_1 \\ \vdots & \qquad \vdots \\ \text{(response)}_k &= \text{(intercept)}_k + \text{(predictors)}_k + \text{(noise)}_k \end{aligned}\]

In symbols, \[\begin{aligned} X_1 &= \beta_{0,1} + \beta_{1,1} Y_1 + \cdots + \beta_{m,1} Y_m + \epsilon_1 \\ \vdots & \qquad \vdots \\ X_k &= \beta_{0,k} + \beta_{1,k} Y_1 + \cdots + \beta_{m,k} Y_m + \epsilon_k \end{aligned}\]

and the noise is \[\begin{aligned} \epsilon &\sim \MVN(0, \Sigma) . \end{aligned}\]

(interlude on the multivariate normal distribution)

Baby birds

https://doi.org/10.1086/662172

Incubation Temperature Affects Growth and Energy Metabolism in Blue Tit Nestlings. Andreas Nord and Jan-Åke Nilsson The American Naturalist 2011 178:5, 639-651

Their data, at ../Datasets/Incubation_temperature/nestling_biometric_data.txt:

nestlings <- read.table("../Datasets/Incubation_temperature/nestling_biometric_data.txt", header=TRUE, stringsAsFactors=TRUE)
nestlings$treat <- factor(nestlings$treat, levels=c("L", "M", "H"))  # enforce a sensible ordering
##            id nest treat mass6 tars6 mass14 tars14 wing14 laydate year
## 1     CK59989  690     H   5.1 11.10   10.8  18.20   41.0      17 2010
## 2     CK59990  690     H   5.6 12.30   11.2  19.00   42.0      17 2010
## 3     CK59991  690     H   5.8 12.00   11.6  18.80   43.0      17 2010
## 4     CK59992  690     H   5.9 12.40   12.0  19.30   42.5      17 2010
## 5     CK59993  690     H   5.6 11.30   12.6  19.40   42.0      17 2010
## 6     CK59994  690     H   5.7 11.20   11.6  18.40   42.0      17 2010
## 7     CK59995  690     H   5.4 11.60   11.1  18.50   42.0      17 2010
## 8     CK59996  690     H   5.7 12.30   11.8  18.90   41.5      17 2010
## 9     CK59997  690     H   5.3 11.20   11.7  18.60   39.5      17 2010
## 10    CK59998  690     H   4.6 10.20   12.2  18.80   41.0      17 2010
## 11    CK59979  623     M   5.3 11.30   11.2  19.20   41.0      17 2010
## 12    CK59980  623     M   4.7 10.50    9.9  18.40   38.5      17 2010
## 13    CK59981  623     M   5.2 11.00    9.4  18.00   38.5      17 2010
## 14    CK59982  623     M   4.6  9.80    9.7  17.40   38.5      17 2010
## 15    CK59983  623     M   5.1 10.90    9.7  18.40   38.0      17 2010
## 16    CK59984  623     M   4.4  9.80    9.3  17.70   37.5      17 2010
## 17    CK59985  623     M   5.3 11.20   11.3  18.80   40.0      17 2010
## 18    CK59810  623     M   3.7  9.20    7.2  15.50   28.5      17 2010
## 19       <NA>  623     M   3.0  8.90     NA     NA     NA      17 2010
## 20    CK59986  867     L   4.8  9.80   10.8  17.70   38.5      17 2010
## 21    CK59987  867     L   6.2 11.10   12.1  18.10   43.5      17 2010
## 22    CK59988  867     L   5.6 10.50   11.2  17.50   41.5      17 2010
## 23    CK59812  867     L   3.8  9.80    9.4  17.00   35.5      17 2010
## 24    CK59813   H8     M   5.3 10.80   10.7  18.40   40.5      18 2010
## 25    CK59814   H8     M   4.4 10.10   10.2  18.60   37.0      18 2010
## 26    CK59815   H8     M   4.0  9.30   10.3  18.00   34.5      18 2010
## 27    CK59816   H8     M   5.0 10.60   10.2  18.50   39.5      18 2010
## 28    CK59817   H8     M   5.1 10.60   10.5  18.30   40.0      18 2010
## 29    CL77375   H8     M   2.6    NA    9.2  17.30   31.0      18 2010
## 30    CK59961  515     M   3.8  9.50   10.6  17.70   34.5      16 2010
## 31    CK59962  515     M   3.8  9.10     NA     NA     NA      16 2010
## 32    CK59963  515     M   4.8 10.20   11.5  18.60   37.5      16 2010
## 33    CK59964  515     M   4.3  9.60   11.4  17.60   35.0      16 2010
## 34       <NA>  515     M   2.6  8.20     NA     NA     NA      16 2010
## 35    CK59806  161     H   5.4 10.90   11.0  18.60   41.0      19 2010
## 36    CK59807  161     H   4.7 10.00   10.1  17.80   39.5      19 2010
## 37    CK59808  161     H   5.6 10.60   11.8  18.80   42.0      19 2010
## 38    CK59809  161     H   4.9 10.00   11.1  18.50   40.0      19 2010
## 39    CL77376  161     H   3.2  8.20   10.2  17.30   32.5      19 2010
## 40       <NA>  161     H   2.9  8.40     NA     NA     NA      19 2010
## 41    CK59999   82     H   3.4  9.30     NA     NA     NA      18 2010
## 42    CK59801   82     H   4.2 10.00     NA     NA     NA      18 2010
## 43    CK59802   82     H   4.1  9.70     NA     NA     NA      18 2010
## 44    CK59803   82     H   4.3 10.10     NA     NA     NA      18 2010
## 45    CK59804   82     H   4.3 10.10     NA     NA     NA      18 2010
## 46    CK59805   82     H   4.2 10.70     NA     NA     NA      18 2010
## 47    CK59818   82     H   3.4  9.30     NA     NA     NA      18 2010
## 48       <NA>   82     H   2.7  8.60     NA     NA     NA      18 2010
## 49       <NA>   82     H   2.3  8.50     NA     NA     NA      18 2010
## 50    CK59819  119     M   5.4 11.40   11.1  18.40   41.0      19 2010
## 51    CK59820  119     M   6.1 11.90   10.8  18.20   41.5      19 2010
## 52    CK59821  119     M   4.7 10.40   10.3  17.50   38.5      19 2010
## 53    CK59822  119     M   6.0 11.60   11.6  19.10   41.5      19 2010
## 54    CK59823  119     M   5.0 10.30   11.4  18.00   39.5      19 2010
## 55    CK59824  119     M   5.7 11.10   11.3  18.00   42.0      19 2010
## 56    CK59825  510     L   6.3 11.30   10.8  18.40   41.5      19 2010
## 57    CK59826  510     L   5.7 10.80   11.1  18.20   41.0      19 2010
## 58    CK59827  510     L   6.1 11.20   10.9  18.40   41.0      19 2010
## 59    CK59828  510     L   5.9 10.60   10.3  17.40   40.0      19 2010
## 60          1  694     M   5.3 10.30   11.6  17.80   43.0      20 2010
## 61          2  694     M   6.4 11.70   11.1  18.10   39.0      20 2010
## 62          3  694     M   5.9 11.20   11.3  18.20   43.0      20 2010
## 63          4  694     M   4.3  9.30   11.1  18.70   41.0      20 2010
## 64          5  694     M   6.2 11.40   11.6  18.40   43.0      20 2010
## 65          6  694     M   6.1 11.30   10.5  18.30   41.5      20 2010
## 66          7  694     M   4.8 10.10   11.2  18.50   41.5      20 2010
## 67          8  694     M   5.8 11.10   11.5  19.30   44.5      20 2010
## 68    CK59850   61     M   3.7  9.00     NA     NA     NA      18 2010
## 69    CK59851   61     M   4.3  9.90     NA     NA     NA      18 2010
## 70    CK59852   61     M   5.1 10.40     NA     NA     NA      18 2010
## 71    CK59853   61     M   4.3  9.70     NA     NA     NA      18 2010
## 72    CK59854   61     M   3.0  9.20     NA     NA     NA      18 2010
## 73    CK59855   61     M   3.9  9.40     NA     NA     NA      18 2010
## 74    CK59856   61     M   4.8 10.30     NA     NA     NA      18 2010
## 75    CK59857  545     H   4.2  9.20   11.3  18.70   37.5      21 2010
## 76    CK59858  545     H   4.5 10.00   11.2  18.30   37.5      21 2010
## 77    CK59859  545     H   4.9 10.40   11.7  18.60   39.0      21 2010
## 78    CK59860  545     H   3.7  9.10   10.7  18.30   36.0      21 2010
## 79    CK59861  545     H   3.7  8.90   10.6  18.70   36.0      21 2010
## 80    CK59862  545     H   4.3  9.30   11.7  18.40   38.0      21 2010
## 81    CK59863  545     H   4.8 10.10   11.0  18.90   39.0      21 2010
## 82    CK59864  223     H   5.2 10.50    9.7  18.10   39.5      18 2010
## 83    CK59865  223     H   5.8 11.80   10.3  18.20   41.0      18 2010
## 84    CK59866  223     H   5.6 10.90    9.5  18.40   40.0      18 2010
## 85    CK59867  223     H   5.8 11.60    9.6  18.20   41.0      18 2010
## 86    CK59868  223     H   6.0 11.50    9.7  18.40   40.0      18 2010
## 87    CK59869  223     H   5.8 11.40   10.3  18.80   40.5      18 2010
## 88    CK59870  223     H   5.7 11.30   10.0  18.30   40.0      18 2010
## 89    CK59871  223     H   4.9  9.60   10.0  18.70   38.5      18 2010
## 90    CK59872  223     H   4.8 10.40    9.7  17.60   39.5      18 2010
## 91    CK59890  806     L   5.4 10.30     NA     NA     NA      20 2010
## 92    CK59891  806     L   4.8  9.80     NA     NA     NA      20 2010
## 93    CK59892  806     L   5.3 10.14     NA     NA     NA      20 2010
## 94    CK59893  806     L   3.3  7.90     NA     NA     NA      20 2010
## 95    CK59894  806     L   5.2 10.90     NA     NA     NA      20 2010
## 96    CK59895  806     L   4.0  9.40     NA     NA     NA      20 2010
## 97    CK59896  806     L   4.9 10.10     NA     NA     NA      20 2010
## 98    CK59897  284     M   6.0 10.90     NA     NA     NA      20 2010
## 99    CK59898  284     M   6.8 12.00     NA     NA     NA      20 2010
## 100   CK59899  284     M   6.0 11.10     NA     NA     NA      20 2010
## 101   CL77300  284     M   5.2 10.20     NA     NA     NA      20 2010
## 102   CL77301  284     M   5.2 10.20     NA     NA     NA      20 2010
## 103   CL77302  284     M   6.4 11.80     NA     NA     NA      20 2010
## 104   CL77303  284     M   5.4 10.60     NA     NA     NA      20 2010
## 105      <NA>  284     M   2.8  7.60     NA     NA     NA      20 2010
## 106   CL77304   68     M   5.1 10.50    9.4  17.60   37.0      17 2010
## 107   CL77305   68     M   4.8  9.80    9.3  17.70   36.0      17 2010
## 108   CL77306   68     M   5.7 10.70   10.8  18.50   40.5      17 2010
## 109   CL77307   68     M   3.1  8.40     NA     NA     NA      17 2010
## 110   CL77308   68     M   5.1 10.60   10.5  17.30   40.0      17 2010
## 111   CL77309   68     M   4.2  9.60     NA     NA     NA      17 2010
## 112   CL77310   68     M   4.1  9.70    9.1  17.50   34.0      17 2010
## 113   CL77311   68     M   5.4 10.30     NA     NA     NA      17 2010
## 114   CL77312   68     M   5.0 10.50   10.7  17.60   37.0      17 2010
## 115   CL77313   68     M   4.6 10.50    9.3  17.40   36.0      17 2010
## 116   CL77314   68     M   5.2 10.80    9.7  17.80   39.0      17 2010
## 117      <NA>   68     M   4.1    NA     NA     NA     NA      17 2010
## 118   CL77315   60     L   4.4  9.80     NA     NA     NA      18 2010
## 119   CL77316   60     L   5.1 10.40     NA     NA     NA      18 2010
## 120   CL77317   60     L   4.7 10.40     NA     NA     NA      18 2010
## 121   CL77318   60     L   4.6 10.20     NA     NA     NA      18 2010
## 122   CL77319   60     L   5.3 10.20     NA     NA     NA      18 2010
## 123   CL77320  637     M   5.2 10.20   10.6  18.40   39.5      21 2010
## 124   CL77321  637     M   6.3 11.80   11.3  18.70   41.0      21 2010
## 125   CL77322  637     M   6.3 11.70   11.1  18.80   40.0      21 2010
## 126   CL77323  637     M   6.6 11.70   11.6  19.00   41.0      21 2010
## 127   CL77324  637     M   5.2 10.20   10.7  18.60   37.5      21 2010
## 128   CL77325  637     M   4.3  8.20   10.3  18.20   37.0      21 2010
## 129   CL77326  637     M   5.6 10.70   10.6  18.10   39.0      21 2010
## 130   CL77327  637     M   5.8 11.20   11.1  18.90   39.0      21 2010
## 131   CL77328  247     L   6.4 11.10   12.6  17.90   44.0      18 2010
## 132   CL77329  247     L   6.1 11.50   11.1  18.20   43.0      18 2010
## 133   CL77330  247     L   5.7 11.30   11.2  18.20   43.5      18 2010
## 134   CL77331  247     L   5.1  9.70   10.5  17.80   40.5      18 2010
## 135   CL77332  247     L   5.3 10.60   10.7  18.40   41.5      18 2010
## 136   CL77338  365     H   5.2 10.60   11.5  17.70   39.5      21 2010
## 137   CL77339  365     H   4.0  9.50   10.5  17.80   35.5      21 2010
## 138   CL77340  365     H   3.8  9.60   10.7  18.30   36.0      21 2010
## 139   CL77341  365     H   4.5  9.80   10.4  17.50   39.0      21 2010
## 140   CL77342  365     H   3.8  9.70   10.3  17.80   35.0      21 2010
## 141   CL77343  365     H   5.3 11.50   11.7  18.80   42.0      21 2010
## 142   CL77344  365     H   5.9 11.40   13.0  18.70   42.5      21 2010
## 143      <NA>  365     H   2.7  8.00     NA     NA     NA      21 2010
## 144   CL77345  571     M   6.4 11.10   10.6  18.10   40.5      20 2010
## 145   CL77346  571     M   5.4 10.10   11.2  18.30   39.5      20 2010
## 146   CL77347  571     M   5.7 10.90   11.3  18.50   41.0      20 2010
## 147   CL77348  571     M   6.3 11.50   11.5  18.50   41.5      20 2010
## 148   CL77349  571     M   4.0  8.90    9.6  17.30   35.0      20 2010
## 149   CL77350  571     M   6.1 11.20   11.2  17.80   41.0      20 2010
## 150   CL77351  571     M   6.1 10.90   10.6  18.00   41.0      20 2010
## 151   CL77352  571     M   4.8 10.00   11.4  18.50   40.0      20 2010
## 152   CL77353  571     M   3.7  8.20    8.0  16.90   29.0      20 2010
## 153   CL77354  571     M   5.8 11.00   10.7  17.60   40.5      20 2010
## 154   CL77355  191     H   6.2 12.20   11.9  19.40   42.5      24 2010
## 155   CL77356  191     H   6.5 12.20   12.6  19.10   42.0      24 2010
## 156   CL77357  191     H   6.3 12.10   12.1  19.10   42.5      24 2010
## 157   CL77358  191     H   4.3  9.40   11.8  18.90   36.5      24 2010
## 158   CL77359  191     H   7.0 12.50   12.9  18.80   44.0      24 2010
## 159   CL77360  191     H   4.4  9.50   10.8  18.20   37.0      24 2010
## 160   CL77361  191     H   6.0 12.00   12.5  19.20   42.0      24 2010
## 161   CL77362  524     L   6.3 11.60   12.4  19.00   42.0      23 2010
## 162   CL77363  524     L   7.0 11.90   11.8  18.00   44.0      23 2010
## 163   CL77364  524     L   7.6 12.60   12.4  18.60   41.5      23 2010
## 164  1ER79930   89     H   7.3 13.30   12.3  19.10   44.5      27 2010
## 165  1ER79931   89     H   7.9 13.30   12.8  19.40   46.0      27 2010
## 166  1ER79932   89     H   7.0 12.50   12.2  19.10   44.0      27 2010
## 167  1ER79933   89     H   6.4 12.30   12.5  19.00   41.5      27 2010
## 168  1ER79934   89     H   5.6 11.10   10.8  18.40   40.0      27 2010
## 169  1ER79935   89     H   6.7 12.40   11.8  18.70   44.0      27 2010
## 170  1ER79936   89     H   6.7 12.60   11.5  19.00   42.0      27 2010
## 171  1ER79937   89     H   7.8 13.00   11.8  18.70   45.0      27 2010
## 172  1ER79938   89     H   6.6 12.80   11.3  18.90   43.0      27 2010
## 173      <NA>   89     H   2.9  7.50     NA     NA     NA      27 2010
## 174  1ER79939  729     M   6.9 12.40   10.8  19.40   43.0      24 2010
## 175  1ER79940  729     M   7.8 13.20   12.2  19.30   45.0      24 2010
## 176  1ER79941  729     M   7.2 13.20   11.1  19.10   44.0      24 2010
## 177  1ER79942  729     M   5.3 11.40    9.6  18.40   41.0      24 2010
## 178  1ER79943  729     M   7.3 13.00   11.0  18.80   45.0      24 2010
## 179  1ER79944  729     M   7.1 13.00   10.7  18.70   41.0      24 2010
## 180  1ER79945  729     M   7.8 12.90   11.9  18.70   43.5      24 2010
## 181  1ER79946  729     M   5.0 11.20   10.5  18.70   39.5      24 2010
## 182  1ER79947  729     M   5.8 11.80   10.6  18.30   40.0      24 2010
## 183  1ER79948  729     M   7.1 13.30   11.9  19.10   44.5      24 2010
## 184  1ER79949  729     M   6.7 12.70   10.7  18.70   41.0      24 2010
## 185  1ER79951  176     H   5.9 11.20   11.8  19.00   40.0      28 2010
## 186  1ER79952  176     H   6.6 12.20   12.4  19.50   42.0      28 2010
## 187  1ER79953  176     H   6.6 12.50   12.8  19.80   43.0      28 2010
## 188   CL77379  176     H   6.0 12.10   11.8  18.80   35.5      28 2010
## 189  1ER79955  176     H   5.6 11.30   11.7  19.30   42.5      28 2010
## 190  1ER79956  176     H   5.1 10.40   12.0  19.40   39.0      28 2010
## 191  1ER79957  176     H   5.2 10.40   12.5  19.00   39.5      28 2010
## 192      <NA>  176     H   3.8  8.50     NA     NA     NA      28 2010
## 193  1ER79959  695     L   4.6  9.80   11.7  18.60   40.0      27 2010
## 194  1ER79960  695     L   7.8 12.80   11.8  18.50   46.0      27 2010
## 195  1ER79961  695     L   7.4 12.50   12.7  18.30   45.5      27 2010
## 196  1ER79962  695     L   7.3 12.90   11.4  18.80   46.0      27 2010
## 197  1ER79963  695     L   6.2 11.50   12.3  19.40   43.0      27 2010
## 198  1ER79826  695     L   3.7  8.60   11.6  18.50   37.0      27 2010
## 199  1ER79827  695     L   2.1    NA    9.8  17.30   25.5      27 2010
## 200  1ER79966  215     M   6.3 12.10     NA     NA     NA      25 2010
## 201  1ER79967  215     M   6.6 12.30     NA     NA     NA      25 2010
## 202  1ER79968  215     M   5.2 10.40     NA     NA     NA      25 2010
## 203  1ER79969  215     M   6.4 12.00     NA     NA     NA      25 2010
## 204  1ER79970  215     M   6.6 11.90     NA     NA     NA      25 2010
## 205  1ER79971  215     M   6.6 11.70     NA     NA     NA      25 2010
## 206  1ER79972  215     M   6.3 11.80     NA     NA     NA      25 2010
## 207  1ER79973  215     M   5.8 11.50     NA     NA     NA      25 2010
## 208  1ER79974  215     M   5.7 10.20     NA     NA     NA      25 2010
## 209      <NA>  215     M   3.3  8.70     NA     NA     NA      25 2010
## 210  1ER79975   96     L   6.5 11.20   10.4  17.00   43.5      21 2010
## 211  1ER79976   96     L   5.2 10.60   11.6  18.00   40.5      21 2010
## 212  1ER79977   96     L   6.2 11.70   11.6  18.00   45.0      21 2010
## 213  1ER79979   96     L   5.4 11.10    9.8  16.60   43.0      21 2010
## 214  1ER79980   96     L   5.6 10.80   10.2  16.80   42.5      21 2010
## 215  1ER79981   96     L   4.7  9.50     NA     NA     NA      21 2010
## 216  1ER79982  654     H   6.3 12.20   10.9  19.10   41.5      26 2010
## 217  1ER79983  654     H   5.6 11.00   10.8  18.10   40.5      26 2010
## 218  1ER79984  654     H   5.0 10.10   11.0  18.70   38.5      26 2010
## 219  1ER79985  654     H   6.1 11.60   11.2  19.20   41.0      26 2010
## 220  1ER79986  654     H   5.6 11.70   11.0  18.30   40.0      26 2010
## 221  1ER79987  654     H   6.1 11.70   10.3  18.50   41.0      26 2010
## 222  1ER79988  654     H   6.6 12.30   11.1  19.20   43.0      26 2010
## 223  1ER79989  654     H   5.5 11.30   11.2  19.20   40.5      26 2010
## 224  1ER79830  654     H   2.9  8.20    9.3  17.10   31.0      26 2010
## 225  1ER79990  211     M   6.1 11.30   11.6  18.80   42.0      28 2010
## 226  1ER79991  211     M   6.7 12.20   11.4  18.40   44.0      28 2010
## 227  1ER79992  211     M   6.1 11.50   11.7  19.00   42.5      28 2010
## 228  1ER79993  211     M   5.4 10.90   11.4  18.70   42.5      28 2010
## 229  1ER79994  211     M   6.4 12.00   11.5  19.40   44.5      28 2010
## 230  1ER79995  211     M   7.3 13.20   11.9  19.30   45.0      28 2010
## 231  1ER79996  211     M   6.6 11.50   11.6  18.40   42.0      28 2010
## 232  1ER79997  211     M   5.4 10.80   11.6  18.90   41.5      28 2010
## 233  1ER79837  211     M   3.6  8.70   10.8  17.50   35.5      28 2010
## 234  1ER79838  211     M   2.7  7.30   10.0  17.60   30.0      28 2010
## 235  1ER79700  625     M   6.7 12.60    9.2  18.40   41.5      24 2010
## 236  1ER79701  625     M   6.3 12.40   10.9  18.40   41.5      24 2010
## 237  1ER79702  625     M   6.7 12.40   11.4  18.20   41.0      24 2010
## 238  1ER79703  625     M   4.5  9.80     NA     NA     NA      24 2010
## 239  1ER79704  625     M   6.2 11.70   10.3  17.80   42.0      24 2010
## 240  1ER79705  625     M   5.6 11.20   10.1  17.70   39.5      24 2010
## 241  1ER79706  625     M   5.2 10.80   10.6  18.50   38.0      24 2010
## 242  1ER79707  625     M   6.7 12.40   10.5  18.40   43.0      24 2010
## 243  1ER79708  625     M   6.1 11.50   10.6  18.00   43.0      24 2010
## 244  1ER79709  625     M   5.7 10.60   11.1  18.20   41.0      24 2010
## 245  1ER79710  625     M   7.6 12.70   11.6  18.30   43.0      24 2010
## 246  1ER79713  213     H   6.5 12.20   12.1  19.00   42.0      28 2010
## 247  1ER79714  213     H   6.6 12.40   11.9  19.40   42.0      28 2010
## 248  1ER79715  213     H   6.6 12.40   11.7  18.90   40.5      28 2010
## 249  1ER79716  213     H   7.1 12.80   12.3  19.20   41.0      28 2010
## 250  1ER79717  213     H   6.2  9.70   10.8  18.10   38.0      28 2010
## 251  1ER79718  213     H   7.5 13.10   12.4  19.90   43.5      28 2010
## 252  1ER79719  213     H   6.7 12.10   11.7  18.30   42.0      28 2010
## 253  1ER79720  213     H   6.1 12.00   11.6  18.80   39.5      28 2010
## 254  1ER79856  213     H   3.2  8.20   10.5  18.10   33.0      28 2010
## 255  1ER79726  779     H   7.1 12.60   11.7  18.00   45.0      29 2010
## 256  1ER79727  779     H   4.6  9.30   12.5  18.20   44.5      29 2010
## 257  1ER79728  779     H   6.9 12.50   12.5  18.20   44.5      29 2010
## 258  1ER79729  779     H   6.2 11.60   11.5  18.20   41.5      29 2010
## 259  1ER79730  779     H   7.4 12.80   11.8  19.10   44.0      29 2010
## 260  1ER79731  779     H   7.3 12.50   12.3  18.40   42.5      29 2010
## 261  1ER79732  779     H   7.3 12.60   12.4  18.70   42.0      29 2010
## 262  1ER79733  779     H   5.6 11.40   10.9  18.10   38.0      29 2010
## 263  1ER79734  779     H   5.8 10.70   11.5  17.70   39.0      29 2010
## 264  1ER79735  779     H   6.7 12.50   11.5  18.10   45.0      29 2010
## 265  1ER79737  112     H   5.6 11.60     NA     NA     NA      29 2010
## 266  1ER79738  112     H   4.4 10.00     NA     NA     NA      29 2010
## 267  1ER79739  112     H   5.6 11.40     NA     NA     NA      29 2010
## 268  1ER79740  112     H   4.7 10.70     NA     NA     NA      29 2010
## 269  1ER79741  112     H   6.5 12.10     NA     NA     NA      29 2010
## 270  1ER79858  112     H   3.7  9.00     NA     NA     NA      29 2010
## 271  1ER79743  112     H   4.3  9.40     NA     NA     NA      29 2010
## 272  1ER79744  112     H   5.4 10.70     NA     NA     NA      29 2010
## 273  1ER79745  112     H   6.4 12.10     NA     NA     NA      29 2010
## 274  1ER79746  112     H   6.2 12.00     NA     NA     NA      29 2010
## 275      <NA>  112     H   2.3    NA     NA     NA     NA      29 2010
## 276  1ER79747  727     L   5.3 10.40   12.2  18.50   42.0      26 2010
## 277  1ER79748  727     L   5.3 10.40   11.3  18.00   40.0      26 2010
## 278  1ER79749  727     L   4.4  9.20   10.8  18.40   37.0      26 2010
## 279  1ER79750  727     L   5.6 10.50   11.4  18.50   41.5      26 2010
## 280  1ER79752   29     L   7.2 12.40     NA     NA     NA      28 2010
## 281  1ER79753   29     L   6.3 11.60     NA     NA     NA      28 2010
## 282  1ER79754   29     L   4.1  9.20     NA     NA     NA      28 2010
## 283  1ER79756   29     L   6.6 12.10     NA     NA     NA      28 2010
## 284  1ER79757   29     L   6.4 11.90     NA     NA     NA      28 2010
## 285  1ER79857   29     L   3.9  8.80     NA     NA     NA      28 2010
## 286  1ER79759  608     M   5.9 11.10     NA     NA     NA      29 2010
## 287  1ER79760  608     M   5.5 11.40     NA     NA     NA      29 2010
## 288  1ER79761  608     M   6.1 11.20     NA     NA     NA      29 2010
## 289  1ER79762  608     M   5.2 10.40     NA     NA     NA      29 2010
## 290  1ER79763  608     M   4.3  9.60     NA     NA     NA      29 2010
## 291  1ER79764  608     M   6.7 12.00     NA     NA     NA      29 2010
## 292  1ER79765  608     M   5.7 11.40     NA     NA     NA      29 2010
## 293  1ER79766  860     L   5.3 10.30   12.1  18.10   42.0      28 2010
## 294  1ER79767  860     L   5.8  9.80   12.5  18.50   40.5      28 2010
## 295  1ER79768  860     L   5.4 10.60   12.2  18.50   42.5      28 2010
## 296  1ER79769  860     L   5.6 10.80   12.4  18.90   42.5      28 2010
## 297  1ER79770  860     L   5.5 10.40   12.7  18.20   43.0      28 2010
## 298  1ER79771  860     L   5.8 10.80   11.7  17.50   42.5      28 2010
## 299  1ER79772  860     L   5.9 10.90   12.1  18.00   43.0      28 2010
## 300  1ER79773   94     M   6.6 13.20   11.8  19.10   45.0      26 2010
## 301  1ER79774   94     M   6.7 12.70   11.5  18.70   44.5      26 2010
## 302  1ER79775   94     M   6.4 12.80   11.2  19.10   44.0      26 2010
## 303  1ER79776   94     M   6.4 13.00   11.0  19.60   43.0      26 2010
## 304  1ER79777   94     M   6.5 12.40   11.1  18.50   39.5      26 2010
## 305  1ER79778   94     M   6.6 12.30   11.0  18.60   42.5      26 2010
## 306  1ER79779   94     M   6.3 12.40   10.1  18.50   45.0      26 2010
## 307  1ER79780   94     M   5.4 11.80   10.6  18.70   43.0      26 2010
## 308  1ER79781   94     M   6.0 12.50   10.3  18.70   44.0      26 2010
## 309  1ER79782   94     M   6.0 12.60   10.0  18.10   44.0      26 2010
## 310  1ER79783  648     M   6.2 11.60   11.5  18.70   39.5      30 2010
## 311  1ER79784  648     M   6.3 11.40   11.8  18.60   41.0      30 2010
## 312  1ER79785  648     M   6.3 11.50   10.9  17.60   39.5      30 2010
## 313  1ER79786  648     M   6.0 10.70   11.5  17.80   39.0      30 2010
## 314  1ER79787  648     M   4.7  9.70     NA     NA     NA      30 2010
## 315  1ER79788  648     M   6.1 11.40   11.3  17.90   39.5      30 2010
## 316  1ER79789  648     M   6.1 11.70   11.7  18.60   40.0      30 2010
## 317  1ER79790  648     M   4.8  9.80   10.3  17.40   36.0      30 2010
## 318      <NA>  648     M   3.0  7.80     NA     NA     NA      30 2010
## 319  1ER79791  689     M   5.3 10.70   11.5  19.40   40.0      32 2010
## 320  1ER79792  689     M   4.2  9.80   10.2  18.10   36.5      32 2010
## 321  1ER79793  689     M   6.3 12.50   11.7  19.60   43.0      32 2010
## 322  1ER79794  689     M   5.3 11.40   10.6  19.50   40.0      32 2010
## 323  1ER79795  689     M   5.4 11.20   10.8  19.30   40.5      32 2010
## 324  1ER79796  689     M   4.8 10.40    9.7  18.70   38.0      32 2010
## 325  1ER79797  689     M   3.9  9.10    9.6  17.90   35.5      32 2010
## 326  1ER79798  689     M   4.6 10.20   10.5  18.80   38.5      32 2010
## 327  1ER79808  364     H   7.1 12.40   11.4  18.80   42.5      30 2010
## 328  1ER79809  364     H   6.9 12.30   12.0  19.10   43.5      30 2010
## 329  1ER79810  364     H   7.0 12.30   11.8  19.50   43.5      30 2010
## 330  1ER79811  364     H   6.2 11.70   11.6  19.50   42.0      30 2010
## 331  1ER79812  364     H   4.2  9.40     NA     NA     NA      30 2010
## 332  1ER79813  364     H   7.1 12.10   12.7  19.20   43.5      30 2010
## 333  1ER79814  364     H   6.9 12.20   11.1  18.20   43.0      30 2010
## 334  1ER79815  364     H   7.3 12.10   11.5  19.20   44.0      30 2010
## 335  1ER79816  364     H   6.8 12.30   12.3  19.60   41.0      30 2010
## 336  1ER79817  364     H   6.9 12.30   11.6  19.00   44.0      30 2010
## 337  1ER79818  506     L   5.5 11.50   12.3  19.50   43.5      31 2010
## 338  1ER79819  506     L   6.1 12.00   12.8  19.80   43.0      31 2010
## 339  1ER79820  506     L   6.7 12.70   12.3  19.40   45.0      31 2010
## 340  1ER79821  506     L   6.2 11.90   12.7  19.70   44.0      31 2010
## 341  1ER79822  506     L   6.2 12.30   11.8  19.20   44.0      31 2010
## 342  1ER79823  506     L   5.2 11.50   11.3  18.30   41.0      31 2010
## 343  1ER17601  838     H   6.4 11.60   11.2  18.60   41.0      16 2009
## 344  1ER17602  838     H   4.9  9.80   10.3  17.70   39.5      16 2009
## 345  1ER17603  838     H   6.6 11.30   11.2  18.20   40.5      16 2009
## 346  1ER17604  838     H   6.9 11.30   11.1  17.90   43.0      16 2009
## 347  1ER17605  838     H   6.2 10.30   10.1  17.20   39.0      16 2009
## 348  1ER17606  838     H   4.1  8.90    9.4  17.60   29.0      16 2009
## 349  1ER17607  838     H   6.1 10.70   11.0  17.90   41.0      16 2009
## 350  1ER17608  838     H   6.4 11.70   10.6  17.60   43.0      16 2009
## 351      <NA>  838     H   3.2  8.20     NA     NA     NA      16 2009
## 352  1EP51290  617     H   7.1 11.50   12.2  18.70   45.0      17 2009
## 353  1EP51291  617     H   6.3 10.50   12.9  18.60   42.5      17 2009
## 354  1EP51292  617     H   6.0 10.20   12.3  18.60   42.0      17 2009
## 355  1EP51294  247     H   6.0 10.80   11.4  18.10   41.0      15 2009
## 356  1EP51295  247     H   5.8 10.90   11.3  18.40   43.0      15 2009
## 357  1EP51296  247     H   5.7 10.30   11.1  18.50   42.0      15 2009
## 358  1EP51297  247     H   5.8  9.80   12.6  18.60   40.0      15 2009
## 359  1EP51298  247     H   4.6  9.10   10.8  18.40   39.0      15 2009
## 360  1EP51299  247     H   5.9 10.30   11.2  17.80   41.0      15 2009
## 361  1ER17600  247     H   5.8 10.70   11.0  17.80   42.0      15 2009
## 362  1ER17610   21     M   6.2 11.30   12.1  19.10   44.5      16 2009
## 363  1ER17611   21     M   6.1 10.80   12.1  19.00   41.5      16 2009
## 364  1ER17612   21     M   5.8 10.60   11.5  18.40   44.0      16 2009
## 365  1ER17613   21     M   5.0  9.90   11.8  18.70   42.5      16 2009
## 366  1ER17614   21     M   6.4 11.60   12.4  18.70   44.0      16 2009
## 367  1ER17615   21     M   6.8 11.60   12.4  18.80   44.0      16 2009
## 368  1ER17616   21     M   6.9 11.40   12.2  18.50   44.0      16 2009
## 369  1ER17617   21     M   4.6  8.50   11.7  17.70   39.5      16 2009
## 370  1ER17618   21     M   5.4 10.40   11.8  18.10   43.0      16 2009
## 371  1ER17620 7001     L   5.2 10.10   13.1  18.90   40.5      18 2009
## 372  1ER17621 7001     L   4.3  9.40   12.5  18.30   34.0      18 2009
## 373  1ER17622 7001     L   5.0  9.60   12.4  18.40   39.0      18 2009
## 374  1ER17623 7001     L   4.9  9.20   12.9  18.20   39.0      18 2009
## 375  1ER17624 7001     L   4.6  9.30   12.3  18.10   38.5      18 2009
## 376      <NA> 7001     L   4.0  8.70     NA     NA     NA      18 2009
## 377      <NA> 7001     L   3.7  8.40     NA     NA     NA      18 2009
## 378      <NA> 7001     L   1.3    NA     NA     NA     NA      18 2009
## 379  1ER17627  875     L   6.6 11.10   10.4  18.30   42.5      14 2009
## 380  1ER17628  875     L   5.7 10.80   10.2  18.40   40.5      14 2009
## 381  1ER17629  875     L   5.6  9.80   10.9  18.90   39.5      14 2009
## 382  1ER17630  875     L   5.2  9.60   10.5  17.70   38.5      14 2009
## 383  1ER17631  875     L   5.0  9.60   10.4  18.30   40.5      14 2009
## 384  1ER17632  875     L   4.6  9.30   10.3  17.90   38.0      14 2009
## 385  1ER17633  875     L   5.6  9.80   10.3  17.80   40.0      14 2009
## 386  1ER17636  674     H   3.9  8.40     NA     NA     NA      18 2009
## 387  1ER17637  674     H   6.1 10.90     NA     NA     NA      18 2009
## 388  1ER17638  674     H   6.0 10.60     NA     NA     NA      18 2009
## 389  1ER17937  674     H   3.8  8.40     NA     NA     NA      18 2009
## 390  1ER17640  674     H   5.8 10.70     NA     NA     NA      18 2009
## 391  1ER17641  674     H   6.3 11.10     NA     NA     NA      18 2009
## 392  1ER17642  674     H   4.8  9.40     NA     NA     NA      18 2009
## 393  1ER17643  674     H   5.4  9.90     NA     NA     NA      18 2009
## 394  1ER17644  674     H   5.7 10.70     NA     NA     NA      18 2009
## 395  1ER17645  674     H   5.3 10.20     NA     NA     NA      18 2009
## 396  1ER17646  674     H   6.2 10.80     NA     NA     NA      18 2009
## 397  1ER17649    4     M   6.2 10.90   11.3  18.30   42.5      18 2009
## 398  1ER17650    4     M   5.9 10.70   11.3  18.00   42.0      18 2009
## 399  1ER17651    4     M   7.2 11.00   12.7  18.30   42.5      18 2009
## 400  1ER17652    4     M   5.8 10.30   11.7  18.40   40.0      18 2009
## 401  1ER17653    4     M   6.1 11.10   12.4  19.30   43.0      18 2009
## 402  1ER17654    4     M   5.8 10.10   12.2  18.10   42.0      18 2009
## 403  1ER17655    4     M   6.3 11.00   12.2  19.20   42.0      18 2009
## 404  1ER17936    4     M   3.5  7.90   11.0  17.20   35.0      18 2009
## 405  1ER17671  836     L   7.2 12.20   11.9  18.55   44.5      17 2009
## 406  1ER17672  836     L   6.0 11.60   11.7  18.72   46.0      17 2009
## 407  1ER17673  836     L   6.7 12.50   11.5  18.52   47.0      17 2009
## 408  1ER17674  836     L   6.7 12.30   11.8  18.49   45.5      17 2009
## 409  1ER17675  836     L   7.3 12.10   12.8  18.81   47.0      17 2009
## 410  1ER17676  836     L   7.0 12.70   11.7  18.97   45.0      17 2009
## 411  1ER17677  836     L   5.7 10.90   11.6  18.93   42.0      17 2009
## 412  1ER17678  836     L   6.8 12.20   11.9  19.20   46.5      17 2009
## 413  1ER17679  836     L   5.1 10.10   11.4  19.25   43.0      17 2009
## 414  1ER17200  753     M   5.6 10.20     NA     NA     NA      18 2009
## 415  1ER17201  753     M   5.9 10.30     NA     NA     NA      18 2009
## 416  1ER17202  753     M   6.7 10.90     NA     NA     NA      18 2009
## 417  1ER17203  753     M   6.3 11.00     NA     NA     NA      18 2009
## 418  1ER17691  548     M   6.9 11.90   11.6  18.40   45.5      19 2009
## 419  1ER17692  548     M   5.6 10.80   11.4  17.70   43.0      19 2009
## 420  1ER17693  548     M   4.8  9.80   11.4  17.70   40.0      19 2009
## 421  1ER17694  548     M   6.3 11.90   11.3  18.30   44.0      19 2009
## 422  1ER17695  548     M   5.7 11.20   11.3  17.40   43.0      19 2009
## 423  1ER17696  548     M   6.8 12.00   12.5  17.90   45.0      19 2009
## 424  1ER17697  548     M   6.0 11.40   12.6  18.60   45.5      19 2009
## 425  1ER17698  548     M   5.7 10.90   12.1  18.40   43.5      19 2009
## 426  1ER17699  548     M   6.4 11.50   11.7  17.70   45.5      19 2009
## 427  1ER17969  548     M   4.2  8.00   11.2  17.10   39.5      19 2009
## 428  1ER17214  277     L   5.6 10.70   10.3  18.00   39.5      17 2009
## 429  1ER17215  277     L   5.5 10.60   10.7  17.90   39.0      17 2009
## 430  1ER17216  277     L   6.2 11.50   11.0  18.40   42.0      17 2009
## 431  1ER17217  277     L   6.7 12.00   11.7  18.60   42.0      17 2009
## 432  1ER17218  277     L   6.4 11.40   11.2  17.90   42.0      17 2009
## 433  1ER17219  277     L   5.3  9.30   10.3  17.40   37.0      17 2009
## 434  1ER17220  277     L   6.1 11.10   10.9  18.20   41.5      17 2009
## 435  1ER17970  277     L   3.2  7.70    6.7  15.10   26.5      17 2009
## 436  1ER17682  207     M   4.6  9.60   12.5  19.30   43.0      19 2009
## 437  1ER17683  207     M   6.2 11.30   11.5  18.40   43.0      19 2009
## 438  1ER17684  207     M   5.7 10.40   12.4  19.00   42.5      19 2009
## 439  1ER17685  207     M   6.2 11.90   11.6  18.40   43.5      19 2009
## 440  1ER17686  207     M   5.7 10.70   12.0  19.00   44.0      19 2009
## 441  1ER17687  207     M   6.7 11.80   12.5  18.80   45.0      19 2009
## 442  1ER17688  207     M   4.2  9.00   10.8  18.10   37.0      19 2009
## 443  1ER17689  207     M   7.2 12.30   12.5  19.30   45.5      19 2009
## 444  1ER17690  207     M   6.7 12.20   11.7  19.00   45.0      19 2009
## 445  1ER17206 7009     H   7.3 11.50   12.4  18.20   45.0      18 2009
## 446  1ER17207 7009     H   7.7 12.00   12.6  18.80   45.5      18 2009
## 447  1ER17208 7009     H   5.9 11.30   11.9  18.50   43.0      18 2009
## 448  1ER17209 7009     H   5.3  9.40   11.9  17.80   41.0      18 2009
## 449  1ER17210 7009     H   7.8 11.80   12.2  18.20   44.0      18 2009
## 450  1ER17211 7009     H   6.3  8.80   12.6  18.30   42.0      18 2009
## 451  1ER17212 7009     H   6.9 10.90   11.4  17.30   42.5      18 2009
## 452  1ER17663   H8     L   5.0 10.30     NA     NA     NA      17 2009
## 453  1ER17664   H8     L   5.6 10.60   10.6  18.60   44.0      17 2009
## 454  1ER17665   H8     L   6.6 12.20   11.9  19.40   45.0      17 2009
## 455  1ER17666   H8     L   5.0  9.80   11.5  18.70   42.0      17 2009
## 456  1ER17667   H8     L   6.2 11.40   11.2  19.10   43.0      17 2009
## 457  1ER17668   H8     L   6.4 12.20   11.1  19.20   45.0      17 2009
## 458  1ER17669   H8     L   5.4 10.50   11.1  18.40   41.5      17 2009
## 459  1ER17968   H8     L   3.8  7.90   10.4  16.50   38.0      17 2009
## 460  1ER17221  730     H   5.3 10.10   11.4  18.10   40.0      16 2009
## 461  1ER17222  730     H   5.7 10.50   12.4  18.30   42.5      16 2009
## 462  1ER17223  730     H   5.5 10.50   10.7  17.90   41.0      16 2009
## 463  1ER17227  693     L   7.7 12.20   13.4  19.20   47.5      20 2009
## 464  1ER17228  693     L   6.3 11.20   11.8  18.30   42.5      20 2009
## 465  1ER17229  693     L   5.6  9.80   12.3  17.70   42.0      20 2009
## 466  1ER17230  693     L   4.7  9.10   12.3  18.80   40.0      20 2009
## 467  1ER17231  693     L   5.8 10.90   10.8  18.70   42.0      20 2009
## 468  1ER17232  693     L   6.8 11.70   11.8  18.60   45.5      20 2009
## 469  1ER17233  693     L   7.1 12.00   11.9  18.50   46.0      20 2009
## 470  1ER17234  693     L   4.8  9.20   11.7  17.50   40.5      20 2009
## 471  1ER17235  606     L   5.7 10.60   11.2  18.10   43.0      19 2009
## 472  1ER17236  606     L   6.1 11.30   10.6  17.90   45.0      19 2009
## 473  1ER17237  606     L   4.2  8.90   11.5  18.50   39.5      19 2009
## 474  1ER17238  606     L   5.6 10.70   11.9  18.60   44.0      19 2009
## 475  1ER17239  606     L   5.6 11.10   11.5  18.70   45.0      19 2009
## 476  1ER17240  606     L   4.7 10.00   11.2  18.70   42.0      19 2009
## 477  1ER17241  606     L   6.4 11.40   12.0  18.50   44.5      19 2009
## 478  1ER17242  606     L   6.1 11.20   11.4  18.40   45.0      19 2009
## 479  1ER17243  606     L   6.0 11.30   11.6  18.10   44.5      19 2009
## 480  1ER17244  606     L   5.6 10.60   10.3  18.20   43.0      19 2009
## 481  1ER17245  606     L   5.4 10.70   11.1  18.30   43.5      19 2009
## 482  1ER17246  533     L   4.3  9.10   10.7  17.10   38.0      18 2009
## 483  1ER17247  533     L   5.9 11.00   12.0  17.70   42.0      18 2009
## 484  1ER17248  533     L   5.4 10.80   11.9  17.60   41.5      18 2009
## 485  1ER17249  533     L   6.4 11.70   12.6  17.80   43.0      18 2009
## 486  1ER17250  533     L   6.7 12.10   12.4  18.00   45.0      18 2009
## 487  1ER17251  533     L   5.7 11.20   11.8  18.00   43.0      18 2009
## 488  1ER17252  533     L   6.6 11.80   12.6  18.10   44.0      18 2009
## 489  1ER17253  533     L   7.4 12.50   12.9  18.50   46.0      18 2009
## 490  1ER17254  533     L   6.2 11.70   12.1  18.50   43.5      18 2009
## 491  1ER17255  533     L   6.7 12.30   13.1  18.90   45.0      18 2009
## 492      <NA>  533     L   1.9    NA     NA     NA     NA      18 2009
## 493  1ER17256 7004     M   7.2 12.10   11.6  17.80   42.0      19 2009
## 494  1ER17257 7004     M   6.9 12.30   11.8  17.90   40.5      19 2009
## 495  1ER17258 7004     M   6.8 12.30   12.0  18.60   42.0      19 2009
## 496  1ER17259 7004     M   5.4 10.80   10.8  18.10   38.0      19 2009
## 497  1ER17260 7004     M   6.1 11.70   10.8  17.90   38.5      19 2009
## 498  1ER17261 7004     M   7.4 12.70   12.2  18.80   43.0      19 2009
## 499  1ER17262 7004     M   6.7 12.10   12.2  18.40   43.0      19 2009
## 500  1ER17275  863     H   7.3 12.80   12.7  19.60   46.5      22 2009
## 501  1ER17276  863     H   6.3 12.00   12.7  19.30   45.5      22 2009
## 502  1ER17277  863     H   7.1 12.20   11.3  18.50   42.5      22 2009
## 503  1ER17278  863     H   7.3 12.80   12.2  19.00   47.0      22 2009
## 504  1ER17279  863     H   6.9 12.10   13.3  20.00   43.0      22 2009
## 505  1ER17280  863     H   7.4 12.50   12.2  19.70   46.0      22 2009
## 506  1ER17281  863     H   5.6 10.00   12.5  18.80   40.0      22 2009
## 507  1ER17282  863     H   5.1  9.80   11.4  18.50   40.0      22 2009
## 508  1ER17283  863     H   7.1 12.30   12.1  18.70   44.0      22 2009
## 509  1ER17264   78     M   4.6 10.30     NA     NA     NA      21 2009
## 510  1ER17265   78     M   5.3 11.10     NA     NA     NA      21 2009
## 511  1ER17266   78     M   5.2 10.70     NA     NA     NA      21 2009
## 512  1ER17267   78     M   6.6 11.70     NA     NA     NA      21 2009
## 513  1ER17268   78     M   5.6 11.60     NA     NA     NA      21 2009
## 514  1ER17269   78     M   6.2 11.60     NA     NA     NA      21 2009
## 515  1ER17270   78     M   6.3 11.70     NA     NA     NA      21 2009
## 516  1ER17271   78     M   5.8 11.90     NA     NA     NA      21 2009
## 517  1ER17272   78     M   5.8 11.40     NA     NA     NA      21 2009
## 518      <NA>   78     M   3.2  8.50     NA     NA     NA      21 2009
## 519  1ER17431  771     L   7.0 12.50   12.9  18.70   45.0      24 2009
## 520  1ER17432  771     L   6.3 11.40   12.6  18.30   44.0      24 2009
## 521  1ER17433  771     L   6.3 11.70   11.7  18.00   42.0      24 2009
## 522  1ER17434  771     L   5.7 10.70    9.2  16.40   36.5      24 2009
## 523  1ER17435  771     L   6.9 12.90   12.1  18.60   46.0      24 2009
## 524  1ER17436  771     L   7.4 12.70   13.2  18.70   45.5      24 2009
## 525  1ER17437  771     L   4.6  9.70   13.3  18.40   45.0      24 2009
## 526  1ER79505  771     L   2.4    NA   10.3  16.70   29.0      24 2009
## 527  1ER17407  274     H   5.6 10.20     NA     NA     NA      23 2009
## 528  1ER17408  274     H   5.2  9.90     NA     NA     NA      23 2009
## 529  1ER17409  274     H   6.2 10.90     NA     NA     NA      23 2009
## 530  1ER17410  274     H   5.3 10.00     NA     NA     NA      23 2009
## 531  1ER17411  274     H   5.6 10.10     NA     NA     NA      23 2009
## 532  1ER17412  274     H   5.8 10.30     NA     NA     NA      23 2009
## 533  1ER17413  274     H   6.3 10.80     NA     NA     NA      23 2009
## 534  1ER17414  274     H   4.7  9.30     NA     NA     NA      23 2009
## 535      <NA>  274     H   3.5  7.90     NA     NA     NA      23 2009
## 536  1ER17296  203     M   5.8 11.10   12.4  18.40   40.0      23 2009
## 537  1ER17297  203     M   4.1  8.70   11.2  16.50   35.0      23 2009
## 538  1ER17298  203     M   6.1 11.20   11.8  18.20   43.0      23 2009
## 539  1ER17299  203     M   5.6 10.40   11.7  17.00   41.0      23 2009
## 540  1ER17400  203     M   5.3 10.00   12.2  17.60   42.0      23 2009
## 541  1ER17401  203     M   6.7 11.20   13.2  18.00   44.0      23 2009
## 542  1ER17402  203     M   6.7 11.10   13.6  17.90   44.0      23 2009
## 543  1ER17403  203     M   6.3 11.50   12.7  18.10   43.0      23 2009
## 544  1ER17404  203     M   6.7 11.50   12.7  17.80   43.0      23 2009
## 545  1ER79506  203     M   2.3    NA   10.3  16.30   25.0      23 2009
## 546  1ER17417   97     M   6.2 11.20   12.5  18.40   42.0      24 2009
## 547  1ER17418   97     M   6.6 11.40   13.5  19.20   44.0      24 2009
## 548  1ER17419   97     M   6.5 11.50   12.7  18.50   43.0      24 2009
## 549  1ER17420   97     M   4.8  9.10   12.0  18.30   39.5      24 2009
## 550  1ER17421   97     M   7.0 11.60   13.9  18.90   44.0      24 2009
## 551  1ER17422   97     M   5.3 10.00   13.1  18.80   40.5      24 2009
## 552  1ER17423   97     M   6.4 10.80   13.3  19.10   42.0      24 2009
## 553  1ER17424   97     M   4.0  8.30   12.6  17.70   36.5      24 2009
## 554  1ER79504   97     M   1.9    NA   11.2  16.30   27.5      24 2009
## 555  1ER17425   61     H   7.2 12.90   12.7  19.80   44.0      25 2009
## 556  1ER17426   61     H   5.6 11.70   12.7  20.10   43.5      25 2009
## 557  1ER17427   61     H   6.3 11.50   11.8  18.90   42.0      25 2009
## 558  1ER17428   61     H   7.1 12.30   12.1  19.50   44.5      25 2009
## 559  1ER17429   61     H   4.3  8.70   12.3  19.30   38.0      25 2009
## 560  1ER17430   61     H   6.9 12.00   13.7  20.00   44.0      25 2009
## 561  1ER17287   26     M   5.7 11.30   12.3  18.60   44.5      22 2009
## 562  1ER17288   26     M   5.2 10.30   12.0  18.50   39.5      22 2009
## 563  1ER17289   26     M   6.2 10.90   12.3  18.60   42.0      22 2009
## 564  1ER17290   26     M   6.3 11.30   13.3  19.00   44.0      22 2009
## 565  1ER17291   26     M   4.9  9.50   11.2  17.80   39.0      22 2009
## 566  1ER17292   26     M   5.3 10.50   12.0  18.30   42.0      22 2009
## 567  1ER17293   26     M   3.8  8.10   12.1  17.50   36.0      22 2009
## 568  1ER17294   26     M   5.2 10.30   12.2  18.60   39.0      22 2009
## 569  1ER17295   26     M   5.1 10.30   11.9  18.20   39.5      22 2009
## 570  1ER17462  729     H   6.2 12.40   10.8  19.10   43.0      24 2009
## 571  1ER17463  729     H   6.2 12.80   11.2  19.60   44.5      24 2009
## 572  1ER17464  729     H   5.0 10.50   11.1  19.00   40.0      24 2009
## 573  1ER17465  729     H   5.6 11.30   11.3  19.40   44.0      24 2009
## 574  1ER17466  729     H   6.2 12.30   11.5  19.30   45.0      24 2009
## 575  1ER17467  729     H   6.0 12.00   11.0  19.10   42.5      24 2009
## 576  1ER17468  729     H   6.5 13.10   10.3  19.40   45.0      24 2009
## 577  1ER17469  729     H   6.5 12.80   10.7  19.00   44.0      24 2009
## 578  1ER17452  629     H   5.8 10.60   12.0  18.50   40.5      23 2009
## 579  1ER17453  629     H   5.7 11.60   11.7  19.00   40.0      23 2009
## 580  1ER17454  629     H   5.2 10.30   11.9  18.70   39.0      23 2009
## 581  1ER17455  629     H   3.8  8.80    9.8  17.70   33.5      23 2009
## 582  1ER17456  629     H   5.1 10.00   11.2  18.40   37.0      23 2009
## 583  1ER17457  629     H   5.8 10.90   12.0  19.30   39.0      23 2009
## 584  1ER17458  629     H   4.8  9.90   11.1  18.30   39.0      23 2009
## 585  1ER17459  629     H   5.8 11.20   11.9  19.40   39.0      23 2009
## 586  1ER17460  629     H   5.1 10.50   11.7  18.80   38.0      23 2009
## 587  1ER17461  629     H   3.7  8.50   10.2  17.40   31.5      23 2009
## 588      <NA>  629     H   1.4    NA     NA     NA     NA      23 2009
## 589  1ER17443   62     M   5.3 10.60     NA     NA     NA      24 2009
## 590  1ER17444   62     M   5.1 10.60     NA     NA     NA      24 2009
## 591  1ER17445   62     M   5.0 10.40     NA     NA     NA      24 2009
## 592  1ER17446   62     M   5.7 11.30     NA     NA     NA      24 2009
## 593  1ER17447   62     M   5.7 11.30     NA     NA     NA      24 2009
## 594  1ER17448   62     M   5.4 11.20     NA     NA     NA      24 2009
## 595  1ER17449   62     M   5.9 10.90     NA     NA     NA      24 2009
## 596  1ER17450   62     M   5.8 11.70     NA     NA     NA      24 2009
## 597  1ER17451   62     M   3.2  8.20     NA     NA     NA      24 2009
## 598  1ER17470 7013     M   5.2 10.10   11.2  17.70   39.0      24 2009
## 599  1ER17471 7013     M   6.8 11.50   12.1  18.60   44.5      24 2009
## 600  1ER17472 7013     M   6.2 10.80   11.7  18.10   44.0      24 2009
## 601  1ER17473 7013     M   5.6 10.40   11.8  18.40   43.0      24 2009
## 602  1ER17474 7013     M   6.1 11.10   12.7  18.20   44.0      24 2009
## 603  1ER17475 7013     M   6.3 11.50   12.5  18.60   45.0      24 2009
## 604  1ER17476 7013     M   5.8 11.20   11.3  18.00   43.0      24 2009
## 605  1ER17478 7013     M   4.7  9.60   12.2  18.30   44.0      24 2009
## 606  1ER17479 7013     M   6.3 11.00   11.7  18.00   41.5      24 2009
## 607  1ER79507 7013     M   3.3  8.00   10.8  17.40   36.0      24 2009
## 608  1ER17481  257     M   6.2 10.90   13.6  18.80   44.5      23 2009
## 609  1ER17482  257     M   6.6 11.80   12.3  18.90   45.0      23 2009
## 610  1ER17483  257     M   5.7 10.80   12.2  18.20   42.5      23 2009
## 611  1ER17484  257     M   4.6  9.60   12.2  18.20   39.0      23 2009
## 612  1ER17485  257     M   5.6 10.50   12.6  18.80   41.5      23 2009
## 613  1ER17486  257     M   5.8 11.20   12.0  18.20   43.0      23 2009
## 614  1ER17487  257     M   4.4  9.20   11.2  17.20   38.0      23 2009
## 615  1ER17488  257     M   6.0 10.80   12.7  18.90   42.5      23 2009
## 616  1ER17489  257     M   4.8  9.90   11.7  18.90   39.5      23 2009
## 617  1ER17490  257     M   4.0  8.50   10.6  17.70   37.5      23 2009
## 618  1ER17491  257     M   5.3 10.50   12.0  18.70   42.0      23 2009
## 619  1ER79528  257     M   2.8    NA   11.2  18.40   34.5      23 2009
## 620  1ER17492  508     L   5.9 10.80   10.9  17.50   43.0      21 2009
## 621  1ER17493  508     L   6.6 11.80   12.2  18.30   44.5      21 2009
## 622  1ER17494  508     L   5.0 10.30   11.3  17.70   40.0      21 2009
## 623  1ER17495  508     L   3.8  8.80   10.6  17.40   37.5      21 2009
## 624  1ER17496  508     L   7.0 12.00   12.3  18.30   44.0      21 2009
## 625  1ER17497  508     L   6.6 12.10   11.8  18.60   45.0      21 2009
## 626  1ER17498  508     L   6.7 12.50   11.7  18.90   45.0      21 2009
## 627  1ER17499  508     L   7.0 12.20   11.9  18.90   44.0      21 2009
## 628  1ER17900  508     L   6.6 12.60   12.0  18.70   45.0      21 2009
## 629  1ER17901  508     L   6.4 11.40   11.7  18.40   45.0      21 2009
## 630  1ER17902  508     L   5.5 10.60   10.7  17.10   41.0      21 2009
## 631  1ER17925  299     H   4.8  9.60     NA     NA     NA      25 2009
## 632  1ER17926  299     H   6.2 11.30     NA     NA     NA      25 2009
## 633  1ER17927  299     H   6.8 12.20     NA     NA     NA      25 2009
## 634  1ER17928  299     H   6.0 11.60     NA     NA     NA      25 2009
## 635  1ER17929  299     H   6.6 11.10     NA     NA     NA      25 2009
## 636  1ER17930  299     H   5.4 10.40     NA     NA     NA      25 2009
## 637  1ER17931  299     H   7.1 12.20     NA     NA     NA      25 2009
## 638  1ER17932  299     H   6.5 11.90     NA     NA     NA      25 2009
## 639  1ER17933  299     H   6.6 10.90     NA     NA     NA      25 2009
## 640  1ER17934  299     H   6.3 11.90     NA     NA     NA      25 2009
## 641  1ER17905  756     L   4.6  9.60   11.7  18.30   38.0      21 2009
## 642  1ER17906  756     L   7.1 11.80   12.3  19.10   43.5      21 2009
## 643  1ER17907  756     L   5.1 10.60   12.0  18.00   41.0      21 2009
## 644  1ER17908  756     L   5.2 10.50   10.6  18.30   41.0      21 2009
## 645  1ER17909  756     L   6.5 11.60   12.2  18.20   44.0      21 2009
## 646  1ER17910  756     L   6.5 11.90   11.4  18.60   42.0      21 2009
## 647  1ER17911  756     L   4.8 10.00   11.8  18.20   40.5      21 2009
## 648  1ER17912  756     L   5.4 10.60   10.7  17.90   43.5      21 2009
## 649  1ER79539  756     L   2.6    NA   11.5  18.00   34.5      21 2009
## 650  1ER17913  569     H   6.6 12.10     NA     NA     NA      24 2009
## 651  1ER17914  569     H   5.8 11.20     NA     NA     NA      24 2009
## 652  1ER17915  569     H   6.0 12.00     NA     NA     NA      24 2009
## 653  1ER17916  569     H   4.6  9.60     NA     NA     NA      24 2009
## 654  1ER17917  569     H   5.6 10.90     NA     NA     NA      24 2009
## 655  1ER17918  569     H   6.6 12.40     NA     NA     NA      24 2009
## 656  1ER17920  569     H   7.1 12.50     NA     NA     NA      24 2009
## 657  1ER17921  569     H   6.3 11.50     NA     NA     NA      24 2009
## 658  1ER17922  569     H   6.3 12.10     NA     NA     NA      24 2009
## 659  1ER17923  569     H   6.4 12.40     NA     NA     NA      24 2009
## 660      <NA>  569     H   2.9    NA     NA     NA     NA      24 2009
## 661      <NA>  569     H   3.1    NA     NA     NA     NA      24 2009
## 662  1ER17958   58     L   5.8 11.50     NA     NA     NA      26 2009
## 663  1ER17959   58     L   5.3 10.80     NA     NA     NA      26 2009
## 664  1ER17961   58     L   5.6 11.20     NA     NA     NA      26 2009
## 665  1ER17962   58     L   6.4 11.40     NA     NA     NA      26 2009
## 666  1ER17963   58     L   5.6 10.50     NA     NA     NA      26 2009
## 667  1ER17964   58     L   5.1  9.40     NA     NA     NA      26 2009
## 668  1ER17965   58     L   6.8 11.80     NA     NA     NA      26 2009
## 669  1ER17966   58     L   5.6 10.90     NA     NA     NA      26 2009
## 670  1ER17967   58     L   4.1  8.80     NA     NA     NA      26 2009
## 671  1ER79540   58     L   2.9    NA     NA     NA     NA      26 2009
## 672  1ER17945  811     M   6.6 12.40   11.5  18.80   45.0      26 2009
## 673  1ER17946  811     M   4.6  9.90   11.3  18.50   40.0      26 2009
## 674  1ER17947  811     M   5.5 11.20   11.9  18.90   46.0      26 2009
## 675  1ER17948  811     M   7.0 13.00   11.8  18.40   45.0      26 2009
## 676  1ER17949  811     M   7.0 12.70   12.2  18.50   44.5      26 2009
## 677  1ER17950  811     M   6.3 11.60   13.4  18.90   45.0      26 2009
## 678  1ER17951  811     M   6.4 11.40   12.5  19.00   44.5      26 2009
## 679  1ER17952  811     M   5.3 10.80   11.1  18.50   42.5      26 2009
## 680  1ER17953  811     M   7.1 12.90   12.6  19.30   46.0      26 2009
## 681  1ER17954  811     M   7.1 12.60   11.8  19.10   44.0      26 2009
## 682  1ER17955  811     M   6.3 11.80   11.9  18.20   44.0      26 2009
## 683  1ER17938  860     L   5.4  9.90   12.3  16.80   43.0      26 2009
## 684  1ER17939  860     L   6.3 11.30   11.5  17.50   44.0      26 2009
## 685  1ER17940  860     L   5.5 10.20   11.9  17.90   43.0      26 2009
## 686  1ER17941  860     L   6.4 11.60   12.1  18.10   45.0      26 2009
## 687  1ER17942  860     L   6.8 10.70   12.8  18.30   45.0      26 2009
## 688  1ER17971  367     L   6.1 11.40   11.4  18.20   44.0      23 2009
## 689  1ER17972  367     L   6.1 11.00   12.1  18.80   44.0      23 2009
## 690  1ER17973  367     L   6.1 11.80   12.0  18.30   44.5      23 2009
## 691  1ER17974  367     L   5.1  9.90   10.8  18.20   41.0      23 2009
## 692  1ER17975  367     L   6.2 11.10   12.5  18.70   45.0      23 2009
## 693  1ER17976  367     L   6.5 12.00   12.2  18.70   45.0      23 2009
## 694  1ER17977  367     L   6.7 12.10   11.6  18.40   44.0      23 2009
## 695  1ER17978  367     L   6.4 12.00   12.6  19.10   45.5      23 2009
## 696  1ER17979  367     L   6.8 12.00   11.2  18.10   45.0      23 2009
## 697  1ER17980  367     L   5.6 11.00   10.9  17.80   38.5      23 2009
## 698  1ER17981  367     L   3.9  9.10   10.9  17.80   38.5      23 2009
## 699   CJ32961  298     L   5.7 11.40   11.3  18.20   43.0      26 2009
## 700   CJ32962  298     L   4.8  9.80   12.2  18.30   42.0      26 2009
## 701   CJ32963  298     L   6.0 11.50   11.6  18.80   42.5      26 2009
## 702   CJ32964  298     L   6.6 11.70   10.8  17.90   41.0      26 2009
## 703   CJ32966  298     L   6.0 11.50   10.4  17.60   44.5      26 2009
## 704   CJ32967  298     L   6.1 12.50   11.3  18.90   42.5      26 2009
## 705   CJ32968  298     L   6.2 11.80   10.8  17.70   43.5      26 2009
## 706   CJ32969  298     L   6.6 12.00   12.0  18.20   45.0      26 2009
## 707   CJ32954  269     M   6.9 11.30     NA     NA     NA      27 2009
## 708   CJ32955  269     M   6.3 11.30     NA     NA     NA      27 2009
## 709   CJ32956  269     M   6.4 11.60     NA     NA     NA      27 2009
## 710   CJ32957  269     M   6.5 11.90     NA     NA     NA      27 2009
## 711   CJ32958  269     M   7.1 11.70     NA     NA     NA      27 2009
## 712   CJ32959  269     M   6.5 11.20     NA     NA     NA      27 2009
## 713   CJ32960  269     M   6.5 11.40     NA     NA     NA      27 2009
## 714  1ER17984   31     M   6.9 12.00   11.1  17.60   44.5      28 2009
## 715  1ER17985   31     M   6.8 12.20   11.3  18.20   42.0      28 2009
## 716  1ER17986   31     M   7.9 13.20   10.9  18.30   44.5      28 2009
## 717  1ER17987   31     M   7.2 13.20   11.2  18.30   45.0      28 2009
## 718  1ER17988   31     M   6.8 12.60   10.1  18.30   42.5      28 2009
## 719  1ER17989   31     M   6.5 11.90   11.0  17.50   42.0      28 2009
## 720  1ER17990   31     M   5.0 10.10    8.3  17.00   35.0      28 2009
## 721  1ER17991   31     M   7.7 13.30   10.6  18.50   44.0      28 2009
## 722  1ER17992   31     M   6.1 12.00    9.9  17.80   42.0      28 2009
## 723  1ER17993   31     M   6.4 12.40   10.1  17.80   43.0      28 2009
## 724  1ER79525  375     L   5.8 10.70     NA     NA     NA      29 2009
## 725  1ER17526  375     L   5.8 10.50     NA     NA     NA      29 2009
## 726   CJ32970  375     L   4.0  9.10     NA     NA     NA      29 2009
## 727   CJ32971  375     L   4.8  9.50     NA     NA     NA      29 2009
## 728  1EP51709   17     H   5.8 11.70   11.1  18.30   41.5      20 2008
## 729  1EP51710   17     H   6.4 12.10   11.3  18.70   42.0      20 2008
## 730  1EP51711   17     H   5.6 10.60   11.4  17.70   38.0      20 2008
## 731  1EP51712   17     H   5.4  9.90   11.7  18.30   40.0      20 2008
## 732  1EP51714   17     H   4.9 10.40   10.3  17.90   37.5      20 2008
## 733  1EP51715   17     H   6.2 11.70   11.2  18.20   41.5      20 2008
## 734  1EP51716   17     H   6.2 11.60   11.4  18.20   41.0      20 2008
## 735  1EP52183   17     H   4.1  9.80    9.8  17.60   35.0      20 2008
## 736  1EP51717   17     H   6.5 11.80   11.6  18.00   40.0      20 2008
## 737  1EP51751   67     H   6.4 11.60   11.8  18.60   43.0      19 2008
## 738  1EP51752   67     H   5.9 11.00   11.4  18.70   43.5      19 2008
## 739  1EP51753   67     H   5.8 10.80   10.9  18.60   41.5      19 2008
## 740  1EP51754   67     H   5.7 11.00   11.2  18.80   42.0      19 2008
## 741  1EP51755   67     H   6.0 11.70   10.4  18.10   43.5      19 2008
## 742  1EP51756   67     H   5.8 11.70   10.6  18.10   41.0      19 2008
## 743  1EP51757   67     H   6.6 12.30   11.3  18.30   44.5      19 2008
## 744  1EP51758   67     H   6.2 12.10   11.7  19.40   45.0      19 2008
## 745  1EP51759   67     H   5.4 11.00   10.5  19.00   41.0      19 2008
## 746  1EP51760   67     H   6.1 11.60   11.2  18.70   42.5      19 2008
## 747  1EP51761  189     H   6.8 12.30   12.1  18.70   44.5      22 2008
## 748  1EP51762  189     H   4.1  8.60   10.6  18.00   36.0      22 2008
## 749  1EP51763  189     H   6.3 11.40   12.2  18.80   41.0      22 2008
## 750  1EP51764  189     H   6.8 12.10   11.9  19.20   45.0      22 2008
## 751  1EP51765  189     H   6.1 12.20   11.4  18.80   44.0      22 2008
## 752  1EP51766  189     H   6.6 12.10   11.8  18.60   43.0      22 2008
## 753  1EP51767  189     H   4.8  9.40   12.2  18.70   39.0      22 2008
## 754  1EP51768  189     H   5.8 10.60   11.6  18.70   42.5      22 2008
## 755  1EP51769  823     H   5.9 11.40   11.9  18.90   42.0      23 2008
## 756  1EP51770  823     H   6.5 11.80   11.1  18.20   44.0      23 2008
## 757  1EP51771  823     H   4.7  9.60   11.6  18.10   40.0      23 2008
## 758  1EP51772  823     H   6.3 12.20   11.7  18.50   44.0      23 2008
## 759  1EP51773  823     H   6.3 12.30   11.3  18.40   43.0      23 2008
## 760  1EP51774  823     H   5.8 11.80   12.1  19.40   43.0      23 2008
## 761  1EP51775  823     H   6.8 12.00   11.7  18.40   45.0      23 2008
## 762  1EP51776  823     H   6.9 12.10   12.3  18.80   46.0      23 2008
## 763  1EP51777  823     H   6.8 11.30   11.9  18.80   43.5      23 2008
## 764  1EP51778  613     H   6.3 11.20   12.2  18.80   43.0      23 2008
## 765  1EP51779  613     H   6.8 11.80   11.5  17.80   42.5      23 2008
## 766  1EP51780  613     H   7.7 13.30   12.1  19.00   46.0      23 2008
## 767  1EP51781  613     H   6.6 12.10   11.3  18.50   44.5      23 2008
## 768  1EP51782  613     H   6.7 11.70   12.0  18.10   44.0      23 2008
## 769  1EP51783  613     H   6.7 12.10   11.3  18.30   45.0      23 2008
## 770  1EP51784  613     H   7.7 13.30   11.6  19.10   47.0      23 2008
## 771  1EP51785  613     H   4.8  9.50   11.3  18.20   40.5      23 2008
## 772  1EP51786  613     H   7.6 12.60   11.8  18.40   45.5      23 2008
## 773  1EP51787  613     H   5.9 10.60   12.8  18.90   43.0      23 2008
## 774  1EP51788  687     H   6.9 12.00   11.2  18.50   44.0      23 2008
## 775  1EP51789  687     H   7.2 12.40   12.1  18.60   45.5      23 2008
## 776  1EP51790  687     H   6.7 12.70   11.2  18.30   44.0      23 2008
## 777  1EP51791  687     H   6.7 12.40   11.2  18.10   45.0      23 2008
## 778  1EP51792  687     H   5.6 11.70   12.1  18.90   42.5      23 2008
## 779  1EP51793  687     H   6.1 11.60   12.2  18.30   43.0      23 2008
## 780  1EP51794  687     H   6.6 12.30   11.8  17.90   44.0      23 2008
## 781  1EP51795  687     H   5.6 10.20   11.3  18.40   41.0      23 2008
## 782  1EP51798  204     H   9.1 13.70   11.7  19.10   46.0      22 2008
## 783  1EP51799  204     H   9.6 14.40   13.1  19.40   47.0      22 2008
## 784  1EP52300  204     H   9.6 10.20   12.6  19.30   40.5      22 2008
## 785  1EP52301  204     H   9.6 14.50   11.6  19.20   47.5      22 2008
## 786  1EP52302  204     H   8.9 13.50   12.1  18.90   48.0      22 2008
## 787  1EP52303  204     H   6.6 10.50   11.3  18.00   39.0      22 2008
## 788  1EP52304  204     H   8.4 12.80   11.5  18.60   46.0      22 2008
## 789  1EP52336  207     H   7.0 11.20   12.2  18.40   44.0      23 2008
## 790  1EP52337  207     H   5.5 11.80   11.2  18.90   44.5      23 2008
## 791  1EP52338  207     H   6.1 11.40   11.6  18.30   44.0      23 2008
## 792  1EP52339  207     H   6.7 11.40   12.0  18.70   44.5      23 2008
## 793  1EP52340  207     H   4.2  8.50   11.6  18.00   39.5      23 2008
## 794  1EP52341  207     H   5.0 10.00   12.8  19.40   41.5      23 2008
## 795  1EP52342  207     H   5.8 10.40   10.4  17.70   42.5      23 2008
## 796  1EP52343  207     H   6.0 10.60   11.1  17.50   42.0      23 2008
## 797  1EP52344  207     H   5.4  9.70   11.5  17.90   42.5      23 2008
## 798  1EP52373  566     H   6.2 11.50     NA     NA     NA      25 2008
## 799  1EP52374  566     H   6.6 11.60     NA     NA     NA      25 2008
## 800  1EP52375  566     H   5.8 10.60     NA     NA     NA      25 2008
## 801  1EP52376  566     H   6.2 10.40     NA     NA     NA      25 2008
## 802  1EP52378  566     H   6.1 11.10     NA     NA     NA      25 2008
## 803  1EP52379  566     H   6.3 10.80     NA     NA     NA      25 2008
## 804  1EP52380  566     H   7.1 11.80     NA     NA     NA      25 2008
## 805  1EP52385  632     H   5.1  9.00     NA     NA     NA      26 2008
## 806  1EP52386  632     H   4.9  9.80     NA     NA     NA      26 2008
## 807  1EP52387  632     H   6.1 10.90     NA     NA     NA      26 2008
## 808  1EP52388  632     H   4.7  9.60     NA     NA     NA      26 2008
## 809  1EP52389  632     H   6.0 11.00     NA     NA     NA      26 2008
## 810  1EP52391  632     H   5.4  9.80     NA     NA     NA      26 2008
## 811  1EP52392  632     H   6.0 10.50     NA     NA     NA      26 2008
## 812  1EP52403  808     H   6.3 11.30   11.7  18.30   43.5      24 2008
## 813  1EP52404  808     H   7.1 12.20   11.6  18.80   44.5      24 2008
## 814  1EP52405  808     H   6.7 11.50   11.8  17.90   45.0      24 2008
## 815  1EP52406  808     H   5.4 10.20   11.3  18.00   43.0      24 2008
## 816  1EP52407  808     H   6.9 11.20   12.5  18.70   45.0      24 2008
## 817  1EP52408  808     H   4.5  8.60   11.7  17.40   39.5      24 2008
## 818  1EP52409  699     H   6.2 10.80   12.1  18.30   44.0      26 2008
## 819  1EP52410  699     H   6.3 11.30   12.2  18.60   43.5      26 2008
## 820  1EP52411  699     H   4.6  8.20   12.1  17.30   39.0      26 2008
## 821  1EP52412  699     H   7.1 11.10   12.4  18.30   45.0      26 2008
## 822  1EP52413  699     H   5.3  9.80   12.2  18.40   42.5      26 2008
## 823  1EP52414  699     H   4.6  8.30   11.8  18.00   39.0      26 2008
## 824  1EP52415  699     H   6.4 11.40   11.6  18.10   45.0      26 2008
## 825  1EP52416  699     H   6.6 11.20   11.4  18.10   45.0      26 2008
## 826  1EP52417  699     H   6.1 10.10   11.1  17.90   43.0      26 2008
## 827  1EP52418  869     H   7.1 11.50   11.6  19.20   43.5      25 2008
## 828  1EP52419  869     H   5.8  9.00   11.5  18.40   40.0      25 2008
## 829  1EP52420  869     H   7.5 11.10   11.9  18.50   44.0      25 2008
## 830  1EP52421  869     H   7.1 11.40   10.4  18.60   42.0      25 2008
## 831  1EP52437  185     H   6.0  9.90   12.1  17.80   44.0      25 2008
## 832  1EP52438  185     H   7.3 12.50   11.8  19.00   46.0      25 2008
## 833  1EP52439  185     H   7.5 12.40   12.1  18.60   44.5      25 2008
## 834  1EP52440  185     H   7.5 11.90   11.8  18.70   46.0      25 2008
## 835  1EP52441  185     H   5.4  9.90   11.7  18.60   45.0      25 2008
## 836  1EP52442  185     H   5.7 10.80   10.6  18.30   42.5      25 2008
## 837  1EP52443  185     H   6.6 12.40   12.2  18.70   44.5      25 2008
## 838  1EP52444  185     H   6.0 10.70   11.6  19.40   43.0      25 2008
## 839  1EP52445  185     H   6.8 11.40   11.7  19.00   46.0      25 2008
## 840  1EP52446  185     H   7.0 12.00   11.9  18.50   44.0      25 2008
## 841  1EP52460  V53     H   6.2 10.70   13.2  19.00   44.0      27 2008
## 842  1EP52461  V53     H   5.7 10.40   13.2  19.50   44.0      27 2008
## 843  1EP52462  V53     H   5.8 10.40   13.7  19.30   43.0      27 2008
## 844  1EP52463  V53     H   5.3  9.80   12.4  17.90   44.5      27 2008
## 845  1EP52464  V53     H   6.0 10.70   11.6  18.30   44.0      27 2008
## 846  1EP52465  V53     H   6.1 10.40   13.1  19.50   43.0      27 2008
## 847  1EP52466  V53     H   5.7 10.50   12.7  18.70   43.5      27 2008
## 848  1EP52467  V53     H   5.5 10.00   12.8  18.80   42.5      27 2008
## 849  1EP51827  V53     H   3.1    NA   11.5  18.10   36.0      27 2008
## 850  1EP52119  279     H   5.6 10.90   11.2  18.40   41.0      26 2008
## 851  1EP52120  279     H   5.8 10.40   10.8  17.90   41.0      26 2008
## 852  1EP52121  279     H   5.2  9.50   10.4  17.70   38.5      26 2008
## 853  1EP52122  279     H   6.2 10.80   11.8  19.10   43.0      26 2008
## 854  1EP52123  279     H   6.3 10.70   11.7  18.60   43.5      26 2008
## 855  1EP52124  279     H   5.9 10.40   10.8  17.80   42.5      26 2008
## 856  1EP52125  279     H   5.9  8.90   10.9  17.50   41.0      26 2008
## 857  1EP52126  279     H   4.9  9.40   10.0  17.80   39.0      26 2008
## 858  1EP51828  279     H   4.6  8.40   11.5  18.00   39.0      26 2008
## 859  1EP52127   46     H   8.2 12.90   13.2  19.50   48.0      29 2008
## 860  1EP52128   46     H   7.6 12.20   12.3  18.40   47.0      29 2008
## 861  1EP52129   46     H   5.2 10.10   11.6  18.60   43.5      29 2008
## 862  1EP52130   46     H   6.3 10.70   12.3  18.30   44.5      29 2008
## 863  1EP52131   46     H   6.8 11.90   11.8  19.00   47.5      29 2008
## 864  1EP52133   46     H   4.3  9.00   11.1  18.00   39.5      29 2008
## 865  1EP52134   46     H   7.3 12.70   11.8  19.40   46.0      29 2008
## 866  1EP52135   46     H   7.8 12.80   12.7  18.90   48.0      29 2008
## 867  1EP52136   46     H   7.6 13.30   11.9  19.70   48.0      29 2008
## 868  1EP52137   46     H   7.6 13.30   12.3  19.60   47.0      29 2008
## 869  1EP52145  826     H   8.0 13.60   12.2  19.80   48.5      27 2008
## 870  1EP52146  826     H   7.5 13.20   13.3  19.20   49.0      27 2008
## 871  1EP52147  826     H   5.8 11.10   12.1  18.60   43.5      27 2008
## 872  1EP52148  826     H   7.2 13.00   12.8  18.60   48.0      27 2008
## 873  1EP52149  826     H   6.8 13.20   11.3  19.20   48.0      27 2008
## 874  1EP52150  826     H   6.7 12.10   12.6  18.90   46.0      27 2008
## 875  1EP52151  826     H   7.5 13.30   12.3  19.00   48.0      27 2008
## 876  1EP52152  826     H   7.7 13.70   13.2  19.70   49.5      27 2008
## 877      <NA>  826     H   1.6    NA     NA     NA     NA      27 2008
## 878      <NA>  826     H   2.9    NA     NA     NA     NA      27 2008
## 879  1EP52163  574     H   6.1 10.20   12.7  19.20   43.0      29 2008
## 880  1EP52164  574     H   6.2 10.00   12.2  18.10   41.5      29 2008
## 881  1EP52165  574     H   6.4 10.50   11.8  18.70   43.5      29 2008
## 882  1EP52166  574     H   6.6 10.90   11.1  18.60   44.0      29 2008
## 883  1EP52167  574     H   6.2 10.30   11.6  18.00   42.5      29 2008
## 884  1EP52168  574     H   5.8 10.40   12.0  18.50   42.5      29 2008
## 885  1EP52169  574     H   5.7  9.90   11.3  17.50   43.0      29 2008
## 886  1EP52170  574     H   7.1 10.50   13.0  18.70   44.5      29 2008
## 887  1EP52171  574     H   4.3  7.10   12.1  18.00   38.0      29 2008
## 888  1EP51726   28     L   5.6 10.30   12.7  18.20   41.0      19 2008
## 889  1EP51727   28     L   5.3 10.40   11.8  18.10   40.0      19 2008
## 890  1EP51728   28     L   5.6 10.20   13.2  18.50   43.5      19 2008
## 891  1EP52196   28     L   3.6  8.20   12.4  17.60   34.0      19 2008
## 892  1EP51742  372     L   5.8 11.00     NA     NA     NA      20 2008
## 893  1EP51743  372     L   6.1 11.70     NA     NA     NA      20 2008
## 894  1EP51744  372     L   5.1 10.40     NA     NA     NA      20 2008
## 895  1EP51745  372     L   5.3 10.00     NA     NA     NA      20 2008
## 896  1EP51746  372     L   6.0 10.40     NA     NA     NA      20 2008
## 897  1EP51747  372     L   5.9 10.70     NA     NA     NA      20 2008
## 898  1EP51748  372     L   4.4  8.80     NA     NA     NA      20 2008
## 899  1EP51749  372     L   5.1  9.50     NA     NA     NA      20 2008
## 900  1EP51750  372     L   5.1 10.00     NA     NA     NA      20 2008
## 901  1EP51796  693     L   5.7 10.50   10.2  16.60   39.5      20 2008
## 902  1EP51797  693     L   5.6 10.20   10.3  17.30   40.0      20 2008
## 903  1EP52331  211     L   5.3  9.60   11.5  18.00   43.5      20 2008
## 904  1EP52332  211     L   5.4    NA   12.2  18.00   42.5      20 2008
## 905  1EP52333  211     L   5.3 10.20   11.9  18.00   43.5      20 2008
## 906  1EP52334  211     L   4.7  9.40   10.8  17.40   41.5      20 2008
## 907  1EP52335  211     L   5.6 10.00   11.6  17.60   42.0      20 2008
## 908  1EP52369  773     L   5.0 10.00   11.7  17.30   43.0      20 2008
## 909  1EP52370  773     L   4.6  8.60   10.9  16.50   38.0      20 2008
## 910  1EP52371  773     L   5.3  9.20   12.0  17.20   42.0      20 2008
## 911  1EP52372  773     L   4.7  8.70   11.1  17.30   40.5      20 2008
## 912  1EP51818  773     L   4.2  7.30   11.4  15.70   38.0      20 2008
## 913  1EP52393  608     L   7.3 12.30   12.2  18.40   47.0      22 2008
## 914  1EP52394  608     L   7.2 11.80   12.3  18.00   47.5      22 2008
## 915  1EP52395  608     L   7.1 12.30   12.5  18.30   46.0      22 2008
## 916  1EP52396  608     L   5.6 11.30   11.1  18.40   44.0      22 2008
## 917  1EP52397  608     L   6.3 11.80   11.7  18.40   45.5      22 2008
## 918  1EP52398  608     L   7.0 11.90   11.8  18.00   46.5      22 2008
## 919  1EP52399  608     L   6.5 10.80   11.7  17.70   44.0      22 2008
## 920  1EP52400  608     L   7.0 11.90   12.2  18.80   47.5      22 2008
## 921  1EP52401  608     L   6.9 11.70   11.7  18.10   46.5      22 2008
## 922  1EP52402  608     L   6.6 12.20   11.3  17.70   45.5      22 2008
## 923  1EP51820  608     L   3.3    NA   12.3  17.60   36.5      22 2008
## 924  1EP52447  736     L   6.9 12.20   11.9  18.50   45.0      24 2008
## 925  1EP52448  736     L   6.6 10.80   12.4  18.30   45.5      24 2008
## 926  1EP52449  736     L   7.2 12.60   11.8  18.90   47.5      24 2008
## 927  1EP52450  736     L   6.3 11.20   11.8  18.80   42.0      24 2008
## 928  1EP52451  736     L   6.9 11.80   11.8  18.30   44.0      24 2008
## 929  1EP52452  736     L   7.3 12.10   12.1  18.70   47.0      24 2008
## 930  1EP52453  736     L   7.3 11.30   11.2  17.40   44.0      24 2008
## 931  1EP52454  736     L   7.6 11.40   12.5  18.40   46.0      24 2008
## 932  1EP52455  845     L   4.5  9.00    9.7  16.90   39.0      25 2008
## 933  1EP52456  845     L   6.2 10.70   10.2  17.50   43.5      25 2008
## 934  1EP52457  845     L   6.5    NA    9.3  17.10   44.0      25 2008
## 935  1EP52458  845     L    NA    NA   11.3  18.10   45.0      25 2008
## 936  1EP52476  727     L   7.1 11.20   12.8  19.10   44.5      25 2008
## 937  1EP52477  727     L   6.1 10.90   12.1  18.60   43.5      25 2008
## 938  1EP52478  727     L   6.7 10.80   12.6  18.80   44.0      25 2008
## 939  1EP52479  727     L   5.8  9.50   12.6  18.80   41.5      25 2008
## 940  1EP52480  727     L   6.2 10.70   12.2  17.80   43.0      25 2008
## 941  1EP52486  835     L   7.2 11.50   12.3  18.00   45.0      26 2008
## 942  1EP52487  835     L   7.6 12.10   12.1  17.50   46.5      26 2008
## 943  1EP52488  835     L   7.1 11.50   11.4  17.80   43.5      26 2008
## 944  1EP52489  835     L   7.9 12.70   13.3  18.70   47.0      26 2008
## 945  1EP52490  835     L   7.9 12.60   11.2  17.90   45.5      26 2008
## 946  1EP52491  835     L   7.2 12.50   11.6  17.30   46.0      26 2008
## 947  1EP52492  835     L   5.6 10.40   10.8  17.90   42.0      26 2008
## 948  1EP52493  835     L   7.2 12.30   11.4  17.60   46.0      26 2008
## 949  1EP52115  603     L   6.1  9.90   11.3  17.40   43.5      26 2008
## 950  1EP52116  603     L   7.1 10.80   11.8  17.80   45.0      26 2008
## 951  1EP52117  603     L   7.1 11.30   11.9  18.30   45.5      26 2008
## 952  1EP52118  603     L   5.2    NA   12.7  17.50   41.0      26 2008
## 953   CG53578  678     L   4.4  8.70   11.1  17.50   40.5      28 2008
## 954   CG53579  678     L   4.2  8.10   10.4  16.80   37.0      28 2008
## 955   CG53580  678     L   4.7  8.80   10.9  17.70   42.5      28 2008
## 956  1EP51832  678     L   2.7    NA    9.5  17.00   33.5      28 2008
## 957  1EP52153  158     L   6.7 11.00   12.3  18.40   45.0      23 2008
## 958  1EP52154  158     L   5.5  8.70   12.2  18.00   42.5      23 2008
## 959  1EP52155  158     L   4.6  8.80   11.1  17.00   38.0      23 2008
## 960  1EP52156  158     L   5.7 10.80   11.6  18.30   43.0      23 2008
## 961  1EP52157  158     L   6.2 11.00   11.9  17.70   43.5      23 2008
## 962  1EP52158  158     L   5.0  9.70   11.6  18.70   39.5      23 2008
## 963  1EP52159  158     L   6.0  9.40   11.7  17.50   43.0      23 2008
## 964  1EP52160  158     L   6.1 10.00   11.8  17.40   44.0      23 2008
## 965  1EP52161  158     L   6.6 11.10   12.6  17.80   43.0      23 2008
## 966  1EP52172  606     L   6.1 10.70   11.2  17.80   45.0      28 2008
## 967  1EP52173  606     L   6.2 11.00   11.6  18.40   45.5      28 2008
## 968  1EP52174  606     L   4.3  8.10   11.2  17.40   40.0      28 2008
## 969  1EP51834  606     L   3.8  8.10   11.6  18.20   38.0      28 2008
## 970  1EP52177  606     L   7.1 11.60   12.1  18.70   48.0      28 2008
## 971  1EP52178  606     L   6.8 11.70   11.4  18.30   47.0      28 2008
## 972  1EP52179  606     L   6.7 11.20   11.6  18.20   46.5      28 2008
## 973  1EP52180  606     L   7.2 11.60   11.6  18.70   47.0      28 2008
## 974  1EP52198  159     L   5.7 10.00   11.9  17.50   44.0      28 2008
## 975  1EP52199  159     L   4.1  8.80   12.1  18.20   41.0      28 2008
## 976  1EP51800  159     L   6.0 11.10   12.6  18.60   45.0      28 2008
## 977  1EP51801  159     L   6.7 12.10   13.2  18.90   46.0      28 2008
## 978  1EP51845  159     L   3.8  8.20   10.2  17.00   38.5      28 2008
## 979  1EP51718  839     M   5.7 10.60    9.1  17.50   38.0      19 2008
## 980  1EP51719  839     M   5.9 11.10    9.8  17.80   38.5      19 2008
## 981  1EP51720  839     M   5.5 10.80   10.1  18.30   37.0      19 2008
## 982  1EP51721  839     M   5.1 10.10   10.3  18.00   38.0      19 2008
## 983  1EP51722  839     M   6.4 11.50   10.2  18.60   39.0      19 2008
## 984  1EP51723  839     M   5.3 10.50    9.8  17.60   38.5      19 2008
## 985  1EP51724  839     M   5.7 11.10   10.7  18.20   37.0      19 2008
## 986  1EP51725  839     M   4.8  9.90    8.6  16.80   34.5      19 2008
## 987  1EP52182  839     M   2.8    NA    6.3  15.70   27.5      19 2008
## 988  1EP52181  839     M   3.2    NA    7.4  16.40   32.0      19 2008
## 989      <NA>  839     M   2.2    NA     NA     NA     NA      19 2008
## 990  1EP51729  723     M   6.7 13.00   11.3  19.30   46.0      22 2008
## 991  1EP51730  723     M   6.7 12.60   11.4  19.90   44.5      22 2008
## 992  1EP51731  723     M   6.3 12.40   12.7  19.70   44.0      22 2008
## 993  1EP51732  723     M   5.1 10.10   12.6  19.60   42.0      22 2008
## 994  1EP51733  723     M   7.2 13.20   12.0  20.00   45.5      22 2008
## 995  1EP51734  723     M   7.2 13.00   12.4  19.90   45.0      22 2008
## 996  1EP51735  817     M   5.1 10.90   11.8  18.80   40.5      23 2008
## 997  1EP51736  817     M   5.2 10.90   12.4  18.60   42.0      23 2008
## 998  1EP51737  817     M   6.2 12.60   12.0  18.70   45.5      23 2008
## 999  1EP51738  817     M   6.1 12.30   11.6  18.30   43.0      23 2008
## 1000 1EP51739  817     M   6.7 13.00   12.5  19.40   44.0      23 2008
## 1001 1EP51740  817     M   5.8 11.60   12.3  18.80   44.0      23 2008
## 1002 1EP51741  817     M   6.5 12.30   12.9  19.10   44.0      23 2008
## 1003 1EP52305  642     M   7.4 12.30   11.6  19.20   47.5      21 2008
## 1004 1EP52306  642     M   5.9 11.00   10.7  17.60   43.5      21 2008
## 1005 1EP52307  642     M   6.6 11.80   10.7  17.80   46.0      21 2008
## 1006 1EP52308  642     M   6.6 11.50   11.2  18.40   46.0      21 2008
## 1007 1EP52309  642     M   7.6 12.10   12.2  19.00   48.5      21 2008
## 1008 1EP52310  642     M   7.1 11.60   11.9  18.60   46.5      21 2008
## 1009 1EP52311  642     M   4.8  8.80   10.9  17.70   40.5      21 2008
## 1010 1EP52312  642     M   5.5  9.30   12.0  18.20   42.0      21 2008
## 1011 1EP51802  642     M   3.6  8.10   11.2  18.10   38.0      21 2008
## 1012 1EP52313  637     M   6.3 11.40   12.3  18.80   46.0      23 2008
## 1013 1EP52314  637     M   5.7 10.40   11.1  17.60   39.5      23 2008
## 1014 1EP52315  637     M   5.8 11.00   12.3  18.90   43.0      23 2008
## 1015 1EP52316  637     M   4.4  8.80   12.3  18.70   40.0      23 2008
## 1016 1EP52317  637     M   5.7 11.10   11.8  18.60   44.5      23 2008
## 1017 1EP52318  637     M   6.0 10.40   12.4  18.10   43.5      23 2008
## 1018 1EP52319  637     M   5.5  9.50   10.6  17.60   42.5      23 2008
## 1019 1EP52320  637     M   5.3  9.70   11.3  18.00   41.5      23 2008
## 1020 1EP52321  637     M   5.3  9.30   11.2  17.10   43.0      23 2008
## 1021 1EP51817  637     M   4.2  7.70   11.8  17.00   38.5      23 2008
## 1022 1EP52322  689     M   5.8 10.70   10.0  18.00   42.5      25 2008
## 1023 1EP52323  689     M   6.2 11.10    9.7  18.30   43.0      25 2008
## 1024 1EP52324  689     M   4.3  8.10    9.3  17.50   38.0      25 2008
## 1025 1EP52325  689     M   5.1  9.50    9.6  17.90   40.0      25 2008
## 1026 1EP52326  689     M   5.6 10.20    9.1  17.40   42.0      25 2008
## 1027 1EP52327  689     M   5.6  9.70    9.5  17.50   40.5      25 2008
## 1028 1EP52328  689     M   6.1 10.30   10.1  18.30   43.0      25 2008
## 1029 1EP52329  689     M   5.8 10.30   11.0  18.60   44.0      25 2008
## 1030 1EP52330  689     M   6.4 10.50   10.4  18.60   44.0      25 2008
## 1031 1EP52345  763     M   6.2 11.40   11.5  18.10   44.0      23 2008
## 1032 1EP52346  763     M   6.6 11.60   12.2  19.10   45.0      23 2008
## 1033 1EP52347  763     M   5.2  9.90   12.6  18.30   42.0      23 2008
## 1034 1EP52348  763     M   6.2 11.30   12.3  19.00   45.0      23 2008
## 1035 1EP52349  763     M   6.7 11.80   12.2  18.90   46.0      23 2008
## 1036 1EP52350  763     M   6.7 11.70   11.9  18.80   45.0      23 2008
## 1037 1EP52351  763     M   4.6  8.90   12.0  18.80   45.0      23 2008
## 1038 1EP52352  651     M   4.3  8.60   11.9  18.20   38.5      24 2008
## 1039 1EP52353  651     M   5.7 10.10   12.3  18.20   43.5      24 2008
## 1040 1EP52354  651     M   5.8 10.10   12.1  18.60   43.5      24 2008
## 1041 1EP52355  651     M   5.8 10.70   11.6  17.80   43.0      24 2008
## 1042 1EP52356  651     M   7.0 11.60   12.4  19.10   45.0      24 2008
## 1043 1EP52357  651     M   6.8 11.70   11.6  18.20   43.0      24 2008
## 1044 1EP52358  651     M   6.3 10.50   12.0  18.40   43.0      24 2008
## 1045 1EP52359  651     M   6.2 10.60   11.3  17.60   44.0      24 2008
## 1046 1EP52360  651     M   6.3 11.30   12.3  19.00   45.5      24 2008
## 1047 1EP52361  801     M   5.6  9.70     NA     NA     NA      25 2008
## 1048 1EP52362  801     M   5.8 10.20     NA     NA     NA      25 2008
## 1049 1EP52363  801     M   5.5  9.60     NA     NA     NA      25 2008
## 1050 1EP52364  801     M   5.1  9.00     NA     NA     NA      25 2008
## 1051 1EP52365  801     M   5.3  9.30     NA     NA     NA      25 2008
## 1052 1EP52366    6     M   5.9 10.40     NA     NA     NA      23 2008
## 1053 1EP52367    6     M   5.9 10.60     NA     NA     NA      23 2008
## 1054 1EP52368    6     M   4.5  8.90     NA     NA     NA      23 2008
## 1055 1EP52381  692     M   7.4 12.10   11.9  18.70   44.5      26 2008
## 1056 1EP52382  692     M   6.9 11.30   12.1  18.20   43.0      26 2008
## 1057 1EP52383  692     M   4.4  7.70   12.2  18.00   38.0      26 2008
## 1058 1EP52384  692     M   6.7 11.00   13.0  18.80   44.0      26 2008
## 1059 1EP52819  692     M   2.6    NA   11.7  17.40   36.0      26 2008
## 1060 1EP52468 V211     M   6.2 10.80   13.0  19.10   43.0      27 2008
## 1061 1EP52469 V211     M   6.6 11.10   12.2  18.30   43.0      27 2008
## 1062 1EP52470 V211     M   7.1 11.40   13.3  19.10   46.0      27 2008
## 1063 1EP52471 V211     M   6.1 10.70   12.6  18.30   45.0      27 2008
## 1064 1EP52472 V211     M   6.7 11.40   12.7  18.50   44.0      27 2008
## 1065 1EP52473 V211     M   5.6  9.70   13.4  19.10   42.0      27 2008
## 1066 1EP52474 V211     M   7.0 11.40   13.5  19.10   44.5      27 2008
## 1067 1EP52475 V211     M   4.7  8.70   12.6  18.20   39.5      27 2008
## 1068 1EP51826 V211     M   3.7    NA   12.8  18.40   44.0      27 2008
## 1069 1EP52481  824     M   7.7 12.70   13.5  19.30   47.0      24 2008
## 1070 1EP52482  824     M   7.7 12.70   12.5  19.40   46.0      24 2008
## 1071 1EP52483  824     M   8.4 12.40   12.6  18.50   44.0      24 2008
## 1072 1EP52484  824     M   7.3 12.30   12.9  19.30   45.5      24 2008
## 1073 1EP52485  824     M   7.4 12.00   12.7  19.20   45.0      24 2008
## 1074 1EP52494  627     M   6.0 10.50   11.2  18.60   44.0      24 2008
## 1075 1EP52495  627     M   4.6  8.90   11.1  19.10   44.0      24 2008
## 1076 1EP52496  627     M   6.2 11.70   11.0  19.60   44.0      24 2008
## 1077 1EP52497  627     M   6.9 11.20   11.2  18.98   45.0      24 2008
## 1078 1EP52498  627     M   7.7 12.50   11.8  19.70   47.0      24 2008
## 1079 1EP52499  627     M   6.3 11.30   11.9  19.60   45.5      24 2008
## 1080 1EP52100  627     M   6.3 10.60   11.8  18.30   44.0      24 2008
## 1081 1EP52101  627     M   6.3 10.60     NA     NA     NA      24 2008
## 1082 1EP52102  627     M   6.5 11.00   12.2  19.80   44.0      24 2008
## 1083 1EP52103  627     M   5.7 10.20   12.2  19.40   44.5      24 2008
## 1084 1EP52104   11     M   6.6 11.10   12.6  19.60   44.5      25 2008
## 1085 1EP52105   11     M   6.8 11.70   12.6  19.50   44.5      25 2008
## 1086 1EP52106   11     M   6.7 11.40   11.7  19.40   43.5      25 2008
## 1087 1EP52107   11     M   6.5 11.40   12.1  19.50   44.0      25 2008
## 1088 1EP52108   11     M   6.6 10.90   12.1  19.60   44.0      25 2008
## 1089 1EP52109   11     M   6.3 10.70   12.6  19.20   44.0      25 2008
## 1090 1EP52110   11     M   6.3 11.20   10.3  19.10   42.5      25 2008
## 1091 1EP52111   11     M   7.1 12.20   12.7  19.70   45.5      25 2008
## 1092 1EP52112   11     M   7.2 11.80   12.0  19.40   45.5      25 2008
## 1093 1EP52113   11     M   5.7 10.50   11.4  19.00   41.0      25 2008
## 1094 1EP52114   11     M   6.1 10.60   12.2  19.30   42.5      25 2008
## 1095 1EP52138  675     M   7.0 10.50     NA     NA     NA      29 2008
## 1096 1EP52139  675     M   5.9 10.40     NA     NA     NA      29 2008
## 1097 1EP52140  675     M   7.1 11.50     NA     NA     NA      29 2008
## 1098 1EP52141  675     M   7.3 12.00     NA     NA     NA      29 2008
## 1099 1EP52142  675     M   6.4 11.00     NA     NA     NA      29 2008
## 1100 1EP52143  675     M   5.7 10.50     NA     NA     NA      29 2008
## 1101 1EP52144  675     M   5.2  9.60     NA     NA     NA      29 2008
## 1102     <NA>  675     M   3.8  8.40     NA     NA     NA      29 2008
## 1103 1EP52185  V32     M   5.9 11.00     NA     NA     NA      27 2008
## 1104 1EP52186  V32     M   7.1 12.30     NA     NA     NA      27 2008
## 1105 1EP52187  V32     M   7.4 12.50     NA     NA     NA      27 2008
## 1106 1EP52188  V32     M   6.7 11.80     NA     NA     NA      27 2008
## 1107 1EP52189  V32     M   6.8 12.20     NA     NA     NA      27 2008
## 1108 1EP52190  V32     M   7.2 12.50     NA     NA     NA      27 2008
## 1109 1EP52191  V32     M   8.1 13.20     NA     NA     NA      27 2008
## 1110 1EP52192  V32     M   6.7 12.30     NA     NA     NA      27 2008
## 1111 1EP52193  V32     M   4.2  8.90     NA     NA     NA      27 2008
## 1112 1EP52194  V32     M   5.8 11.30     NA     NA     NA      27 2008
## 1113 1EP52195  V32     M   7.7 12.80     NA     NA     NA      27 2008

Sample sizes:

There are 1113 nestlings from 127 distinct nests; of these, 291 were incubated at the L temperature, 434 at the M temperature, and 388 at the H temperature.

with(nestlings, table(nest, treat))
##       treat
## nest    L  M  H
##   11    0 11  0
##   112   0  0 11
##   119   0  6  0
##   158   9  0  0
##   159   5  0  0
##   161   0  0  6
##   17    0  0  9
##   176   0  0  8
##   185   0  0 10
##   189   0  0  8
##   191   0  0  7
##   203   0 10  0
##   204   0  0  7
##   207   0  9  9
##   21    0  9  0
##   211   5 10  0
##   213   0  0  9
##   215   0 10  0
##   223   0  0  9
##   247   5  0  7
##   257   0 12  0
##   26    0  9  0
##   269   0  7  0
##   274   0  0  9
##   277   8  0  0
##   279   0  0  9
##   28    4  0  0
##   284   0  8  0
##   29    6  0  0
##   298   8  0  0
##   299   0  0 10
##   31    0 10  0
##   364   0  0 10
##   365   0  0  8
##   367  11  0  0
##   372   9  0  0
##   375   4  0  0
##   4     0  8  0
##   46    0  0 10
##   506   6  0  0
##   508  11  0  0
##   510   4  0  0
##   515   0  5  0
##   524   3  0  0
##   533  11  0  0
##   545   0  0  7
##   548   0 10  0
##   566   0  0  7
##   569   0  0 12
##   571   0 10  0
##   574   0  0  9
##   58   10  0  0
##   6     0  3  0
##   60    5  0  0
##   603   4  0  0
##   606  19  0  0
##   608  11  7  0
##   61    0  7  6
##   613   0  0 10
##   617   0  0  3
##   62    0  9  0
##   623   0  9  0
##   625   0 11  0
##   627   0 10  0
##   629   0  0 11
##   632   0  0  7
##   637   0 18  0
##   642   0  9  0
##   648   0  9  0
##   651   0  9  0
##   654   0  0  9
##   67    0  0 10
##   674   0  0 11
##   675   0  8  0
##   678   4  0  0
##   68    0 12  0
##   687   0  0  8
##   689   0 17  0
##   690   0  0 10
##   692   0  5  0
##   693  10  0  0
##   694   0  8  0
##   695   7  0  0
##   699   0  0  9
##   7001  8  0  0
##   7004  0  7  0
##   7009  0  0  7
##   7013  0 10  0
##   723   0  6  0
##   727   9  0  0
##   729   0 11  8
##   730   0  0  3
##   736   8  0  0
##   753   0  4  0
##   756   9  0  0
##   763   0  7  0
##   771   8  0  0
##   773   5  0  0
##   779   0  0 10
##   78    0 10  0
##   801   0  5  0
##   806   7  0  0
##   808   0  0  6
##   811   0 11  0
##   817   0  7  0
##   82    0  0  9
##   823   0  0  9
##   824   0  5  0
##   826   0  0 10
##   835   8  0  0
##   836   9  0  0
##   838   0  0  9
##   839   0 11  0
##   845   4  0  0
##   860  12  0  0
##   863   0  0  9
##   867   4  0  0
##   869   0  0  4
##   875   7  0  0
##   89    0  0 10
##   94    0 10  0
##   96    6  0  0
##   97    0  9  0
##   H8    8  6  0
##   V211  0  9  0
##   V32   0 11  0
##   V53   0  0  9

Response variables

We’ll look at mass and tarsus length on days 6 and 14. Color = treatment: plot of chunk plot_pairs

Questions

  1. How much are the four measures of size affected by incubation temperature?
  2. How does this effect (if any) compare to other sources of variation?
  3. Is there an effect of “nest”, and is it correlated between the four measures?
  4. How correlated are the four measurements, after accounting for other variables?

A multivariate model

We’d like to fit \[\begin{aligned} \text{mass6}[i] &= \beta_{1,\text{treat[i]}} + \epsilon_{1,i} \\ \text{tarsus6}[i] &= \beta_{2,\text{treat[i]}} + \epsilon_{2,i} \\ \text{mass14}[i] &= \beta_{3,\text{treat[i]}} + \epsilon_{3,i} \\ \text{tarsus14}[i] &= \beta_{4,\text{treat[i]}} + \epsilon_{4,i} \\ \epsilon_i &\sim \MVN(0, \Sigma) . \end{aligned}\]

  • \(\beta\) is the mean effect of each treatment
  • \(\epsilon\) are the residuals
  • \(\Sigma\) is the covariance matrix of the residuals
bf1 <- brmsformula( mvbind(mass6, tars6, mass14, tars14) ~ treat )
priors <- set_prior("normal(0, 5)", class='b', resp=c("mass6", "tars6", "mass14", "tars14"))
bm1 <- brm( bf1 + set_rescor(TRUE), data=nestlings,
            prior=priors)
## Warning: Rows containing NAs were excluded from the model.
## Compiling Stan program...
## Start sampling

Notes:

  1. set_rescor(TRUE) enables residual correlations.

  2. With set_rescor(FALSE) we’d just have four independent linear models.

  3. This:

brmsformula(
    mvbind(mass6, tars6, mass14, tars14) ~ treat
)

is just shorthand for this:

( brmsformula( mass6 ~ treat ) +
   brmsformula( tars6 ~ treat ) +
   brmsformula( mass14 ~ treat ) +
   brmsformula( tars14 ~ treat ) )

Summary:

summary(bm1)
##  Family: MV(gaussian, gaussian, gaussian, gaussian) 
##   Links: mu = identity; sigma = identity
##          mu = identity; sigma = identity
##          mu = identity; sigma = identity
##          mu = identity; sigma = identity 
## Formula: mass6 ~ treat 
##          tars6 ~ treat 
##          mass14 ~ treat 
##          tars14 ~ treat 
##    Data: nestlings (Number of observations: 862) 
## Samples: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
##          total post-warmup samples = 4000
## 
## Population-Level Effects: 
##                  Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## mass6_Intercept      5.93      0.07     5.80     6.05 1.00     2169     2037
## tars6_Intercept     10.87      0.08    10.71    11.02 1.00     2245     2755
## mass14_Intercept    11.61      0.06    11.49    11.73 1.00     2866     2837
## tars14_Intercept    18.16      0.04    18.07    18.23 1.00     2795     2321
## mass6_treatM        -0.05      0.09    -0.22     0.12 1.00     2289     2677
## mass6_treatH         0.17      0.09    -0.01     0.33 1.00     2195     2452
## tars6_treatM         0.13      0.11    -0.08     0.34 1.00     2350     2833
## tars6_treatH         0.34      0.11     0.13     0.54 1.00     2188     2599
## mass14_treatM       -0.18      0.08    -0.33    -0.02 1.00     2936     3113
## mass14_treatH        0.00      0.08    -0.15     0.16 1.00     2881     3302
## tars14_treatM        0.28      0.05     0.17     0.38 1.00     3040     3226
## tars14_treatH        0.44      0.06     0.34     0.55 1.00     2784     2874
## 
## Family Specific Parameters: 
##              Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## sigma_mass6      1.00      0.02     0.95     1.05 1.00     2721     2855
## sigma_tars6      1.22      0.03     1.16     1.28 1.00     2514     2772
## sigma_mass14     0.91      0.02     0.87     0.95 1.00     3840     3243
## sigma_tars14     0.63      0.02     0.61     0.67 1.00     3650     3014
## 
## Residual Correlations: 
##                       Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## rescor(mass6,tars6)       0.85      0.01     0.83     0.87 1.00     2940     2985
## rescor(mass6,mass14)      0.43      0.03     0.37     0.48 1.00     3842     3043
## rescor(tars6,mass14)      0.27      0.03     0.21     0.34 1.00     3260     2787
## rescor(mass6,tars14)      0.45      0.03     0.39     0.50 1.00     2886     2622
## rescor(tars6,tars14)      0.53      0.02     0.48     0.57 1.00     3139     2995
## rescor(mass14,tars14)     0.54      0.02     0.49     0.59 1.00     3959     3460
## 
## Samples were drawn using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).

Let’s compare the residual correlations with the actual correlations in the data:

cor(nestlings[,4:7], use='pairwise')
##         mass6  tars6 mass14 tars14
## mass6  1.0000 0.8684 0.4433 0.4894
## tars6  0.8684 1.0000 0.2776 0.5347
## mass14 0.4433 0.2776 1.0000 0.5409
## tars14 0.4894 0.5347 0.5409 1.0000

Conditional effects:

plot of chunk ce1

But what about the “nest” effect?

We assumed that the residuals are independent between nestlings.

But, some are from the same nest?

Seems bad.

Also, recall that

  1. Is there an effect of “nest”, and is it correlated between the four measures?

Let’s add random effects by nest: \[\begin{aligned} \text{mass6}[i] &= \beta_{1,\text{treat[i]}} + \gamma_{1,\text{nest[i]}} + \epsilon_{1,i} \\ \text{tarsus6}[i] &= \beta_{2,\text{treat[i]}} + \gamma_{2,\text{nest[i]}} + \epsilon_{2,i} \\ \text{mass14}[i] &= \beta_{3,\text{treat[i]}} + \gamma_{3,\text{nest[i]}} + \epsilon_{3,i} \\ \text{tarsus14}[i] &= \beta_{4,\text{treat[i]}} + \gamma_{4,\text{nest[i]}} + \epsilon_{4,i} \\ \epsilon_i &\sim \MVN(0, \Sigma) \\ \gamma_n &\sim \MVN(0, S) \end{aligned}\]

Here:

  • \(\beta\) is the mean effect of each treatment
  • \(\epsilon\) are the residuals
  • \(\Sigma\) is the covariance matrix of the residuals
  • \(S\) is the covariance matrix of nest effects
bf2 <- brmsformula( mvbind(mass6, tars6, mass14, tars14) ~ treat + (1|n|nest) )
bm2 <- brm( bf2 + set_rescor(TRUE), data=nestlings, prior=priors )
## Warning: Rows containing NAs were excluded from the model.
## Compiling Stan program...
## Start sampling

New feature:

mvbind(mass6, tars6, mass14, tars14) ~ treat + (1|n|nest)

Here (1|n|nest) means that the random effects will be correlated with each other (and all others with label n).

Summary:

summary(bm2)
##  Family: MV(gaussian, gaussian, gaussian, gaussian) 
##   Links: mu = identity; sigma = identity
##          mu = identity; sigma = identity
##          mu = identity; sigma = identity
##          mu = identity; sigma = identity 
## Formula: mass6 ~ treat + (1 | n | nest) 
##          tars6 ~ treat + (1 | n | nest) 
##          mass14 ~ treat + (1 | n | nest) 
##          tars14 ~ treat + (1 | n | nest) 
##    Data: nestlings (Number of observations: 862) 
## Samples: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
##          total post-warmup samples = 4000
## 
## Group-Level Effects: 
## ~nest (Number of levels: 103) 
##                                        Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## sd(mass6_Intercept)                        0.61      0.06     0.51     0.73 1.00     1112     1626
## sd(tars6_Intercept)                        0.72      0.07     0.60     0.86 1.01      458      728
## sd(mass14_Intercept)                       0.67      0.05     0.57     0.79 1.01      779     1547
## sd(tars14_Intercept)                       0.40      0.04     0.34     0.48 1.00      725     1437
## cor(mass6_Intercept,tars6_Intercept)       0.72      0.06     0.60     0.82 1.00      791     1563
## cor(mass6_Intercept,mass14_Intercept)      0.35      0.10     0.13     0.54 1.01      455      722
## cor(tars6_Intercept,mass14_Intercept)      0.12      0.12    -0.11     0.36 1.00      512      899
## cor(mass6_Intercept,tars14_Intercept)      0.37      0.11     0.14     0.55 1.01      689     1255
## cor(tars6_Intercept,tars14_Intercept)      0.48      0.09     0.28     0.64 1.01      638     1300
## cor(mass14_Intercept,tars14_Intercept)     0.45      0.09     0.26     0.62 1.01      702     1290
## 
## Population-Level Effects: 
##                  Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## mass6_Intercept      5.92      0.11     5.71     6.14 1.00     1089     1656
## tars6_Intercept     10.82      0.13    10.57    11.08 1.00      893     1410
## mass14_Intercept    11.62      0.11    11.40    11.84 1.00      934     1497
## tars14_Intercept    18.12      0.07    17.98    18.27 1.01      953     1670
## mass6_treatM         0.02      0.14    -0.26     0.29 1.00     1019     1880
## mass6_treatH         0.10      0.14    -0.18     0.38 1.00     1018     1586
## tars6_treatM         0.26      0.17    -0.08     0.59 1.00      922     1625
## tars6_treatH         0.26      0.17    -0.08     0.59 1.00     1025     1562
## mass14_treatM       -0.12      0.14    -0.38     0.15 1.00      943     1752
## mass14_treatH       -0.01      0.14    -0.29     0.26 1.00      895     1475
## tars14_treatM        0.33      0.09     0.15     0.52 1.00     1085     1720
## tars14_treatH        0.44      0.10     0.24     0.61 1.00     1117     1893
## 
## Family Specific Parameters: 
##              Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## sigma_mass6      0.82      0.02     0.78     0.86 1.00     1944     2784
## sigma_tars6      1.01      0.03     0.96     1.06 1.00     1800     2629
## sigma_mass14     0.64      0.02     0.61     0.68 1.00     2459     2702
## sigma_tars14     0.51      0.01     0.49     0.54 1.00     2143     2778
## 
## Residual Correlations: 
##                       Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## rescor(mass6,tars6)       0.90      0.01     0.89     0.92 1.00     2276     2748
## rescor(mass6,mass14)      0.48      0.03     0.43     0.54 1.00     2177     3121
## rescor(tars6,mass14)      0.39      0.03     0.33     0.45 1.00     2138     2853
## rescor(mass6,tars14)      0.48      0.03     0.43     0.54 1.00     2089     2898
## rescor(tars6,tars14)      0.54      0.03     0.49     0.59 1.00     2177     2508
## rescor(mass14,tars14)     0.64      0.02     0.60     0.68 1.00     2291     2405
## 
## Samples were drawn using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).

Conditional effects:

plot of chunk ce2

Exercise:

Without fitting a model, how would you estimate the strength of correlation between (the nest effect on tarsus length at 6 days) and (the nest effect on mass at 6 days)?

Back to the questions

1. How much are the four measures of size affected by incubation temperature?

plot of chunk ce3

2. How does this effect (if any) compare to other sources of variation?

plot of chunk plot_effects

3. Is there an effect of “nest”

Well,* plot of chunk plot_hists

…and is it correlated between the four measures?

plot of chunk plot_effects_cor

4. How correlated are the four measurements, after accounting for other variables?

plot of chunk plot_effects_cor2

Model comparison

Did adding nest improve the model?

Comparing the models comes down to asking: (how much) does knowing measurements of others in the same nest help predict a nestling’s measurements?

this is a job for

crossvalidation

Fancy-pants crossvalidation

bm1_kfold <- kfold(bm1, K=5, save_fits=TRUE)
## Fitting model 1 out of 5
## Fitting model 2 out of 5
## Fitting model 3 out of 5
## Fitting model 4 out of 5
## Fitting model 5 out of 5
## Start sampling
## Start sampling
## Start sampling
## Start sampling
## Start sampling
bm2_kfold <- kfold(bm2, K=5, folds='stratified', group='nest', save_fits=TRUE)
## Fitting model 1 out of 5
## Fitting model 2 out of 5
## Fitting model 3 out of 5
## Fitting model 4 out of 5
## Fitting model 5 out of 5
## Start sampling
## Start sampling
## Start sampling
## Start sampling
## Start sampling

Extract out-of-sample predictions

## DOESN'T WORK?!?!
kpred1 <- kfold_predict(bm1_kfold)
kpred2 <- kfold_predict(bm2_kfold)
get_rmse <- function (kf) {
    varnames <- c("mass6", "tars6", "mass14", "tars14")
    rmse <- expand.grid(variable=varnames,
                        fold=1:5)
    rmse$rmse <- NA
    for (k in 1:5) {
        oos_indices <- kf$fits[[5+k]]
        stopifnot(all(oos_indices <= nrow(kf$data)))
        obs <- kf$data[oos_indices,]
        oos_data <- data.frame(lapply(kf$data, "[", oos_indices))
        exp <- predict(kf$fits[[k]], newdata=oos_data)[,"Estimate",]
        for (vn in varnames) {
            which_row <- which(rmse$variable == vn & rmse$fold == k)
            stopifnot(length(which_row) == 1)
            rmse$rmse[which_row] <- sqrt(mean((obs[,vn] - exp[,vn])^2))
        }
    }
    return(rmse)
}
rmse <- rbind(cbind(get_rmse(bm1_kfold), model=1),
              cbind(get_rmse(bm2_kfold), model=2))

Compute root-mean-squared-error

plot of chunk get_rmse

Conclusions?

// reveal.js plugins