TROM (version 1.0)

select.associated.orthologs: Select the associated genes for biological samples of a species among the genes with orthologs in the other species

Description

select.associated.orthologs() finds the ids of associated genes of a species within orthologous genes based on criterion: Z-scores $\ge$ z_thre and FPKM $\ge$ FPKM_thre and saves the results to an .xlsx file.

Usage

select.associated.orthologs(sp_gene_expr, sp1_sp2_orthologs, 
z_thre = 1.5, FPKM_thre = 1, i, save = TRUE, plot_distribution = FALSE)

Arguments

sp_gene_expr

a data frame containing gene expression estimates of the species; rows correspond to genes; columns (from the second to the last) correspond to samples, with the first column as gene IDs.

sp1_sp2_orthologs

a data frame containing ortholog gene pairs between species 1 and 2.

z_thre

a numeric value specifiying the Z-score threshold used to select associated genes, whose Z-scores $\ge$ z_thre. Defaults to 1.5.

FPKM_thre

a numeric value specifying FPKM threshold used to select associated genes, whose FPKM $\ge$ FPKM_thre. Defaults to 1.

i

an integer specifying which column of sp1_sp2_orthologs the species corresponds to. 1 for the first column and 2 for the second column.

save

a Boolean value specifying whether to save the associated orthologs to an Excel file. If save = TRUE, this function saves the results to an Excel file named "associated genes within ortholog genes.xlsx".

plot_distribution

a Boolean value specifying whether to output the distribution of the number of associated orthologs across different samples. If plot_distribution = TRUE, this function outputs a barplot of the number of associated orthologs for each sample i

Value

• A data frame containing the associated genes within orthologous genes of each samples from the specified species. Every column in the data frame stores the associated gene IDs for the corresponding sample.

References

Li JJ., Huang H, Bickel PJ, & Brenner SE (2014). Comparison of D. melanogaster and C. elegans developmental stages, tissues, and cells by modENCODE RNA-seq data. Genome Research, 24(7), 1086-1101.

See Also

select.associated.genes

Examples

## Find the associaed genes of D. melanogaster samples
## among the D. melanogaster genes having orthologs in C.elegans

## dm_gene_expr.rda and dm_ce_orthologs.rda can be downloaded and unzipped from
## http://www.stat.ucla.edu/~jingyi.li/packages/TROM/TROM_Rdata.zip.
load("dm_gene_expr.rda")
load("dm_ce_orthologs.rda")
dm_associated_orthologs <- select.associated.orthologs(sp_gene_expr = dm_gene_expr, 
                           sp1_sp2_orthologs = dm_ce_orthologs, z_thre = 1.5, 
                           FPKM_thre = 1, 1, save = TRUE, plot_distribution = TRUE)