<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.WordSection1
{page:WordSection1;}
--></style></head><body lang=ES link=blue vlink="#954F72"><div class=WordSection1><p class=MsoNormal>Estimado Carlos y Francisco</p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>El lowess no es spline <a href="http://www.inside-r.org/r-doc/stats/lowess">http://www.inside-r.org/r-doc/stats/lowess</a>.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Por la pregunta yo pensé en algo como muchos sujetos con medidas de peso corporal a lo largo del tiempo relacionado a dietas distintas, esa fue mi analogía mental para mi respuesta anterior, donde spline se encuentra dentro de las funciones de ajuste posible.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Javier Rubén Marcuzzi</p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><o:p> </o:p></span></p><div style='mso-element:para-border-div;border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal style='border:none;padding:0cm'><b>De: </b><a href="mailto:cof@qualityexcellence.es">Carlos Ortega</a><br><b>Enviado: </b>viernes, 4 de marzo de 2016 3:27<br><b>Para: </b><a href="mailto:iterador10@hotmail.com">Francisco Javier</a><br><b>CC: </b><a href="mailto:r-help-es@r-project.org">r-help-es@r-project.org</a><br><b>Asunto: </b>Re: [R-es] Representar datos longitudinales mediante splines</p></div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><o:p> </o:p></span></p><div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>Hola,</span><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>El suavizado con "lowess()" te da error:<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>> for(i in unique(df$id)) {<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>+ lines(lowess(df$x[df$id == i], df$y[df$id == i], col = df$id)) }<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>Error in lowess(df$x[df$id == i], df$y[df$id == i], col = df$id) : <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'> unused argument (col = df$id)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>Y el mensaje de error es bastante explícito. Incluyes el parámetro "col" dentro de lowess que no lo considera.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>Es un simple problema de dónde has colocado un paréntesis de cierre.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>Así sí que funciona (destaco en rojo la posición del paréntesis de cierre):<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>#-----------------<o:p></o:p></span></p></div><div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>for(i in unique(df$id)) {<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'> lines(lowess(df$x[df$id == i], df$y[df$id == i]<b><span style='color:red'>)</span></b>, col = df$id) }<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'> points(x = df$x, y = df$y, pch=16, col = df$id) <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'># Los datos asociados a cada id con un color de referencia<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>title("Scatter Plot of value vs. time")<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>#-----------------<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><img border=0 width=472 height=248 style='width:4.9166in;height:2.5833in' id="Imagen_x0020_3" src="cid:image003.png@01D175EB.D4C241C0" alt="Imágenes integradas 1"></span><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>Saludos,<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>Carlos Ortega<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><a href="http://www.qualityexcellence.es">www.qualityexcellence.es</a><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'><o:p> </o:p></span></p></div></div></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><o:p> </o:p></span></p><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>El 3 de marzo de 2016, 21:35, Francisco Javier <<a href="mailto:iterador10@hotmail.com" target="_blank">iterador10@hotmail.com</a>> escribió:<o:p></o:p></span></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm'><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>Buenas noches a todos,<br><br><br>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.<br><br><br>A modo de ejemplo, tengo el siguiente data frame:<br><br><br>df <- data.frame(<br><br> id = factor(rep(c(1,2,3), c(3,4,3))),<br> x = c( 1, 14, 22, 2, 9, 20, 25, 4, 15, 22),<br> y = c(35, 28, 52, 79, 64, 70, 95, 45, 102, 105))<br><br><br>E intento (sin conseguir resultado):<br><br><br>plot(df$x, df$y, xlab = "time", ylab="value", type="n", xlim = c(0,25), ylim=c(20,120))<br>for(i in unique(df$id)) {<br> lines(lowess(df$x[df$id == i], df$y[df$id == i], col = df$id)) }<br>points(x = df$x, y = df$y, pch=16, col = df$id) # Los datos asociados a cada id con un color de referencia<br>title("Scatter Plot of value vs. time")<br><br><br>¡Muchas gracias por cualquier ayuda al respecto!<br><br><br>Francisco<br><br> [[alternative HTML version deleted]]<br><br><br>_______________________________________________<br>R-help-es mailing list<br><a href="mailto:R-help-es@r-project.org">R-help-es@r-project.org</a><br><a href="https://stat.ethz.ch/mailman/listinfo/r-help-es" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-help-es</a><o:p></o:p></span></p></blockquote></div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><br><br clear=all><o:p></o:p></span></p><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><o:p> </o:p></span></p></div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>-- <o:p></o:p></span></p></div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Verdana",sans-serif'>Saludos,<br>Carlos Ortega<br><a href="http://www.qualityexcellence.es" target="_blank">www.qualityexcellence.es</a></span><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>