Subset of eggs_and_nests data the contains only Dome, Cup, and ground surface nests

Subset of eggs_and_nests data the contains only Dome, Cup, and ground surface nests

eggs_and_nests2

Format

A data frame

family

Taxonomic family

order

Taxonomic order

spp

Genus and species

asymmetry

Index of how asymmetric an egg is

ellipticity

Index of how elliptical an egg is

ave.length.com

Mean length of egg

nest

nest type of species

source

author

bird.ellipticity

ellipticity of bird in grams

References

....

Examples

library(ggplot2)
library(ggpubr)

ggerrorplot(eggs_and_nests2,
x = "nest",
y = "ellipticity",
desc_stat = "mean_ci")

## 1-way ANOVA

### null model
model.null <- lm(ellipticity ~ 1, data = eggs_and_nests2)

### model of interest
model.alt <- lm(ellipticity ~ nest, data = eggs_and_nests2)

### compare models
anova(model.null, model.alt)
#> Analysis of Variance Table
#> 
#> Model 1: ellipticity ~ 1
#> Model 2: ellipticity ~ nest
#>   Res.Df    RSS Df Sum of Sq      F    Pr(>F)    
#> 1    700 3.2337                                  
#> 2    698 3.0981  2   0.13554 15.268 3.235e-07 ***
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1


## Tukey test

### re-fit model with aov()
model.alt.aov <- aov(ellipticity ~ nest, data = eggs_and_nests2)

### TukeyHSD() on model from aov()
TukeyHSD(model.alt.aov)
#>   Tukey multiple comparisons of means
#>     95% family-wise confidence level
#> 
#> Fit: aov(formula = ellipticity ~ nest, data = eggs_and_nests2)
#> 
#> $nest
#>                           diff         lwr        upr     p adj
#> Dome-Cup            0.01768136 0.003768425 0.03159430 0.0082470
#> ground surface-Cup  0.07301529 0.038175706 0.10785487 0.0000032
#> ground surface-Dome 0.05533393 0.019116251 0.09155161 0.0010362
#> 

### Plot effect sizes
plotTukeysHSD(TukeyHSD(model.alt.aov))