R/plotStats.R

In akopich/covcp: R package for change-point detection in covariance structure

plot.CPDResult = function(statsAndCritValue) {
  stats = statsAndCritValue$statistics$window2statistics
  crits = statsAndCritValue$criticalValue
  
  for (index in 1:length(stats)) {
    forWindow = stats[[index]]
    crit = crits[index]
    plot(forWindow$centralPoints, 
         forWindow$distances, 
         col = ifelse(forWindow$distances > crit,"red", "black"), 
         xlab = "Central point",
         ylab = "Parameter distance", 
         ylim = c(0, max(crit, max(forWindow$distances))))
    title(paste("Window size = ", forWindow$windowSize))
    abline(h = crit, col = "red")
  }
}