Nothing
testOutliers.profilet.metaplus <- function(object,R=999, cores) {
meta.fun <- function(data) {
isreg <- (dim(data)[2]>2)
if (isreg) betadata.meta <- metaplus(yi=data[,1], sei=data[,2], mods=data[,3:dim(data)[2],drop=FALSE],justfit=TRUE,plotci=FALSE,slab=NULL)
else betadata.meta <- metaplus(yi=data[,1], sei=data[,2],justfit=TRUE,plotci=FALSE,slab=NULL)
if (isreg) betadata.meta2 <- metaplus(yi=data[,1], sei=data[,2], mods=data[,3:dim(data)[2],drop=FALSE],justfit=TRUE,plotci=FALSE,slab=NULL,
random="t-dist")
else betadata.meta2 <- metaplus(yi=data[,1], sei=data[,2],justfit=TRUE,plotci=FALSE,slab=NULL,random="t-dist")
if (abs((logLik(betadata.meta2)-logLik(betadata.meta))/logLik(betadata.meta))<1.0e-4) sim.chisq <- 0
else sim.chisq <- 2*(logLik(betadata.meta2)-logLik(betadata.meta))
if (sim.chisq < 0) {
warning(paste("t-dist chisq< 0","is",sim.chisq))
sim.chisq <- 0
}
return(sim.chisq)
}
meta.rg <- function(data, mle) {
isreg <- (dim(data)[2]>2)
data[,2] <- sample(data[,2], length(data[,2]), replace = TRUE)
if (isreg) mods <- data[,3:dim(data)[2],drop=FALSE]
if (isreg) data[,1] <- rnorm(dim(data)[1],mean=mle$b[1,1],sd=sqrt(mle$tau2))+
as.matrix(mods) %*% mle$b[2:dim(mle$b)[1],1]+
rnorm(dim(data)[1],mean=0,sd=data[,2])
else data[,1] <- rnorm(dim(data)[1],mean=mle$b[1,1],sd=sqrt(mle$tau2))+
rnorm(dim(data)[1],mean=0,sd=data[,2])
return(data)
}
if (object$justfit) stop("Cannot use with objects fitted with justfit=TRUE")
if (is.null(object$mods)) meta.ml <- rma(yi=object$yi, sei=object$sei, method="ML")
else meta.ml <- rma(yi=object$yi, sei=object$sei, mods=object$mods, method="ML")
if (cores>1) {
if(.Platform$OS.type=="unix") parallel <- "multicore"
else parallel <- "snow"
} else parallel <- "no"
meta.boot <- boot(data=cbind(object$yi,object$sei,object$mods), meta.fun,
R = R, sim = "parametric", ran.gen = meta.rg, mle=meta.ml, parallel=parallel,
ncpus = cores)
pvalue <- (sum(meta.boot$t[,1] >= meta.boot$t0[1],na.rm=TRUE)+1)/(1+sum(!is.nan(meta.boot$t[,1])))
return(list(pvalue=pvalue,observed=meta.boot$t0[1],sims=meta.boot$t[,1]))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.