powered by
Fit survival function, used for Bayesian simulation
fitsurv(parm, Dm, eta)
the initial parameter value
the D matrix
equals to X*bea