HatchingSuccess.MHmcmc_p: Generates set of parameters to be used with...

In embryogrowth: Tools to Analyze the Thermal Reaction Norm of Embryo Growth

Generates set of parameters to be used with HatchingSuccess.MHmcmc()

Description

Interactive script used to generate set of parameters to be used with HatchingSuccess.MHmcmc().

Usage

HatchingSuccess.MHmcmc_p(
  result = NULL,
  parameters = NULL,
  fixed.parameters = NULL,
  accept = FALSE
)

Arguments

result	An object obtained after a HatchingSuccess.fit() fit
parameters	A set of parameters. Replace the one from result
fixed.parameters	A set of fixed parameters. Replace the one from result
accept	If TRUE, the script does not wait user information

Details

HatchingSuccess.MHmcmc_p generates set of parameters to be used with HatchingSuccess.MHmcmc()

Value

A matrix with the parameters

Author(s)

Marc Girondot

See Also

Other Hatching success: HatchingSuccess.MHmcmc(), HatchingSuccess.fit(), HatchingSuccess.lnL(), HatchingSuccess.model(), logLik.HatchingSuccess(), nobs.HatchingSuccess(), predict.HatchingSuccess()

Examples

## 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)
         
g <- HatchingSuccess.fit(par=par, data=totalIncubation_Cc)
pMCMC <- HatchingSuccess.MHmcmc_p(g, accept=TRUE)
mcmc <- HatchingSuccess.MHmcmc(result=g, parameters = pMCMC, 
                  adaptive=TRUE, n.iter=100000, trace=1000)

## End(Not run)

embryogrowth documentation built on Sept. 11, 2024, 8:16 p.m.