R - por favor, comparta sus experiencias - página 6

 
RandomWorker:

He introducido EURUSD en R. ¿Cómo puedo dibujar el gráfico y alinearlo con el kotir?


Líneas de la tramapuntos

Por cierto, acabas de estimar el modelo. La predicción se realiza mediante la función predict.

n <- 100
n.ahead <- 15
ts <- sin(0.05 * pi * (1 : n)) + rnorm(n = n, mean = 0, sd = 0.1)
ts.model <- ar(ts, method = 'burg')
ts.pred <- predict(ts.model, n.ahead = n.ahead)
plot(
  1 : (n + n.ahead),
  c(ts, rep(NA, n.ahead)),
  t = 'l'
  )
lines(
  (n + 1) : (n + n.ahead),
  as.numeric(ts.pred$pred),
  col = 'gray'
  )
lines(
  (n + 1) : (n + n.ahead),
  as.numeric(ts.pred$pred) + as.numeric(ts.pred$se),
  col = 'gray',
  lty = 'dashed'
  )
lines(
  (n + 1) : (n + n.ahead),
  as.numeric(ts.pred$pred) - as.numeric(ts.pred$se),
  col = 'gray',
  lty = 'dashed'
  )
 
anonymous:


?trazar ?líneas ?puntos

Por cierto, sólo ha estimado el modelo. La predicción se realiza mediante la función predict.

Gracias. Sin embargo, es difícil para mí. Lo intentaré.

¿Lo que escribiste se basará en la historia?

 
anonymous:


?trazar ?líneas ?puntos

Por cierto, sólo has estimado el modelo. La predicción se realiza mediante la función predict.

Me he dado cuenta un poco.

No me interesa la predicción.

Tengo un modelo que considero como indicador algo así como MA. ¿Cómo lo dibujo en el historial donde se calcula el modelo?

 
RandomWorker:

¿Cómo se dibuja en la historia sobre la que se calculó el modelo?


Líneas de la parcela

library(TTR)

x <- cumsum(rnorm(100))
x.ma <- EMA(x, 10)

plot(x, t = 'l')
lines(x.ma, col = 'red')
 
anonymous:


Líneas de la trama

> x<-ar(eur[1:256],method="mle")
> x

llamar:
ar(x = eur[1:256],method="mle")

Coeficientes:
1 2 3
0.9420 0.1955 -0.1644

Orden seleccionado 3 sigma^2 estimado como 2,73e-06
> plot(x, t = 'l')
Error en xy.coords(x, y, xlabel, ylabel, log) :

'x' es una lista, pero no tiene componentes 'x' e 'y'

Lo siento, pero debo ser extraordinariamente estúpido.

 
RandomWorker:
> plot(x, t = 'l')


Pruebe esto: plot(eur[1:256], t='l')

 
anonymous:


Inténtelo de esta manera: plot(eur[1:256], t='l')

Sin preguntas. Lo he dibujado. Pero es un vector.

x<-ar(eur[1:256],method="mle")

Aquí x no es un vector. Tiene una fórmula en su interior con la que se pueden hacer cálculos como en la predicción, sólo que sobre la historia

 
RandomWorker:

Sin preguntas. Lo he dibujado. Pero es un vector.

x<-ar(eur[1:256],method="mle")

Aquí x no es un vector. Tiene una fórmula en su interior con la que se pueden hacer cálculos como en el prefijo, sólo que en la historia


El comando ?ar nos dice que x es una lista, no una fórmula. También lo hace class(x).

El comando str(x) muestra el contenido del objeto.

El fitted(x) da NULL, pero x tiene un componente de residuo (residuos para los datos originales) - los valores del modelo todavía pueden ser calculados: (eur[1:256]-x$resid).

Plot(eur[1:256], t='l')

Añadir una línea de valores del modelo: lines(eur[1:256]-x$resid, col='red')

 
anonymous:


El comando ?ar nos dice que x es una lista, no una fórmula. También lo hace class(x).

El comando str(x) muestra el contenido del objeto.

fitted(x) da NULL, pero x tiene un componente de residuo (residuos para los datos brutos) - los valores del modelo aún pueden ser calculados: (eur[1:256]-x$resid).

Plot(eur[1:256], t='l')

Añadir una línea de valores del modelo: lines(eur[1:256]-x$resid, col='red')

¡Lo tengo! Gracias.

Pero de alguna manera es sorprendente que nosotros, teniendo una fórmula, hayamos calculado indirectamente su valor. El predicado está ahí, el resultado del ajuste no.

 

Algún tipo de ajuste extraño.

> x<-ar.ols(eur[1:256], order.max = 20, demean = TRUE)

> x


llamar:

ar.ols(x = eur[1:256], order.max = 20, demean = TRUE)


Coeficientes:

1 2 3

0.9425 0.1967 -0.1647

Orden.max = 20, pero sólo tres coeficientes. Yo lo entiendo así:

eur = 0,9425 + 0,1967*eur(-1) + (-0,1647)*eur(-2)

¿Se supone que hay 20 miembros?

Aparte de eso he cambiado diferentes parámetros ar siempre tres coeficientes y son muy ligeramente diferentes.

Por favor, comenten esto.

Razón de la queja: