View source: R/fitRMU_MHmcmc_p.R
fitRMU_MHmcmc_p | R Documentation |
Interactive or automatic script used to generate set of parameters to be
used with fitRMU_MHmcmc().
If density="dgamma" is used, a uniform distribution is used for r, as
r can be negative.
fitRMU_MHmcmc_p(
result = stop("An output from fitRMU() must be provided"),
density = "dunif",
accept = FALSE
)
result |
An object obtained after a fitRMU() fit |
density |
Preset of density; can be dnorm, dunif, or dgamma |
accept |
If TRUE, does not wait for user interaction |
fitRMU_MHmcmc_p generates set of parameters to be used with fitRMU_MHmcmc()
A matrix with the parameters
Marc Girondot
Other Fill gaps in RMU:
CI.RMU()
,
fitRMU()
,
fitRMU_MHmcmc()
,
logLik.fitRMU()
,
plot.fitRMU()
## Not run:
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)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.