R/psm.R

Defines functions psm

psm <- function(fit) {
  X <- as.matrix(fit$jagsfit$mcmc)
  Bait <- fit$Data$Bait
  multi <- fit$Data$multiBait
  if (multi) {
    dim(X) <- c(nrow(X), dim(Bait)[1], dim(Bait)[3])
    pm <- cbind(apply(X, 2:3, median), apply(X[,,2]-X[,,1], 2, median))
    ps <- cbind(apply(X, 2:3, mad), apply(X[,,2]-X[,,1], 2, mad))
    b <- dimnames(Bait)[[3]]
    colnames(pm) <- colnames(ps) <- c(b, paste(b[2], b[1], sep="_"))
    rownames(pm) <- rownames(ps) <- dimnames(Bait)[[1]]
  } else {
    pm <- apply(X, 2, median)
    ps <- apply(X, 2, mad)
    names(pm) <- names(ps) <- dimnames(Bait)[[1]]
  }
  structure(list(Data=fit$Data, pm=pm, ps=ps),
            class = "psm.deepn")
}
pbreheny/deepn documentation built on May 11, 2020, 3:49 a.m.