[R-es] Eventos en una serie de tiempo

Mie Ene 7 18:23:07 CET 2015

Estimado Carlos Andres Perez Angarita

Entiendo lo que busca, o creo entenderlo, yo uso algo de esa tecnología, 
podría ser adaptable a sus requerimientos.

En primer lugar es más simple trabajar con los datos ordenados de otra 
forma, lo que usted dice columnas deberían ser filas, no es algo 
indispensable pero en el caso de miles de datos y muchas mediciónes a 
distintos tiempos, ..., facilita.

Facilita para ser leído e ingresado a lme4, nlme, etc. Porque una 
algernativa es usar modelos de regresiónes aleatorias.

Básicamente, lo que se usar, es realizar un ajuste de curva con 
legendre, polinomio, cuadrática, spline, ..., habría que ver en su caso 
que dicen las publicaciónes.

Cuándo el modelo realiza un ajuste de curva, hacer otro modelo pero 
pidiendo para cada individuo.

Con esos datos puede realizar la predicción para el día que desee (dado 
por el ajuste de curva para cada individuo), calcular entre fechas 
(derivadas, integrales y todas esas cosas muy lindas para los matemáticos).

Yo en su caso, leería alguna publicación sobre su área de la ciencia y 
buscaría una librería de R que tenga el método estadístico que 
corresponde o documentado, luego ordeno los datos para ser leídos en 
forma sencilla para esa herramienta específica dentro de R,  obtengo los 
valores del ajuste de curva con mis datos, y a partir de ahí busco 
extraer o calcular toda la inormación que me interese.

La forma de sus datos es semejante a las hojas de cálculo, lo cuál no es 
malo, pero pensar en R se aproxima a volúmenes de datos que no se 
almacenan en hojas de cálculos, lo cuál implica que los desarrolladores 
tengan otra forma de acomodar los datos, porque si usted tendría 1000 
columnas (fechas) donde no todas las fechas tienen datos para todos los 
individuos, es una complicación que fué resuelta al ordenar distinto 
(fila en columnas).

Yo pienso que si ordena sus datos en una forma más amigable con R podría 
encontrar muchos ejemplos. Pero posiblemente esto es muy influenciado 
por mi forma de trabajar con R, siendo ilógico para su caso en particular.

Javier Marcuzzi

El 07/01/2015 a las 10:53 a.m., Carlos Andres Perez Angarita escibió:
> Muchas gracias javier
>
> Si javier cada columna representa un individuo y todos los individuos 
> tienen curvas con comportamientos muy parecidas entre si y lo q me 
> referia de los eventos es poder dividir cada curva en 4 secciones 
> donde pueda hallar la posicion o el indice de cada fila los maximos o 
> minimos de esta seccion de curva. Es como dividir cada columna en 4 
> partes (no iguales), y como son 160 columnas seria muy dispendioso, 
> pense q con un algoritmo podria hallar el maximo o  minimo de cada columna
>
> CARLOS ANDRES
>
>
> El 7 de enero de 2015, 6:00, <r-help-es-request en r-project.org 
> <mailto:r-help-es-request en r-project.org>> escribió:
>
>     Envíe los mensajes para la lista R-help-es a
>     r-help-es en r-project.org <mailto:r-help-es en r-project.org>
>
>     Para subscribirse o anular su subscripción a través de la WEB
>     https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>     O por correo electrónico, enviando un mensaje con el texto "help" en
>     el asunto (subject) o en el cuerpo a:
>     r-help-es-request en r-project.org
>     <mailto:r-help-es-request en r-project.org>
>
>     Puede contactar con el responsable de la lista escribiendo a:
>     r-help-es-owner en r-project.org <mailto:r-help-es-owner en r-project.org>
>
>     Si responde a algún contenido de este mensaje, por favor, edite la
>     linea del asunto (subject) para que el texto sea mas especifico que:
>     "Re: Contents of R-help-es digest...". Además, por favor, incluya en
>     la respuesta sólo aquellas partes del mensaje a las que está
>     respondiendo.
>
>
>     Asuntos del día:
>
>        1. Re: Eventos de una serie de tiempo (Marcuzzi, Javier Rubén)
>
>
>     ----------------------------------------------------------------------
>
>     Message: 1
>     Date: Tue, 06 Jan 2015 11:43:54 -0300
>     From: "Marcuzzi, Javier Rubén"  <javier.ruben.marcuzzi en gmail.com
>     <mailto:javier.ruben.marcuzzi en gmail.com>>
>     To: r-help-es en r-project.org <mailto:r-help-es en r-project.org>
>     Subject: Re: [R-es] Eventos de una serie de tiempo
>     Message-ID: <54ABF4AA.4010107 en gmail.com
>     <mailto:54ABF4AA.4010107 en gmail.com>>
>     Content-Type: text/plain; charset="UTF-8"
>
>     Estimado Carlos
>
>     Creo que hay un error o por la simpleza lógica del correo se resume
>     queda a interpretación.
>
>     Compendí que tienes algo como lo siguiente:
>     dato, tiempo1, tiempo2, tiempo3, tiempo4 ...
>
>     Deseas que el resultado sea algo como
>     dato1, tiempo mínimo columna tiempo1, tiempo máximo columna tiempo4
>     dato2, tiempo mínimo columna tiempo1, tiempo máximo columna tiempo3
>     ...
>
>     Pero también pides los máximos y mínimos para esas "curvas", y es ahí
>     donde puede haber un error, porque: ¿los datos de sus tablas son los
>     observados o tiene la curva con los parámetros calculada con R?
>
>     Porque si son obervados puede buscarlos en la tabla.
>     Si es una curva debe buscarlos en la curva (en todos los puntos que
>     dibujan la curva), no en la tabla, en la preducción del modelo
>     para cada
>     fila (individio, dato, ...).
>
>     Javier Marcuzzi
>
>     El 06/01/2015 a las 01:28 a.m., Jorge I Velez escibió:
>     > Hola Carlos,
>     >
>     > No conozco una libreria, pero puedes hacer lo que necesitas
>     utilizando
>     >
>     > t(apply(series, 2, function(x) c(which.min(x), which.max(x))))
>     >
>     > donde "series" es el objeto que contiene las series de tiempo
>     con las que
>     > estas trabajando.
>     >
>     > En cuanto a los "eventos", no estoy seguro a que te refieres. 
>     Sin embargo,
>     > puede que el paquete strucchange en
>     > http://cran.r-project.org/web/packages/strucchange/index.html te
>     sirva de
>     > ayuda.
>     >
>     > Saludos cordiales,
>     > Jorge.-
>     >
>     >
>     > 2015-01-06 14:40 GMT+11:00 Carlos Andres Perez Angarita <
>     > caaperezan en gmail.com <mailto:caaperezan en gmail.com>>:
>     >
>     >> Saludo estimados compa?eros y compa?eras
>     >>
>     >> Tengo una matriz de datos de 51 filas por 160 columnas, Cada
>     columna es una
>     >> serie de tiempo.
>     >>
>     >> Existe alguna libreria q me calcule la posicion de la fila
>     donde se hallan
>     >> estos maximos y minimos o los eventos de estas curvas
>     >>
>     >> Agradezco la atencion
>     >>
>     >> CARLOS ANDRES
>     >>
>     >>          [[alternative HTML version deleted]]
>     >>
>     >> _______________________________________________
>     >> R-help-es mailing list
>     >> R-help-es en r-project.org <mailto:R-help-es en r-project.org>
>     >> https://stat.ethz.ch/mailman/listinfo/r-help-es
>     >>
>     >       [[alternative HTML version deleted]]
>     >
>     >
>     >
>     > _______________________________________________
>     > R-help-es mailing list
>     > R-help-es en r-project.org <mailto:R-help-es en r-project.org>
>     > https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
>             [[alternative HTML version deleted]]
>
>
>
>     ------------------------------
>
>     Subject: Pié de página del digest
>
>     _______________________________________________
>     R-help-es mailing list
>     R-help-es en r-project.org <mailto:R-help-es en r-project.org>
>     https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
>     ------------------------------
>
>     Fin de Resumen de R-help-es, Vol 71, Envío 7
>     ********************************************
>
>

	[[alternative HTML version deleted]]