Fits general rate model. The formula takes `reSurv` from `reReg`.

ssfit(
  formula,
  data,
  shp.ind = FALSE,
  bIni = NULL,
  rIni = NULL,
  optMethod = NULL,
  smoothFirst = NULL,
  interval = c(0, 2 * pi),
  h = NULL,
  weights = NULL
)

Arguments

formula

a formula object, with the response on the left of a '~' operator.

data

an optional argument for data.

shp.ind

a character string indicating either to proceed with shape independence assumption.

bIni

initial value for estimating beta

rIni

initial value for estimating gamma

smoothFirst

estimate the smooth estimator first and use it to guide? Default: FALSE when p = 1, TRUE when p > 1

opt

how many optimization methods to try? 1: optim 2: 1 + spg Default: 1 when p = 1, 2 when p > 1