Model Summary

Model Summary

# S3 method for slouch
summary(object, ...)

Arguments

object	An object of class 'slouch'
...	Additional arguments, unused.

Examples

data(artiodactyla)
data(neocortex)

neocortex <- neocortex[match(artiodactyla$tip.label, neocortex$species), ]

m0 <- slouch.fit(phy = artiodactyla,
                 hl_values = seq(0.001, 4, length.out = 10),
                 vy_values = seq(0.001, 0.05, length.out = 10),
                 species = neocortex$species,
                 response = neocortex$neocortex_area_mm2_log_mean,
                 mv.response = neocortex$neocortex_se_squared,
                 random.cov = neocortex$brain_mass_g_log_mean,
                 mv.random.cov = neocortex$brain_se_squared,
                 fixed.fact = neocortex$diet,
                 hillclimb = FALSE)

summary(m0)
#> Important - Always inspect the likelihood surface of the model parameters with
#>           grid search before evaluating model fit & results.
#> 
#>   Maximum-likelihood estimates
#>                        Estimate
#> Phylogenetic half-life 2.667   
#> Stationary variance    0.001   
#> 
#>   Stochastic predictor(s)
#>                                 Phylogenetic mean Diffusion variance
#> neocortex$brain_mass_g_log_mean          5.139644         0.02922457
#> 
#>   Inferred maximum-likelihood parameters
#>                                     Value       
#> Mean phylogenetic correction factor 0.8586099   
#> Rate of adaptation                  0.2598977   
#> Diffusion variance                  0.0005197954
#> 
#>   Interval of parameters in 3d plot (Sensitive to grid mesh, grid size and local ML estimate)
#>                        Minimum    Maximum
#> Phylogenetic half-life   0.001 4.00000000
#> Stationary variance      0.001 0.03911111
#> 
#>   Regime optima
#>    Estimates Std. error
#> Br  5.323085  0.2175527
#> Gr  5.465051  0.2661451
#> MF  5.449717  0.2314937
#> 
#>   Optimal regression slope
#>                                 Estimates Std. error
#> neocortex$brain_mass_g_log_mean 0.8444137  0.0429896
#> 
#>   Evolutionary regression slope
#>                                 Predictions Std. error
#> neocortex$brain_mass_g_log_mean   0.8444137  0.0429896
#> 
#>   Attenuation factor. Linear model coefficients (above) are not corrected for bias.
#>                                              
#> Br                              1 0 0 -0.0328
#> Gr                              0 1 0 -0.0354
#> MF                              0 0 1 -0.0342
#> neocortex$brain_mass_g_log_mean 0 0 0  1.0100
#> 
#>   Model fit summary
#>             Values
#> Support     14.100
#> AIC        -16.200
#> AICc       -13.800
#> SIC         -5.600
#> R squared    0.927
#> SST        550.000
#> SSE         40.000
#> N (params)   6.000

plot(m0, theta = 150)