https://github.com/cran/ensembleBMA
Tip revision: 65bc4d8d72e5060987d37deb8279042bbcd033ff authored by Chris Fraley on 08 May 2008, 00:00:00 UTC
version 3.0-2
version 3.0-2
Tip revision: 65bc4d8
matchEnsembleMembers.fitBMAgamma0.R
`matchEnsembleMembers.fitBMAgamma0` <-
function( fit, ensembleData)
{
# match ensemble members in ensembleData and fit
if (!is.null(dim(fit$weights))) {
fitMems <- dimnames(fit$weights)[[1]]
}
else {
fitMems <- names(fit$weights)
}
ensMems <- ensembleMemberLabels(ensembleData)
if (!is.null(fitMems) && !is.null(ensMems)
&& length(fitMems) > length(ensMems))
stop("model fit has more ensemble members than ensemble data")
WARN <- rep(FALSE,3)
WARN[1] <- is.null(fitMems) && !is.null(ensMems)
WARN[2] <- !is.null(fitMems) && is.null(ensMems)
WARN[3] <- is.null(fitMems) && is.null(ensMems)
if (any(WARN) && length(fitMems) != length(ensMems))
stop("model fit and ensemble data differ in ensemble size")
if (any(WARN))
warning("cannot check correspondence between model fit and ensemble data members")
M <- match(fitMems, ensMems, nomatch = 0)
if (any(!M)) stop("ensembleData is missing a member used in fit")
}