Survival analysis approach consists of the analysis of the time it takes for events to occur. The data consist in a set of individuals that are followed until the event of interest occurs. This repository shows how to analyse seed germination (event) data through a time-to-event approach, using a data set of seeds of Ipomoea hederacea (Convolvulaceae). The analysis is part of an assessment of inbreeding depression in this species carried out in a population in the State of Veracruz (Mexico) during 2012 and 2013 (Delgado-Dávila and Martén-Rodríguez, 2021).
Because the study design include only a short interval of time to evaluate if the event occurs, the data analysed here corresponds to censored data, i.e. not all seeds had germinated by the time the monitoring period was over.
The objective is to evaluate the number of days seeds took to germinate through a survival analysis (Clark et al., 2003; Onofri et al., 2010), by fitting cure models to test for differences in germination times between pollination treatments per year (Onofri et al., 2011).
The complete data set of the study is available through Zenodo. For this repository only the file Germination_ID.csv is used. The file includes the identity of the plants and the fruits, the seeds weight, and the day of germination. The seeds were obtained after the application of the following pollination experiments in 38 plants:
-
Hand-self pollination (HSP): flowers were manually pollinated with pollen from the two longest anthers of the flower.
-
Hand-cross pollination (HCP): emasculated flowers were manually pollinated with pollen obtained from the two longest anthers of two-three donor plants.
The file CureModel.Rmd
shows how to fit a parametric Cure Model following Onofri et al. 2011, using the R package flexsurvcure
(Amdahl, 2017). The file also shows how to plot the survival curves to graphical comparison between observed and predicted values.