0

我几乎找不到关于这个主题的信息，但我希望将 Pearson I 曲线拟合到生育力数据。下面的数据和解决方案来自一本书（生物学家应用人口学，Carey 1993）。我希望使用最大似然来自动化这个过程，然后应用拟合优度测试。目前，它是一个纯粹的视觉/试错练习，将曲线拟合到数据中。我在使用 PearsonDS 软件包或其他方法（例如 fitdist）方面没有取得成功。

#par(mar=c(5,5,3,4),cex=1.2)
###taken from Carey 1993,egg/fem/day
egg=c(0,0,0,2.32,17.82,36.34,80.77,68.49,75.00,82.12,75.16,65.20,56.04,43.28,53.52,49.50,27.85,37.62,28.53,30.34,29.57,48.17,27.25,25.37,28.66,13.00,17.13,8,12.70,10.25,4)


##Vector for days,x
d=seq(from=0,to=length(med)-1)

###Parameters for the Pearson I function
a1=.55
a2=34
m1=1.46
m2=3

##The Pearson function  
y=2.8
pearI=function(pearson)
    {
    y*((1+d/a1)^m1)*((1-d/a2)^m2)
    }   

pr2=(pearI(egg))    
##
plot(egg,main=expression(bold(
           paste("Pearson Type I function fitted to ", M[x]))), 
     xlab="Age (days)",ylab="Eggs/fem", pch=16,cex.lab=1.4)
#lines(egg,col="red")
lines(pr2,type="b",col="blue",pch=16, cex=1.3)
4

0 回答 0