<HTML><HEAD></HEAD>
<BODY dir=ltr>
<DIV dir=ltr>
<DIV style="FONT-FAMILY: 'Calibri'; COLOR: #000000; FONT-SIZE: 12pt">
<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>
<DIV> </DIV>
<DIV>Javier Marcuzzi</DIV></DIV></DIV></BODY></HTML>