Carlos Andres Perez Angarita
2015-Jan-07 13:53 UTC
[R-es] Resumen de R-help-es, Vol 71, Envío 7
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> escribió:> Envíe los mensajes para la lista R-help-es a > 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 > > Puede contactar con el responsable de la lista escribiendo a: > 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> > To: r-help-es en r-project.org > Subject: Re: [R-es] Eventos de una serie de tiempo > Message-ID: <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>: > > > >> 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 > >> 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 > > 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 > 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]]
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]]