Pessoal,
quero encontrar o horizonte de previsões, usando o expoente de Lyapunov.
O exemplo proposto no R, a serie de lorenz.ts funciona. Porém, quando troco
pela minha variável que estou predizendo, da erro. Gostaria de entender
qual parâmetro estou colocando errado. Segue o código mínimo reproduzível.
#=====================================
set.seed(22)
y_obs <- 1:100 + rnorm(100,mean=0, sd=5) # supondo essa ser a variavel
observada (100 observações temporais)
y_estimado <- 1:80 # supondo essa ser a variavel estimada, até o tempo
80
y_pred <- 81:100 + rnorm(20,mean=0, sd=3) # suponha que usando uma rede
recorrente chegue nesse resultado, como codigo minimo reproduzido,
#simplifiquemos com uma rnorm com tendencia linear, para predição usando 80
tempos de predição
plot(y_obs)
points(1:80,1:80, type = "l")
points(81:100,y_pred, type = "l", col="blue")
legend("topleft", c("Observado","Estimado", "Predito"),
col=c("black","black", "blue"), pch = c(1,NA, NA), lty = c(NA,1,
1))
# Horizonte de predicao usando o expoente de Lyapunov
require("tseriesChaos")
y_pred = as.ts(y_pred)
output <-lyap_k(y_pred, m=3, d=2, s=200, t=2, ref=2, k=2, eps=4)
par = lyap(output, 0.73, 2.47)
lambda = par[2]
plot(output)
#===================================
O erro é esse:
> output <-lyap_k(y_pred, m=3, d=2, s=200, t=2, ref=2, k=2, eps=4)
Finding nearests
Error in find_knearests(series, m = m, d = d, t = t, ref = length(ref), :
cannot allocate memory block of size 134217728 Tb
OBS: No meu caso, já havia testado e a dimensão m e d, são exatamente 3 e
2. Porém, não sei como controlar os demais parâmetros.
Fico no aguardo, se alguém conseguir me explicar o problema.
Agradeço e obrigado.
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
