data.prep: Pre-process and re-structurize the input data

In JingXieMIZZOU/hmmASE: Detection of Differentially Allelic Expressed Regions (DAERs) by a Bayesian Hidden Markov Model

Description

A function to pre-process the raw input data by applying the Haldane-Anscombe correction, and then convert the raw data into a structure which contains necessary information for analysis, i.e., Chromosome, GeneID, Gene_name, genomic position, and the logit transformation of the ASE ratios for normal (O1) and abnormal (O2) groups. This function can be easily modified to apply to other similar real data structures.

Usage

data.prep(norm.M, norm.P, abnorm.M, abnorm.P, rep)

Arguments

norm.M	A dataframe contains the read counts for maternal allele in normal group. The dimension of this dataframe can vary with the number of available biological replicates but the first 4 columns should be Chromosome, GeneID, Gene_name, and Position.
norm.P	A dataframe contains the read counts for paternal allele in normal group, with similar structure with norm.M.
abnorm.M	A dataframe contains the read counts for maternal allele in abnormal group, with similar structure with norm.M and norm.P.
abnorm.P	A dataframe contains the read counts for paternal allele in abnormal group, with similar structure with other three input datasets.
rep	A vector indicates the number of biological replicates for normal and abnormal groups.

Value

A data.frame in the structure that would be needed for analysis by hmmASE method.

JingXieMIZZOU/hmmASE documentation built on Dec. 31, 2019, 12:52 a.m.