Title: | Easily Performs Differential Coexpression Analysis |
---|---|
Description: | A function that reads in the GEO code of a list of probesets or gene symbols, a gene expression dataset GEO accession code, the name of the dataset feature discriminating the two conditions for the differential coexpression, and the values of the two different conditions for the differential coexpression, and returns the significant pairs of genes/probesets with highest differential coexpression (p-value < 0.005). If the input gene list is made of gene symbols, this package associates the probesets to these gene symbols, if found. Platforms available: GPL80, GPL8300, GPL80, GPL96, GPL570, GPL571, GPL20115, GPL1293, GPL6102, GPL6104, GPL6883, GPL6884, GPL13497, GPL14550, GPL17077, GPL6480. GEO: Gene Expression Omnibus. ID: identifier code. The GEO datasets are downloaded from the URL <https://ftp.ncbi.nlm.nih.gov/geo/series/>. |
Authors: | Davide Chicco [aut, cre] |
Maintainer: | Davide Chicco <[email protected]> |
License: | GPL-3 |
Version: | 1.4 |
Built: | 2024-12-02 07:01:57 UTC |
Source: | CRAN |
Function that computes the differential coexpression of a list of probesets in a specific dataset and returns the most significant pairs
easyDifferentialGeneCoexpression( list_of_probesets_to_select, GSE_code, featureNameToDiscriminateConditions, firstConditionName, secondConditionName, batchCorrection = TRUE, verbose = FALSE )
easyDifferentialGeneCoexpression( list_of_probesets_to_select, GSE_code, featureNameToDiscriminateConditions, firstConditionName, secondConditionName, batchCorrection = TRUE, verbose = FALSE )
list_of_probesets_to_select |
list of probesets for which the differential coexpression should be computed |
GSE_code |
GEO accession code of the dataset to analyze |
featureNameToDiscriminateConditions |
name of the feature of the dataset that contains the two conditions to investigate |
firstConditionName |
name of the first condition in the feature to discriminate (for example, "healthy") |
secondConditionName |
name of the second condition in the feature to discriminate (for example, "cancer") |
batchCorrection |
says if the script should perform the batch correction with limma::removeBatchEffect() or not |
verbose |
prints all the intermediate message to standard output or not |
a dataframe containing the significantly differentially co-expressed pairs of genes
probesetList <- c("200738_s_at", "217356_s_at", "206686_at") verboseFlag <- "TRUE" batchCorrection <- "TRUE" signDiffCoexpressGenePairs <- easyDifferentialGeneCoexpression(probesetList, "GSE3268", "description", "Normal", "Tumor", verboseFlag)
probesetList <- c("200738_s_at", "217356_s_at", "206686_at") verboseFlag <- "TRUE" batchCorrection <- "TRUE" signDiffCoexpressGenePairs <- easyDifferentialGeneCoexpression(probesetList, "GSE3268", "description", "Normal", "Tumor", verboseFlag)
Function that associates a gene symbol to a probeset for some Affymetrix platforms
fromProbesetToGeneSymbol( thisProbeset, thisPlatform, this_platform_ann_df, verbose = FALSE )
fromProbesetToGeneSymbol( thisProbeset, thisPlatform, this_platform_ann_df, verbose = FALSE )
thisProbeset |
probeset in input |
thisPlatform |
GEO platform accession code |
this_platform_ann_df |
annotation dataframe of the platform |
verbose |
prints all the intermediate message to standard output or not |
a gene symbol as string
Function that downloads gene expression data from GEO, after checking the connection
geoDataDownload(GSE_code, verbose = FALSE)
geoDataDownload(GSE_code, verbose = FALSE)
GSE_code |
GEO code dataset |
verbose |
prints all the intermediate message to standard output or not |
a gene set gene expression AnnotationDataFrame
Function that downloads the annotations of a GEO platform
geoPlatformAnnotationsDownload(platformID, verbose = FALSE)
geoPlatformAnnotationsDownload(platformID, verbose = FALSE)
platformID |
GEO platform ID |
verbose |
prints all the intermediate message to standard output or not |
a dataframe containing the annotations of the GEO platform
Function that reads a CSV file of probesets or gene symbols and, in the latter case, it retrieves the original probesets
probesetRetrieval( probesets_or_gene_symbols, csv_file_name, platformCode, verbose = FALSE )
probesetRetrieval( probesets_or_gene_symbols, csv_file_name, platformCode, verbose = FALSE )
probesets_or_gene_symbols |
flag saying if we're reading probesets or gene symbols |
csv_file_name |
complete name of CSV file containing the probesets or the gene symbols |
platformCode |
code of the microarray platform for which the probeset-gene symbol mapping should be done |
verbose |
prints all the intermediate message to standard output or not |
a vector of probesets