femininity <- mutate_branch(
1 * (MasFem > 6), 1 * (MasFem > mean(MasFem))
)
intensity <- mutate_branch(
Minpressure_Updated_2014,
Category,
NDAM,
HighestWindSpeed
)
model <- formula_branch(
log(alldeaths + 1) ~ femininity,
log(alldeaths + 1) ~ femininity * intensity
)
mv <- mverse(hurricane) %>%
add_mutate_branch(femininity) %>%
add_mutate_branch(intensity) %>%
add_formula_branch(model) %>%
lm_mverse()
spec_summary(mv, "femininity")
Run the code above in your browser using DataLab