0

このテーマに関する情報はほとんど見つかりませんが、ピアソン 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