View source: R/HatchingSuccess.logLik.R
logLik.HatchingSuccess | R Documentation |
Set of functions to study the hatching success.
## S3 method for class 'HatchingSuccess'
logLik(object, ...)
object |
The return of a fit done with fitHS. |
... |
Not used |
logLik.HatchingSuccess returns -log L of a fit
Return -log L of a fit
Marc Girondot
Other Hatching success:
HatchingSuccess.MHmcmc()
,
HatchingSuccess.MHmcmc_p()
,
HatchingSuccess.fit()
,
HatchingSuccess.lnL()
,
HatchingSuccess.model()
,
nobs.HatchingSuccess()
,
predict.HatchingSuccess()
## 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=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)
t <- seq(from=20, to=40, by=0.1)
CIq <- predict(g, temperature=t)
par(mar=c(4, 4, 1, 1), +0.4)
plot(g)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.