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 <iterado...@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@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
