r - R中绘图的最佳拟合曲线

我在一个名为 ph 的图中有一个概率密度函数，我通过 stackoverflow 用户的帮助从两个数据样本中导出，以这种方式

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

我想拟合 ph 曲线的最佳曲线，但找不到工作方法。
我怎样才能做到这一点？我必须从 ph 中提取 vaule 然后对它们进行处理？或者有相同的功能

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

直接地？

最佳答案

假设你的意思是你想对 ph 中的数据执行曲线拟合，那么沿着nls(FUN, cbind(ph$counts, ph$mids),...) 可能有效。您需要知道您认为直方图数据应该适合哪种函数“FUN”，例如正态分布。阅读 nls() 上的帮助文件，了解如何为 FUN 中的系数设置起始“猜测”值。

如果你只是想在直方图上叠加一条曲线，那么 smoo<-spline(ph$mids,ph$counts);lines(smoo$x,smoo$y)
将接近做到这一点。您可能需要调整 x 和/或 y 缩放比例。

关于r - R中绘图的最佳拟合曲线，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/7502027/