[R-es] Representar datos longitudinales mediante splines

Carlos Ortega cof en qualityexcellence.es
Vie Mar 4 07:26:59 CET 2016


Hola,

El suavizado con "lowess()" te da error:

> for(i in unique(df$id)) {
+   lines(lowess(df$x[df$id == i], df$y[df$id == i], col = df$id))  }
Error in lowess(df$x[df$id == i], df$y[df$id == i], col = df$id) :
  unused argument (col = df$id)

Y el mensaje de error es bastante explícito. Incluyes el parámetro "col"
dentro de lowess que no lo considera.

Es un simple problema de dónde has colocado un paréntesis de cierre.

Así sí que funciona (destaco en rojo la posición del paréntesis de cierre):

#-----------------
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")
#-----------------

[image: Imágenes integradas 1]


Saludos,
Carlos Ortega
www.qualityexcellence.es



El 3 de marzo de 2016, 21:35, Francisco Javier <iterador10 en hotmail.com>
escribió:

> 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]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20160304/1d5def66/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 25195 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20160304/1d5def66/attachment-0001.png>


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