Tenho uma função de densidade de probabilidade em um gráfico chamado ph, derivado de duas amostras de dados, com a ajuda de um usuário do stackoverflow, dessa maneira

 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")

Gostaria de ajustar a melhor curva do gráfico de ph, mas não consigo encontrar um método de trabalho. como posso fazer isso? Eu tenho que extrair o vaule de ph e depois trabalha com eles? ou existe a mesma função que funciona em

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

diretamente