select.associated.orthologs: Select the associated genes for biological samples of a...
In TROM: Transcriptome Overlap Measure

Description Usage Arguments Value Author(s) References See Also Examples

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

1
2
3

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

`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 ≥ `z_thre`. Defaults to 1.5. This can be specified by users or calculated using `choose.z()`.
`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 in a pdf file named "number of sample associated orthologous genes.pdf".

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.

Jingyi Jessica Li, Wei Vivian Li

Li WV, Chen Y and Li JJ (2016). TROM: A Testing-Based Method for Finding Transcriptomic Similarity of Biological Samples. Statistics in Biosciences. DOI: 10.1007/s12561-016-9163-y

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.

select.associated.genes

## 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.
## Not run: 
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, 
                           i = 1, save = TRUE, plot_distribution = TRUE)
## End(Not run)