inst/doc/HDTD.R

In HDTD: Statistical Inference about the Mean Matrix and the Covariance Matrices in High-Dimensional Transposable Data (HDTD)

## ---- echo=FALSE--------------------------------------------------------------
    knitr::opts_chunk$set(
    tidy = TRUE,
    collapse = TRUE, 
    comment = "#>")

## -----------------------------------------------------------------------------
library("HDTD")
data("VEGFmouse")
dim(VEGFmouse)
rownames(VEGFmouse)

## -----------------------------------------------------------------------------
colnames(VEGFmouse)[1:18]

## -----------------------------------------------------------------------------
meanmat.ts(VEGFmouse, N = 40, group.sizes = 9)

## -----------------------------------------------------------------------------
colnames(VEGFmouse)[1:9]
columnorder <- c(1, 4, 5, 8, 2, 3, 6, 7, 9)
VEGForder <- orderdata(VEGFmouse, N = 40, order.cols = columnorder) 
colnames(VEGForder)[1:9]

## -----------------------------------------------------------------------------
meanmat.ts(VEGForder, N = 40, group.sizes = c(4, 1, 1, 1, 1, 1))

## -----------------------------------------------------------------------------
sample_mean <- meanmat.hat(VEGFmouse, N = 40)
sample_mean

## -----------------------------------------------------------------------------
head(round(sample_mean$estmeanmat, 4), n = 10)

## -----------------------------------------------------------------------------
est_cov_mat <- covmat.hat(datamat = VEGFmouse, N = 40)
est_cov_mat

## -----------------------------------------------------------------------------
round(est_cov_mat$cols.covmat, 3)

## -----------------------------------------------------------------------------
covmat.ts(datamat = VEGFmouse, N = 40, voi = "columns")

## -----------------------------------------------------------------------------
print(citation("HDTD"), bibtex = TRUE)