Automated build.

.Rbuildignore

@@ -11,6 +11,7 @@
 ^\.git.*$
 ^\.library$
 ^\.setup$
+^\.sif$
 ^\.sim$
 ^detritus$
 ^index\.qmd$

.setup/latex/pdf/.gitignore

@@ -0,0 +1,5 @@
+*
+*/
+!*.pdf
+!.gitignore
+!bib.bib

NAMESPACE

@@ -13,6 +13,7 @@ S3method(print,ctmedmed)
 S3method(summary,ctmeddelta)
 S3method(summary,ctmedmc)
 S3method(summary,ctmedmed)
+S3method(summary,ctmedposteriorphi)
 export(DeltaBeta)
 export(DeltaIndirectCentral)
 export(DeltaMed)

R/cTMed-methods.R

@@ -508,6 +508,63 @@ summary.ctmedmed <- function(object,
   }
 }
 
+#' Summary Method for Object of Class `ctmedposteriorphi`
+#'
+#' @author Ivan Jacob Agaloos Pesigan
+#' @param object an object of class `ctmedposteriorphi`.
+#' @param ... further arguments.
+#' @return Returns a list of the posterior means
+#'   (in matrix form) and covariance matrix.
+#'
+#' @examples
+#' set.seed(42)
+#' phi <- matrix(
+#'   data = c(
+#'     -0.357, 0.771, -0.450,
+#'     0.0, -0.511, 0.729,
+#'     0, 0, -0.693
+#'   ),
+#'   nrow = 3
+#' )
+#' colnames(phi) <- rownames(phi) <- c("x", "m", "y")
+#' mc <- MCPhi(
+#'   phi = phi,
+#'   vcov_phi_vec = 0.1 * diag(9),
+#'   R = 100L # use a large value for R in actual research
+#' )$output
+#' summary(mc)
+#'
+#' @keywords methods
+#' @export
+summary.ctmedposteriorphi <- function(object,
+                                      ...) {
+  varnames <- colnames(object[[1]])
+  p <- dim(object[[1]])[1]
+  x <- do.call(
+    what = "rbind",
+    args = lapply(
+      X = object,
+      FUN = function(x) {
+        dim(x) <- NULL
+        return(x)
+      }
+    )
+  )
+  phi_vec <- colMeans(x)
+  phi <- matrix(
+    data = phi_vec,
+    nrow = p
+  )
+  colnames(phi) <- rownames(phi) <- varnames
+  vcov_phi_vec <- stats::var(x)
+  return(
+    list(
+      phi = phi,
+      vcov_phi_vec = vcov_phi_vec
+    )
+  )
+}
+
 #' Summary Method for an Object of Class `ctmeddelta`
 #'
 #' @author Ivan Jacob Agaloos Pesigan

R/cTMed-posterior-phi.R

@@ -21,14 +21,17 @@ PosteriorPhi <- function(object) {
   )
   varnames <- object$ctstanmodelbase$latentNames
   posterior <- ctsem::ctExtract(object)$pop_DRIFT
-  return(
-    lapply(
-      X = seq_len(dim(posterior)[1]),
-      FUN = function(i) {
-        phi <- posterior[i, , ]
-        colnames(phi) <- rownames(phi) <- varnames
-        return(phi)
-      }
-    )
+  out <- lapply(
+    X = seq_len(dim(posterior)[1]),
+    FUN = function(i) {
+      phi <- posterior[i, , ]
+      colnames(phi) <- rownames(phi) <- varnames
+      return(phi)
+    }
+  )
+  class(out) <- c(
+    "ctmedposteriorphi",
+    class(out)
   )
+  return(out)
 }

README.md

@@ -1,7 +1,7 @@
 cTMed
 ================
 Ivan Jacob Agaloos Pesigan
-2024-05-03
+2024-05-10
 
 <!-- README.md is generated from README.Rmd. Please edit that file -->
 <!-- badges: start -->