Find initial estimates for model parameters using smoothing splines.
init_mle(ModelInfo, Dat)Model information object from set_model_info().
Data frame containg x and y data.
Names vector containing initial parameters estimates.
## Simulate data and get initial model parameter estimates
Models <- set_models (mean_fun="gaussian", err_dist="zip")
ModelInfo <- set_model_info (Models[1,])
Pars <- create_default_par_list (models=Models)
Data <- create_simulated_datasets (Pars, seed=12345)
theta <- init_mle (ModelInfo, Dat=data.frame(x=Data$x, y=Data$gaussian_zip))
print (rbind (Pars[[1]]$theta, theta))
#> H m s pi
#> 60.00000 50.00000 10.00000 0.3000000
#> theta 45.44242 49.23082 17.99716 0.1191357