Nothing
#' Fitted Values of IAR model
#'
#' Fit an IAR model to an irregularly observed time series.
#'
#' @param phi Estimated phi parameter by the iAR model.
#' @param y Array with the time series observations.
#' @param st Array with the irregular observational times.
#' @param standardized logical; if TRUE, the array y is standardized; if FALSE, y contains the raw time series
#' @param zero.mean logical; if TRUE, the array y has zero mean; if FALSE, y has a mean different from zero.
#'
#' @return Fitted values of the iAR model
#' @export
#' @references
#' \insertRef{Eyheramendy_2018}{iAR}
#'
#' @seealso
#' \code{\link{gentime}}, \code{\link{IARsample}}, \code{\link{IARloglik}}, \code{\link{IARkalman}}
#' @examples
#'
#' set.seed(6714)
#' st<-gentime(n=100)
#' y<-IARsample(phi=0.99,st=st,n=100)
#' y<-y$series
#' phi=IARloglik(y=y,st=st)$phi
#' fit=IARfit(phi=phi,y=y,st=st)
IARfit<-function(phi,y,st,standardized=TRUE,zero.mean=TRUE)
{
sigma = 1
mu = 0
if (standardized == FALSE)
sigma = var(y)
if (zero.mean == TRUE)
mu = mean(y)
if (zero.mean == TRUE)
y = y-mu
if (standardized == FALSE)
y=y/sqrt(sigma)
delta=diff(st)
y1=y[-c(length(y))]
fit=c(0,(phi**delta)*y1)
fit=(fit*sqrt(sigma)+mu)
return(fitted=fit)
}
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.