Foro de discusión y soporte para usuarios de R

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 !!!!