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@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]]

_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Responder a