```r
f1 <- fit_cyclomort(gazelle$T, n.seasons = 1)
f1
```

```
## Multi-seasonal hazard function fit with 1 seasons with periodicity 1.
## 
## $meanhazard
##    estimate    CI.low   CI.high        se
## 1 0.4139975 0.2256053 0.6023897 0.0941961
## 
## $`season 1`
##   parameter season  estimate       CI.low  CI.high
## 1      peak      1 0.1980790 -0.769357793 1.165516
## 2  duration      1 0.4687249  0.003530531 0.500000
## 
## Log-likelihood: -37.5215; AIC: 81.043
```

```r
plot(f1)
```

![](GazelleAnalysis_files/figure-latex/unnamed-chunk-9-1.pdf)<!-- -->