Расстояние между временным рядом и представлением (кусочно-совокупная аппроксимация - PAA) (R) - PullRequest
Новая
       13

Расстояние между временным рядом и представлением (кусочно-совокупная аппроксимация - PAA) (R)

0 голосов
/ 06 сентября 2018

Я пытаюсь найти способ измерения расстояния между временным рядом и его представлением, чтобы увидеть, насколько точно представление описывает исходный ряд. Я чувствую, что есть простой способ сделать это, но я не могу понять это. График, сгенерированный приведенным ниже кодом, покажет красные линии для среднего значения каждого окна. Это расстояние между этой линией и каждой точкой в ​​этом окне, которое я пытаюсь вычислить.

Заранее спасибо, 

################################################################################
install.packages("fitdistrplus")
install.packages("seewave")
install.packages("jmotif")
install.packages("TSdist")


library(fitdistrplus)
library(seewave)
library(jmotif)
library(TSdist)

#####################################

# set size of the SAX alphabet 
alphabet_size <- 5

# set number of elements of the PAA (Number of windows)
PAA_number <- 10

#Generate data
 x <- rlnorm(100)
 time <- seq(1,100)

#SAX (seewave) ###
  # SAX - breakpoints = sample quantiles
  sax_2 <- SAX(x,  alphabet=alphabet_size, PAA=PAA_number, breakpoints="quantiles")

  ras <- 1/alphabet_size

  #PAA lenghts for graph
  y_paa <- paa(x, PAA_number)
  x_paa <- seq(time[1], time[length(time)], length.out=PAA_number+1)


  ### 8. Plot SAX (jmotif) / sample quantiles ###
  sample_bp <- quantile(x, probs = seq(ras, 1-ras, ras))
  title2 <- "SAX sample quantiles)"
  dev.new()
  plot(time, x, type="l", col="blue", ylab="", main=title2)
  segments(x_paa[-length(x_paa)], y_paa, x_paa[-1], y_paa, col="red", lwd=2)
  abline(h=sample_bp, col="green", lty=2)
  text(x_paa[-length(x_paa)]+(x_paa[2]-x_paa[1])/2, y_paa, sax_2, pos=3, cex=0.7)
  #dev.print(pdf, paste(x_name, "SAX sample.pdf"))
...