[R-es] Representar datos longitudinales mediante splines

Javier Marcuzzi javier.ruben.marcuzzi en gmail.com
Jue Mar 3 22:06:35 CET 2016


Estimado Francisco Javier

Puede realizar eso con xy, o la parte gráfica que prefiera (hay que intentar y se llegaría), pero en su código no veo alguna función de ajuste de curva, spline permite ajustar curvas y fijar algunos puntos, hay varias posibilidades al respecto, debería seleccionar una librería con la palabra clave spline que cree apropiada. 
 Yo cuándo uso splines al data.frame le agrego las referencias al spline (una cantidad de columnas de acuerdo al modelo), y al realizar la gráfica la predicción correspondiente al ajunte de curvas. 

Javier Rubén Marcuzzi

De: Francisco Javier
Enviado: jueves, 3 de marzo de 2016 17:36
Para: r-help-es en r-project.org
Asunto: [R-es] Representar datos longitudinales mediante splines

Buenas noches a todos,


Me dirijo a vosotros porque estoy trabajando con medidas repetidas sobre un grupo de sujetos y quisiera graficar la evoluci�n temporal de cada uno de ellos pero suavizada mediante un spline (o en su defecto alg�n tipo de l�nea suavizadora que no necesariamente tenga que pasar por los puntos de cada sujeto). Adem�s, quisiera hacerlo en R b�sico, esto es, sin recurrir al uso de paquetes adicionales tipo lattice,  ggplot2.


A modo de ejemplo, tengo el siguiente data frame:


df <- data.frame(

                id = factor(rep(c(1,2,3), c(3,4,3))),
                x = c( 1, 14, 22,  2, 9,  20, 25,  4,  15,  22),
                y = c(35, 28, 52, 79, 64, 70, 95, 45, 102, 105))


E intento (sin conseguir resultado):


plot(df$x, df$y, xlab = "time", ylab="value", type="n", xlim = c(0,25), ylim=c(20,120))
for(i in unique(df$id)) {
 lines(lowess(df$x[df$id == i], df$y[df$id == i], col = df$id))  }
points(x = df$x, y = df$y, pch=16, col = df$id)              # Los datos asociados a cada id con un color de referencia
title("Scatter Plot of value vs. time")


�Muchas gracias por cualquier ayuda al respecto!


Francisco

	[[alternative HTML version deleted]]



	[[alternative HTML version deleted]]



Más información sobre la lista de distribución R-help-es