Abordagem Longitudinal para Análise da Altura de Plantas em População F2 de Café - Dados sem plantas de replante

Introdução

Descreva aqui seus dados, seu experimento, os objetivos, etc … A análise de medidas repeidas precisa considerar a correlação entre medidas tomadas no mesmo indivíduo em tempos distintos, para isso, vamos especificar um modelo de variância do erro para o conjunto de dados. O ASReml possui uma gama de modelos de variâncias distintas. Vamos comparar nesse relatório, 5 modelos comumente utilizados na análise de medidas repetidas, utilizando o BIC como parâmetro.

Bibliotecas

library(ggplot2)
library(asreml)
## Online License checked out Sun Dec 15 22:46:52 2024
library(ggpubr)
library(agriutilities)
library(tidyr)
library(dplyr)
library(tibble)

Dados Fenotípicos

data <- xlsx::read.xlsx("dadosbruna.xlsx", header = TRUE, sheetIndex = 10,startRow = 2, colIndex = 6:9)
data$Trat[data$Trat == '23'|data$Trat == '22'] <- 'F2'
data$Pl <- rep(c(1:208),each = 6)

data <- within(data, {
  check <- ifelse(data$Trat == "T1"|data$Trat == "T2", 1, 0)
})

str(data)
## 'data.frame':    1248 obs. of  6 variables:
##  $ Trat  : chr  "T1" "T1" "T1" "T1" ...
##  $ Bl    : num  19 19 19 19 19 19 19 19 19 19 ...
##  $ Ano   : num  2018 2020 2021 2022 2023 ...
##  $ Altura: num  100 144 167 177 185 195 90 123 128 127 ...
##  $ Pl    : int  1 1 1 1 1 1 2 2 2 2 ...
##  $ check : num  1 1 1 1 1 1 1 1 1 1 ...
factors <- c(colnames(data)[c(1:3,5,6)]); factors
## [1] "Trat"  "Bl"    "Ano"   "Pl"    "check"
data[,factors] <- lapply(data[,factors], as.factor)

table(data$Trat)
## 
##   F2   T1   T2 
## 1116   72   60
table(data$Pl)
## 
##   1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19  20 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
##  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38  39  40 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
##  41  42  43  44  45  46  47  48  49  50  51  52  53  54  55  56  57  58  59  60 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
##  61  62  63  64  65  66  67  68  69  70  71  72  73  74  75  76  77  78  79  80 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
##  81  82  83  84  85  86  87  88  89  90  91  92  93  94  95  96  97  98  99 100 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
## 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
## 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
## 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
## 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
## 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 
##   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6 
## 201 202 203 204 205 206 207 208 
##   6   6   6   6   6   6   6   6
table(data$Bl)
## 
##  19  20  21  22 
## 282 312 324 330
table(data$Ano)
## 
## 2018 2020 2021 2022 2023 2024 
##  208  208  208  208  208  208
table(data$check)
## 
##    0    1 
## 1116  132
table(data$Trat,data$Pl)
##     
##      1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
##   F2 0 0 0 0 0 6 6 6 6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6
##   T1 6 6 6 0 0 0 0 0 0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
##   T2 0 0 0 6 6 0 0 0 0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
##     
##      29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53
##   F2  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  0  0  0  0  0  6
##   T1  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  6  6  6  0  0  0
##   T2  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  6  6  0
##     
##      54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78
##   F2  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6
##   T1  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
##   T2  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
##     
##      79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102
##   F2  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6   0   0   0
##   T1  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0   6   6   6
##   T2  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0   0   0   0
##     
##      103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120
##   F2   0   0   0   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6
##   T1   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##   T2   6   6   6   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##     
##      121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138
##   F2   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6
##   T1   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##   T2   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##     
##      139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156
##   F2   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   0   0   0
##   T1   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   6   6   6
##   T2   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##     
##      157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174
##   F2   0   0   0   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6
##   T1   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##   T2   6   6   6   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##     
##      175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192
##   F2   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6
##   T1   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##   T2   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##     
##      193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208
##   F2   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6   6
##   T1   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
##   T2   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0

Visualização dos Dados

ggplot(data, aes(x=Ano,y=Altura, group = Pl, colour=Pl))+
  geom_line()+
  facet_grid(.~Trat) +
  theme_minimal(base_size = 15) +
  theme(legend.position = 'none')

Modelo Uniforme

Assume que as correlações entre medidas de altura de um mesmo indivíduo são constantes ao longo do tempo

ans.uniform <- asreml(Altura ~ Trat + Ano + Trat:Ano,
                      random = ~ Bl + at(check,"0"):Pl:Ano,
                      residual = ~ Pl:corv(Ano),
                      data = data)
## ASReml Version 4.2 15/12/2024 22:47:02
##           LogLik        Sigma2     DF     wall
##  1     -3324.279           1.0    877   22:47:02
##  2     -3265.429           1.0    877   22:47:02
##  3     -3201.083           1.0    877   22:47:02
##  4     -3164.754           1.0    877   22:47:02
##  5     -3152.764           1.0    877   22:47:02
##  6     -3151.165           1.0    877   22:47:02
##  7     -3151.134           1.0    877   22:47:02
##  8     -3151.134           1.0    877   22:47:02
summary(ans.uniform)$varcomp
##                         component   std.error    z.ratio bound %ch
## Bl                     12.6045702 22.28425720  0.5656267     P   0
## at(check, '0'):Pl:Ano  96.7415662 41.13313392  2.3519133     P   0
## Pl:Ano!R                1.0000000          NA         NA     F   0
## Pl:Ano!Ano!cor          0.6865543  0.04589865 14.9580508     U   0
## Pl:Ano!Ano!var        704.4987837 72.75255262  9.6834923     P   0

Modelo Exponencial

Assume que as correlações decrescem de forma linear à medida em que o tempo entre medidas aumenta

ans.power <- asreml(Altura ~ Trat + Ano + Trat:Ano,
                    random = ~ Bl + at(check,"0"):Pl:Ano,
                    residual = ~ Pl:exp(Ano),
                    data = data)
## ASReml Version 4.2 15/12/2024 22:47:02
##           LogLik        Sigma2     DF     wall
##  1     -3337.447      638.3977    877   22:47:02
##  2     -3096.236      391.1353    877   22:47:02  (  1 restrained)
##  3     -2997.100      744.3443    877   22:47:02
##  4     -2995.507      750.7055    877   22:47:02
##  5     -2995.492      747.0486    877   22:47:02
##  6     -2995.492      746.5519    877   22:47:02
summary(ans.power)$varcomp
##                         component   std.error    z.ratio bound %ch
## Bl                     10.5367877 20.58583278  0.5118466     P 0.0
## at(check, '0'):Pl:Ano  34.6412883 12.86987389  2.6916572     P 0.1
## Pl:Ano!R              746.5518505 68.32983932 10.9257077     P 0.0
## Pl:Ano!Ano!pow          0.8623158  0.01856505 46.4483345     U 0.0
plot(data$Ano,ans.power$residuals)

O gráfico dos resíduos mostra que as variâncias residuais mudam com o tempo, por isso, vamos incluir uma matriz heterogênea para ajustar os dados

Modelo Exponencial Heterogêneo

ans.hetpower <- asreml(Altura ~ Trat + Ano + Trat:Ano,
                       random = ~ Bl + at(check,"0"):Pl:Ano,
                       residual = ~ Pl:exph(Ano),
                       data = data)
## ASReml Version 4.2 15/12/2024 22:47:03
##           LogLik        Sigma2     DF     wall
##  1     -3299.890           1.0    877   22:47:03  (  1 restrained)
##  2     -3179.320           1.0    877   22:47:03
##  3     -3046.865           1.0    877   22:47:03
##  4     -2975.090           1.0    877   22:47:03
##  5     -2942.566           1.0    877   22:47:03
##  6     -2935.098           1.0    877   22:47:03
##  7     -2934.318           1.0    877   22:47:03
##  8     -2934.272           1.0    877   22:47:03
##  9     -2934.269           1.0    877   22:47:03
summary(ans.hetpower)$varcomp
##                          component  std.error    z.ratio bound %ch
## Bl                       3.5917211   8.413649  0.4268922     P 0.3
## at(check, '0'):Pl:Ano   68.2504458  11.256844  6.0630178     P 0.1
## Pl:Ano!R                 1.0000000         NA         NA     F 0.0
## Pl:Ano!Ano!pow           0.8949425   0.015279 58.5733541     U 0.0
## Pl:Ano!Ano_2018        196.3640111  30.536670  6.4304331     P 0.1
## Pl:Ano!Ano_2020        414.8703818  59.125065  7.0168275     P 0.0
## Pl:Ano!Ano_2021        571.6775375  76.683063  7.4550692     P 0.1
## Pl:Ano!Ano_2022        798.9098647  98.231104  8.1329623     P 0.1
## Pl:Ano!Ano_2023       1116.3329070 133.287977  8.3753459     P 0.1
## Pl:Ano!Ano_2024       1206.7965881 141.517830  8.5275232     P 0.1

Modelo de Antedependência de Ordem 1

Permite que as correlações mudem ao longo do tempo

ans.ante <- asreml(Altura ~ Trat + Ano + Trat:Ano,
                   random = ~ Bl + at(check,"0"):Pl:Ano,
                   residual = ~ Pl:ante(Ano),
                   data = data)
## ASReml Version 4.2 15/12/2024 22:47:03
##           LogLik        Sigma2     DF     wall
##  1     -3029.248           1.0    877   22:47:03  (  1 restrained)
##  2     -2981.320           1.0    877   22:47:03
##  3     -2934.871           1.0    877   22:47:03
##  4     -2914.343           1.0    877   22:47:03
##  5     -2910.234           1.0    877   22:47:03
##  6     -2910.160           1.0    877   22:47:03
##  7     -2910.159           1.0    877   22:47:03
summary(ans.ante)$varcomp
##                          component    std.error     z.ratio bound %ch
## Bl                     5.612751634 8.751416e+00   0.6413535     P 0.2
## at(check, '0'):Pl:Ano 52.730480580 1.134128e+01   4.6494295     P 0.0
## Pl:Ano!R               1.000000000           NA          NA     F 0.0
## Pl:Ano!Ano_2018:2018   0.005900376 9.026076e-04   6.5370332     U 0.0
## Pl:Ano!Ano_2020:2018  -0.823540238 1.336249e-01  -6.1630742     U 0.0
## Pl:Ano!Ano_2020:2020   0.003765399 6.031628e-04   6.2427568     U 0.0
## Pl:Ano!Ano_2021:2020  -0.993856690 8.095966e-02 -12.2759491     U 0.0
## Pl:Ano!Ano_2021:2021   0.005151530 9.626788e-04   5.3512444     U 0.0
## Pl:Ano!Ano_2022:2021  -1.127516010 6.011906e-02 -18.7547177     U 0.0
## Pl:Ano!Ano_2022:2022   0.006191031 1.253463e-03   4.9391400     U 0.0
## Pl:Ano!Ano_2023:2022  -1.156298500 5.068391e-02 -22.8139171     U 0.0
## Pl:Ano!Ano_2023:2023   0.005552667 1.064946e-03   5.2140338     U 0.0
## Pl:Ano!Ano_2024:2023  -1.001284319 3.724974e-02 -26.8803057     U 0.0
## Pl:Ano!Ano_2024:2024   0.006946649 1.381945e-03   5.0267192     U 0.0

Modelo Não Estruturado

Permite que todas as correlações entre tempos sejam únicas, mas o número de parâmetros é muito superior

ans.uns <- asreml(Altura ~ Trat + Ano + Trat:Ano,
                  random = ~ Bl + at(check,"0"):Pl:Ano,
                  residual = ~ Pl:us(Ano),
                  data = data)
## ASReml Version 4.2 15/12/2024 22:47:04
##           LogLik        Sigma2     DF     wall
##  1     -2945.954           1.0    877   22:47:04  (  1 restrained)
##  2     -2923.667           1.0    877   22:47:04
##  3     -2907.354           1.0    877   22:47:04
##  4     -2901.272           1.0    877   22:47:04
##  5     -2899.944           1.0    877   22:47:04
##  6     -2899.895           1.0    877   22:47:04
##  7     -2899.895           1.0    877   22:47:04
summary(ans.uns)$varcomp
##                         component std.error  z.ratio bound %ch
## Bl                       8.968396  11.61207 0.772334     P 0.0
## at(check, '0'):Pl:Ano   65.432592  11.39871 5.740351     P 0.1
## Pl:Ano!R                 1.000000        NA       NA     F 0.0
## Pl:Ano!Ano_2018:2018   163.612651  25.90351 6.316234     P 0.0
## Pl:Ano!Ano_2020:2018   140.601855  27.51253 5.110467     P 0.0
## Pl:Ano!Ano_2020:2020   361.157514  49.76830 7.256778     P 0.0
## Pl:Ano!Ano_2021:2018   104.327662  31.84927 3.275669     P 0.0
## Pl:Ano!Ano_2021:2020   383.047640  52.63002 7.278121     P 0.0
## Pl:Ano!Ano_2021:2021   561.863789  72.84794 7.712830     P 0.0
## Pl:Ano!Ano_2022:2018   125.909480  39.06715 3.222899     P 0.0
## Pl:Ano!Ano_2022:2020   440.002017  62.89773 6.995515     P 0.0
## Pl:Ano!Ano_2022:2021   642.600446  82.37923 7.800516     P 0.0
## Pl:Ano!Ano_2022:2022   874.386412 109.34151 7.996839     P 0.0
## Pl:Ano!Ano_2023:2018   108.412884  47.47613 2.283524     P 0.0
## Pl:Ano!Ano_2023:2020   488.303616  75.30962 6.483948     P 0.0
## Pl:Ano!Ano_2023:2021   765.892347 100.41965 7.626917     P 0.0
## Pl:Ano!Ano_2023:2022  1008.583394 127.06653 7.937444     P 0.0
## Pl:Ano!Ano_2023:2023  1340.120999 167.14635 8.017650     P 0.0
## Pl:Ano!Ano_2024:2018   130.446416  50.68974 2.573428     P 0.0
## Pl:Ano!Ano_2024:2020   514.857627  80.09740 6.427894     P 0.0
## Pl:Ano!Ano_2024:2021   755.483388 103.65761 7.288257     P 0.0
## Pl:Ano!Ano_2024:2022  1053.134684 134.45005 7.832907     P 0.0
## Pl:Ano!Ano_2024:2023  1368.931441 171.05748 8.002757     P 0.0
## Pl:Ano!Ano_2024:2024  1532.645627 190.79238 8.033055     P 0.0

Escolha do modelo mais parcimonioso

Critério: BIC

summary(ans.uniform)$bic
## [1] 6329.374
## attr(,"parameters")
## [1] 4
summary(ans.power)$bic
## [1] 6018.09
## attr(,"parameters")
## [1] 4
summary(ans.hetpower)$bic
## [1] 5929.526
## attr(,"parameters")
## [1] 9
summary(ans.ante)$bic
## [1] 5908.413
## attr(,"parameters")
## [1] 13
summary(ans.uns)$bic
## [1] 5955.649
## attr(,"parameters")
## [1] 23

A comparação entre os diferentes valores de BIC mostra que o Modelo Antedependência de Ordem 1 foi o que melhor se ajustou aos dados, a partir dele faremos as estimações e predições.

Teste Wald

Vamos usar o teste Wald para testar a significância dos efeitos fixos do modelo ans.ante

wald(ans.ante)
## Wald tests for fixed effects.
## Response: Altura
## Terms added sequentially; adjusted for those above.
## 
##               Df Sum of Sq Wald statistic Pr(Chisq)    
## (Intercept)    1    3456.7         3456.7 < 2.2e-16 ***
## Trat           2      19.9           19.9 4.673e-05 ***
## Ano            5    1439.1         1439.1 < 2.2e-16 ***
## Trat:Ano      10       9.2            9.2    0.5175    
## residual (MS)          1.0                             
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Interpretações: Os tratamentos (T1, T2 e F2) possuem estimativas de alturas diferentes entre si (p-valor = 5.916e-05), as estimativas médias de alturas nos diferentes anos diferem entre si (p-valor = < 2.2e-16). A interação entre tratamentos e anos é não-significativa (p-valor = 0.3934), ou seja, não há evidências de que os tratamentos crescem em padrões diferentes ao longo do tempo (ex:F2 cresce a uma mesma taxa que o T1 e que o T2).

Estimativas de efeitos fixos

Podemos obter as estimativas (BLUEs) para cada efeito fixo

BLUES.trat    <- predict(ans.ante, classify = "Trat")$pvals;BLUES.trat
## ASReml Version 4.2 15/12/2024 22:47:04
##           LogLik        Sigma2     DF     wall
##  1     -2910.159           1.0    877   22:47:04
##  2     -2910.159           1.0    877   22:47:04
##  3     -2910.159           1.0    877   22:47:04
## 
## Notes:
## - The predictions are obtained by averaging across the hypertable
##   calculated from model terms constructed solely from factors in
##   the averaging and classify sets.
## - Use 'average' to move ignored factors into the averaging set.
## - The simple averaging set: Ano
## - The ignored set: Bl,at(check, '0'),Pl
## 
##   Trat predicted.value std.error    status
## 1   F2        150.5048  2.413144 Estimable
## 2   T1        153.1114  7.758199 Estimable
## 3   T2        146.2257  9.986474 Estimable
BLUES.ano     <- predict(ans.ante, classify = "Ano")$pvals;BLUES.ano
## ASReml Version 4.2 15/12/2024 22:47:05
##           LogLik        Sigma2     DF     wall
##  1     -2910.159           1.0    877   22:47:05
##  2     -2910.159           1.0    877   22:47:05
##  3     -2910.159           1.0    877   22:47:05
## 
## Notes:
## - The predictions are obtained by averaging across the hypertable
##   calculated from model terms constructed solely from factors in
##   the averaging and classify sets.
## - Use 'average' to move ignored factors into the averaging set.
## - The simple averaging set: Trat
## - The ignored set: Bl,at(check, '0'),Pl
## 
##    Ano predicted.value std.error    status
## 1 2018        82.28388  2.576507 Estimable
## 2 2020       135.17223  3.614587 Estimable
## 3 2021       145.73927  4.340272 Estimable
## 4 2022       164.25772  5.335286 Estimable
## 5 2023       182.78414  6.561749 Estimable
## 6 2024       189.44641  6.895718 Estimable
BLUES.TratAno <- predict(ans.ante, classify = "Trat:Ano")$pvals;BLUES.TratAno
## ASReml Version 4.2 15/12/2024 22:47:05
##           LogLik        Sigma2     DF     wall
##  1     -2910.159           1.0    877   22:47:05
##  2     -2910.159           1.0    877   22:47:05
##  3     -2910.159           1.0    877   22:47:05
## 
## Notes:
## - The predictions are obtained by averaging across the hypertable
##   calculated from model terms constructed solely from factors in
##   the averaging and classify sets.
## - Use 'average' to move ignored factors into the averaging set.
## - The ignored set: Bl,at(check, '0'),Pl
## 
##    Trat  Ano predicted.value std.error    status
## 1    F2 2018        88.37568  1.742400 Estimable
## 2    F2 2020       140.30740  2.143272 Estimable
## 3    F2 2021       146.84186  2.453328 Estimable
## 4    F2 2022       162.06387  2.892330 Estimable
## 5    F2 2023       177.87647  3.464520 Estimable
## 6    F2 2024       187.56327  3.625541 Estimable
## 7    T1 2018        79.94472  4.286868 Estimable
## 8    T1 2020       136.84472  6.283438 Estimable
## 9    T1 2021       150.84472  7.643956 Estimable
## 10   T1 2022       170.34472  9.489060 Estimable
## 11   T1 2023       186.44472 11.744119 Estimable
## 12   T1 2024       194.24472 12.355758 Estimable
## 13   T2 2018        78.53124  5.481086 Estimable
## 14   T2 2020       128.36457  8.075657 Estimable
## 15   T2 2021       139.53124  9.838548 Estimable
## 16   T2 2022       160.36457 12.226366 Estimable
## 17   T2 2023       184.03124 15.142241 Estimable
## 18   T2 2024       186.53124 15.932822 Estimable
xlsx::write.xlsx2(BLUES.trat, paste0("genetic_parameters.xlsx"), sheetName = "BLUES.trat")
xlsx::write.xlsx2(BLUES.ano, paste0("genetic_parameters.xlsx"), sheetName = "BLUES.ano", append = T)
xlsx::write.xlsx2(BLUES.TratAno, paste0("genetic_parameters.xlsx"), sheetName = "BLUES.TratAno", append = T)

Teste Razão de Verossimilhança - LRT

Vamos testar a significância dos componentes de variância de planta. Para isso vamos obter o modelo reduzido para posterior comparação.

ans.ante.red <- asreml(Altura ~ Trat + Ano + Trat:Ano,
                   random = ~ Bl, 
                   residual = ~ Pl:ante(Ano),
                   data = data)
## ASReml Version 4.2 15/12/2024 22:47:06
##           LogLik        Sigma2     DF     wall
##  1     -3130.391           1.0    877   22:47:06  (  1 restrained)
##  2     -3043.319           1.0    877   22:47:06
##  3     -2960.374           1.0    877   22:47:06
##  4     -2926.423           1.0    877   22:47:06
##  5     -2919.966           1.0    877   22:47:06
##  6     -2919.767           1.0    877   22:47:06
##  7     -2919.764           1.0    877   22:47:06
lrt.asreml(ans.ante,ans.ante.red)
## Likelihood ratio test(s) assuming nested random models.
## (See Self & Liang, 1987)
## 
##                       df LR-statistic Pr(Chisq)    
## ans.ante/ans.ante.red  1        19.21 5.854e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Interpretação: Há variação de altura entre plantas na população F2 nos diferentes anos (p-valor = 5.883e-05).

Predição de valores genéticos

Podemos obter as predições de cada genótipo (BLUPs) dentro de cada ano.

pred <- predict(ans.ante, classify = "check:Pl:Ano", levels = list(check = '0'))$pvals; pred
## ASReml Version 4.2 15/12/2024 22:47:06
##           LogLik        Sigma2     DF     wall
##  1     -2910.159           1.0    877   22:47:06
##  2     -2910.159           1.0    877   22:47:06
##  3     -2910.159           1.0    877   22:47:06
## 
## Notes:
## - The predictions are obtained by averaging across the hypertable
##   calculated from model terms constructed solely from factors in
##   the averaging and classify sets.
## - Use 'average' to move ignored factors into the averaging set.
## - at(check, '0') is evaluated at 1.00000
## - The simple averaging set: Trat
## - The ignored set: Bl
## 
##      check  Pl  Ano predicted.value std.error    status
## 1        0   1 2018        82.28388  7.704892 Estimable
## 2        0   1 2020       135.17223  8.111240 Estimable
## 3        0   1 2021       145.73927  8.459606 Estimable
## 4        0   1 2022       164.25772  9.010674 Estimable
## 5        0   1 2023       182.78414  9.786905 Estimable
## 6        0   1 2024       189.44641 10.013885 Estimable
## 7        0   2 2018        82.28388  7.704892 Estimable
## 8        0   2 2020       135.17223  8.111240 Estimable
## 9        0   2 2021       145.73927  8.459606 Estimable
## 10       0   2 2022       164.25772  9.010674 Estimable
## 11       0   2 2023       182.78414  9.786905 Estimable
## 12       0   2 2024       189.44641 10.013885 Estimable
## 13       0   3 2018        82.28388  7.704892 Estimable
## 14       0   3 2020       135.17223  8.111240 Estimable
## 15       0   3 2021       145.73927  8.459606 Estimable
## 16       0   3 2022       164.25772  9.010674 Estimable
## 17       0   3 2023       182.78414  9.786905 Estimable
## 18       0   3 2024       189.44641 10.013885 Estimable
## 19       0   4 2018        82.28388  7.704892 Estimable
## 20       0   4 2020       135.17223  8.111240 Estimable
## 21       0   4 2021       145.73927  8.459606 Estimable
## 22       0   4 2022       164.25772  9.010674 Estimable
## 23       0   4 2023       182.78414  9.786905 Estimable
## 24       0   4 2024       189.44641 10.013885 Estimable
## 25       0   5 2018        82.28388  7.704892 Estimable
## 26       0   5 2020       135.17223  8.111240 Estimable
## 27       0   5 2021       145.73927  8.459606 Estimable
## 28       0   5 2022       164.25772  9.010674 Estimable
## 29       0   5 2023       182.78414  9.786905 Estimable
## 30       0   5 2024       189.44641 10.013885 Estimable
## 31       0   6 2018        77.86870  6.591210 Estimable
## 32       0   6 2020       135.54125  7.105057 Estimable
## 33       0   6 2021       147.65275  7.254742 Estimable
## 34       0   6 2022       159.32271  7.880627 Estimable
## 35       0   6 2023       184.57346  8.848074 Estimable
## 36       0   6 2024       188.95366  9.408065 Estimable
## 37       0   7 2018        82.28388  7.704892 Estimable
## 38       0   7 2020       135.17223  8.111240 Estimable
## 39       0   7 2021       145.73927  8.459606 Estimable
## 40       0   7 2022       164.25772  9.010674 Estimable
## 41       0   7 2023       182.78414  9.786905 Estimable
## 42       0   7 2024       189.44641 10.013885 Estimable
## 43       0   8 2018        81.60504  6.591210 Estimable
## 44       0   8 2020       133.23192  7.105057 Estimable
## 45       0   8 2021       144.11482  7.254742 Estimable
## 46       0   8 2022       166.56251  7.880627 Estimable
## 47       0   8 2023       184.75208  8.848074 Estimable
## 48       0   8 2024       187.37759  9.408065 Estimable
## 49       0   9 2018        79.01389  6.591210 Estimable
## 50       0   9 2020       132.13001  7.105057 Estimable
## 51       0   9 2021       149.85145  7.254742 Estimable
## 52       0   9 2022       160.67722  7.880627 Estimable
## 53       0   9 2023       181.56095  8.848074 Estimable
## 54       0   9 2024       192.95285  9.408065 Estimable
## 55       0  10 2018        77.66690  6.591210 Estimable
## 56       0  10 2020       132.94500  7.105057 Estimable
## 57       0  10 2021       150.19687  7.254742 Estimable
## 58       0  10 2022       161.74599  7.880627 Estimable
## 59       0  10 2023       193.69220  8.848074 Estimable
## 60       0  10 2024       180.38275  9.408065 Estimable
## 61       0  11 2018        77.59351  6.591210 Estimable
## 62       0  11 2020       128.92835  7.105057 Estimable
## 63       0  11 2021       148.78595  7.254742 Estimable
## 64       0  11 2022       161.95025  7.880627 Estimable
## 65       0  11 2023       184.13690  8.848074 Estimable
## 66       0  11 2024       192.56047  9.408065 Estimable
## 67       0  12 2018        78.87216  6.591210 Estimable
## 68       0  12 2020       130.88562  7.105057 Estimable
## 69       0  12 2021       152.95791  7.254742 Estimable
## 70       0  12 2022       165.22885  7.880627 Estimable
## 71       0  12 2023       176.63463  8.848074 Estimable
## 72       0  12 2024       195.89814  9.408065 Estimable
## 73       0  13 2018        82.28388  7.704892 Estimable
## 74       0  13 2020       135.17223  8.111240 Estimable
## 75       0  13 2021       145.73927  8.459606 Estimable
## 76       0  13 2022       164.25772  9.010674 Estimable
## 77       0  13 2023       182.78414  9.786905 Estimable
## 78       0  13 2024       189.44641 10.013885 Estimable
## 79       0  14 2018        78.62090  6.591210 Estimable
## 80       0  14 2020       134.38130  7.105057 Estimable
## 81       0  14 2021       151.52887  7.254742 Estimable
## 82       0  14 2022       168.24062  7.880627 Estimable
## 83       0  14 2023       172.88731  8.848074 Estimable
## 84       0  14 2024       192.25380  9.408065 Estimable
## 85       0  15 2018        82.60163  6.591210 Estimable
## 86       0  15 2020       129.59194  7.105057 Estimable
## 87       0  15 2021       149.04893  7.254742 Estimable
## 88       0  15 2022       159.32335  7.880627 Estimable
## 89       0  15 2023       186.30181  8.848074 Estimable
## 90       0  15 2024       190.46756  9.408065 Estimable
## 91       0  16 2018        79.91118  6.591210 Estimable
## 92       0  16 2020       134.79125  7.105057 Estimable
## 93       0  16 2021       146.99758  7.254742 Estimable
## 94       0  16 2022       161.74351  7.880627 Estimable
## 95       0  16 2023       178.55898  8.848074 Estimable
## 96       0  16 2024       193.77277  9.408065 Estimable
## 97       0  17 2018        82.98674  6.591210 Estimable
## 98       0  17 2020       127.99745  7.105057 Estimable
## 99       0  17 2021       151.79820  7.254742 Estimable
## 100      0  17 2022       161.82781  7.880627 Estimable
## 101      0  17 2023       183.52427  8.848074 Estimable
## 102      0  17 2024       189.72879  9.408065 Estimable
## 103      0  18 2018        75.50926  6.591210 Estimable
## 104      0  18 2020       140.13846  7.105057 Estimable
## 105      0  18 2021       158.04661  7.254742 Estimable
## 106      0  18 2022       148.99939  7.880627 Estimable
## 107      0  18 2023       183.93639  8.848074 Estimable
## 108      0  18 2024       192.51679  9.408065 Estimable
## 109      0  19 2018        82.28388  7.704892 Estimable
## 110      0  19 2020       135.17223  8.111240 Estimable
## 111      0  19 2021       145.73927  8.459606 Estimable
## 112      0  19 2022       164.25772  9.010674 Estimable
## 113      0  19 2023       182.78414  9.786905 Estimable
## 114      0  19 2024       189.44641 10.013885 Estimable
## 115      0  20 2018        83.63690  6.591210 Estimable
## 116      0  20 2020       133.20867  7.105057 Estimable
## 117      0  20 2021       146.66809  7.254742 Estimable
## 118      0  20 2022       163.50329  7.880627 Estimable
## 119      0  20 2023       184.81549  8.848074 Estimable
## 120      0  20 2024       188.72823  9.408065 Estimable
## 121      0  21 2018        82.28388  7.704892 Estimable
## 122      0  21 2020       135.17223  8.111240 Estimable
## 123      0  21 2021       145.73927  8.459606 Estimable
## 124      0  21 2022       164.25772  9.010674 Estimable
## 125      0  21 2023       182.78414  9.786905 Estimable
## 126      0  21 2024       189.44641 10.013885 Estimable
## 127      0  22 2018        82.28388  7.704892 Estimable
## 128      0  22 2020       135.17223  8.111240 Estimable
## 129      0  22 2021       145.73927  8.459606 Estimable
## 130      0  22 2022       164.25772  9.010674 Estimable
## 131      0  22 2023       182.78414  9.786905 Estimable
## 132      0  22 2024       189.44641 10.013885 Estimable
## 133      0  23 2018        82.28388  7.704892 Estimable
## 134      0  23 2020       135.17223  8.111240 Estimable
## 135      0  23 2021       145.73927  8.459606 Estimable
## 136      0  23 2022       164.25772  9.010674 Estimable
## 137      0  23 2023       182.78414  9.786905 Estimable
## 138      0  23 2024       189.44641 10.013885 Estimable
## 139      0  24 2018        86.55444  6.591210 Estimable
## 140      0  24 2020       137.38511  7.105057 Estimable
## 141      0  24 2021       151.24021  7.254742 Estimable
## 142      0  24 2022       160.14276  7.880627 Estimable
## 143      0  24 2023       182.64629  8.848074 Estimable
## 144      0  24 2024       188.15444  9.408065 Estimable
## 145      0  25 2018        82.28388  7.704892 Estimable
## 146      0  25 2020       135.17223  8.111240 Estimable
## 147      0  25 2021       145.73927  8.459606 Estimable
## 148      0  25 2022       164.25772  9.010674 Estimable
## 149      0  25 2023       182.78414  9.786905 Estimable
## 150      0  25 2024       189.44641 10.013885 Estimable
## 151      0  26 2018        91.02171  6.591219 Estimable
## 152      0  26 2020       121.69216  7.106143 Estimable
## 153      0  26 2021       146.79284  7.284732 Estimable
## 154      0  26 2022       163.68362  8.407764 Estimable
## 155      0  26 2023       182.78414  9.786905 Estimable
## 156      0  26 2024       189.44641 10.013885 Estimable
## 157      0  27 2018        83.99499  6.834141 Estimable
## 158      0  27 2020       135.17223  8.111240 Estimable
## 159      0  27 2021       145.73927  8.459606 Estimable
## 160      0  27 2022       164.25772  9.010674 Estimable
## 161      0  27 2023       182.78414  9.786905 Estimable
## 162      0  27 2024       189.44641 10.013885 Estimable
## 163      0  28 2018        82.28388  7.704892 Estimable
## 164      0  28 2020       135.17223  8.111240 Estimable
## 165      0  28 2021       145.73927  8.459606 Estimable
## 166      0  28 2022       164.25772  9.010674 Estimable
## 167      0  28 2023       182.78414  9.786905 Estimable
## 168      0  28 2024       189.44641 10.013885 Estimable
## 169      0  29 2018        82.28388  7.704892 Estimable
## 170      0  29 2020       135.17223  8.111240 Estimable
## 171      0  29 2021       145.73927  8.459606 Estimable
## 172      0  29 2022       164.25772  9.010674 Estimable
## 173      0  29 2023       182.78414  9.786905 Estimable
## 174      0  29 2024       189.44641 10.013885 Estimable
## 175      0  30 2018        82.28388  7.704892 Estimable
## 176      0  30 2020       135.17223  8.111240 Estimable
## 177      0  30 2021       145.73927  8.459606 Estimable
## 178      0  30 2022       164.25772  9.010674 Estimable
## 179      0  30 2023       182.78414  9.786905 Estimable
## 180      0  30 2024       189.44641 10.013885 Estimable
## 181      0  31 2018        82.28388  7.704892 Estimable
## 182      0  31 2020       135.17223  8.111240 Estimable
## 183      0  31 2021       145.73927  8.459606 Estimable
## 184      0  31 2022       164.25772  9.010674 Estimable
## 185      0  31 2023       182.78414  9.786905 Estimable
## 186      0  31 2024       189.44641 10.013885 Estimable
## 187      0  32 2018        82.28388  7.704892 Estimable
## 188      0  32 2020       135.17223  8.111240 Estimable
## 189      0  32 2021       145.73927  8.459606 Estimable
## 190      0  32 2022       164.25772  9.010674 Estimable
## 191      0  32 2023       182.78414  9.786905 Estimable
## 192      0  32 2024       189.44641 10.013885 Estimable
## 193      0  33 2018        78.63988  6.591210 Estimable
## 194      0  33 2020       137.91679  7.105057 Estimable
## 195      0  33 2021       147.35820  7.254742 Estimable
## 196      0  33 2022       160.40759  7.880627 Estimable
## 197      0  33 2023       186.09062  8.848074 Estimable
## 198      0  33 2024       187.73519  9.408065 Estimable
## 199      0  34 2018        82.28388  7.704892 Estimable
## 200      0  34 2020       135.17223  8.111240 Estimable
## 201      0  34 2021       145.73927  8.459606 Estimable
## 202      0  34 2022       164.25772  9.010674 Estimable
## 203      0  34 2023       182.78414  9.786905 Estimable
## 204      0  34 2024       189.44641 10.013885 Estimable
## 205      0  35 2018        82.28388  7.704892 Estimable
## 206      0  35 2020       135.17223  8.111240 Estimable
## 207      0  35 2021       145.73927  8.459606 Estimable
## 208      0  35 2022       164.25772  9.010674 Estimable
## 209      0  35 2023       182.78414  9.786905 Estimable
## 210      0  35 2024       189.44641 10.013885 Estimable
## 211      0  36 2018        89.18614  6.591210 Estimable
## 212      0  36 2020       138.48683  7.105057 Estimable
## 213      0  36 2021       146.44269  7.254742 Estimable
## 214      0  36 2022       158.43053  7.880627 Estimable
## 215      0  36 2023       185.66902  8.848074 Estimable
## 216      0  36 2024       190.28757  9.408065 Estimable
## 217      0  37 2018        87.18645  6.591210 Estimable
## 218      0  37 2020       135.33280  7.105057 Estimable
## 219      0  37 2021       147.57549  7.254742 Estimable
## 220      0  37 2022       163.21502  7.880627 Estimable
## 221      0  37 2023       182.01718  8.848074 Estimable
## 222      0  37 2024       187.99572  9.408065 Estimable
## 223      0  38 2018        84.80920  6.591210 Estimable
## 224      0  38 2020       132.94304  7.105057 Estimable
## 225      0  38 2021       144.35496  7.254742 Estimable
## 226      0  38 2022       166.56753  7.880627 Estimable
## 227      0  38 2023       183.24541  8.848074 Estimable
## 228      0  38 2024       189.64380  9.408065 Estimable
## 229      0  39 2018        82.16037  6.591210 Estimable
## 230      0  39 2020       130.96336  7.105057 Estimable
## 231      0  39 2021       151.53921  7.254742 Estimable
## 232      0  39 2022       158.67334  7.880627 Estimable
## 233      0  39 2023       188.23435  8.848074 Estimable
## 234      0  39 2024       186.96683  9.408065 Estimable
## 235      0  40 2018        87.07529  6.591210 Estimable
## 236      0  40 2020       128.52918  7.105057 Estimable
## 237      0  40 2021       151.06700  7.254742 Estimable
## 238      0  40 2022       170.35841  7.880627 Estimable
## 239      0  40 2023       175.42202  8.848074 Estimable
## 240      0  40 2024       190.24494  9.408065 Estimable
## 241      0  41 2018        83.65546  6.591210 Estimable
## 242      0  41 2020       135.02287  7.105057 Estimable
## 243      0  41 2021       145.80329  7.254742 Estimable
## 244      0  41 2022       172.75715  7.880627 Estimable
## 245      0  41 2023       178.66380  8.848074 Estimable
## 246      0  41 2024       191.07113  9.408065 Estimable
## 247      0  42 2018        82.28388  7.704892 Estimable
## 248      0  42 2020       135.17223  8.111240 Estimable
## 249      0  42 2021       145.73927  8.459606 Estimable
## 250      0  42 2022       164.25772  9.010674 Estimable
## 251      0  42 2023       182.78414  9.786905 Estimable
## 252      0  42 2024       189.44641 10.013885 Estimable
## 253      0  43 2018        81.69053  6.591210 Estimable
## 254      0  43 2020       134.89385  7.105057 Estimable
## 255      0  43 2021       141.96427  7.254742 Estimable
## 256      0  43 2022       161.65597  7.880627 Estimable
## 257      0  43 2023       183.14583  8.848074 Estimable
## 258      0  43 2024       192.31643  9.408065 Estimable
## 259      0  44 2018        88.99331  6.591210 Estimable
## 260      0  44 2020       137.68747  7.105057 Estimable
## 261      0  44 2021       137.08655  7.254742 Estimable
## 262      0  44 2022       168.11631  7.880627 Estimable
## 263      0  44 2023       186.59528  8.848074 Estimable
## 264      0  44 2024       187.85205  9.408065 Estimable
## 265      0  45 2018        81.25629  6.591210 Estimable
## 266      0  45 2020       134.88328  7.105057 Estimable
## 267      0  45 2021       144.01729  7.254742 Estimable
## 268      0  45 2022       168.75664  7.880627 Estimable
## 269      0  45 2023       182.84718  8.848074 Estimable
## 270      0  45 2024       189.52677  9.408065 Estimable
## 271      0  46 2018        78.51794  6.591210 Estimable
## 272      0  46 2020       135.83846  7.105057 Estimable
## 273      0  46 2021       147.86818  7.254742 Estimable
## 274      0  46 2022       163.88087  7.880627 Estimable
## 275      0  46 2023       183.74202  8.848074 Estimable
## 276      0  46 2024       188.44179  9.408065 Estimable
## 277      0  47 2018        81.10474  6.591210 Estimable
## 278      0  47 2020       134.35484  7.105057 Estimable
## 279      0  47 2021       142.70825  7.254742 Estimable
## 280      0  47 2022       166.70666  7.880627 Estimable
## 281      0  47 2023       184.71411  8.848074 Estimable
## 282      0  47 2024       188.69933  9.408065 Estimable
## 283      0  48 2018        82.28388  7.704892 Estimable
## 284      0  48 2020       135.17223  8.111240 Estimable
## 285      0  48 2021       145.73927  8.459606 Estimable
## 286      0  48 2022       164.25772  9.010674 Estimable
## 287      0  48 2023       182.78414  9.786905 Estimable
## 288      0  48 2024       189.44641 10.013885 Estimable
## 289      0  49 2018        82.28388  7.704892 Estimable
## 290      0  49 2020       135.17223  8.111240 Estimable
## 291      0  49 2021       145.73927  8.459606 Estimable
## 292      0  49 2022       164.25772  9.010674 Estimable
## 293      0  49 2023       182.78414  9.786905 Estimable
## 294      0  49 2024       189.44641 10.013885 Estimable
## 295      0  50 2018        82.28388  7.704892 Estimable
## 296      0  50 2020       135.17223  8.111240 Estimable
## 297      0  50 2021       145.73927  8.459606 Estimable
## 298      0  50 2022       164.25772  9.010674 Estimable
## 299      0  50 2023       182.78414  9.786905 Estimable
## 300      0  50 2024       189.44641 10.013885 Estimable
## 301      0  51 2018        82.28388  7.704892 Estimable
## 302      0  51 2020       135.17223  8.111240 Estimable
## 303      0  51 2021       145.73927  8.459606 Estimable
## 304      0  51 2022       164.25772  9.010674 Estimable
## 305      0  51 2023       182.78414  9.786905 Estimable
## 306      0  51 2024       189.44641 10.013885 Estimable
## 307      0  52 2018        82.28388  7.704892 Estimable
## 308      0  52 2020       135.17223  8.111240 Estimable
## 309      0  52 2021       145.73927  8.459606 Estimable
## 310      0  52 2022       164.25772  9.010674 Estimable
## 311      0  52 2023       182.78414  9.786905 Estimable
## 312      0  52 2024       189.44641 10.013885 Estimable
## 313      0  53 2018        82.28388  7.704892 Estimable
## 314      0  53 2020       135.17223  8.111240 Estimable
## 315      0  53 2021       145.73927  8.459606 Estimable
## 316      0  53 2022       164.25772  9.010674 Estimable
## 317      0  53 2023       182.78414  9.786905 Estimable
## 318      0  53 2024       189.44641 10.013885 Estimable
## 319      0  54 2018        82.28388  7.704892 Estimable
## 320      0  54 2020       135.17223  8.111240 Estimable
## 321      0  54 2021       145.73927  8.459606 Estimable
## 322      0  54 2022       164.25772  9.010674 Estimable
## 323      0  54 2023       182.78414  9.786905 Estimable
## 324      0  54 2024       189.44641 10.013885 Estimable
## 325      0  55 2018        79.71961  6.596646 Estimable
## 326      0  55 2020       130.89800  7.106743 Estimable
## 327      0  55 2021       146.80095  7.256008 Estimable
## 328      0  55 2022       164.35507  7.880928 Estimable
## 329      0  55 2023       183.99006  8.847289 Estimable
## 330      0  55 2024       189.85035  9.407862 Estimable
## 331      0  56 2018        80.35129  6.596646 Estimable
## 332      0  56 2020       138.95464  7.106743 Estimable
## 333      0  56 2021       151.73790  7.256008 Estimable
## 334      0  56 2022       156.33493  7.880928 Estimable
## 335      0  56 2023       181.33527  8.847289 Estimable
## 336      0  56 2024       194.48576  9.407862 Estimable
## 337      0  57 2018        79.53113  6.596646 Estimable
## 338      0  57 2020       137.59895  7.106743 Estimable
## 339      0  57 2021       146.30974  7.256008 Estimable
## 340      0  57 2022       160.05632  7.880928 Estimable
## 341      0  57 2023       186.58557  8.847289 Estimable
## 342      0  57 2024       187.86963  9.407862 Estimable
## 343      0  58 2018        77.01038  6.596646 Estimable
## 344      0  58 2020       137.29898  7.106743 Estimable
## 345      0  58 2021       148.37893  7.256008 Estimable
## 346      0  58 2022       163.21814  7.880928 Estimable
## 347      0  58 2023       181.15043  8.847289 Estimable
## 348      0  58 2024       189.10336  9.407862 Estimable
## 349      0  59 2018        76.33577  6.596646 Estimable
## 350      0  59 2020       136.42750  7.106743 Estimable
## 351      0  59 2021       148.80550  7.256008 Estimable
## 352      0  59 2022       161.58363  7.880928 Estimable
## 353      0  59 2023       185.18440  8.847289 Estimable
## 354      0  59 2024       188.82378  9.407862 Estimable
## 355      0  60 2018        82.28388  7.704892 Estimable
## 356      0  60 2020       135.17223  8.111240 Estimable
## 357      0  60 2021       145.73927  8.459606 Estimable
## 358      0  60 2022       164.25772  9.010674 Estimable
## 359      0  60 2023       182.78414  9.786905 Estimable
## 360      0  60 2024       189.44641 10.013885 Estimable
## 361      0  61 2018        80.91963  6.596646 Estimable
## 362      0  61 2020       134.96275  7.106743 Estimable
## 363      0  61 2021       140.70176  7.256008 Estimable
## 364      0  61 2022       165.29905  7.880928 Estimable
## 365      0  61 2023       183.98144  8.847289 Estimable
## 366      0  61 2024       188.52793  9.407862 Estimable
## 367      0  62 2018        79.59319  6.596646 Estimable
## 368      0  62 2020       135.87418  7.106743 Estimable
## 369      0  62 2021       141.94108  7.256008 Estimable
## 370      0  62 2022       163.23079  7.880928 Estimable
## 371      0  62 2023       186.77097  8.847289 Estimable
## 372      0  62 2024       189.25133  9.407862 Estimable
## 373      0  63 2018        80.75495  6.596646 Estimable
## 374      0  63 2020       134.51533  7.106743 Estimable
## 375      0  63 2021       146.42336  7.256008 Estimable
## 376      0  63 2022       161.98369  7.880928 Estimable
## 377      0  63 2023       184.70260  8.847289 Estimable
## 378      0  63 2024       190.03822  9.407862 Estimable
## 379      0  64 2018        79.44371  6.596646 Estimable
## 380      0  64 2020       136.82588  7.106743 Estimable
## 381      0  64 2021       147.16200  7.256008 Estimable
## 382      0  64 2022       165.20159  7.880928 Estimable
## 383      0  64 2023       180.99856  8.847289 Estimable
## 384      0  64 2024       189.05583  9.407862 Estimable
## 385      0  65 2018        81.46162  6.596646 Estimable
## 386      0  65 2020       135.86747  7.106743 Estimable
## 387      0  65 2021       145.17164  7.256008 Estimable
## 388      0  65 2022       167.71468  7.880928 Estimable
## 389      0  65 2023       185.62742  8.847289 Estimable
## 390      0  65 2024       187.58203  9.407862 Estimable
## 391      0  66 2018        79.03569  6.596652 Estimable
## 392      0  66 2020       138.38825  7.107788 Estimable
## 393      0  66 2021       144.48734  7.285771 Estimable
## 394      0  66 2022       159.52238  8.407083 Estimable
## 395      0  66 2023       182.78414  9.786905 Estimable
## 396      0  66 2024       189.44641 10.013885 Estimable
## 397      0  67 2018        82.28388  7.704892 Estimable
## 398      0  67 2020       135.17223  8.111240 Estimable
## 399      0  67 2021       145.73927  8.459606 Estimable
## 400      0  67 2022       164.25772  9.010674 Estimable
## 401      0  67 2023       182.78414  9.786905 Estimable
## 402      0  67 2024       189.44641 10.013885 Estimable
## 403      0  68 2018        79.24620  6.596646 Estimable
## 404      0  68 2020       133.98513  7.106743 Estimable
## 405      0  68 2021       145.17490  7.256008 Estimable
## 406      0  68 2022       158.70486  7.880928 Estimable
## 407      0  68 2023       189.75573  8.847289 Estimable
## 408      0  68 2024       188.71041  9.407862 Estimable
## 409      0  69 2018        80.94377  6.596646 Estimable
## 410      0  69 2020       135.91156  7.106743 Estimable
## 411      0  69 2021       143.47815  7.256008 Estimable
## 412      0  69 2022       165.24453  7.880928 Estimable
## 413      0  69 2023       185.10877  8.847289 Estimable
## 414      0  69 2024       187.47324  9.407862 Estimable
## 415      0  70 2018        84.30818  6.596646 Estimable
## 416      0  70 2020       132.30434  7.106743 Estimable
## 417      0  70 2021       144.87198  7.256008 Estimable
## 418      0  70 2022       163.46338  7.880928 Estimable
## 419      0  70 2023       184.73317  8.847289 Estimable
## 420      0  70 2024       188.71788  9.407862 Estimable
## 421      0  71 2018        83.22136  6.596646 Estimable
## 422      0  71 2020       136.15538  7.106743 Estimable
## 423      0  71 2021       146.19245  7.256008 Estimable
## 424      0  71 2022       156.94299  7.880928 Estimable
## 425      0  71 2023       190.23533  8.847289 Estimable
## 426      0  71 2024       188.82900  9.407862 Estimable
## 427      0  72 2018        83.90483  6.596646 Estimable
## 428      0  72 2020       137.84053  7.106743 Estimable
## 429      0  72 2021       142.81727  7.256008 Estimable
## 430      0  72 2022       167.71428  7.880928 Estimable
## 431      0  72 2023       184.62311  8.847289 Estimable
## 432      0  72 2024       187.32598  9.407862 Estimable
## 433      0  73 2018        82.28388  7.704892 Estimable
## 434      0  73 2020       135.17223  8.111240 Estimable
## 435      0  73 2021       145.73927  8.459606 Estimable
## 436      0  73 2022       164.25772  9.010674 Estimable
## 437      0  73 2023       182.78414  9.786905 Estimable
## 438      0  73 2024       189.44641 10.013885 Estimable
## 439      0  74 2018        99.78662  6.603497 Estimable
## 440      0  74 2020       124.48978  7.586290 Estimable
## 441      0  74 2021       145.73927  8.459606 Estimable
## 442      0  74 2022       164.25772  9.010674 Estimable
## 443      0  74 2023       182.78414  9.786905 Estimable
## 444      0  74 2024       189.44641 10.013885 Estimable
## 445      0  75 2018        82.28388  7.704892 Estimable
## 446      0  75 2020       135.17223  8.111240 Estimable
## 447      0  75 2021       145.73927  8.459606 Estimable
## 448      0  75 2022       164.25772  9.010674 Estimable
## 449      0  75 2023       182.78414  9.786905 Estimable
## 450      0  75 2024       189.44641 10.013885 Estimable
## 451      0  76 2018        80.31299  6.596646 Estimable
## 452      0  76 2020       141.52183  7.106743 Estimable
## 453      0  76 2021       147.12267  7.256008 Estimable
## 454      0  76 2022       160.39865  7.880928 Estimable
## 455      0  76 2023       183.13232  8.847289 Estimable
## 456      0  76 2024       189.62011  9.407862 Estimable
## 457      0  77 2018        84.72979  6.596646 Estimable
## 458      0  77 2020       134.49448  7.106743 Estimable
## 459      0  77 2021       147.19374  7.256008 Estimable
## 460      0  77 2022       165.72970  7.880928 Estimable
## 461      0  77 2023       182.11517  8.847289 Estimable
## 462      0  77 2024       188.01179  9.407862 Estimable
## 463      0  78 2018        86.85091  6.596646 Estimable
## 464      0  78 2020       132.25168  7.106743 Estimable
## 465      0  78 2021       145.03800  7.256008 Estimable
## 466      0  78 2022       163.72975  7.880928 Estimable
## 467      0  78 2023       185.14545  8.847289 Estimable
## 468      0  78 2024       187.48985  9.407862 Estimable
## 469      0  79 2018        86.83938  6.596646 Estimable
## 470      0  79 2020       138.42319  7.106743 Estimable
## 471      0  79 2021       143.62140  7.256008 Estimable
## 472      0  79 2022       163.62199  7.880928 Estimable
## 473      0  79 2023       179.37429  8.847289 Estimable
## 474      0  79 2024       189.98052  9.407862 Estimable
## 475      0  80 2018        85.67044  6.596646 Estimable
## 476      0  80 2020       137.63335  7.106743 Estimable
## 477      0  80 2021       150.30896  7.256008 Estimable
## 478      0  80 2022       160.55821  7.880928 Estimable
## 479      0  80 2023       192.19586  8.847289 Estimable
## 480      0  80 2024       177.32392  9.407862 Estimable
## 481      0  81 2018        77.80202  6.596646 Estimable
## 482      0  81 2020       132.94072  7.106743 Estimable
## 483      0  81 2021       144.20490  7.256008 Estimable
## 484      0  81 2022       164.31769  7.880928 Estimable
## 485      0  81 2023       186.21284  8.847289 Estimable
## 486      0  81 2024       190.43007  9.407862 Estimable
## 487      0  82 2018        84.18263  6.596646 Estimable
## 488      0  82 2020       132.56180  7.106743 Estimable
## 489      0  82 2021       148.43541  7.256008 Estimable
## 490      0  82 2022       165.96953  7.880928 Estimable
## 491      0  82 2023       184.10290  8.847289 Estimable
## 492      0  82 2024       185.86456  9.407862 Estimable
## 493      0  83 2018        79.04243  6.596652 Estimable
## 494      0  83 2020       135.08910  7.107788 Estimable
## 495      0  83 2021       147.70906  7.285771 Estimable
## 496      0  83 2022       161.62234  8.407083 Estimable
## 497      0  83 2023       182.78414  9.786905 Estimable
## 498      0  83 2024       189.44641 10.013885 Estimable
## 499      0  84 2018        83.83204  6.596646 Estimable
## 500      0  84 2020       135.73506  7.106743 Estimable
## 501      0  84 2021       145.93226  7.256008 Estimable
## 502      0  84 2022       161.09974  7.880928 Estimable
## 503      0  84 2023       186.52084  8.847289 Estimable
## 504      0  84 2024       187.84888  9.407862 Estimable
## 505      0  85 2018        77.49385  6.596646 Estimable
## 506      0  85 2020       135.66848  7.106743 Estimable
## 507      0  85 2021       151.00233  7.256008 Estimable
## 508      0  85 2022       156.52492  7.880928 Estimable
## 509      0  85 2023       187.69701  8.847289 Estimable
## 510      0  85 2024       188.16458  9.407862 Estimable
## 511      0  86 2018        85.47115  6.596646 Estimable
## 512      0  86 2020       139.60633  7.106743 Estimable
## 513      0  86 2021       148.59783  7.256008 Estimable
## 514      0  86 2022       162.71591  7.880928 Estimable
## 515      0  86 2023       183.54509  8.847289 Estimable
## 516      0  86 2024       187.04677  9.407862 Estimable
## 517      0  87 2018        89.92892  6.596646 Estimable
## 518      0  87 2020       138.67658  7.106743 Estimable
## 519      0  87 2021       147.76599  7.256008 Estimable
## 520      0  87 2022       166.92838  7.880928 Estimable
## 521      0  87 2023       184.50648  8.847289 Estimable
## 522      0  87 2024       183.28199  9.407862 Estimable
## 523      0  88 2018        82.28388  7.704892 Estimable
## 524      0  88 2020       135.17223  8.111240 Estimable
## 525      0  88 2021       145.73927  8.459606 Estimable
## 526      0  88 2022       164.25772  9.010674 Estimable
## 527      0  88 2023       182.78414  9.786905 Estimable
## 528      0  88 2024       189.44641 10.013885 Estimable
## 529      0  89 2018        81.53217  6.596646 Estimable
## 530      0  89 2020       135.92656  7.106743 Estimable
## 531      0  89 2021       145.04987  7.256008 Estimable
## 532      0  89 2022       165.25113  7.880928 Estimable
## 533      0  89 2023       176.42245  8.847289 Estimable
## 534      0  89 2024       191.87573  9.407862 Estimable
## 535      0  90 2018        82.28388  7.704892 Estimable
## 536      0  90 2020       135.17223  8.111240 Estimable
## 537      0  90 2021       145.73927  8.459606 Estimable
## 538      0  90 2022       164.25772  9.010674 Estimable
## 539      0  90 2023       182.78414  9.786905 Estimable
## 540      0  90 2024       189.44641 10.013885 Estimable
## 541      0  91 2018        81.37229  6.596646 Estimable
## 542      0  91 2020       140.06869  7.106743 Estimable
## 543      0  91 2021       144.62523  7.256008 Estimable
## 544      0  91 2022       155.72055  7.880928 Estimable
## 545      0  91 2023       172.18508  8.847289 Estimable
## 546      0  91 2024       205.45381  9.407862 Estimable
## 547      0  92 2018        79.13164  6.596646 Estimable
## 548      0  92 2020       136.69861  7.106743 Estimable
## 549      0  92 2021       152.65379  7.256008 Estimable
## 550      0  92 2022       161.04794  7.880928 Estimable
## 551      0  92 2023       182.90168  8.847289 Estimable
## 552      0  92 2024       185.56073  9.407862 Estimable
## 553      0  93 2018        89.50105  6.596646 Estimable
## 554      0  93 2020       139.34587  7.106743 Estimable
## 555      0  93 2021       140.17992  7.256008 Estimable
## 556      0  93 2022       173.32595  7.880928 Estimable
## 557      0  93 2023       174.76083  8.847289 Estimable
## 558      0  93 2024       192.72785  9.407862 Estimable
## 559      0  94 2018        82.28388  7.704892 Estimable
## 560      0  94 2020       135.17223  8.111240 Estimable
## 561      0  94 2021       145.73927  8.459606 Estimable
## 562      0  94 2022       164.25772  9.010674 Estimable
## 563      0  94 2023       182.78414  9.786905 Estimable
## 564      0  94 2024       189.44641 10.013885 Estimable
## 565      0  95 2018        87.55719  6.596652 Estimable
## 566      0  95 2020       139.87381  7.107788 Estimable
## 567      0  95 2021       137.89803  7.285771 Estimable
## 568      0  95 2022       168.11643  8.407083 Estimable
## 569      0  95 2023       182.78414  9.786905 Estimable
## 570      0  95 2024       189.44641 10.013885 Estimable
## 571      0  96 2018        82.28388  7.704892 Estimable
## 572      0  96 2020       135.17223  8.111240 Estimable
## 573      0  96 2021       145.73927  8.459606 Estimable
## 574      0  96 2022       164.25772  9.010674 Estimable
## 575      0  96 2023       182.78414  9.786905 Estimable
## 576      0  96 2024       189.44641 10.013885 Estimable
## 577      0  97 2018        78.79602  6.596646 Estimable
## 578      0  97 2020       137.01310  7.106743 Estimable
## 579      0  97 2021       144.57797  7.256008 Estimable
## 580      0  97 2022       166.26809  7.880928 Estimable
## 581      0  97 2023       185.03622  8.847289 Estimable
## 582      0  97 2024       188.77049  9.407862 Estimable
## 583      0  98 2018        83.56533  6.596646 Estimable
## 584      0  98 2020       131.19715  7.106743 Estimable
## 585      0  98 2021       152.87674  7.256008 Estimable
## 586      0  98 2022       166.03388  7.880928 Estimable
## 587      0  98 2023       182.67269  8.847289 Estimable
## 588      0  98 2024       186.80922  9.407862 Estimable
## 589      0  99 2018        77.05842  6.596646 Estimable
## 590      0  99 2020       133.72378  7.106743 Estimable
## 591      0  99 2021       145.84027  7.256008 Estimable
## 592      0  99 2022       162.61833  7.880928 Estimable
## 593      0  99 2023       183.12127  8.847289 Estimable
## 594      0  99 2024       192.29453  9.407862 Estimable
## 595      0 100 2018        82.28388  7.704892 Estimable
## 596      0 100 2020       135.17223  8.111240 Estimable
## 597      0 100 2021       145.73927  8.459606 Estimable
## 598      0 100 2022       164.25772  9.010674 Estimable
## 599      0 100 2023       182.78414  9.786905 Estimable
## 600      0 100 2024       189.44641 10.013885 Estimable
## 601      0 101 2018        82.28388  7.704892 Estimable
## 602      0 101 2020       135.17223  8.111240 Estimable
## 603      0 101 2021       145.73927  8.459606 Estimable
## 604      0 101 2022       164.25772  9.010674 Estimable
## 605      0 101 2023       182.78414  9.786905 Estimable
## 606      0 101 2024       189.44641 10.013885 Estimable
## 607      0 102 2018        82.28388  7.704892 Estimable
## 608      0 102 2020       135.17223  8.111240 Estimable
## 609      0 102 2021       145.73927  8.459606 Estimable
## 610      0 102 2022       164.25772  9.010674 Estimable
## 611      0 102 2023       182.78414  9.786905 Estimable
## 612      0 102 2024       189.44641 10.013885 Estimable
## 613      0 103 2018        82.28388  7.704892 Estimable
## 614      0 103 2020       135.17223  8.111240 Estimable
## 615      0 103 2021       145.73927  8.459606 Estimable
## 616      0 103 2022       164.25772  9.010674 Estimable
## 617      0 103 2023       182.78414  9.786905 Estimable
## 618      0 103 2024       189.44641 10.013885 Estimable
## 619      0 104 2018        82.28388  7.704892 Estimable
## 620      0 104 2020       135.17223  8.111240 Estimable
## 621      0 104 2021       145.73927  8.459606 Estimable
## 622      0 104 2022       164.25772  9.010674 Estimable
## 623      0 104 2023       182.78414  9.786905 Estimable
## 624      0 104 2024       189.44641 10.013885 Estimable
## 625      0 105 2018        82.28388  7.704892 Estimable
## 626      0 105 2020       135.17223  8.111240 Estimable
## 627      0 105 2021       145.73927  8.459606 Estimable
## 628      0 105 2022       164.25772  9.010674 Estimable
## 629      0 105 2023       182.78414  9.786905 Estimable
## 630      0 105 2024       189.44641 10.013885 Estimable
## 631      0 106 2018        81.37356  6.598485 Estimable
## 632      0 106 2020       135.54057  7.107270 Estimable
## 633      0 106 2021       150.29637  7.256399 Estimable
## 634      0 106 2022       175.00150  7.881017 Estimable
## 635      0 106 2023       166.10027  8.847059 Estimable
## 636      0 106 2024       194.45024  9.407803 Estimable
## 637      0 107 2018        82.11739  6.598485 Estimable
## 638      0 107 2020       135.39180  7.107270 Estimable
## 639      0 107 2021       151.56045  7.256399 Estimable
## 640      0 107 2022       162.29680  7.881017 Estimable
## 641      0 107 2023       177.83490  8.847059 Estimable
## 642      0 107 2024       192.22832  9.407803 Estimable
## 643      0 108 2018        83.81570  6.598485 Estimable
## 644      0 108 2020       135.05122  7.107270 Estimable
## 645      0 108 2021       145.58574  7.256399 Estimable
## 646      0 108 2022       167.22790  7.881017 Estimable
## 647      0 108 2023       182.26486  8.847059 Estimable
## 648      0 108 2024       188.04741  9.407803 Estimable
## 649      0 109 2018        82.89634  6.598485 Estimable
## 650      0 109 2020       134.55362  7.107270 Estimable
## 651      0 109 2021       140.78647  7.256399 Estimable
## 652      0 109 2022       170.85371  7.881017 Estimable
## 653      0 109 2023       181.68373  8.847059 Estimable
## 654      0 109 2024       189.22843  9.407803 Estimable
## 655      0 110 2018        90.14870  6.605338 Estimable
## 656      0 110 2020       134.34542  7.586823 Estimable
## 657      0 110 2021       145.73927  8.459606 Estimable
## 658      0 110 2022       164.25772  9.010674 Estimable
## 659      0 110 2023       182.78414  9.786905 Estimable
## 660      0 110 2024       189.44641 10.013885 Estimable
## 661      0 111 2018        81.01335  6.598485 Estimable
## 662      0 111 2020       132.44043  7.107270 Estimable
## 663      0 111 2021       147.28669  7.256399 Estimable
## 664      0 111 2022       162.52940  7.881017 Estimable
## 665      0 111 2023       186.58433  8.847059 Estimable
## 666      0 111 2024       187.86306  9.407803 Estimable
## 667      0 112 2018        87.19510  6.598485 Estimable
## 668      0 112 2020       136.49134  7.107270 Estimable
## 669      0 112 2021       147.13246  7.256399 Estimable
## 670      0 112 2022       161.39383  7.881017 Estimable
## 671      0 112 2023       182.02817  8.847059 Estimable
## 672      0 112 2024       190.66633  9.407803 Estimable
## 673      0 113 2018        82.28388  7.704892 Estimable
## 674      0 113 2020       135.17223  8.111240 Estimable
## 675      0 113 2021       145.73927  8.459606 Estimable
## 676      0 113 2022       164.25772  9.010674 Estimable
## 677      0 113 2023       182.78414  9.786905 Estimable
## 678      0 113 2024       189.44641 10.013885 Estimable
## 679      0 114 2018        86.53235  6.598485 Estimable
## 680      0 114 2020       130.62169  7.107270 Estimable
## 681      0 114 2021       145.85704  7.256399 Estimable
## 682      0 114 2022       165.83229  7.881017 Estimable
## 683      0 114 2023       184.63029  8.847059 Estimable
## 684      0 114 2024       187.34533  9.407803 Estimable
## 685      0 115 2018        82.90338  6.598485 Estimable
## 686      0 115 2020       136.43643  7.107270 Estimable
## 687      0 115 2021       148.53668  7.256399 Estimable
## 688      0 115 2022       160.07222  7.881017 Estimable
## 689      0 115 2023       187.05655  8.847059 Estimable
## 690      0 115 2024       186.65112  9.407803 Estimable
## 691      0 116 2018        85.80781  6.598485 Estimable
## 692      0 116 2020       129.83800  7.107270 Estimable
## 693      0 116 2021       142.06147  7.256399 Estimable
## 694      0 116 2022       170.60801  7.881017 Estimable
## 695      0 116 2023       173.74138  8.847059 Estimable
## 696      0 116 2024       193.83737  9.407803 Estimable
## 697      0 117 2018        82.28388  7.704892 Estimable
## 698      0 117 2020       135.17223  8.111240 Estimable
## 699      0 117 2021       145.73927  8.459606 Estimable
## 700      0 117 2022       164.25772  9.010674 Estimable
## 701      0 117 2023       182.78414  9.786905 Estimable
## 702      0 117 2024       189.44641 10.013885 Estimable
## 703      0 118 2018        82.28388  7.704892 Estimable
## 704      0 118 2020       135.17223  8.111240 Estimable
## 705      0 118 2021       145.73927  8.459606 Estimable
## 706      0 118 2022       164.25772  9.010674 Estimable
## 707      0 118 2023       182.78414  9.786905 Estimable
## 708      0 118 2024       189.44641 10.013885 Estimable
## 709      0 119 2018        83.49940  6.598485 Estimable
## 710      0 119 2020       134.43359  7.107270 Estimable
## 711      0 119 2021       145.98025  7.256399 Estimable
## 712      0 119 2022       157.41532  7.881017 Estimable
## 713      0 119 2023       187.65625  8.847059 Estimable
## 714      0 119 2024       190.82647  9.407803 Estimable
## 715      0 120 2018        82.03740  6.598485 Estimable
## 716      0 120 2020       132.86383  7.107270 Estimable
## 717      0 120 2021       143.34875  7.256399 Estimable
## 718      0 120 2022       165.84031  7.881017 Estimable
## 719      0 120 2023       184.49640  8.847059 Estimable
## 720      0 120 2024       189.97832  9.407803 Estimable
## 721      0 121 2018        81.63755  6.598485 Estimable
## 722      0 121 2020       135.50734  7.107270 Estimable
## 723      0 121 2021       139.59686  7.256399 Estimable
## 724      0 121 2022       166.55049  7.881017 Estimable
## 725      0 121 2023       184.36505  8.847059 Estimable
## 726      0 121 2024       189.94982  9.407803 Estimable
## 727      0 122 2018        81.76040  6.598485 Estimable
## 728      0 122 2020       127.49072  7.107270 Estimable
## 729      0 122 2021       143.34340  7.256399 Estimable
## 730      0 122 2022       161.88978  7.881017 Estimable
## 731      0 122 2023       188.69343  8.847059 Estimable
## 732      0 122 2024       188.44269  9.407803 Estimable
## 733      0 123 2018        79.79468  6.598485 Estimable
## 734      0 123 2020       133.92839  7.107270 Estimable
## 735      0 123 2021       145.99341  7.256399 Estimable
## 736      0 123 2022       163.47348  7.881017 Estimable
## 737      0 123 2023       183.62733  8.847059 Estimable
## 738      0 123 2024       188.42833  9.407803 Estimable
## 739      0 124 2018        79.88955  6.598485 Estimable
## 740      0 124 2020       136.13495  7.107270 Estimable
## 741      0 124 2021       145.01434  7.256399 Estimable
## 742      0 124 2022       161.71571  7.881017 Estimable
## 743      0 124 2023       186.19675  8.847059 Estimable
## 744      0 124 2024       187.76917  9.407803 Estimable
## 745      0 125 2018        86.09393  6.598485 Estimable
## 746      0 125 2020       134.93301  7.107270 Estimable
## 747      0 125 2021       147.51480  7.256399 Estimable
## 748      0 125 2022       160.63424  7.881017 Estimable
## 749      0 125 2023       186.82105  8.847059 Estimable
## 750      0 125 2024       187.92153  9.407803 Estimable
## 751      0 126 2018        81.49877  6.598485 Estimable
## 752      0 126 2020       131.73305  7.107270 Estimable
## 753      0 126 2021       141.79553  7.256399 Estimable
## 754      0 126 2022       172.87194  7.881017 Estimable
## 755      0 126 2023       181.75409  8.847059 Estimable
## 756      0 126 2024       187.91031  9.407803 Estimable
## 757      0 127 2018        82.28388  7.704892 Estimable
## 758      0 127 2020       135.17223  8.111240 Estimable
## 759      0 127 2021       145.73927  8.459606 Estimable
## 760      0 127 2022       164.25772  9.010674 Estimable
## 761      0 127 2023       182.78414  9.786905 Estimable
## 762      0 127 2024       189.44641 10.013885 Estimable
## 763      0 128 2018        79.15703  6.598485 Estimable
## 764      0 128 2020       128.74252  7.107270 Estimable
## 765      0 128 2021       145.63606  7.256399 Estimable
## 766      0 128 2022       166.66285  7.881017 Estimable
## 767      0 128 2023       182.87955  8.847059 Estimable
## 768      0 128 2024       190.89147  9.407803 Estimable
## 769      0 129 2018        82.84645  6.598485 Estimable
## 770      0 129 2020       132.29179  7.107270 Estimable
## 771      0 129 2021       147.00269  7.256399 Estimable
## 772      0 129 2022       165.06787  7.881017 Estimable
## 773      0 129 2023       182.60017  8.847059 Estimable
## 774      0 129 2024       189.47948  9.407803 Estimable
## 775      0 130 2018        78.05851  6.598485 Estimable
## 776      0 130 2020       136.12241  7.107270 Estimable
## 777      0 130 2021       145.40359  7.256399 Estimable
## 778      0 130 2022       166.39443  7.881017 Estimable
## 779      0 130 2023       182.80870  8.847059 Estimable
## 780      0 130 2024       188.19677  9.407803 Estimable
## 781      0 131 2018        82.40860  6.598485 Estimable
## 782      0 131 2020       135.51408  7.107270 Estimable
## 783      0 131 2021       143.07181  7.256399 Estimable
## 784      0 131 2022       165.83505  7.881017 Estimable
## 785      0 131 2023       183.84486  8.847059 Estimable
## 786      0 131 2024       189.80344  9.407803 Estimable
## 787      0 132 2018        79.60248  6.598485 Estimable
## 788      0 132 2020       139.23765  7.107270 Estimable
## 789      0 132 2021       142.63631  7.256399 Estimable
## 790      0 132 2022       163.69677  7.881017 Estimable
## 791      0 132 2023       182.88690  8.847059 Estimable
## 792      0 132 2024       190.89345  9.407803 Estimable
## 793      0 133 2018        80.01144  6.598485 Estimable
## 794      0 133 2020       134.21289  7.107270 Estimable
## 795      0 133 2021       144.68438  7.256399 Estimable
## 796      0 133 2022       166.08569  7.881017 Estimable
## 797      0 133 2023       180.23752  8.847059 Estimable
## 798      0 133 2024       191.53115  9.407803 Estimable
## 799      0 134 2018        81.22116  6.598485 Estimable
## 800      0 134 2020       133.37221  7.107270 Estimable
## 801      0 134 2021       147.35327  7.256399 Estimable
## 802      0 134 2022       165.83108  7.881017 Estimable
## 803      0 134 2023       182.32922  8.847059 Estimable
## 804      0 134 2024       189.40169  9.407803 Estimable
## 805      0 135 2018        74.23015  6.598485 Estimable
## 806      0 135 2020       150.99459  7.107270 Estimable
## 807      0 135 2021       139.37995  7.256399 Estimable
## 808      0 135 2022       161.41110  7.881017 Estimable
## 809      0 135 2023       185.33485  8.847059 Estimable
## 810      0 135 2024       187.53613  9.407803 Estimable
## 811      0 136 2018        82.68190  6.598485 Estimable
## 812      0 136 2020       133.56319  7.107270 Estimable
## 813      0 136 2021       146.40521  7.256399 Estimable
## 814      0 136 2022       167.45565  7.881017 Estimable
## 815      0 136 2023       182.60026  8.847059 Estimable
## 816      0 136 2024       186.80551  9.407803 Estimable
## 817      0 137 2018        86.68055  6.598485 Estimable
## 818      0 137 2020       134.29105  7.107270 Estimable
## 819      0 137 2021       136.28247  7.256399 Estimable
## 820      0 137 2022       166.26962  7.881017 Estimable
## 821      0 137 2023       184.24127  8.847059 Estimable
## 822      0 137 2024       189.93350  9.407803 Estimable
## 823      0 138 2018        82.28388  7.704892 Estimable
## 824      0 138 2020       135.17223  8.111240 Estimable
## 825      0 138 2021       145.73927  8.459606 Estimable
## 826      0 138 2022       164.25772  9.010674 Estimable
## 827      0 138 2023       182.78414  9.786905 Estimable
## 828      0 138 2024       189.44641 10.013885 Estimable
## 829      0 139 2018        89.09403  6.598485 Estimable
## 830      0 139 2020       138.19836  7.107270 Estimable
## 831      0 139 2021       137.58020  7.256399 Estimable
## 832      0 139 2022       167.28378  7.881017 Estimable
## 833      0 139 2023       183.58269  8.847059 Estimable
## 834      0 139 2024       189.73645  9.407803 Estimable
## 835      0 140 2018        83.66140  6.598485 Estimable
## 836      0 140 2020       135.59259  7.107270 Estimable
## 837      0 140 2021       146.76223  7.256399 Estimable
## 838      0 140 2022       161.60294  7.881017 Estimable
## 839      0 140 2023       178.99300  8.847059 Estimable
## 840      0 140 2024       193.87955  9.407803 Estimable
## 841      0 141 2018        82.28388  7.704892 Estimable
## 842      0 141 2020       135.17223  8.111240 Estimable
## 843      0 141 2021       145.73927  8.459606 Estimable
## 844      0 141 2022       164.25772  9.010674 Estimable
## 845      0 141 2023       182.78414  9.786905 Estimable
## 846      0 141 2024       189.44641 10.013885 Estimable
## 847      0 142 2018        84.27880  6.598485 Estimable
## 848      0 142 2020       129.95815  7.107270 Estimable
## 849      0 142 2021       146.39452  7.256399 Estimable
## 850      0 142 2022       164.24162  7.881017 Estimable
## 851      0 142 2023       182.50291  8.847059 Estimable
## 852      0 142 2024       190.79714  9.407803 Estimable
## 853      0 143 2018        83.78548  6.598485 Estimable
## 854      0 143 2020       135.41627  7.107270 Estimable
## 855      0 143 2021       146.04065  7.256399 Estimable
## 856      0 143 2022       158.68451  7.881017 Estimable
## 857      0 143 2023       185.40959  8.847059 Estimable
## 858      0 143 2024       190.23695  9.407803 Estimable
## 859      0 144 2018        82.08509  6.598485 Estimable
## 860      0 144 2020       138.28559  7.107270 Estimable
## 861      0 144 2021       137.13394  7.256399 Estimable
## 862      0 144 2022       169.41874  7.881017 Estimable
## 863      0 144 2023       178.22852  8.847059 Estimable
## 864      0 144 2024       193.67098  9.407803 Estimable
## 865      0 145 2018        75.54661  6.598485 Estimable
## 866      0 145 2020       138.17423  7.107270 Estimable
## 867      0 145 2021       144.20563  7.256399 Estimable
## 868      0 145 2022       165.45040  7.881017 Estimable
## 869      0 145 2023       184.11671  8.847059 Estimable
## 870      0 145 2024       187.21255  9.407803 Estimable
## 871      0 146 2018        82.28388  7.704892 Estimable
## 872      0 146 2020       135.17223  8.111240 Estimable
## 873      0 146 2021       145.73927  8.459606 Estimable
## 874      0 146 2022       164.25772  9.010674 Estimable
## 875      0 146 2023       182.78414  9.786905 Estimable
## 876      0 146 2024       189.44641 10.013885 Estimable
## 877      0 147 2018        76.62608  6.598485 Estimable
## 878      0 147 2020       131.17302  7.107270 Estimable
## 879      0 147 2021       145.34257  7.256399 Estimable
## 880      0 147 2022       169.06297  7.881017 Estimable
## 881      0 147 2023       182.66551  8.847059 Estimable
## 882      0 147 2024       188.16002  9.407803 Estimable
## 883      0 148 2018        82.39724  6.598485 Estimable
## 884      0 148 2020       134.59296  7.107270 Estimable
## 885      0 148 2021       146.82054  7.256399 Estimable
## 886      0 148 2022       165.41560  7.881017 Estimable
## 887      0 148 2023       181.73703  8.847059 Estimable
## 888      0 148 2024       189.24950  9.407803 Estimable
## 889      0 149 2018        82.78217  6.598485 Estimable
## 890      0 149 2020       137.18520  7.107270 Estimable
## 891      0 149 2021       146.20241  7.256399 Estimable
## 892      0 149 2022       163.37156  7.881017 Estimable
## 893      0 149 2023       180.22745  8.847059 Estimable
## 894      0 149 2024       192.85531  9.407803 Estimable
## 895      0 150 2018        80.57507  6.598485 Estimable
## 896      0 150 2020       135.83828  7.107270 Estimable
## 897      0 150 2021       144.32770  7.256399 Estimable
## 898      0 150 2022       166.13941  7.881017 Estimable
## 899      0 150 2023       181.32261  8.847059 Estimable
## 900      0 150 2024       190.47695  9.407803 Estimable
## 901      0 151 2018        81.20624  6.598485 Estimable
## 902      0 151 2020       129.77855  7.107270 Estimable
## 903      0 151 2021       145.80535  7.256399 Estimable
## 904      0 151 2022       168.47305  7.881017 Estimable
## 905      0 151 2023       180.53056  8.847059 Estimable
## 906      0 151 2024       190.26942  9.407803 Estimable
## 907      0 152 2018        82.28388  7.704892 Estimable
## 908      0 152 2020       135.17223  8.111240 Estimable
## 909      0 152 2021       145.73927  8.459606 Estimable
## 910      0 152 2022       164.25772  9.010674 Estimable
## 911      0 152 2023       182.78414  9.786905 Estimable
## 912      0 152 2024       189.44641 10.013885 Estimable
## 913      0 153 2018        81.91779  6.598485 Estimable
## 914      0 153 2020       132.43843  7.107270 Estimable
## 915      0 153 2021       146.56263  7.256399 Estimable
## 916      0 153 2022       164.46158  7.881017 Estimable
## 917      0 153 2023       183.76460  8.847059 Estimable
## 918      0 153 2024       188.45672  9.407803 Estimable
## 919      0 154 2018        82.28388  7.704892 Estimable
## 920      0 154 2020       135.17223  8.111240 Estimable
## 921      0 154 2021       145.73927  8.459606 Estimable
## 922      0 154 2022       164.25772  9.010674 Estimable
## 923      0 154 2023       182.78414  9.786905 Estimable
## 924      0 154 2024       189.44641 10.013885 Estimable
## 925      0 155 2018        82.28388  7.704892 Estimable
## 926      0 155 2020       135.17223  8.111240 Estimable
## 927      0 155 2021       145.73927  8.459606 Estimable
## 928      0 155 2022       164.25772  9.010674 Estimable
## 929      0 155 2023       182.78414  9.786905 Estimable
## 930      0 155 2024       189.44641 10.013885 Estimable
## 931      0 156 2018        82.28388  7.704892 Estimable
## 932      0 156 2020       135.17223  8.111240 Estimable
## 933      0 156 2021       145.73927  8.459606 Estimable
## 934      0 156 2022       164.25772  9.010674 Estimable
## 935      0 156 2023       182.78414  9.786905 Estimable
## 936      0 156 2024       189.44641 10.013885 Estimable
## 937      0 157 2018        82.28388  7.704892 Estimable
## 938      0 157 2020       135.17223  8.111240 Estimable
## 939      0 157 2021       145.73927  8.459606 Estimable
## 940      0 157 2022       164.25772  9.010674 Estimable
## 941      0 157 2023       182.78414  9.786905 Estimable
## 942      0 157 2024       189.44641 10.013885 Estimable
## 943      0 158 2018        82.28388  7.704892 Estimable
## 944      0 158 2020       135.17223  8.111240 Estimable
## 945      0 158 2021       145.73927  8.459606 Estimable
## 946      0 158 2022       164.25772  9.010674 Estimable
## 947      0 158 2023       182.78414  9.786905 Estimable
## 948      0 158 2024       189.44641 10.013885 Estimable
## 949      0 159 2018        82.28388  7.704892 Estimable
## 950      0 159 2020       135.17223  8.111240 Estimable
## 951      0 159 2021       145.73927  8.459606 Estimable
## 952      0 159 2022       164.25772  9.010674 Estimable
## 953      0 159 2023       182.78414  9.786905 Estimable
## 954      0 159 2024       189.44641 10.013885 Estimable
## 955      0 160 2018        84.68329  6.602888 Estimable
## 956      0 160 2020       144.18680  7.108154 Estimable
## 957      0 160 2021       131.40303  7.257038 Estimable
## 958      0 160 2022       170.81395  7.881157 Estimable
## 959      0 160 2023       179.01615  8.846729 Estimable
## 960      0 160 2024       191.21248  9.407713 Estimable
## 961      0 161 2018        88.70166  6.602888 Estimable
## 962      0 161 2020       147.36256  7.108154 Estimable
## 963      0 161 2021       127.52879  7.257038 Estimable
## 964      0 161 2022       166.26856  7.881157 Estimable
## 965      0 161 2023       175.76225  8.846729 Estimable
## 966      0 161 2024       195.72427  9.407713 Estimable
## 967      0 162 2018        86.06257  6.602888 Estimable
## 968      0 162 2020       143.20439  7.108154 Estimable
## 969      0 162 2021       132.69302  7.257038 Estimable
## 970      0 162 2022       172.37508  7.881157 Estimable
## 971      0 162 2023       177.11809  8.846729 Estimable
## 972      0 162 2024       190.70978  9.407713 Estimable
## 973      0 163 2018        87.49173  6.602888 Estimable
## 974      0 163 2020       137.94549  7.108154 Estimable
## 975      0 163 2021       145.56196  7.257038 Estimable
## 976      0 163 2022       161.89578  7.881157 Estimable
## 977      0 163 2023       181.76689  8.846729 Estimable
## 978      0 163 2024       191.92715  9.407713 Estimable
## 979      0 164 2018        82.62911  6.602888 Estimable
## 980      0 164 2020       141.20985  7.108154 Estimable
## 981      0 164 2021       147.33125  7.257038 Estimable
## 982      0 164 2022       163.23596  7.881157 Estimable
## 983      0 164 2023       178.09637  8.846729 Estimable
## 984      0 164 2024       193.61932  9.407713 Estimable
## 985      0 165 2018        81.43594  6.602888 Estimable
## 986      0 165 2020       143.64969  7.108154 Estimable
## 987      0 165 2021       139.04784  7.257038 Estimable
## 988      0 165 2022       162.96417  7.881157 Estimable
## 989      0 165 2023       184.65878  8.846729 Estimable
## 990      0 165 2024       190.02600  9.407713 Estimable
## 991      0 166 2018        84.32844  6.602888 Estimable
## 992      0 166 2020       137.96352  7.108154 Estimable
## 993      0 166 2021       144.37232  7.257038 Estimable
## 994      0 166 2022       173.96403  7.881157 Estimable
## 995      0 166 2023       179.00932  8.846729 Estimable
## 996      0 166 2024       185.83391  9.407713 Estimable
## 997      0 167 2018        84.21389  6.602888 Estimable
## 998      0 167 2020       129.12109  7.108154 Estimable
## 999      0 167 2021       147.46368  7.257038 Estimable
## 1000     0 167 2022       165.18585  7.881157 Estimable
## 1001     0 167 2023       184.33965  8.846729 Estimable
## 1002     0 167 2024       188.60334  9.407713 Estimable
## 1003     0 168 2018        82.28388  7.704892 Estimable
## 1004     0 168 2020       135.17223  8.111240 Estimable
## 1005     0 168 2021       145.73927  8.459606 Estimable
## 1006     0 168 2022       164.25772  9.010674 Estimable
## 1007     0 168 2023       182.78414  9.786905 Estimable
## 1008     0 168 2024       189.44641 10.013885 Estimable
## 1009     0 169 2018        75.21411  6.602888 Estimable
## 1010     0 169 2020       135.24444  7.108154 Estimable
## 1011     0 169 2021       147.48311  7.257038 Estimable
## 1012     0 169 2022       169.67138  7.881157 Estimable
## 1013     0 169 2023       180.97612  8.846729 Estimable
## 1014     0 169 2024       187.70052  9.407713 Estimable
## 1015     0 170 2018        82.70363  6.602888 Estimable
## 1016     0 170 2020       136.31286  7.108154 Estimable
## 1017     0 170 2021       143.70581  7.257038 Estimable
## 1018     0 170 2022       168.76040  7.881157 Estimable
## 1019     0 170 2023       181.71996  8.846729 Estimable
## 1020     0 170 2024       187.90356  9.407713 Estimable
## 1021     0 171 2018        82.28388  7.704892 Estimable
## 1022     0 171 2020       135.17223  8.111240 Estimable
## 1023     0 171 2021       145.73927  8.459606 Estimable
## 1024     0 171 2022       164.25772  9.010674 Estimable
## 1025     0 171 2023       182.78414  9.786905 Estimable
## 1026     0 171 2024       189.44641 10.013885 Estimable
## 1027     0 172 2018        81.10312  6.602888 Estimable
## 1028     0 172 2020       136.43302  7.108154 Estimable
## 1029     0 172 2021       152.53518  7.257038 Estimable
## 1030     0 172 2022       159.44562  7.881157 Estimable
## 1031     0 172 2023       183.38822  8.846729 Estimable
## 1032     0 172 2024       189.68158  9.407713 Estimable
## 1033     0 173 2018        80.65184  6.602888 Estimable
## 1034     0 173 2020       134.43346  7.108154 Estimable
## 1035     0 173 2021       145.70230  7.257038 Estimable
## 1036     0 173 2022       165.95070  7.881157 Estimable
## 1037     0 173 2023       182.20255  8.846729 Estimable
## 1038     0 173 2024       189.37347  9.407713 Estimable
## 1039     0 174 2018        82.28388  7.704892 Estimable
## 1040     0 174 2020       135.17223  8.111240 Estimable
## 1041     0 174 2021       145.73927  8.459606 Estimable
## 1042     0 174 2022       164.25772  9.010674 Estimable
## 1043     0 174 2023       182.78414  9.786905 Estimable
## 1044     0 174 2024       189.44641 10.013885 Estimable
## 1045     0 175 2018        79.28286  6.602888 Estimable
## 1046     0 175 2020       136.60106  7.108154 Estimable
## 1047     0 175 2021       146.51685  7.257038 Estimable
## 1048     0 175 2022       163.84961  7.881157 Estimable
## 1049     0 175 2023       182.32821  8.846729 Estimable
## 1050     0 175 2024       190.73913  9.407713 Estimable
## 1051     0 176 2018        81.03462  6.602888 Estimable
## 1052     0 176 2020       134.35665  7.108154 Estimable
## 1053     0 176 2021       145.02336  7.257038 Estimable
## 1054     0 176 2022       164.57373  7.881157 Estimable
## 1055     0 176 2023       182.19849  8.846729 Estimable
## 1056     0 176 2024       188.05566  9.407713 Estimable
## 1057     0 177 2018        81.86021  6.602888 Estimable
## 1058     0 177 2020       137.67505  7.108154 Estimable
## 1059     0 177 2021       145.03690  7.257038 Estimable
## 1060     0 177 2022       164.15523  7.881157 Estimable
## 1061     0 177 2023       183.68472  8.846729 Estimable
## 1062     0 177 2024       189.75779  9.407713 Estimable
## 1063     0 178 2018        78.17766  6.602888 Estimable
## 1064     0 178 2020       136.45120  7.108154 Estimable
## 1065     0 178 2021       140.47858  7.257038 Estimable
## 1066     0 178 2022       166.95150  7.881157 Estimable
## 1067     0 178 2023       184.50635  8.846729 Estimable
## 1068     0 178 2024       187.32292  9.407713 Estimable
## 1069     0 179 2018        76.88393  6.602888 Estimable
## 1070     0 179 2020       136.73722  7.108154 Estimable
## 1071     0 179 2021       144.64673  7.257038 Estimable
## 1072     0 179 2022       167.86534  7.881157 Estimable
## 1073     0 179 2023       180.09582  8.846729 Estimable
## 1074     0 179 2024       188.82152  9.407713 Estimable
## 1075     0 180 2018        82.28388  7.704892 Estimable
## 1076     0 180 2020       135.17223  8.111240 Estimable
## 1077     0 180 2021       145.73927  8.459606 Estimable
## 1078     0 180 2022       164.25772  9.010674 Estimable
## 1079     0 180 2023       182.78414  9.786905 Estimable
## 1080     0 180 2024       189.44641 10.013885 Estimable
## 1081     0 181 2018        83.91990  6.602888 Estimable
## 1082     0 181 2020       133.15364  7.108154 Estimable
## 1083     0 181 2021       143.63082  7.257038 Estimable
## 1084     0 181 2022       161.82082  7.881157 Estimable
## 1085     0 181 2023       184.56481  8.846729 Estimable
## 1086     0 181 2024       190.01768  9.407713 Estimable
## 1087     0 182 2018        82.28388  7.704892 Estimable
## 1088     0 182 2020       135.17223  8.111240 Estimable
## 1089     0 182 2021       145.73927  8.459606 Estimable
## 1090     0 182 2022       164.25772  9.010674 Estimable
## 1091     0 182 2023       182.78414  9.786905 Estimable
## 1092     0 182 2024       189.44641 10.013885 Estimable
## 1093     0 183 2018        82.39596  6.602888 Estimable
## 1094     0 183 2020       138.67741  7.108154 Estimable
## 1095     0 183 2021       141.71160  7.257038 Estimable
## 1096     0 183 2022       159.91026  7.881157 Estimable
## 1097     0 183 2023       188.31792  8.846729 Estimable
## 1098     0 183 2024       187.00223  9.407713 Estimable
## 1099     0 184 2018        82.28388  7.704892 Estimable
## 1100     0 184 2020       135.17223  8.111240 Estimable
## 1101     0 184 2021       145.73927  8.459606 Estimable
## 1102     0 184 2022       164.25772  9.010674 Estimable
## 1103     0 184 2023       182.78414  9.786905 Estimable
## 1104     0 184 2024       189.44641 10.013885 Estimable
## 1105     0 185 2018        82.28388  7.704892 Estimable
## 1106     0 185 2020       135.17223  8.111240 Estimable
## 1107     0 185 2021       145.73927  8.459606 Estimable
## 1108     0 185 2022       164.25772  9.010674 Estimable
## 1109     0 185 2023       182.78414  9.786905 Estimable
## 1110     0 185 2024       189.44641 10.013885 Estimable
## 1111     0 186 2018        80.25914  6.602888 Estimable
## 1112     0 186 2020       134.04610  7.108154 Estimable
## 1113     0 186 2021       144.26496  7.257038 Estimable
## 1114     0 186 2022       165.81836  7.881157 Estimable
## 1115     0 186 2023       181.10322  8.846729 Estimable
## 1116     0 186 2024       190.42892  9.407713 Estimable
## 1117     0 187 2018        82.28388  7.704892 Estimable
## 1118     0 187 2020       135.17223  8.111240 Estimable
## 1119     0 187 2021       145.73927  8.459606 Estimable
## 1120     0 187 2022       164.25772  9.010674 Estimable
## 1121     0 187 2023       182.78414  9.786905 Estimable
## 1122     0 187 2024       189.44641 10.013885 Estimable
## 1123     0 188 2018        83.08489  6.602888 Estimable
## 1124     0 188 2020       133.76964  7.108154 Estimable
## 1125     0 188 2021       143.64916  7.257038 Estimable
## 1126     0 188 2022       169.66615  7.881157 Estimable
## 1127     0 188 2023       180.69530  8.846729 Estimable
## 1128     0 188 2024       187.64205  9.407713 Estimable
## 1129     0 189 2018        82.28388  7.704892 Estimable
## 1130     0 189 2020       135.17223  8.111240 Estimable
## 1131     0 189 2021       145.73927  8.459606 Estimable
## 1132     0 189 2022       164.25772  9.010674 Estimable
## 1133     0 189 2023       182.78414  9.786905 Estimable
## 1134     0 189 2024       189.44641 10.013885 Estimable
## 1135     0 190 2018        82.28388  7.704892 Estimable
## 1136     0 190 2020       135.17223  8.111240 Estimable
## 1137     0 190 2021       145.73927  8.459606 Estimable
## 1138     0 190 2022       164.25772  9.010674 Estimable
## 1139     0 190 2023       182.78414  9.786905 Estimable
## 1140     0 190 2024       189.44641 10.013885 Estimable
## 1141     0 191 2018        82.28388  7.704892 Estimable
## 1142     0 191 2020       135.17223  8.111240 Estimable
## 1143     0 191 2021       145.73927  8.459606 Estimable
## 1144     0 191 2022       164.25772  9.010674 Estimable
## 1145     0 191 2023       182.78414  9.786905 Estimable
## 1146     0 191 2024       189.44641 10.013885 Estimable
## 1147     0 192 2018        83.11805  6.602888 Estimable
## 1148     0 192 2020       140.79822  7.108154 Estimable
## 1149     0 192 2021       144.08029  7.257038 Estimable
## 1150     0 192 2022       167.44247  7.881157 Estimable
## 1151     0 192 2023       176.76762  8.846729 Estimable
## 1152     0 192 2024       190.60894  9.407713 Estimable
## 1153     0 193 2018        84.86577  6.602888 Estimable
## 1154     0 193 2020       129.25339  7.108154 Estimable
## 1155     0 193 2021       151.98141  7.257038 Estimable
## 1156     0 193 2022       158.67737  7.881157 Estimable
## 1157     0 193 2023       183.31127  8.846729 Estimable
## 1158     0 193 2024       192.34507  9.407713 Estimable
## 1159     0 194 2018        82.28388  7.704892 Estimable
## 1160     0 194 2020       135.17223  8.111240 Estimable
## 1161     0 194 2021       145.73927  8.459606 Estimable
## 1162     0 194 2022       164.25772  9.010674 Estimable
## 1163     0 194 2023       182.78414  9.786905 Estimable
## 1164     0 194 2024       189.44641 10.013885 Estimable
## 1165     0 195 2018        76.34582  6.602888 Estimable
## 1166     0 195 2020       135.25216  7.108154 Estimable
## 1167     0 195 2021       149.13285  7.257038 Estimable
## 1168     0 195 2022       165.87308  7.881157 Estimable
## 1169     0 195 2023       181.45908  8.846729 Estimable
## 1170     0 195 2024       187.84029  9.407713 Estimable
## 1171     0 196 2018        81.82676  6.602888 Estimable
## 1172     0 196 2020       134.73810  7.108154 Estimable
## 1173     0 196 2021       152.32166  7.257038 Estimable
## 1174     0 196 2022       172.94923  7.881157 Estimable
## 1175     0 196 2023       174.76044  8.846729 Estimable
## 1176     0 196 2024       186.05243  9.407713 Estimable
## 1177     0 197 2018        77.49777  6.602888 Estimable
## 1178     0 197 2020       135.44777  7.108154 Estimable
## 1179     0 197 2021       141.57368  7.257038 Estimable
## 1180     0 197 2022       170.44922  7.881157 Estimable
## 1181     0 197 2023       183.16657  8.846729 Estimable
## 1182     0 197 2024       189.61195  9.407713 Estimable
## 1183     0 198 2018        81.38467  6.602888 Estimable
## 1184     0 198 2020       137.00805  7.108154 Estimable
## 1185     0 198 2021       146.56554  7.257038 Estimable
## 1186     0 198 2022       168.74888  7.881157 Estimable
## 1187     0 198 2023       178.40751  8.846729 Estimable
## 1188     0 198 2024       188.36666  9.407713 Estimable
## 1189     0 199 2018        82.28388  7.704892 Estimable
## 1190     0 199 2020       130.55384  7.190187 Estimable
## 1191     0 199 2021       144.84431  7.259587 Estimable
## 1192     0 199 2022       165.44621  7.881311 Estimable
## 1193     0 199 2023       180.52745  8.846740 Estimable
## 1194     0 199 2024       192.94668  9.407716 Estimable
## 1195     0 200 2018        83.84066  6.602888 Estimable
## 1196     0 200 2020       135.36764  7.108154 Estimable
## 1197     0 200 2021       146.90360  7.257038 Estimable
## 1198     0 200 2022       162.69429  7.881157 Estimable
## 1199     0 200 2023       184.44777  8.846729 Estimable
## 1200     0 200 2024       187.30212  9.407713 Estimable
## 1201     0 201 2018        79.87436  6.847080 Estimable
## 1202     0 201 2020       135.17223  8.111240 Estimable
## 1203     0 201 2021       145.73927  8.459606 Estimable
## 1204     0 201 2022       164.25772  9.010674 Estimable
## 1205     0 201 2023       182.78414  9.786905 Estimable
## 1206     0 201 2024       189.44641 10.013885 Estimable
## 1207     0 202 2018        84.14113  6.602888 Estimable
## 1208     0 202 2020       138.21303  7.108154 Estimable
## 1209     0 202 2021       143.98287  7.257038 Estimable
## 1210     0 202 2022       160.37713  7.881157 Estimable
## 1211     0 202 2023       179.15792  8.846729 Estimable
## 1212     0 202 2024       196.59177  9.407713 Estimable
## 1213     0 203 2018        84.72685  6.602888 Estimable
## 1214     0 203 2020       132.29845  7.108154 Estimable
## 1215     0 203 2021       147.11212  7.257038 Estimable
## 1216     0 203 2022       159.91901  7.881157 Estimable
## 1217     0 203 2023       199.79293  8.846729 Estimable
## 1218     0 203 2024       176.65651  9.407713 Estimable
## 1219     0 204 2018        82.28388  7.704892 Estimable
## 1220     0 204 2020       135.17223  8.111240 Estimable
## 1221     0 204 2021       145.73927  8.459606 Estimable
## 1222     0 204 2022       164.25772  9.010674 Estimable
## 1223     0 204 2023       182.78414  9.786905 Estimable
## 1224     0 204 2024       189.44641 10.013885 Estimable
## 1225     0 205 2018        82.28388  7.704892 Estimable
## 1226     0 205 2020       135.17223  8.111240 Estimable
## 1227     0 205 2021       145.73927  8.459606 Estimable
## 1228     0 205 2022       164.25772  9.010674 Estimable
## 1229     0 205 2023       182.78414  9.786905 Estimable
## 1230     0 205 2024       189.44641 10.013885 Estimable
## 1231     0 206 2018        84.51706  6.602890 Estimable
## 1232     0 206 2020       136.80200  7.109162 Estimable
## 1233     0 206 2021       161.36986  7.286607 Estimable
## 1234     0 206 2022       146.38943  8.406588 Estimable
## 1235     0 206 2023       182.78414  9.786905 Estimable
## 1236     0 206 2024       189.44641 10.013885 Estimable
## 1237     0 207 2018        85.79259  6.602888 Estimable
## 1238     0 207 2020       129.99630  7.108154 Estimable
## 1239     0 207 2021       152.12682  7.257038 Estimable
## 1240     0 207 2022       167.25561  7.881157 Estimable
## 1241     0 207 2023       178.58870  8.846729 Estimable
## 1242     0 207 2024       188.41360  9.407713 Estimable
## 1243     0 208 2018        82.28388  7.704892 Estimable
## 1244     0 208 2020       135.17223  8.111240 Estimable
## 1245     0 208 2021       145.73927  8.459606 Estimable
## 1246     0 208 2022       164.25772  9.010674 Estimable
## 1247     0 208 2023       182.78414  9.786905 Estimable
## 1248     0 208 2024       189.44641 10.013885 Estimable
xlsx::write.xlsx2(pred, paste0("genetic_parameters.xlsx"), sheetName = "BLUPs", append = T)

Visualização dos valores genéticos preditos

ggplot(pred, aes(x=Ano,y=predicted.value, group = Pl, colour=Pl))+
  geom_line()+
  facet_grid() +
  theme_minimal(base_size = 15) +
  theme(legend.position = 'none')

Análise de Marca Individual

Vamos agora investigar a associação entre os marcadores para altura e a pop