В коде где-то ошибка - PHP - Ответ 9774493

library(quantmod)
library(fractaldim)
 
# https://Asirikuy.com
# http://mechanicalforex.com
 
getSymbols("SPY",src="yahoo", from="1995-01-01")
 
# Predict the last 100 bars
endingIndex <-length(mainData$data)-101
mainData <- SPY$SPY.Adjusted
colnames(mainData) <- c("data")
TEST <- mainData[1:endingIndex]
TEST
total_error <- 0
error_per_prediction <- c()
 
#These are the fractal dimension calculation parameters
#see the fractaldim library reference for more info
 
method <- "rodogram"
 
Sm <- as.data.frame(TEST, row.names = NULL)
 
delta <- c()
 
# calculate delta between consecutive Sm values to use as guesses
for(j in 2:length(Sm$data)){
  delta <- rbind(delta, (Sm$data[j]-Sm$data[j-1])/Sm$data[j-1])
}
 
Sm_guesses <- delta
 
#do 100 predictions of next values in Sm
for(i in 1:100){
  
  #update fractal dimension used as reference
  V_Reference <- fd.estimate(Sm$data, method=method)$fd
   minDifference = 1000000
  
  # check the fractal dimension of Sm plus each different guess and
  # choose the value with the least difference with the reference
  for(j in 1:length(Sm_guesses)){
    
    new_Sm <- rbind(Sm, Sm_guesses[j]*Sm$data[length(Sm$data)]+Sm$data[length(Sm$data)])
    new_V_Reference <- fd.estimate(new_Sm$data, method=method)$fd
    
    if (abs(new_V_Reference - V_Reference) < minDifference ){******
      Sm_prediction <- Sm$data[length(Sm$data)]+Sm_guesses[j]*Sm$data[length(Sm$data)]
      minDifference = abs(new_V_Reference - V_Reference)
    }
  }
  
  print(i)
  #add prediction to Sm
  Sm <- rbind(Sm, Sm_prediction)
  Sm_real <- as.numeric(mainData$data[endingIndex+i])
  error_per_prediction <- rbind(error_per_prediction, (Sm_prediction-Sm_real )/Sm_real )
  total_error <- total_error + ((Sm_prediction-Sm_real )/Sm_real )^2
}
 
total_error <- sqrt(total_error)
print(total_error)
plot(error_per_prediction*100, xlab="Prediction Index", ylab="Error (%)")
plot(Sm$data[endingIndex:(endingIndex+100)], type="l", xlab="Value Index", ylab="Adjusted Close")
lines(as.data.frame(mainData$data[endingIndex:(endingIndex+100)], row.names = NULL), col="blue")