# NOT RUN {
# fit factorial model
data(seasonalsex)
seasonalsex.factorfit <- factorfit_cyclomort(event ~ sex, data = seasonalsex, n.seasons = 1)
# summary
summary(seasonalsex.factorfit, coefs = TRUE)
# }
# NOT RUN {
plot(seasonalsex.factorfit)
# }
Run the code above in your browser using DataLab