fitted.rma <- function(object, ...) {
mstyle <- .get.mstyle()
.chkclass(class(object), must="rma")
na.act <- getOption("na.action")
if (!is.element(na.act, c("na.omit", "na.exclude", "na.fail", "na.pass")))
stop(mstyle$stop("Unknown 'na.action' specified under options()."))
### note: fitted values can be calculated for all studies including those that
### have NA on yi/vi (and with "na.pass" these will be provided); but if there
### is an NA in the X's, then the fitted value will also be NA
out <- c(object$X.f %*% object$beta)
names(out) <- object$slab
#not.na <- !is.na(out)
if (na.act == "na.omit")
out <- out[object$not.na]
if (na.act == "na.exclude")
out[!object$not.na] <- NA_real_
if (na.act == "na.fail" && any(!object$not.na))
stop(mstyle$stop("Missing values in results."))
if (inherits(object, "rma.ls")) {
out <- list(location = out)
out$scale <- c(object$Z.f %*% object$alpha)
names(out$scale) <- object$slab
#not.na <- !is.na(out$scale)
if (na.act == "na.omit")
out$scale <- out$scale[object$not.na]
if (na.act == "na.exclude")
out$scale[!object$not.na] <- NA_real_
if (na.act == "na.fail" && any(!object$not.na))
stop(mstyle$stop("Missing values in results."))
}
return(out)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.