# NOT RUN {
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=logit(0.8))
HatchingSuccesss.lnL(par=par, data=totalIncubation_Cc)
g <- HatchingSuccesss.fit(par=par, data=totalIncubation_Cc)
HatchingSuccesss.lnL(par=g$par, data=totalIncubation_Cc)
plot(g)
# }
Run the code above in your browser using DataLab