R/slupacf.R

In statmanrobin/ts343AJ: Time Series Functions and Data for Stat 343

#'Compute the autocorrelation of a time series
#'
#' @param series A time series
#' @param maxlag Number of partialautocorrelations to return (Null default picks a reasoanble number)
#' @param plot Plot of the PACF with bounds? (default is TRUE)
#' @param ndiff Number of regular differences (default is 0)
#' @param sdiff Number of seasonal differences (default is zero)
#' @param out Output the numerical ACF? (default is TRUE)
#' @param ci.type Type of confidence bounds ("ma" defualt or "white" for white noise)
#'
#' @return An partial autocorrelation object
#'
#' @export
#'

slupacf=function(series,maxlag=NULL,plot=TRUE,ndiff=0,sdiff=0,out=TRUE,ci.type="ma"){
  sname=deparse(substitute(series))
  if(ndiff>0){series=diff(series,differences=ndiff)
  sname=paste(sname,"d=",ndiff)}
  if(sdiff>0){series=diff(series,differenes=sdiff,lag=frequency(series))
  sname=paste(sname,"D=",sdiff)}
  seriespacf=pacf(series,lag.max=maxlag,plot=F)
  seriespacf$series=sname
  if(frequency(series)>1){seriespacf$lag=seriespacf$lag*frequency(series)}
  if(plot){
    marold=par("mar")    #adjust margins to handle title better
    par(mar=c(5.1,4.1,4.1,2.1))
    plot(seriespacf)
    par(mar=marold)}
  if(out){return(seriespacf)}
}