Tengo una función de densidad de probabilidad en un gráfico llamado ph que deriva de dos muestras de datos, con la ayuda de un usuario de stackoverflow, de esta manera

 few <-read.table('outcome.dat',head=TRUE)
 many<-read.table('alldata.dat',head=TRUE)
 mh <- hist(many$G,breaks=seq(0,1.,by=0.03), plot=FALSE)
 fh <- hist(few$G, breaks=mh$breaks, plot=FALSE)
 ph <- fh
 ph$density <- fh$counts/(mh$counts+0.001)
 plot(ph,freq=FALSE,col="blue")

Me gustaría ajustar la mejor curva de la gráfica de ph, pero no puedo encontrar un método de trabajo. ¿Cómo puedo hacer esto? Tengo que extraer el vaule de ph y luego trabaja en ellos? o hay la misma función que funciona en

 plot(ph,freq=FALSE,col="blue")

¿directamente