d <- read.csv("Raw data of Table1 and2_migration_type_and_wetland_use_area.csv") result <- glm(home_range_in_wetland~Migration_type+capture_site+offset(log(annual_home_range)),data=d) summary(result) anova(result, test="Chi")