Print the gamem
object in two ways. By default, the results are shown
in the R console. The results can also be exported to the directory.
# S3 method for gamem print(x, export = FALSE, file.name = NULL, digits = 4, ...)
x | An object fitted with the function |
---|---|
export | A logical argument. If |
file.name | The name of the file if |
digits | The significant digits to be shown. |
... | Options used by the tibble package to format the output. See
|
Tiago Olivoto tiagoolivoto@gmail.com
# \donttest{ library(metan) alpha <- gamem(data_alpha, gen = GEN, rep = REP, block = BLOCK, resp = YIELD )#>#>#>#>#> --------------------------------------------------------------------------- #> P-values for Likelihood Ratio Test of the analyzed traits #> --------------------------------------------------------------------------- #> model YIELD #> Complete NA #> Genotype 1.18e-06 #> rep:block 3.35e-03 #> --------------------------------------------------------------------------- #> All variables with significant (p < 0.05) genotype effectprint(alpha)#> Variable YIELD #> --------------------------------------------------------------------------- #> Fixed-effect anova table #> --------------------------------------------------------------------------- #> # A tibble: 1 x 7 #> SOURCE `Sum Sq` `Mean Sq` NumDF DenDF `F value` `Pr(>F)` #> <chr> <dbl> <dbl> <int> <dbl> <dbl> <dbl> #> 1 REP 1.381 0.6907 2 11.76 8.463 0.005287 #> --------------------------------------------------------------------------- #> Variance components for random effects #> --------------------------------------------------------------------------- #> # A tibble: 3 x 2 #> Group Variance #> <chr> <dbl> #> 1 GEN 0.1429 #> 2 REP:BLOCK 0.07022 #> 3 Residual 0.08162 #> --------------------------------------------------------------------------- #> Likelihood ratio test for random effects #> --------------------------------------------------------------------------- #> # A tibble: 3 x 7 #> model npar logLik AIC LRT Df `Pr(>Chisq)` #> <chr> <int> <dbl> <dbl> <dbl> <dbl> <dbl> #> 1 Complete 6 -46.60 105.2 NA NA NA #> 2 Genotype 5 -58.41 126.8 23.62 1 0.000001176 #> 3 rep:block 5 -50.90 111.8 8.606 1 0.003351 #> --------------------------------------------------------------------------- #> Details of the analysis #> --------------------------------------------------------------------------- #> # A tibble: 6 x 2 #> Parameters Values #> <chr> <chr> #> 1 Ngen 24 #> 2 OVmean 4.4795 #> 3 Min 2.8873 (G03 in B6 of R3) #> 4 Max 5.8757 (G05 in B1 of R1) #> 5 MinGEN 3.3431 (G03) #> 6 MaxGEN 5.1625 (G01) #> --------------------------------------------------------------------------- #> Genetic parameters #> --------------------------------------------------------------------------- #> # A tibble: 13 x 2 #> Parameters Values #> <chr> <dbl> #> 1 Gen_var 0.1429 #> 2 Gen (%) 48.48 #> 3 rep:block_var 0.07022 #> 4 rep:block (%) 23.82 #> 5 Res_var 0.08162 #> 6 Res (%) 27.69 #> 7 Phen_var 0.2947 #> 8 H2 0.4848 #> 9 h2mg 0.7980 #> 10 Accuracy 0.8933 #> 11 CVg 8.439 #> 12 CVr 6.378 #> 13 CV ratio 1.323 #> #> #># }