The array of starting parameters is built from the time and threshold data in the object, obj
.
Each parameter is assumed to have a possible range given the data.
Each range is constructed as follows; the time points; alpha (cone-rod transition \(\alpha\) minutes) and beta (rod-rod transition \(\beta\) minutes) are assumed to fall in the first and second halves of the time data respectively (obj$time
).
The cone threshold is assumed to be in the upper half of the threshold data (obj$thrs
log units). The cone coefficient (log units) or threshold at time zero is presumed to be positive and the same values are used for the time constants (\(tau\) minutes).
The rate of rod recovery S2, and the combined parameter -(S2 + S3) are ranged between -0.6 and 0 log units/minute.
These ranges of possible values are complied into an array by sampling without replacement from each range for each parameter.