R: Maximum Likelihood Parameter Estimation of a Refractory...
rexpMLE
R Documentation
Maximum Likelihood Parameter Estimation of a Refractory Exponential Model with Possibly
Censored Data
Description
Estimate refractory exponential model parameters by the maximum likelihood
method using possibly censored data.
Usage
rexpMLE(yi, ni = numeric(length(yi)) + 1,
si = numeric(length(yi)) + 1)
Arguments
yi
vector of (possibly binned) observations or a
spikeTrain object.
ni
vector of counts for each value of yi; default: numeric(length(yi))+1.
si
vector of counts of uncensored observations for each
value of yi; default: numeric(length(yi))+1.
Details
The MLE are available in closed form even in the censored case for
this model. The likelihood function cannot be differentiated with
respect to the rp (refractory period) parameter at the
maximum. COnfidence intervals for this parameter are therefore not
available.
Value
A list of class durationFit with the following components:
estimate
the estimated parameters, a named vector.
se
the standard errors, a named vector.
logLik
the log likelihood at maximum.
r
a function returning the log of the relative likelihood function.
mll
a function returning the opposite of the log likelihood
function using the log of the parameters.