Print the ge_reg
object in two ways. By default, the results are shown
in the R console. The results can also be exported to the directory into a
*.txt file.
Usage
# S3 method for ge_reg
print(x, export = FALSE, file.name = NULL, digits = 3, ...)
Arguments
- x
An object of class
ge_reg
.- export
A logical argument. If
TRUE
, a *.txt file is exported to the working directory.- file.name
The name of the file if
export = TRUE
- digits
The significant digits to be shown.
- ...
Options used by the tibble package to format the output. See
tibble::print()
for more details.
Author
Tiago Olivoto tiagoolivoto@gmail.com
Examples
# \donttest{
library(metan)
model <- ge_reg(data_ge2, ENV, GEN, REP, PH)
#> Evaluating trait PH |============================================| 100% 00:00:01
print(model)
#> Variable PH
#> ---------------------------------------------------------------------------
#> Joint-regression Analysis of variance
#> ---------------------------------------------------------------------------
#> # A tibble: 20 × 6
#> SV Df `Sum Sq` `Mean Sq` `F value` `Pr(>F)`
#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1 "Total" 51 15.0 0.294 NA NA
#> 2 "GEN" 12 1.86 0.155 0.870 0.585
#> 3 "ENV + (GEN x ENV)" 39 13.1 0.336 NA NA
#> 4 "ENV (linear)" 1 7.72 7.72 NA NA
#> 5 " GEN x ENV (linear)" 12 0.755 0.0629 0.352 0.969
#> 6 "Pooled deviation" 26 4.64 0.179 NA NA
#> 7 "H1" 2 0.471 0.235 10.5 0.0000756
#> 8 "H10" 2 0.236 0.118 5.25 0.00685
#> 9 "H11" 2 0.211 0.106 4.71 0.0112
#> 10 "H12" 2 0.450 0.225 10.0 0.000110
#> 11 "H13" 2 0.410 0.205 9.15 0.000231
#> 12 "H2" 2 0.534 0.267 11.9 0.0000241
#> 13 "H3" 2 0.584 0.292 13.0 0.00000995
#> 14 "H4" 2 0.354 0.177 7.89 0.000671
#> 15 "H5" 2 0.140 0.0701 3.13 0.0484
#> 16 "H6" 2 0.334 0.167 7.44 0.000994
#> 17 "H7" 2 0.241 0.121 5.38 0.00609
#> 18 "H8" 2 0.505 0.253 11.3 0.0000403
#> 19 "H9" 2 0.172 0.0861 3.84 0.0249
#> 20 "Pooled error" 96 2.15 0.0224 NA NA
#> ---------------------------------------------------------------------------
#> Regression parameters
#> ---------------------------------------------------------------------------
#> # A tibble: 13 × 10
#> GEN b0 b1 `t(b1=1)`[,1] pval_t[,1] s2di F(s2d…¹ pval_f RMSE R2
#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1 H1 2.62 0.806 -0.997 0.321 0.0710 10.5 7.56e-5 0.198 0.451
#> 2 H10 2.31 1.22 1.12 0.265 0.0318 5.25 6.85e-3 0.140 0.789
#> 3 H11 2.39 1.08 0.411 0.682 0.0278 4.71 1.12e-2 0.133 0.766
#> 4 H12 2.44 0.465 -2.75 0.00711 0.0676 10.0 1.10e-4 0.194 0.222
#> 5 H13 2.54 0.306 -3.57 0.000558 0.0609 9.15 2.31e-4 0.185 0.119
#> 6 H2 2.60 0.963 -0.188 0.851 0.0815 11.9 2.41e-5 0.211 0.508
#> 7 H3 2.59 1.35 1.83 0.0711 0.0898 13.0 9.95e-6 0.221 0.651
#> 8 H4 2.58 1.27 1.38 0.171 0.0515 7.89 6.71e-4 0.172 0.730
#> 9 H5 2.57 1.17 0.887 0.378 0.0159 3.13 4.84e-2 0.108 0.853
#> 10 H6 2.56 0.936 -0.330 0.742 0.0481 7.44 9.94e-4 0.167 0.609
#> 11 H7 2.40 0.992 -0.0393 0.969 0.0328 5.38 6.09e-3 0.142 0.708
#> 12 H8 2.33 1.01 0.0720 0.943 0.0767 11.3 4.03e-5 0.205 0.547
#> 13 H9 2.36 1.42 2.18 0.0317 0.0212 3.84 2.49e-2 0.120 0.875
#> # … with abbreviated variable name ¹`F(s2di=0)`
#> ---------------------------------------------------------------------------
#> Variance of b0: 0.001868774
#> Variance of b1: 0.03777002
#>
#>
#>
# }