To fit the grey-zone model, one would need to call the functions in the order of em.func, cov.func, and greyzone.func.
The package also provides a function bestcut2 to fit a 2-group model, that is, it will find an optimal cutoff of the marker to divide patients into high and low 2 risk groups. Plus there is a function genSurvData to generate survival data with a fixed censoring rate.