sim1: A set of simulated data

In yunzhang813/PBtest-R-Package: Highly efficient hypothesis testing methods for regression-type tests with correlated observations and heterogeneous variance structure

Description

Simulated data of 2000 genes (1000 null genes and 1000 differentially expressed genes) and 40 samples. Paired samples are collected from 20 subjects. These samples come from two groups. Subject IDs and weights are available.

Usage

data("sim1")

Format

The object sim1 is a list with the following elements:

simdat

Data matrix with samples in rows and genes in columns. Column names indicate whether the gene is a null gene (e.g. null.G1) or a differentially expressed gene (e.g. sig.G1).

Grp

Vector of group labels.

ww

Vector of weights.

Subj

Vector of subject IDs, which are also the row names of simdat.

References

Zhang et al. (2019) Highly efficient hypothesis testing methods for regression-type tests with correlated observations and heterogeneous variance structure. BMC Bioinformatics, 20:185.

Examples

library(PBtest)
data(sim1)

## attach list elements to the global environment
list2env(sim1,globalenv())

## the PB-transformed t-test, with Sigma estimation
pvec.PB <- PBtest(YY=simdat, xx=Grp, test="t", id=Subj, weights=ww)$p.value

yunzhang813/PBtest-R-Package documentation built on March 18, 2020, 5:29 p.m.