Nothing
#' Maximum Likelihood Estimation of the IAR-T model
#'
#' Maximum Likelihood Estimation of the IAR-T model.
#'
#' @param y Array with the time series observations
#' @param st Array with the irregular observational times
#' @param nu degrees of freedom
#'
#' @return A list with the following components:
#' \itemize{
#' \item{phi}{ MLE of the phi parameter of the IAR-T model.}
#' \item{sigma}{ MLE of the sigma parameter of the IAR-T model.}
#' \item{ll}{ Value of the negative log likelihood evaluated in phi and sigma.}
#' }
#' @export
#' @references
#' \insertRef{Eyheramendy_2018}{iAR}
#'
#' @seealso
#'
#' \code{\link{gentime}}, \code{\link{IARtsample}}, \code{\link{IARphit}}
#'
#' @examples
#' n=300
#' set.seed(6714)
#' st<-gentime(n)
#' y<-IARtsample(n,0.9,st,sigma2=1,nu=3)
#' model<-IARt(y$y, st=st)
#' phi=model$phi
#' sigmaest=model$sigma
IARt<-function (y, st,nu=3)
{
aux<-1e10
value<-1e10
br<-0
for(i in 1:20)
{
phi=runif(1)
sigma=var(y)*runif(1)
optim<-nlminb(start=c(phi,sigma),objective=IARphit,y=y,st=st,nu=nu,yest=0,lower=c(0,0.0001),upper=c(0.9999,2*var(y)))
value<-optim$objective
if(aux>value)
{
par<-optim$par
aux<-value
br<-br+1
}
if(aux<=value & br>10 & i>15)
break;
}
if(aux==1e10)
par<-c(0,0)
return(list(phi=par[1],sigma=par[2],ll=aux))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.