Hola.
Tengo dos bases de datos que valoran las veces que rumia una vaca con calor (Q1) y sin calor (Q4).
Genero el mismo modelo sinoidal:
mod.RumiQ1 <-nls(Rumination ~ a + b*sin(c*Hour-d)+e*cos(f*Hour-g),data=datosQ1,start=list(a =23.5, b = 6, c=3.1, d=9, e=-2.5, f=0.8,g=1))
co<-coef(mod.RumiQ1 )
co
#co
# a b c d e f g
#23.13338 -4.33017 0.33730 97.53357 -5.14952 0.50374 -1.82722
Y lo dibujo con :
ggplot(datosQ1, aes(x=Hour, y=Rumination))+
# geom_point()+
geom_smooth(method="nls", se=FALSE,
method.args=list(formula=y~ a + b*sin(c*x-d)+e*cos(f*x-g),
start=c(a=23.13338,b= -4.33017,c= 0.33730,d= 97.53357,e= -5.14952, f=0.50374, g= -1.82722 )) )
Ahora utilizo los datos Q4 :
mod.RumiQ4 <-nls(Rumination ~ a + b*sin(c*Hour-d)+e*cos(f*Hour-g),data=datosQ4,start=list(a =22, b = 5, c=3, d=9, e=-3.3, f=1,g=1.5))
co<-coef(mod.RumiQ4 )
co
# a b c d e f g
#22.03420 0.88497 1.89659 -11.66763 -3.35122 1.08926 3.51036
Y lo dibujo con:
ggplot(datosQ4,
aes(x=Hour, y=Rumination))+
# geom_point()+
geom_smooth(method="nls", se=FALSE,
method.args=list(formula=y~ a + b*sin(c*x-d)+e*cos(f*x-g),
start=c(a=22.03,b= 0.88,c= 1.89,d= -11,e= -3, f=1.00, g= 3.51)) )
Me gustaría tener las dos curvas en el mismo dibujo para poder comparar las gráficas
¿Cómo lo hago?
Gracias !!!!
Foro de debate
Unir gráficas de dos curvas nls
This forum has a limit to the number of forum postings you can make in a given time period - this is currently set at 5 posting(s) in 1 day