首页 文章

R中曲线的最佳拟合曲线

提问于
浏览 次
3

我在一个名为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")

直？

r curve-fitting

2 回答

  • 2

    假设您的意思是要对ph中的数据执行曲线拟合，那么 nls(FUN, cbind(ph$counts, ph$mids),...) 中的某些内容可能会起作用 . 您需要知道您认为直方图数据应该适合的函数类型，例如：正态分布 . 阅读 nls() 上的帮助文件，了解如何为FUN中的系数设置起始"guess"值 .

    如果您只想将曲线叠加到直方图上，那么 smoo<-spline(ph$mids,ph$counts); lines(smoo$x,smoo$y)

    将接近这样做 . 您可能必须调整x和/或y缩放 .

    回复于
  • 3

    你想要密度函数吗？

    x = rnorm(1000)
hist(x, breaks = 30, freq = FALSE)
lines(density(x), col = "red")
    回复于

相关问题