Dados experimentais
Os dados avaliados referem-se a um planejamento de acasalamento fatorial interpopulacional de eucalipto. A população parental feminina foi composta por 33 árvores selecionadas com base em um teste de progênie e a população parental masculina foi composta por 64 árvores selecionadas em bancos clonais, testes clonais e pomares de sementes, totalizando 286 famílias de irmãos completos. O experimento foi implantado em novembro de 2003, em delineamento de blocos ao acaso, com oito repetições e seis plantas por parcela, dispostas em fileiras, espaçadas de 3 m entre fileiras e 2 m entre plantas. Foi instalado no campo experimental da Empresa Aperam BioEnergia, Itamarandiba, Estado de Minas Gerais, Brasil (17 ° 44 ′ 45 ″ S e 42 ° 45 ′ 11 ″ W; 1000 metros a.s.l.). A caractrerística avaliada foi diâmetro na altura do peito (DAP), medido em centrímetros (cm).
Dados
head(dados,20)
## Fam dap rep
## 1 3 0.00000 1
## 2 3 0.00000 1
## 3 3 0.00000 1
## 4 3 0.00000 1
## 5 3 0.00000 1
## 6 3 0.00000 1
## 7 7 18.78028 1
## 8 7 17.50704 1
## 9 7 14.64225 1
## 10 7 0.00000 1
## 11 7 0.00000 1
## 12 7 0.00000 1
## 13 9 10.37690 1
## 14 9 16.55211 1
## 15 9 14.19662 1
## 16 9 18.52564 1
## 17 9 17.60254 1
## 18 9 0.00000 1
## 19 13 16.77493 1
## 20 13 19.60789 1
\[ Y = Xb + Za + e \]
onde \(Y\) é o vetor de dados fenotípicos, \(b\) é o vetor de efeitos fixos de blocos adicionados à média geral, \(a\) é o vetor de efeitos aditivos (assumidos como aleatórios), onde \(a ~ N (0, \sigma^2 )\), \(e\) é o vetor de resíduos (aleatório), onde \(e ~ N (0, \sigma^2)\). As letras maiúsculas \(X\) e \(Z\) representam a matriz de incidência dos referidos efeitos.
"
DATAFILE
dados_vol.txt
TRAITS
2
FIELDS_PASSED TO OUTPUT
WEIGHT(S)
RESIDUAL_VARIANCE
0.1
EFFECT
3 cross alpha
EFFECT
1 cross alpha
RANDOM
animal
FILE
pedigree.txt
(CO)VARIANCES
0.1
OPTION missing 0
OPTION sol se
OPTION se_covar_function H2 G_2_2_1_1 / (G_2_2_1_1 + R_1_1)
"
## [1] "\nDATAFILE\ndados_vol.txt\nTRAITS\n2\nFIELDS_PASSED TO OUTPUT\n\nWEIGHT(S)\n\nRESIDUAL_VARIANCE\n0.1\nEFFECT\n3 cross alpha\nEFFECT\n1 cross alpha\nRANDOM\nanimal\nFILE\npedigree.txt\n(CO)VARIANCES\n0.1\nOPTION missing 0\nOPTION sol se\nOPTION se_covar_function H2 G_2_2_1_1 / (G_2_2_1_1 + R_1_1)\n"
\(\sigma^2_a\) = 1.63 and \(SE\) = (0.20)
\(\sigma^2_e\) = 14.73 and \(SE\) = (0.19)
\(h^2_a\) = 0.0997