R/quantile.fts.R

In ftsa: Functional Time Series Analysis

quantile.fts <- function(x, probs, ...)
{
  if(missing(probs))
  {
	probs = c(0, 0.25, 0.5, 0.75, 1)
  }
  if (class(x)[1] == "fts"|class(x)[1] == "fds"|class(x)[1] =="sfts"){
      y = x$y
      p = dim(y)[1]
      q = matrix(, p, length(probs))
      for(i in 1:p){
          q[i,] = stats::quantile(y[i,], probs = probs)
      }
      colnames(q) = paste(round(100*probs,0), "%", sep="")
      rownames(q) = x$x
      if (class(x)[1] == "fds"){
          warning("Object is not a functional time series.")
      }
      return(q)
   }
   else {
        stop("Not a functional object.")
   }
}