if (FALSE) {
library(embryogrowth)
totalIncubation_Cc <- subset(DatabaseTSD,
Species=="Caretta caretta" &
Note != "Sinusoidal pattern" &
!is.na(Total) & Total != 0)
par <- c(S.low=0.5, S.high=0.3,
P.low=25, deltaP=10, MaxHS=0.8)
HatchingSuccess.lnL(par=par, data=totalIncubation_Cc)
g <- HatchingSuccess.fit(par=par, data=totalIncubation_Cc)
HatchingSuccess.lnL(par=g$par, data=totalIncubation_Cc)
plot(g)
nobs(g)
}
Run the code above in your browser using DataLab