phenology (version 4.2.4)

fitRMU_MHmcmc_p: Generates set of parameters to be used with fitRMU_MHmcmc()

Description

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

Usage

fitRMU_MHmcmc_p(result = stop("An output from fitRMU() must be provided"),
  accept = FALSE)

Arguments

result

An object obtained after a fitRMU() fit

accept

If TRUE, does not wait for use interaction

Value

• A matrix with the parameters

Details

fitRMU_MHmcmc_p generates set of parameters to be used with fitRMU_MHmcmc()

See Also

Other Fill gaps in RMU: fitRMU_MHmcmc; fitRMU; logLik.fitRMU; plot.fitRMU

Examples

library("phenology")
RMU.name.AtlanticW <- data.frame(mean=c("Yalimapo.French.Guiana",
                                         "Galibi.Suriname",
                                         "Irakumpapy.French.Guiana"),
                                 se=c("se_Yalimapo.French.Guiana",
                                      "se_Galibi.Suriname",
                                      "se_Irakumpapy.French.Guiana"))
data.AtlanticW <- data.frame(Year=c(1990:2000),
      Yalimapo.French.Guiana=c(2076, 2765, 2890, 2678, NA,
                               6542, 5678, 1243, NA, 1566, 1566),
      se_Yalimapo.French.Guiana=c(123.2, 27.7, 62.5, 126, NA,
                                 230, 129, 167, NA, 145, 20),
      Galibi.Suriname=c(276, 275, 290, NA, 267,
                       542, 678, NA, 243, 156, 123),
      se_Galibi.Suriname=c(22.3, 34.2, 23.2, NA, 23.2,
                           4.3, 2.3, NA, 10.3, 10.1, 8.9),
      Irakumpapy.French.Guiana=c(1076, 1765, 1390, 1678, NA,
                               3542, 2678, 243, NA, 566, 566),
      se_Irakumpapy.French.Guiana=c(23.2, 29.7, 22.5, 226, NA,
                                 130, 29, 67, NA, 15, 20))

cst <- fitRMU(data=data.AtlanticW, RMU.name=RMU.name.AtlanticW,
               colname.year="Year", model.trend="Constant",
               model.SD="Zero")
pMCMC <- fitRMU_MHmcmc_p(result=cst, accept=TRUE)