<font><font face="arial,helvetica,sans-serif">Estimado Javier,</font></font><div><font><font face="arial,helvetica,sans-serif"><br></font></font></div><div><font><font face="arial,helvetica,sans-serif">Quizas lo siguiente pueda servirte:</font></font></div>
<div><font><font face="arial,helvetica,sans-serif"><br></font></font></div><div><font><font face="arial,helvetica,sans-serif"><div># algunso datos</div><div>n <- 100</div><div>x <- rnorm(n)</div><div>y <- 1.4 + 1.2*x + 1.3*x^2+ rnorm(n)</div>
<div><br></div><div># lo que tienes</div><div>par(mfrow = c(1, 2))</div><div>plot(x, y)</div><div>fit <- lm(y ~ x + I(x^2))</div><div>points(x, p <- predict(fit), type = 'l', col = 2)</div><div><br></div><div>
# usando order()</div><div>o <- order(x, decreasing = TRUE)</div><div>plot(x, y)</div><div>points(x[o], p[o], type = 'l', col = 2)</div><div><br></div><div>Un saludo,</div><div>Jorge.-</div></font></font><br><div class="gmail_quote">
2012/7/27 Marcuzzi, Javier Rubén <span dir="ltr"><></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div dir="ltr">
<div style="font-size:12pt;font-family:'Calibri'">
<div>Buenas noches:</div>
<div> </div>
<div>Les consulto por lo siguiente, tengo este gráfico</div>
<div><img style="BACKGROUND-IMAGE:none;BORDER-BOTTOM:0px;BORDER-LEFT:0px;MARGIN:0px;PADDING-LEFT:0px;PADDING-RIGHT:0px;DISPLAY:inline;BORDER-TOP:0px;BORDER-RIGHT:0px;PADDING-TOP:0px" title="image" border="0" alt="image" src="cid:0F00BD74C2044E07A025725B24B97CA7@JavierMarcuzzi" width="644" height="325"></div>
<div>Estos son el resultado de modelos mixtos, los datos son almacenados en un
data.frame, resumiendo el código, quitando el x <- data.frame() como el
plot()</div>
<div> </div>
<div>El primero es fitted:</div>
<div>goatMilk$dim, fitted(fitWilm_goat)</div>
<div>plot la línea de arriba</div>
<div> </div>
<div>El segundo es predict:</div>
<div>goatMilk$dim, predict(fitWilm_goat)</div>
<div> </div>
<div>El tercero es plot + lines</div>
<div>goatMilk$dim, goatMilk$milk</div>
<div>plot a la línea de arriba +</div>
<div>lines(y=predict(fitWilm_goat), col='red', x=goatMilk$dim)</div>
<div> </div>
<div>Ahora me concentro el el último:</div>
<div><img style="BACKGROUND-IMAGE:none;BORDER-BOTTOM:0px;BORDER-LEFT:0px;MARGIN:0px;PADDING-LEFT:0px;PADDING-RIGHT:0px;DISPLAY:inline;BORDER-TOP:0px;BORDER-RIGHT:0px;PADDING-TOP:0px" title="image" border="0" alt="image" src="cid:8A48F420E8C14F0CBAD82BCA918A0B6C@JavierMarcuzzi" width="586" height="484"></div>
<div>Como se puede ver, es horrible. Tendrían que ser curvas como las anteriores
(que son una sucesión de puntos).</div>
<div>Pensé en tres cosas, la primera en que como los puntos son de varios
animales y estos están como aleatorios me crea una curva para cada uno, pero
esto me me explica (o no comprendo) porqué esas líneas “horizontales”.</div>
<div>Otro pensamiento fue ordenar los datos por dim:</div>
<div>goatMilk <- goatMilk[order(goatMilk$dim),]</div>
<div>Superponiendo ambos gráficos obtengo un mamarracho</div>
<div><img style="BACKGROUND-IMAGE:none;BORDER-BOTTOM:0px;BORDER-LEFT:0px;MARGIN:0px;PADDING-LEFT:0px;PADDING-RIGHT:0px;DISPLAY:inline;BORDER-TOP:0px;BORDER-RIGHT:0px;PADDING-TOP:0px" title="image" border="0" alt="image" src="cid:6B739B44B36A4EA4BB85C6AFBA7E7515@JavierMarcuzzi" width="554" height="439"></div>
<div>Lo tercero que pensé es que no se debe usar “lines” cuándo tengo factores
aleatorios, sin separarlos en varios gráficos.</div>
<div>¿pero si deseo graficar la curva teórica resultante del modelo?</div>
<div> </div>
<div>Resulta que a partir de la curva del modelo pienso usar integrales para
calcular valores bajo ciertos tramos. Pero me desorientan los resultados en
rojo, </div>
<div>¿es un error mío al usar lines dentro de plot?</div>
<div> </div>
<div>¿que usarían para graficar la línea (curva) sobre los puntos? (aquí puse
pocos datos, si uso todos la superposición de puntos no sirve). </div>
<div> </div>
<div>Ahora lo mismo pero:</div>
<div>Con pocos datos la sucesión de puntos crea una curva, si utilizo datos de
otra fuente, donde la cantidad es mayor pero de no muy buena calidad, no queda
como una curva presentable fruto de una sucesión de puntos. </div>
<div><img style="BACKGROUND-IMAGE:none;BORDER-BOTTOM:0px;BORDER-LEFT:0px;MARGIN:0px;PADDING-LEFT:0px;PADDING-RIGHT:0px;DISPLAY:inline;BORDER-TOP:0px;BORDER-RIGHT:0px;PADDING-TOP:0px" title="image" border="0" alt="image" src="cid:C250BCAE64D649E4B0E770C1E044F15A@JavierMarcuzzi" width="644" height="368"></div>
<div> </div>
<div>Aquí en realidad puedo hacer más de una curva o un modelo mixto con
parámetros comunes, donde un factor sería quien reemplazaría varias curvas. O
mejor dicho, me interesa un modelo que englobe a todos y poder realizar
integrales a la curva resultante de ese modelo.</div>
<div> </div>
<div>¿Cómo puedo graficar varias líneas dentro de esos puntos? ¿O una línea que
grafique el modelo?</div>
<div>Mirando las líneas rojas, al plantear la integral, ¿hay problemas o al
graficar mal no hay ningún inconveniente con algo como integral<- function(x)
{(ajust$coefficients["x"]) * x + ajust$coefficients["(Intercept)"]} ?</div>
<div> </div>
<div>¿Alguna sugerencia, idea, crítica, enseñanza, corrección?</div>
<div> </div>
<div>Eso rojo me desconcierta, el último gráfico si está todo bien me indica un
problema que no es de R.</div><span class="HOEnZb"><font color="#888888">
<div> </div>
<div>Javier Marcuzzi</div></font></span></div></div></div>
<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><br>
<br></blockquote></div><br></div>