quickCrossCorr: quickCrossCorr

In Krutik6/TimiRGeN: Time sensitive microRNA-mRNA integration, analysis and network generation tool

quickCrossCorr

Description

Plots a cross-correlation plot to compare the miRNA and mRNA of a selected pair. This is a useful test of the similarities between the the two time series. It tracks movement of two time series relative to one another to determine how well they match and at which point the best match occurs.

Usage

quickCrossCorr(filt_df, pair, miRNA_exp, mRNA_exp, scale,
Interpolation, timecourse)

Arguments

filt_df	Dataframe from the matrixFilter function.
pair	Interger representing the pair to be explored.
miRNA_exp	miRNA data from using the diffExpressRes function on miRNA data.
mRNA_exp	mRNA data from using the diffExpressRes function on miRNA data
scale	TRUE or FALSE. Should data be scales. Default is FALSE. If the correlation is based on Log2FC values scale should be TRUE.
Interpolation	TRUE or FALSE. Should the whole time course be interpolated over by a smooth spline? Default is FALSE. This is most useful for longer and regular time courses.
timecourse	If Iterpolation is TRUE, how many time points should be interpolated over?

Value

A cross correlation plot.

Examples

library(org.Mm.eg.db)

miR <- mm_miR[1:50,]

mRNA <- mm_mRNA[1:100,]

MAE <- startObject(miR = miR, mRNA = mRNA)

MAE <- getIdsMir(MAE, assay(MAE, 1), orgDB = org.Mm.eg.db, 'mmu')

MAE <- getIdsMrna(MAE, assay(MAE, 2), "useast", 'mmusculus', orgDB = org.Mm.eg.db)

MAE <- diffExpressRes(MAE, df = assay(MAE, 1), dataType = 'Log2FC',
                     genes_ID = assay(MAE, 3),
                      idColumn = 'GENENAME',
                      name = "miRNA_log2fc")

MAE <- diffExpressRes(MAE, df = assay(MAE, 2), dataType = 'Log2FC',
                     genes_ID = assay(MAE, 7),
                     idColumn = 'GENENAME',
                     name = "mRNA_log2fc")

Filt_df <- data.frame(row.names = c("mmu-miR-145a-3p:Adamts15",
                                   "mmu-miR-146a-5p:Acy1"),
                     corr = c(-0.9191653, 0.7826041),
                     miR = c("mmu-miR-145a-3p", "mmu-miR-146a-5p"),
                     mRNA = c("Adamts15", "Acy1"),
                     miR_Entrez = c(387163, NA),
                     mRNA_Entrez = c(235130, 109652),
                     TargetScan = c(1, 0),
                     miRDB = c(0, 0),
                     Predicted_Interactions = c(1, 0),
                     miRTarBase = c(0, 1),
                     Pred_Fun = c(1, 1))

MAE <- matrixFilter(MAE, miningMatrix = Filt_df, negativeOnly = FALSE,
                   threshold = 1, predictedOnly = FALSE)

quickCrossCorr(filt_df=MAE[[11]], pair=1, miRNA_exp=MAE[[9]],
              mRNA_exp=MAE[[10]],scale = FALSE, Interpolation = FALSE)

Krutik6/TimiRGeN documentation built on Jan. 27, 2024, 7:46 p.m.