## ---- 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)