This will be a very minimal analysis until we get some replicates.
I am using mm38_100.
## My load_biomart_annotations() function defaults to human, so that will be quick.
load_biomart_annotations(species="mmusculus") mm_annot <-
## The biomart annotations file already exists, loading from it.
mm_annot[["annotation"]]
mm_annot <-"txid"]] <- paste0(mm_annot[["ensembl_transcript_id"]], ".", mm_annot[["version"]])
mm_annot[[rownames(mm_annot) <- make.names(mm_annot[["ensembl_gene_id"]], unique=TRUE)
mm_annot[, c("txid", "ensembl_gene_id")] tx_gene_map <-
So, I now have a table of mouse annotations.
I am going to write a quick sample sheet in the current working directory called ‘all_samples.xlsx’ and put the names of the count tables in it.
Here I combine the metadata, count data, and annotations.
It is worth noting that the gene IDs from htseq-count probably do not match the annotations retrieved because they are likely exon-based rather than gene based. This is not really a problem, but don’t forget it!
mm_annot
hisat_annot <-rownames(hisat_annot) <- paste0("gene:", rownames(hisat_annot))
create_expt("sample_sheets/all_samples.xlsx",
mm38_hisat <-gene_info=hisat_annot)
## Reading the sample metadata.
## Dropped 1 rows from the sample metadata because they were blank.
## The sample definitions comprises: 23 rows(samples) and 24 columns(metadata fields).
## Reading count tables.
## Reading count files with read.table().
## /mnt/sshfs/cbcbsub01/fs/cbcb-lab/nelsayed/scratch/atb/rnaseq/mmusculus_iprgc_2019/preprocessing/iprgc_01/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows.
## preprocessing/iprgc_02/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/iprgc_03/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/iprgc_04/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/iprgc_05/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/iprgc_06/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/iprgc_07/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/iprgc_08/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/01_wt_retina/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/02_wt_scn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/03_wt_dlgn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/04_ko_retina/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/05_ko_scn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/06_ko_dlgn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/07_het1_retina/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/08_het2_retina/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/09_het3_retina/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/10_het1_scn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/11_het2_scn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/12_het3_scn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/13_het1_dlgn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/14_het2_dlgn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## preprocessing/202009/15_het3_dlgn/outputs/hisat2_mm38_100/r1_ca.count_mm38_100_sno_gene_ID.count.xz contains 25765 rows and merges to 25765 rows.
## Finished reading count data.
## Matched 25419 annotations and counts.
## Bringing together the count matrix and gene information.
## Some annotations were lost in merging, setting them to 'undefined'.
## Saving the expressionset to 'expt.rda'.
## The final expressionset has 25760 rows and 23 columns.
plot_libsize(mm38_hisat)$plot
normalize_expt(mm38_hisat, filter=TRUE, convert="cpm", norm="quant", transform="log2") mm38_first <-
## This function will replace the expt$expressionset slot with:
## log2(cpm(quant(cbcb(data))))
## It will save copies of each step along the way
## in expt$normalized with the corresponding libsizes. Keep libsizes in mind
## when invoking limma. The appropriate libsize is non-log(cpm(normalized)).
## This is most likely kept at:
## 'new_expt$normalized$intermediate_counts$normalization$libsizes'
## A copy of this may also be found at:
## new_expt$best_libsize
## Not correcting the count-data for batch effects. If batch is
## included in EdgerR/limma's model, then this is probably wise; but in extreme
## batch effects this is a good parameter to play with.
## Step 1: performing count filter with option: cbcb
## Removing 11159 low-count genes (14601 remaining).
## Step 2: normalizing the data with quant.
## Step 3: converting the data with cpm.
## Step 4: transforming the data with log2.
## transform_counts: Found 20 values equal to 0, adding 1 to the matrix.
## Step 5: not doing batch correction.
pp(file="pca_norm.png", image=plot_pca(mm38_first)$plot)
## Writing the image to: pca_norm.png and calling dev.off().
normalize_expt(mm38_hisat, filter=TRUE, convert="cpm", batch="svaseq") mm38_norm <-
## This function will replace the expt$expressionset slot with:
## svaseq(cpm(cbcb(data)))
## It will save copies of each step along the way
## in expt$normalized with the corresponding libsizes. Keep libsizes in mind
## when invoking limma. The appropriate libsize is non-log(cpm(normalized)).
## This is most likely kept at:
## 'new_expt$normalized$intermediate_counts$normalization$libsizes'
## A copy of this may also be found at:
## new_expt$best_libsize
## Leaving the data in its current base format, keep in mind that
## some metrics are easier to see when the data is log2 transformed, but
## EdgeR/DESeq do not accept transformed data.
## Leaving the data unnormalized. This is necessary for DESeq, but
## EdgeR/limma might benefit from normalization. Good choices include quantile,
## size-factor, tmm, etc.
## Step 1: performing count filter with option: cbcb
## Removing 11159 low-count genes (14601 remaining).
## Step 2: not normalizing the data.
## Step 3: converting the data with cpm.
## Step 4: not transforming the data.
## Step 5: doing batch correction with svaseq.
## Using the current state of normalization.
## Passing the data to all_adjusters using the svaseq estimate type.
## batch_counts: Before batch/surrogate estimation, 309316 entries are x>1: 92%.
## batch_counts: Before batch/surrogate estimation, 2857 entries are x==0: 1%.
## batch_counts: Before batch/surrogate estimation, 23650 entries are 0<x<1: 7%.
## The be method chose 2 surrogate variables.
## Attempting svaseq estimation with 2 surrogates.
## There are 2108 (1%) elements which are < 0 after batch correction.
## Setting low elements to zero.
normalize_expt(mm38_norm, transform="log2") mm38_norm <-
## This function will replace the expt$expressionset slot with:
## log2(data)
## It will save copies of each step along the way
## in expt$normalized with the corresponding libsizes. Keep libsizes in mind
## when invoking limma. The appropriate libsize is non-log(cpm(normalized)).
## This is most likely kept at:
## 'new_expt$normalized$intermediate_counts$normalization$libsizes'
## A copy of this may also be found at:
## new_expt$best_libsize
## Filter is false, this should likely be set to something, good
## choices include cbcb, kofa, pofa (anything but FALSE). If you want this to
## stay FALSE, keep in mind that if other normalizations are performed, then the
## resulting libsizes are likely to be strange (potentially negative!)
## Leaving the data unconverted. It is often advisable to cpm/rpkm
## the data to normalize for sampling differences, keep in mind though that rpkm
## has some annoying biases, and voom() by default does a cpm (though hpgl_voom()
## will try to detect this).
## Leaving the data unnormalized. This is necessary for DESeq, but
## EdgeR/limma might benefit from normalization. Good choices include quantile,
## size-factor, tmm, etc.
## Not correcting the count-data for batch effects. If batch is
## included in EdgerR/limma's model, then this is probably wise; but in extreme
## batch effects this is a good parameter to play with.
## Step 1: not doing count filtering.
## Step 2: not normalizing the data.
## Step 3: not converting the data.
## Step 4: transforming the data with log2.
## transform_counts: Found 2108 values equal to 0, adding 1 to the matrix.
## Step 5: not doing batch correction.
pp(file="pca_sva.png", image=plot_pca(mm38_norm)$plot)
## Writing the image to: pca_sva.png and calling dev.off().
create_expt("sample_sheets/all_samples.xlsx", tx_gene_map=tx_gene_map,
mm38_salmon <-gene_info=mm_annot, file_column="salmonfile")
mm_annot
mmtx_annot <-rownames(mmtx_annot) <- mm_annot[["txid"]]
sm(create_expt("sample_sheets/all_samples.xlsx",
mm38_saltx <-gene_info=mmtx_annot, file_column="salmonfile"))
In this block I will calculate all the diagnostic plots, but not show them. I will show them next with a little annotation.
I will leave the output for the first of each invocation and silence it for the second.
sm(graph_metrics(mm38_salmon))
mm38_plots_sa <-
normalize_expt(mm38_salmon, norm="quant", convert="cpm",
mm38_norm_sa <-transform="log2", filter=TRUE)
sm(graph_metrics(mm38_norm_sa)) mm38n_plots_sa <-
$legend
mm38_plots_sa$libsize
mm38_plots_sa$nonzero
mm38_plots_sa$density
mm38n_plots_sa$pc_plot mm38n_plots_sa
sm(graph_metrics(mm38_hisat)) mm38_plots_hi <-
normalize_expt(mm38_hisat, norm="quant", convert="cpm",
mm38_norm_hi <-transform="log2", filter=TRUE)
## This function will replace the expt$expressionset slot with:
## log2(cpm(quant(cbcb(data))))
## It will save copies of each step along the way
## in expt$normalized with the corresponding libsizes. Keep libsizes in mind
## when invoking limma. The appropriate libsize is non-log(cpm(normalized)).
## This is most likely kept at:
## 'new_expt$normalized$intermediate_counts$normalization$libsizes'
## A copy of this may also be found at:
## new_expt$best_libsize
## Not correcting the count-data for batch effects. If batch is
## included in EdgerR/limma's model, then this is probably wise; but in extreme
## batch effects this is a good parameter to play with.
## Step 1: performing count filter with option: cbcb
## Removing 11159 low-count genes (14601 remaining).
## Step 2: normalizing the data with quant.
## Step 3: converting the data with cpm.
## Step 4: transforming the data with log2.
## transform_counts: Found 20 values equal to 0, adding 1 to the matrix.
## Step 5: not doing batch correction.
sm(graph_metrics(mm38_norm_hi)) mm38n_plots_hi <-
$libsize mm38_plots_hi
$nonzero mm38_plots_hi
$density mm38n_plots_hi
$pc_plot mm38n_plots_hi
The only interesting DE I see in this is to compare the retinas to the dlgns. I can treat them as replicates and compare.
These differential expression analyses are EXPLICITLY NOT what you care about. However, they are useful for two purposes:
When we receive full replicate sets, this cheater method of encapsulating the data will not longer be required.
set_expt_conditions(mm38_salmon, fact="celltype")
mm_sa <- sm(normalize_expt(mm_sa, convert="rpkm", transform="log2", column="cds_length"))
mm_norm_sa <-plot_pca(mm_norm_sa)$plot
all_pairwise(mm_sa, model_batch=FALSE) mm_de_sa <-
## mm_hi <- set_expt_conditions(mm38_hisat, fact="celltype")
sm(normalize_expt(mm38_hisat, convert="rpkm",
mm_norm_hi <-column="cds_length"))
normalize_expt(mm_norm_hi, filter=TRUE, batch="svaseq") mm_norm_hi <-
## This function will replace the expt$expressionset slot with:
## svaseq(cbcb(data))
## It will save copies of each step along the way
## in expt$normalized with the corresponding libsizes. Keep libsizes in mind
## when invoking limma. The appropriate libsize is non-log(cpm(normalized)).
## This is most likely kept at:
## 'new_expt$normalized$intermediate_counts$normalization$libsizes'
## A copy of this may also be found at:
## new_expt$best_libsize
## Leaving the data in its current base format, keep in mind that
## some metrics are easier to see when the data is log2 transformed, but
## EdgeR/DESeq do not accept transformed data.
## Leaving the data unconverted. It is often advisable to cpm/rpkm
## the data to normalize for sampling differences, keep in mind though that rpkm
## has some annoying biases, and voom() by default does a cpm (though hpgl_voom()
## will try to detect this).
## Leaving the data unnormalized. This is necessary for DESeq, but
## EdgeR/limma might benefit from normalization. Good choices include quantile,
## size-factor, tmm, etc.
## Step 1: performing count filter with option: cbcb
## Removing 11432 low-count genes (14328 remaining).
## Step 2: not normalizing the data.
## Step 3: not converting the data.
## Step 4: not transforming the data.
## Step 5: doing batch correction with svaseq.
## Using the current state of normalization.
## Passing the data to all_adjusters using the svaseq estimate type.
## batch_counts: Before batch/surrogate estimation, 295319 entries are x>1: 90%.
## batch_counts: Before batch/surrogate estimation, 4171 entries are x==0: 1%.
## batch_counts: Before batch/surrogate estimation, 30054 entries are 0<x<1: 9%.
## The be method chose 2 surrogate variables.
## Attempting svaseq estimation with 2 surrogates.
## There are 2651 (1%) elements which are < 0 after batch correction.
## Setting low elements to zero.
normalize_expt(mm_norm_hi, transform="log2") mm_norm_hi <-
## This function will replace the expt$expressionset slot with:
## log2(data)
## It will save copies of each step along the way
## in expt$normalized with the corresponding libsizes. Keep libsizes in mind
## when invoking limma. The appropriate libsize is non-log(cpm(normalized)).
## This is most likely kept at:
## 'new_expt$normalized$intermediate_counts$normalization$libsizes'
## A copy of this may also be found at:
## new_expt$best_libsize
## Filter is false, this should likely be set to something, good
## choices include cbcb, kofa, pofa (anything but FALSE). If you want this to
## stay FALSE, keep in mind that if other normalizations are performed, then the
## resulting libsizes are likely to be strange (potentially negative!)
## Leaving the data unconverted. It is often advisable to cpm/rpkm
## the data to normalize for sampling differences, keep in mind though that rpkm
## has some annoying biases, and voom() by default does a cpm (though hpgl_voom()
## will try to detect this).
## Leaving the data unnormalized. This is necessary for DESeq, but
## EdgeR/limma might benefit from normalization. Good choices include quantile,
## size-factor, tmm, etc.
## Not correcting the count-data for batch effects. If batch is
## included in EdgerR/limma's model, then this is probably wise; but in extreme
## batch effects this is a good parameter to play with.
## Step 1: not doing count filtering.
## Step 2: not normalizing the data.
## Step 3: not converting the data.
## Step 4: transforming the data with log2.
## transform_counts: Found 2651 values equal to 0, adding 1 to the matrix.
## Step 5: not doing batch correction.
plot_pca(mm_norm_hi)$plot
list(
keepers <-"wt_dlgnret" = c("wt_dlgn", "wt_retina"),
"wt_scnret" = c("wt_scn", "wt_retina"),
"wt_dlgnscn" = c("wt_dlgn", "wt_scn"),
"normret" = c("het_retina", "wt_retina"),
"koret" = c("ko_retina", "wt_retina"),
"normscn" = c("het_scn", "wt_scn"),
"koscn" = c("ko_scn", "wt_scn"),
"normdlgn" = c("het_dlgn", "wt_dlgn"),
"kodlgn" = c("ko_dlgn", "wt_dlgn"),
"normdlgn_vs_normret" = c("normdlgn", "normret"),
"normscn_vs_normret" = c("normscn", "normret"),
"kodlgn_vs_koret" = c("kodlgn", "koret"),
"koscn_vs_koret" = c("koscn", "koret"),
"koscn_vs_kodlgn" = c("koscn", "kodlgn"),
"koret_vs_normret" = c("ko_retina", "het_retina"),
"koscn_vs_normscn" = c("ko_scn", "het_scn"),
"kodlgn_vs_normdlgn" = c("ko_dlgn", "het_dlgn"),
"normko_retdlgn" = c("normdlgn_vs_normret", "kodlgn_vs_koret"),
"normko_retscn" = c("normscn_vs_normret", "koscn_vs_koret"))
"normdlgn_vs_normret = (het_dlgn-wt_dlgn)-(het_retina-wt_retina),
extras <- normscn_vs_normret = (het_scn-wt_scn)-(het_retina-wt_retina),
kodlgn_vs_koret = (ko_dlgn-wt_dlgn)-(ko_retina-wt_retina),
koscn_vs_koret = (ko_scn-wt_scn)-(ko_retina-wt_retina),
koscn_vs_kodlgn = (ko_scn-wt_scn)-(ko_dlgn-wt_dlgn),
normko_vs_retdlgn = ((het_dlgn-wt_dlgn)-(het_retina-wt_retina)) - ((ko_dlgn-wt_dlgn)-(ko_retina-wt_retina)),
normko_vs_retscn = ((het_scn-wt_scn)-(het_retina-wt_retina)) - ((ko_scn-wt_scn)-(ko_retina-wt_retina))"
normalize_expt(mm38_hisat, filter=TRUE) mm_filt <-
## This function will replace the expt$expressionset slot with:
## cbcb(data)
## It will save copies of each step along the way
## in expt$normalized with the corresponding libsizes. Keep libsizes in mind
## when invoking limma. The appropriate libsize is non-log(cpm(normalized)).
## This is most likely kept at:
## 'new_expt$normalized$intermediate_counts$normalization$libsizes'
## A copy of this may also be found at:
## new_expt$best_libsize
## Leaving the data in its current base format, keep in mind that
## some metrics are easier to see when the data is log2 transformed, but
## EdgeR/DESeq do not accept transformed data.
## Leaving the data unconverted. It is often advisable to cpm/rpkm
## the data to normalize for sampling differences, keep in mind though that rpkm
## has some annoying biases, and voom() by default does a cpm (though hpgl_voom()
## will try to detect this).
## Leaving the data unnormalized. This is necessary for DESeq, but
## EdgeR/limma might benefit from normalization. Good choices include quantile,
## size-factor, tmm, etc.
## Not correcting the count-data for batch effects. If batch is
## included in EdgerR/limma's model, then this is probably wise; but in extreme
## batch effects this is a good parameter to play with.
## Step 1: performing count filter with option: cbcb
## Removing 11159 low-count genes (14601 remaining).
## Step 2: not normalizing the data.
## Step 3: not converting the data.
## Step 4: not transforming the data.
## Step 5: not doing batch correction.
limma_pairwise(mm_filt, model_batch="svaseq", extra_contrasts=extras) mm_limma <-
## Starting limma pairwise comparison.
## libsize was not specified, this parameter has profound effects on limma's result.
## Using the libsize from expt$best_libsize.
## Limma step 1/6: choosing model.
## Extracting surrogate estimates from svaseq and adding them to the model.
## batch_counts: Before batch/surrogate estimation, 331257 entries are x>1: 99%.
## batch_counts: Before batch/surrogate estimation, 2857 entries are x==0: 1%.
## The be method chose 2 surrogate variables.
## Attempting svaseq estimation with 2 surrogates.
## Choosing the non-intercept containing model.
## Limma step 2/6: running limma::voom(), switch with the argument 'which_voom'.
## Using normalize.method=quantile for voom.
## Limma step 3/6: running lmFit with method: ls.
## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)
## Limma step 5/6: Running eBayes with robust=FALSE and trend=FALSE.
## Limma step 6/6: Writing limma outputs.
## Limma step 6/6: 1/43: Creating table: het_retina_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 2/43: Creating table: het_scn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 3/43: Creating table: ko_dlgn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 4/43: Creating table: ko_retina_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 5/43: Creating table: ko_scn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 6/43: Creating table: wt_dlgn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 7/43: Creating table: wt_retina_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 8/43: Creating table: wt_scn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 9/43: Creating table: het_scn_vs_het_retina. Adjust=BH
## Limma step 6/6: 10/43: Creating table: ko_dlgn_vs_het_retina. Adjust=BH
## Limma step 6/6: 11/43: Creating table: ko_retina_vs_het_retina. Adjust=BH
## Limma step 6/6: 12/43: Creating table: ko_scn_vs_het_retina. Adjust=BH
## Limma step 6/6: 13/43: Creating table: wt_dlgn_vs_het_retina. Adjust=BH
## Limma step 6/6: 14/43: Creating table: wt_retina_vs_het_retina. Adjust=BH
## Limma step 6/6: 15/43: Creating table: wt_scn_vs_het_retina. Adjust=BH
## Limma step 6/6: 16/43: Creating table: ko_dlgn_vs_het_scn. Adjust=BH
## Limma step 6/6: 17/43: Creating table: ko_retina_vs_het_scn. Adjust=BH
## Limma step 6/6: 18/43: Creating table: ko_scn_vs_het_scn. Adjust=BH
## Limma step 6/6: 19/43: Creating table: wt_dlgn_vs_het_scn. Adjust=BH
## Limma step 6/6: 20/43: Creating table: wt_retina_vs_het_scn. Adjust=BH
## Limma step 6/6: 21/43: Creating table: wt_scn_vs_het_scn. Adjust=BH
## Limma step 6/6: 22/43: Creating table: ko_retina_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 23/43: Creating table: ko_scn_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 24/43: Creating table: wt_dlgn_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 25/43: Creating table: wt_retina_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 26/43: Creating table: wt_scn_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 27/43: Creating table: ko_scn_vs_ko_retina. Adjust=BH
## Limma step 6/6: 28/43: Creating table: wt_dlgn_vs_ko_retina. Adjust=BH
## Limma step 6/6: 29/43: Creating table: wt_retina_vs_ko_retina. Adjust=BH
## Limma step 6/6: 30/43: Creating table: wt_scn_vs_ko_retina. Adjust=BH
## Limma step 6/6: 31/43: Creating table: wt_dlgn_vs_ko_scn. Adjust=BH
## Limma step 6/6: 32/43: Creating table: wt_retina_vs_ko_scn. Adjust=BH
## Limma step 6/6: 33/43: Creating table: wt_scn_vs_ko_scn. Adjust=BH
## Limma step 6/6: 34/43: Creating table: wt_retina_vs_wt_dlgn. Adjust=BH
## Limma step 6/6: 35/43: Creating table: wt_scn_vs_wt_dlgn. Adjust=BH
## Limma step 6/6: 36/43: Creating table: wt_scn_vs_wt_retina. Adjust=BH
## Limma step 6/6: 37/43: Creating table: normdlgn_vs_normret. Adjust=BH
## Limma step 6/6: 38/43: Creating table: normscn_vs_normret. Adjust=BH
## Limma step 6/6: 39/43: Creating table: kodlgn_vs_koret. Adjust=BH
## Limma step 6/6: 40/43: Creating table: koscn_vs_koret. Adjust=BH
## Limma step 6/6: 41/43: Creating table: koscn_vs_kodlgn. Adjust=BH
## Limma step 6/6: 42/43: Creating table: normko_vs_retdlgn. Adjust=BH
## Limma step 6/6: 43/43: Creating table: normko_vs_retscn. Adjust=BH
## Limma step 6/6: 1/9: Creating table: het_dlgn. Adjust=BH
## Limma step 6/6: 2/9: Creating table: het_retina. Adjust=BH
## Limma step 6/6: 3/9: Creating table: het_scn. Adjust=BH
## Limma step 6/6: 4/9: Creating table: ko_dlgn. Adjust=BH
## Limma step 6/6: 5/9: Creating table: ko_retina. Adjust=BH
## Limma step 6/6: 6/9: Creating table: ko_scn. Adjust=BH
## Limma step 6/6: 7/9: Creating table: wt_dlgn. Adjust=BH
## Limma step 6/6: 8/9: Creating table: wt_retina. Adjust=BH
## Limma step 6/6: 9/9: Creating table: wt_scn. Adjust=BH
deseq_pairwise(mm_filt, model_batch="svaseq", extra_contrasts=extras) mm_deseq <-
## Hey you, use deseq2 pairwise.
## Starting DESeq2 pairwise comparisons.
## The data should be suitable for EdgeR/DESeq/EBSeq. If they freak out, check the state of the count table and ensure that it is in integer counts.
## Extracting surrogate estimates from svaseq and adding them to the model.
## batch_counts: Before batch/surrogate estimation, 331257 entries are x>1: 99%.
## batch_counts: Before batch/surrogate estimation, 2857 entries are x==0: 1%.
## The be method chose 2 surrogate variables.
## Attempting svaseq estimation with 2 surrogates.
## Choosing the non-intercept containing model.
## DESeq2 step 1/5: Including a matrix of batch estimates in the deseq model.
## converting counts to integer mode
## DESeq2 step 2/5: Estimate size factors.
## DESeq2 step 3/5: Estimate dispersions.
## gene-wise dispersion estimates
## mean-dispersion relationship
## final dispersion estimates
## Using a parametric fitting seems to have worked.
## DESeq2 step 4/5: nbinomWaldTest.
## The contrast normdlgn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast normscn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast kodlgn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast koscn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast koscn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast normko is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast normko is not in the results.
## If this is not an extra contrast, then this is an error.
edger_pairwise(mm_filt, model_batch="svaseq", extra_contrasts=extras) mm_edger <-
## Starting edgeR pairwise comparisons.
## The data should be suitable for EdgeR/DESeq/EBSeq. If they freak out, check the state of the count table and ensure that it is in integer counts.
## Extracting surrogate estimates from svaseq and adding them to the model.
## batch_counts: Before batch/surrogate estimation, 331257 entries are x>1: 99%.
## batch_counts: Before batch/surrogate estimation, 2857 entries are x==0: 1%.
## The be method chose 2 surrogate variables.
## Attempting svaseq estimation with 2 surrogates.
## Choosing the non-intercept containing model.
## EdgeR step 1/9: Importing and normalizing data.
## EdgeR step 2/9: Estimating the common dispersion.
## EdgeR step 3/9: Estimating dispersion across genes.
## EdgeR step 4/9: Estimating GLM Common dispersion.
## EdgeR step 5/9: Estimating GLM Trended dispersion.
## EdgeR step 6/9: Estimating GLM Tagged dispersion.
## EdgeR step 7/9: Running glmFit, switch to glmQLFit by changing the argument 'edger_test'.
## EdgeR step 8/9: Making pairwise contrasts.
basic_pairwise(mm_filt, model_batch="svaseq", extra_contrasts=extras) mm_basic <-
## Starting basic pairwise comparison.
## Basic step 0/3: Normalizing data.
## Basic step 0/3: Converting data.
## Basic step 0/3: Transforming data.
## Basic step 1/3: Creating mean and variance tables.
## Basic step 2/3: Performing 45 comparisons.
## Basic step 3/3: Creating faux DE Tables.
## Basic: Returning tables.
## mm_ebseq <- ebseq_pairwise(mm_filt, model_batch="svaseq", extra_contrasts=extras)
all_pairwise(mm_filt, model_batch="svaseq", parallel=FALSE, extra_contrasts=extras) mm_de_hi <-
## batch_counts: Before batch/surrogate estimation, 331257 entries are x>1: 99%.
## batch_counts: Before batch/surrogate estimation, 2857 entries are x==0: 1%.
## The be method chose 2 surrogate variables.
## Attempting svaseq estimation with 2 surrogates.
## Plotting a PCA before surrogate/batch inclusion.
## Not putting labels on the PC plot.
## Using svaseq to visualize before/after batch inclusion.
## Performing a test normalization with: raw
## This function will replace the expt$expressionset slot with:
## log2(svaseq(cpm(cbcb(data))))
## It will save copies of each step along the way
## in expt$normalized with the corresponding libsizes. Keep libsizes in mind
## when invoking limma. The appropriate libsize is non-log(cpm(normalized)).
## This is most likely kept at:
## 'new_expt$normalized$intermediate_counts$normalization$libsizes'
## A copy of this may also be found at:
## new_expt$best_libsize
## Leaving the data unnormalized. This is necessary for DESeq, but
## EdgeR/limma might benefit from normalization. Good choices include quantile,
## size-factor, tmm, etc.
## Step 1: performing count filter with option: cbcb
## Removing 0 low-count genes (14601 remaining).
## Step 2: not normalizing the data.
## Step 3: converting the data with cpm.
## Step 4: transforming the data with log2.
## transform_counts: Found 2857 values equal to 0, adding 1 to the matrix.
## Step 5: doing batch correction with svaseq.
## Using the current state of normalization.
## Passing the data to all_adjusters using the svaseq estimate type.
## batch_counts: Before batch/surrogate estimation, 309316 entries are x>1: 92%.
## batch_counts: Before batch/surrogate estimation, 2857 entries are x==0: 1%.
## batch_counts: Before batch/surrogate estimation, 23650 entries are 0<x<1: 7%.
## The be method chose 2 surrogate variables.
## Attempting svaseq estimation with 2 surrogates.
## There are 2108 (1%) elements which are < 0 after batch correction.
## Setting low elements to zero.
## Not putting labels on the PC plot.
## Starting basic_pairwise().
## Starting basic pairwise comparison.
## Basic step 0/3: Normalizing data.
## Basic step 0/3: Converting data.
## Basic step 0/3: Transforming data.
## Basic step 1/3: Creating mean and variance tables.
## Basic step 2/3: Performing 45 comparisons.
## Basic step 3/3: Creating faux DE Tables.
## Basic: Returning tables.
## Starting deseq_pairwise().
## Starting DESeq2 pairwise comparisons.
## The data should be suitable for EdgeR/DESeq/EBSeq. If they freak out, check the state of the count table and ensure that it is in integer counts.
## Including batch estimates from sva/ruv/pca in the model.
## Choosing the non-intercept containing model.
## DESeq2 step 1/5: Including a matrix of batch estimates in the deseq model.
## converting counts to integer mode
## DESeq2 step 2/5: Estimate size factors.
## DESeq2 step 3/5: Estimate dispersions.
## gene-wise dispersion estimates
## mean-dispersion relationship
## final dispersion estimates
## Using a parametric fitting seems to have worked.
## DESeq2 step 4/5: nbinomWaldTest.
## The contrast normdlgn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast normscn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast kodlgn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast koscn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast koscn is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast normko is not in the results.
## If this is not an extra contrast, then this is an error.
## The contrast normko is not in the results.
## If this is not an extra contrast, then this is an error.
## Starting edger_pairwise().
## Starting edgeR pairwise comparisons.
## The data should be suitable for EdgeR/DESeq/EBSeq. If they freak out, check the state of the count table and ensure that it is in integer counts.
## Including batch estimates from sva/ruv/pca in the model.
## Choosing the non-intercept containing model.
## EdgeR step 1/9: Importing and normalizing data.
## EdgeR step 2/9: Estimating the common dispersion.
## EdgeR step 3/9: Estimating dispersion across genes.
## EdgeR step 4/9: Estimating GLM Common dispersion.
## EdgeR step 5/9: Estimating GLM Trended dispersion.
## EdgeR step 6/9: Estimating GLM Tagged dispersion.
## EdgeR step 7/9: Running glmFit, switch to glmQLFit by changing the argument 'edger_test'.
## EdgeR step 8/9: Making pairwise contrasts.
## Starting limma_pairwise().
## Starting limma pairwise comparison.
## libsize was not specified, this parameter has profound effects on limma's result.
## Using the libsize from expt$best_libsize.
## Limma step 1/6: choosing model.
## Including batch estimates from sva/ruv/pca in the model.
## Choosing the non-intercept containing model.
## Limma step 2/6: running limma::voom(), switch with the argument 'which_voom'.
## Using normalize.method=quantile for voom.
## Limma step 3/6: running lmFit with method: ls.
## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)
## Limma step 5/6: Running eBayes with robust=FALSE and trend=FALSE.
## Limma step 6/6: Writing limma outputs.
## Limma step 6/6: 1/43: Creating table: het_retina_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 2/43: Creating table: het_scn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 3/43: Creating table: ko_dlgn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 4/43: Creating table: ko_retina_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 5/43: Creating table: ko_scn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 6/43: Creating table: wt_dlgn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 7/43: Creating table: wt_retina_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 8/43: Creating table: wt_scn_vs_het_dlgn. Adjust=BH
## Limma step 6/6: 9/43: Creating table: het_scn_vs_het_retina. Adjust=BH
## Limma step 6/6: 10/43: Creating table: ko_dlgn_vs_het_retina. Adjust=BH
## Limma step 6/6: 11/43: Creating table: ko_retina_vs_het_retina. Adjust=BH
## Limma step 6/6: 12/43: Creating table: ko_scn_vs_het_retina. Adjust=BH
## Limma step 6/6: 13/43: Creating table: wt_dlgn_vs_het_retina. Adjust=BH
## Limma step 6/6: 14/43: Creating table: wt_retina_vs_het_retina. Adjust=BH
## Limma step 6/6: 15/43: Creating table: wt_scn_vs_het_retina. Adjust=BH
## Limma step 6/6: 16/43: Creating table: ko_dlgn_vs_het_scn. Adjust=BH
## Limma step 6/6: 17/43: Creating table: ko_retina_vs_het_scn. Adjust=BH
## Limma step 6/6: 18/43: Creating table: ko_scn_vs_het_scn. Adjust=BH
## Limma step 6/6: 19/43: Creating table: wt_dlgn_vs_het_scn. Adjust=BH
## Limma step 6/6: 20/43: Creating table: wt_retina_vs_het_scn. Adjust=BH
## Limma step 6/6: 21/43: Creating table: wt_scn_vs_het_scn. Adjust=BH
## Limma step 6/6: 22/43: Creating table: ko_retina_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 23/43: Creating table: ko_scn_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 24/43: Creating table: wt_dlgn_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 25/43: Creating table: wt_retina_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 26/43: Creating table: wt_scn_vs_ko_dlgn. Adjust=BH
## Limma step 6/6: 27/43: Creating table: ko_scn_vs_ko_retina. Adjust=BH
## Limma step 6/6: 28/43: Creating table: wt_dlgn_vs_ko_retina. Adjust=BH
## Limma step 6/6: 29/43: Creating table: wt_retina_vs_ko_retina. Adjust=BH
## Limma step 6/6: 30/43: Creating table: wt_scn_vs_ko_retina. Adjust=BH
## Limma step 6/6: 31/43: Creating table: wt_dlgn_vs_ko_scn. Adjust=BH
## Limma step 6/6: 32/43: Creating table: wt_retina_vs_ko_scn. Adjust=BH
## Limma step 6/6: 33/43: Creating table: wt_scn_vs_ko_scn. Adjust=BH
## Limma step 6/6: 34/43: Creating table: wt_retina_vs_wt_dlgn. Adjust=BH
## Limma step 6/6: 35/43: Creating table: wt_scn_vs_wt_dlgn. Adjust=BH
## Limma step 6/6: 36/43: Creating table: wt_scn_vs_wt_retina. Adjust=BH
## Limma step 6/6: 37/43: Creating table: normdlgn_vs_normret. Adjust=BH
## Limma step 6/6: 38/43: Creating table: normscn_vs_normret. Adjust=BH
## Limma step 6/6: 39/43: Creating table: kodlgn_vs_koret. Adjust=BH
## Limma step 6/6: 40/43: Creating table: koscn_vs_koret. Adjust=BH
## Limma step 6/6: 41/43: Creating table: koscn_vs_kodlgn. Adjust=BH
## Limma step 6/6: 42/43: Creating table: normko_vs_retdlgn. Adjust=BH
## Limma step 6/6: 43/43: Creating table: normko_vs_retscn. Adjust=BH
## Limma step 6/6: 1/9: Creating table: het_dlgn. Adjust=BH
## Limma step 6/6: 2/9: Creating table: het_retina. Adjust=BH
## Limma step 6/6: 3/9: Creating table: het_scn. Adjust=BH
## Limma step 6/6: 4/9: Creating table: ko_dlgn. Adjust=BH
## Limma step 6/6: 5/9: Creating table: ko_retina. Adjust=BH
## Limma step 6/6: 6/9: Creating table: ko_scn. Adjust=BH
## Limma step 6/6: 7/9: Creating table: wt_dlgn. Adjust=BH
## Limma step 6/6: 8/9: Creating table: wt_retina. Adjust=BH
## Limma step 6/6: 9/9: Creating table: wt_scn. Adjust=BH
## Comparing analyses.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, het_retina_vs_het_dlgn and deseq,
## het_retina_vs_het_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, het_scn_vs_het_dlgn and deseq,
## het_scn_vs_het_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_dlgn_vs_het_dlgn and deseq,
## ko_dlgn_vs_het_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_retina_vs_het_dlgn and deseq,
## ko_retina_vs_het_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_scn_vs_het_dlgn and deseq,
## ko_scn_vs_het_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_dlgn_vs_het_dlgn and deseq,
## wt_dlgn_vs_het_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_retina_vs_het_dlgn and deseq,
## wt_retina_vs_het_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_scn_vs_het_dlgn and deseq,
## wt_scn_vs_het_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, het_scn_vs_het_retina and deseq,
## het_scn_vs_het_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_dlgn_vs_het_retina and deseq,
## ko_dlgn_vs_het_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_retina_vs_het_retina and deseq,
## ko_retina_vs_het_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_scn_vs_het_retina and deseq,
## ko_scn_vs_het_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_dlgn_vs_het_retina and deseq,
## wt_dlgn_vs_het_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_retina_vs_het_retina and deseq,
## wt_retina_vs_het_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_scn_vs_het_retina and deseq,
## wt_scn_vs_het_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_dlgn_vs_het_scn and deseq,
## ko_dlgn_vs_het_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_retina_vs_het_scn and deseq,
## ko_retina_vs_het_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_scn_vs_het_scn and deseq,
## ko_scn_vs_het_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_dlgn_vs_het_scn and deseq,
## wt_dlgn_vs_het_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_retina_vs_het_scn and deseq,
## wt_retina_vs_het_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_scn_vs_het_scn and deseq,
## wt_scn_vs_het_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_retina_vs_ko_dlgn and deseq,
## ko_retina_vs_ko_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_scn_vs_ko_dlgn and deseq,
## ko_scn_vs_ko_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_dlgn_vs_ko_dlgn and deseq,
## wt_dlgn_vs_ko_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_retina_vs_ko_dlgn and deseq,
## wt_retina_vs_ko_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_scn_vs_ko_dlgn and deseq,
## wt_scn_vs_ko_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, ko_scn_vs_ko_retina and deseq,
## ko_scn_vs_ko_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_dlgn_vs_ko_retina and deseq,
## wt_dlgn_vs_ko_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_retina_vs_ko_retina and deseq,
## wt_retina_vs_ko_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_scn_vs_ko_retina and deseq,
## wt_scn_vs_ko_retina failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_dlgn_vs_ko_scn and deseq,
## wt_dlgn_vs_ko_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_retina_vs_ko_scn and deseq,
## wt_retina_vs_ko_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_scn_vs_ko_scn and deseq,
## wt_scn_vs_ko_scn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_retina_vs_wt_dlgn and deseq,
## wt_retina_vs_wt_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_scn_vs_wt_dlgn and deseq,
## wt_scn_vs_wt_dlgn failed.
## Used reverse contrast for deseq.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of limma, wt_scn_vs_wt_retina and deseq,
## wt_scn_vs_wt_retina failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, het_retina_vs_het_dlgn=het_retina-
## het_dlgn, and edger, het_retina_vs_het_dlgn=het_retina-het_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, het_scn_vs_het_dlgn=het_scn-het_dlgn, and
## edger, het_scn_vs_het_dlgn=het_scn-het_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn, and
## edger, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,
## and edger, ko_retina_vs_het_dlgn=ko_retina-het_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_scn_vs_het_dlgn=ko_scn-het_dlgn, and
## edger, ko_scn_vs_het_dlgn=ko_scn-het_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn, and
## edger, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,
## and edger, wt_retina_vs_het_dlgn=wt_retina-het_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_het_dlgn=wt_scn-het_dlgn, and
## edger, wt_scn_vs_het_dlgn=wt_scn-het_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, het_scn_vs_het_retina=het_scn-het_retina,
## and edger, het_scn_vs_het_retina=het_scn-het_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,
## and edger, ko_dlgn_vs_het_retina=ko_dlgn-het_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, ko_retina_vs_het_retina=ko_retina-
## het_retina, and edger, ko_retina_vs_het_retina=ko_retina-het_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_scn_vs_het_retina=ko_scn-het_retina,
## and edger, ko_scn_vs_het_retina=ko_scn-het_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,
## and edger, wt_dlgn_vs_het_retina=wt_dlgn-het_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_retina_vs_het_retina=wt_retina-
## het_retina, and edger, wt_retina_vs_het_retina=wt_retina-het_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_het_retina=wt_scn-het_retina,
## and edger, wt_scn_vs_het_retina=wt_scn-het_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_dlgn_vs_het_scn=ko_dlgn-het_scn, and
## edger, ko_dlgn_vs_het_scn=ko_dlgn-het_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_retina_vs_het_scn=ko_retina-het_scn,
## and edger, ko_retina_vs_het_scn=ko_retina-het_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, ko_scn_vs_het_scn=ko_scn-het_scn, and
## edger, ko_scn_vs_het_scn=ko_scn-het_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_het_scn=wt_dlgn-het_scn, and
## edger, wt_dlgn_vs_het_scn=wt_dlgn-het_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_het_scn=wt_retina-het_scn,
## and edger, wt_retina_vs_het_scn=wt_retina-het_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_scn_vs_het_scn=wt_scn-het_scn, and
## edger, wt_scn_vs_het_scn=wt_scn-het_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,
## and edger, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn, and
## edger, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn, and
## edger, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,
## and edger, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn, and
## edger, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_scn_vs_ko_retina=ko_scn-ko_retina, and
## edger, ko_scn_vs_ko_retina=ko_scn-ko_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,
## and edger, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_retina_vs_ko_retina=wt_retina-
## ko_retina, and edger, wt_retina_vs_ko_retina=wt_retina-ko_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_ko_retina=wt_scn-ko_retina, and
## edger, wt_scn_vs_ko_retina=wt_scn-ko_retina, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn, and
## edger, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_ko_scn=wt_retina-ko_scn, and
## edger, wt_retina_vs_ko_scn=wt_retina-ko_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_ko_scn=wt_scn-ko_scn, and edger,
## wt_scn_vs_ko_scn=wt_scn-ko_scn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,
## and edger, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn, and
## edger, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn, failed.
## Used reverse contrast for edger.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_wt_retina=wt_scn-wt_retina, and
## edger, wt_scn_vs_wt_retina=wt_scn-wt_retina, failed.
## Used reverse contrast for deseq.
## Used reverse contrast for edger.
## Used reverse contrast for deseq.
## Used reverse contrast for edger.
## Used reverse contrast for deseq.
## Used reverse contrast for edger.
## Used reverse contrast for deseq.
## Used reverse contrast for edger.
## Used reverse contrast for deseq.
## Used reverse contrast for edger.
## Used reverse contrast for deseq.
## Used reverse contrast for edger.
## Used reverse contrast for deseq.
## Used reverse contrast for edger.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, het_retina_vs_het_dlgn=het_retina-
## het_dlgn, and basic, het_retina_vs_het_dlgn=het_retina-het_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, het_scn_vs_het_dlgn=het_scn-het_dlgn, and
## basic, het_scn_vs_het_dlgn=het_scn-het_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn, and
## basic, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,
## and basic, ko_retina_vs_het_dlgn=ko_retina-het_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_scn_vs_het_dlgn=ko_scn-het_dlgn, and
## basic, ko_scn_vs_het_dlgn=ko_scn-het_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn, and
## basic, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,
## and basic, wt_retina_vs_het_dlgn=wt_retina-het_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_het_dlgn=wt_scn-het_dlgn, and
## basic, wt_scn_vs_het_dlgn=wt_scn-het_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, het_scn_vs_het_retina=het_scn-het_retina,
## and basic, het_scn_vs_het_retina=het_scn-het_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,
## and basic, ko_dlgn_vs_het_retina=ko_dlgn-het_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, ko_retina_vs_het_retina=ko_retina-
## het_retina, and basic, ko_retina_vs_het_retina=ko_retina-het_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_scn_vs_het_retina=ko_scn-het_retina,
## and basic, ko_scn_vs_het_retina=ko_scn-het_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,
## and basic, wt_dlgn_vs_het_retina=wt_dlgn-het_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_retina_vs_het_retina=wt_retina-
## het_retina, and basic, wt_retina_vs_het_retina=wt_retina-het_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_het_retina=wt_scn-het_retina,
## and basic, wt_scn_vs_het_retina=wt_scn-het_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_dlgn_vs_het_scn=ko_dlgn-het_scn, and
## basic, ko_dlgn_vs_het_scn=ko_dlgn-het_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_retina_vs_het_scn=ko_retina-het_scn,
## and basic, ko_retina_vs_het_scn=ko_retina-het_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, ko_scn_vs_het_scn=ko_scn-het_scn, and
## basic, ko_scn_vs_het_scn=ko_scn-het_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_het_scn=wt_dlgn-het_scn, and
## basic, wt_dlgn_vs_het_scn=wt_dlgn-het_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_het_scn=wt_retina-het_scn,
## and basic, wt_retina_vs_het_scn=wt_retina-het_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_scn_vs_het_scn=wt_scn-het_scn, and
## basic, wt_scn_vs_het_scn=wt_scn-het_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,
## and basic, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn, and
## basic, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn, and
## basic, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,
## and basic, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn, and
## basic, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, ko_scn_vs_ko_retina=ko_scn-ko_retina, and
## basic, ko_scn_vs_ko_retina=ko_scn-ko_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,
## and basic, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_retina_vs_ko_retina=wt_retina-
## ko_retina, and basic, wt_retina_vs_ko_retina=wt_retina-ko_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_ko_retina=wt_scn-ko_retina, and
## basic, wt_scn_vs_ko_retina=wt_scn-ko_retina, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn, and
## basic, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_ko_scn=wt_retina-ko_scn, and
## basic, wt_retina_vs_ko_scn=wt_retina-ko_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_ko_scn=wt_scn-ko_scn, and basic,
## wt_scn_vs_ko_scn=wt_scn-ko_scn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,
## and basic, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts
## = extra_contrasts): The merge of deseq, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn, and
## basic, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn, failed.
## Used reverse contrast for basic.
## Warning in correlate_de_tables(results, annot_df = annot_df, extra_contrasts =
## extra_contrasts): The merge of deseq, wt_scn_vs_wt_retina=wt_scn-wt_retina, and
## basic, wt_scn_vs_wt_retina=wt_scn-wt_retina, failed.
## Used reverse contrast for deseq.
## Used reverse contrast for basic.
## Used reverse contrast for deseq.
## Used reverse contrast for basic.
## Used reverse contrast for deseq.
## Used reverse contrast for basic.
## Used reverse contrast for deseq.
## Used reverse contrast for basic.
## Used reverse contrast for deseq.
## Used reverse contrast for basic.
## Used reverse contrast for deseq.
## Used reverse contrast for basic.
## Used reverse contrast for deseq.
## Used reverse contrast for basic.
combine_de_tables(mm_de_hi, excel="excel/testing_202010.xlsx", keepers=keepers) mm_de_tables <-
## Deleting the file excel/testing_202010.xlsx before writing the tables.
## Writing a legend of columns.
## Printing a pca plot before/after surrogates/batch estimation.
## Working on 1/19: wt_dlgnret which is: wt_dlgn/wt_retina.
## Found inverse table with wt_retina_vs_wt_dlgn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 2/19: wt_scnret which is: wt_scn/wt_retina.
## Found table with wt_scn_vs_wt_retina
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 3/19: wt_dlgnscn which is: wt_dlgn/wt_scn.
## Found inverse table with wt_scn_vs_wt_dlgn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 4/19: normret which is: het_retina/wt_retina.
## Found inverse table with wt_retina_vs_het_retina
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 5/19: koret which is: ko_retina/wt_retina.
## Found inverse table with wt_retina_vs_ko_retina
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 6/19: normscn which is: het_scn/wt_scn.
## Found inverse table with wt_scn_vs_het_scn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 7/19: koscn which is: ko_scn/wt_scn.
## Found inverse table with wt_scn_vs_ko_scn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 8/19: normdlgn which is: het_dlgn/wt_dlgn.
## Found inverse table with wt_dlgn_vs_het_dlgn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 9/19: kodlgn which is: ko_dlgn/wt_dlgn.
## Found inverse table with wt_dlgn_vs_ko_dlgn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 10/19: normdlgn_vs_normret which is: normdlgn/normret.
## Found table with normdlgn_vs_normret
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Used the inverse table, might need to -1 the logFC.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The basic table seems to be missing.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 11/19: normscn_vs_normret which is: normscn/normret.
## Found table with normscn_vs_normret
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Used the inverse table, might need to -1 the logFC.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The basic table seems to be missing.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 12/19: kodlgn_vs_koret which is: kodlgn/koret.
## Found table with kodlgn_vs_koret
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Used the inverse table, might need to -1 the logFC.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The basic table seems to be missing.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 13/19: koscn_vs_koret which is: koscn/koret.
## Found table with koscn_vs_koret
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Used the inverse table, might need to -1 the logFC.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The basic table seems to be missing.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 14/19: koscn_vs_kodlgn which is: koscn/kodlgn.
## Found table with koscn_vs_kodlgn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Used the inverse table, might need to -1 the logFC.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The basic table seems to be missing.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 15/19: koret_vs_normret which is: ko_retina/het_retina.
## Found table with ko_retina_vs_het_retina
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 16/19: koscn_vs_normscn which is: ko_scn/het_scn.
## Found table with ko_scn_vs_het_scn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 17/19: kodlgn_vs_normdlgn which is: ko_dlgn/het_dlgn.
## Found table with ko_dlgn_vs_het_dlgn
## Used the inverse table, might need to -1 the logFC and stat.
## Warning in combine_single_de_table(li = limma, ed = edger, eb = ebseq, de =
## deseq, : The deseq table seems to be missing.
## The ebseq table is null.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Working on 18/19: normko_retdlgn which is: normdlgn_vs_normret/kodlgn_vs_koret.
## FOUND NEITHER normdlgn_vs_normret_vs_kodlgn_vs_koret NOR kodlgn_vs_koret_vs_normdlgn_vs_normret!
## Adding venn plots for wt_dlgnret.
## Limma expression coefficients for wt_dlgnret; R^2: 0.748; equation: y = 0.834x + 0.769
## Deseq expression coefficients for wt_dlgnret; R^2: 0.75; equation: y = 0.833x + 1.47
## Edger expression coefficients for wt_dlgnret; R^2: 0.751; equation: y = 0.835x + 1.85
## Adding venn plots for wt_scnret.
## Limma expression coefficients for wt_scnret; R^2: 0.709; equation: y = 0.802x + 0.844
## Deseq expression coefficients for wt_scnret; R^2: 0.719; equation: y = 0.831x + 1.44
## Edger expression coefficients for wt_scnret; R^2: 0.719; equation: y = 0.831x + 1.86
## Adding venn plots for wt_dlgnscn.
## Limma expression coefficients for wt_dlgnscn; R^2: 0.897; equation: y = 0.905x + 0.445
## Deseq expression coefficients for wt_dlgnscn; R^2: 0.898; equation: y = 0.938x + 0.539
## Edger expression coefficients for wt_dlgnscn; R^2: 0.898; equation: y = 0.939x + 0.573
## Adding venn plots for normret.
## Limma expression coefficients for normret; R^2: 0.974; equation: y = 0.997x + 0.0231
## Deseq expression coefficients for normret; R^2: 0.975; equation: y = 0.971x + 0.304
## Edger expression coefficients for normret; R^2: 0.975; equation: y = 0.971x + 0.354
## Adding venn plots for koret.
## Limma expression coefficients for koret; R^2: 0.952; equation: y = 0.98x + 0.0965
## Deseq expression coefficients for koret; R^2: 0.954; equation: y = 0.941x + 0.573
## Edger expression coefficients for koret; R^2: 0.954; equation: y = 0.941x + 0.69
## Adding venn plots for normscn.
## Limma expression coefficients for normscn; R^2: 0.975; equation: y = 1x - 0.0313
## Deseq expression coefficients for normscn; R^2: 0.974; equation: y = 1x - 0.05
## Edger expression coefficients for normscn; R^2: 0.974; equation: y = 1x - 0.0377
## Adding venn plots for koscn.
## Limma expression coefficients for koscn; R^2: 0.966; equation: y = 0.978x + 0.0697
## Deseq expression coefficients for koscn; R^2: 0.967; equation: y = 0.956x + 0.362
## Edger expression coefficients for koscn; R^2: 0.967; equation: y = 0.957x + 0.419
## Adding venn plots for normdlgn.
## Limma expression coefficients for normdlgn; R^2: 0.966; equation: y = 0.997x + 0.0135
## Deseq expression coefficients for normdlgn; R^2: 0.967; equation: y = 0.988x + 0.152
## Edger expression coefficients for normdlgn; R^2: 0.967; equation: y = 0.988x + 0.123
## Adding venn plots for kodlgn.
## Limma expression coefficients for kodlgn; R^2: 0.97; equation: y = 1x - 0.0142
## Deseq expression coefficients for kodlgn; R^2: 0.972; equation: y = 0.968x + 0.329
## Edger expression coefficients for kodlgn; R^2: 0.972; equation: y = 0.968x + 0.313
## Adding venn plots for normdlgn_vs_normret.
## Adding venn plots for normscn_vs_normret.
## Adding venn plots for kodlgn_vs_koret.
## Adding venn plots for koscn_vs_koret.
## Adding venn plots for koscn_vs_kodlgn.
## Adding venn plots for koret_vs_normret.
## Limma expression coefficients for koret_vs_normret; R^2: 0.976; equation: y = 0.99x + 0.0452
## Deseq expression coefficients for koret_vs_normret; R^2: 0.978; equation: y = 1x + 0.00545
## Edger expression coefficients for koret_vs_normret; R^2: 0.977; equation: y = 1x + 0.00654
## Adding venn plots for koscn_vs_normscn.
## Limma expression coefficients for koscn_vs_normscn; R^2: 0.966; equation: y = 0.999x + 0.0197
## Deseq expression coefficients for koscn_vs_normscn; R^2: 0.967; equation: y = 1.01x - 0.0983
## Edger expression coefficients for koscn_vs_normscn; R^2: 0.967; equation: y = 1.01x - 0.112
## Adding venn plots for kodlgn_vs_normdlgn.
## Limma expression coefficients for kodlgn_vs_normdlgn; R^2: 0.986; equation: y = 0.988x + 0.0572
## Deseq expression coefficients for kodlgn_vs_normdlgn; R^2: 0.985; equation: y = 1.01x - 0.106
## Edger expression coefficients for kodlgn_vs_normdlgn; R^2: 0.985; equation: y = 1.01x - 0.122
## Writing summary information, compare_plot is: TRUE.
## Performing save of excel/testing_202010.xlsx.
extract_significant_genes(mm_de_tables, excel="excel/testing_sig_202010.xlsx") mm_de_sig <-
## Writing a legend of columns.
## Did not find the ebseq_logfc, skipping ebseq.
## Printing significant genes to the file: excel/testing_sig_202010.xlsx
## 1/17: Creating significant table up_limma_wt_dlgnret
## 2/17: Creating significant table up_limma_wt_scnret
## 3/17: Creating significant table up_limma_wt_dlgnscn
## The up table normret is empty.
## The down table normret is empty.
## The up table koret is empty.
## The down table koret is empty.
## The up table normscn is empty.
## The down table normscn is empty.
## The up table koscn is empty.
## The down table koscn is empty.
## The up table normdlgn is empty.
## The up table kodlgn is empty.
## The up table normdlgn_vs_normret is empty.
## The down table normdlgn_vs_normret is empty.
## The up table normscn_vs_normret is empty.
## The down table normscn_vs_normret is empty.
## The up table kodlgn_vs_koret is empty.
## The down table kodlgn_vs_koret is empty.
## The up table koscn_vs_koret is empty.
## The down table koscn_vs_koret is empty.
## The up table koscn_vs_kodlgn is empty.
## The down table koscn_vs_kodlgn is empty.
## The up table koret_vs_normret is empty.
## The down table koret_vs_normret is empty.
## The up table koscn_vs_normscn is empty.
## The down table koscn_vs_normscn is empty.
## The up table kodlgn_vs_normdlgn is empty.
## The down table kodlgn_vs_normdlgn is empty.
## Printing significant genes to the file: excel/testing_sig_202010.xlsx
## 1/17: Creating significant table up_edger_wt_dlgnret
## 2/17: Creating significant table up_edger_wt_scnret
## 3/17: Creating significant table up_edger_wt_dlgnscn
## 4/17: Creating significant table up_edger_normret
## The down table normret is empty.
## 5/17: Creating significant table up_edger_koret
## The up table normscn is empty.
## The down table normscn is empty.
## 7/17: Creating significant table up_edger_koscn
## The down table koscn is empty.
## The up table normdlgn is empty.
## 9/17: Creating significant table up_edger_kodlgn
## The up table normdlgn_vs_normret is empty.
## The up table normscn_vs_normret is empty.
## The down table normscn_vs_normret is empty.
## The up table kodlgn_vs_koret is empty.
## The up table koscn_vs_koret is empty.
## The down table koscn_vs_koret is empty.
## 14/17: Creating significant table up_edger_koscn_vs_kodlgn
## The down table koscn_vs_kodlgn is empty.
## 15/17: Creating significant table up_edger_koret_vs_normret
## The down table koret_vs_normret is empty.
## 16/17: Creating significant table up_edger_koscn_vs_normscn
## The up table kodlgn_vs_normdlgn is empty.
## The down table kodlgn_vs_normdlgn is empty.
## Unable to find the table in the set of possible tables.
## The possible tables are: het_retina_vs_het_dlgn=het_retina-het_dlgn,, het_scn_vs_het_dlgn=het_scn-het_dlgn,, ko_dlgn_vs_het_dlgn=ko_dlgn-het_dlgn,, ko_retina_vs_het_dlgn=ko_retina-het_dlgn,, ko_scn_vs_het_dlgn=ko_scn-het_dlgn,, wt_dlgn_vs_het_dlgn=wt_dlgn-het_dlgn,, wt_retina_vs_het_dlgn=wt_retina-het_dlgn,, wt_scn_vs_het_dlgn=wt_scn-het_dlgn,, het_scn_vs_het_retina=het_scn-het_retina,, ko_dlgn_vs_het_retina=ko_dlgn-het_retina,, ko_retina_vs_het_retina=ko_retina-het_retina,, ko_scn_vs_het_retina=ko_scn-het_retina,, wt_dlgn_vs_het_retina=wt_dlgn-het_retina,, wt_retina_vs_het_retina=wt_retina-het_retina,, wt_scn_vs_het_retina=wt_scn-het_retina,, ko_dlgn_vs_het_scn=ko_dlgn-het_scn,, ko_retina_vs_het_scn=ko_retina-het_scn,, ko_scn_vs_het_scn=ko_scn-het_scn,, wt_dlgn_vs_het_scn=wt_dlgn-het_scn,, wt_retina_vs_het_scn=wt_retina-het_scn,, wt_scn_vs_het_scn=wt_scn-het_scn,, ko_retina_vs_ko_dlgn=ko_retina-ko_dlgn,, ko_scn_vs_ko_dlgn=ko_scn-ko_dlgn,, wt_dlgn_vs_ko_dlgn=wt_dlgn-ko_dlgn,, wt_retina_vs_ko_dlgn=wt_retina-ko_dlgn,, wt_scn_vs_ko_dlgn=wt_scn-ko_dlgn,, ko_scn_vs_ko_retina=ko_scn-ko_retina,, wt_dlgn_vs_ko_retina=wt_dlgn-ko_retina,, wt_retina_vs_ko_retina=wt_retina-ko_retina,, wt_scn_vs_ko_retina=wt_scn-ko_retina,, wt_dlgn_vs_ko_scn=wt_dlgn-ko_scn,, wt_retina_vs_ko_scn=wt_retina-ko_scn,, wt_scn_vs_ko_scn=wt_scn-ko_scn,, wt_retina_vs_wt_dlgn=wt_retina-wt_dlgn,, wt_scn_vs_wt_dlgn=wt_scn-wt_dlgn,, wt_scn_vs_wt_retina=wt_scn-wt_retina,
## Error in single_ma[["ma"]]: subscript out of bounds
## I changed the sample sheet so that the column is renamed to 'location'
## alt_model <- "~ location + location:genotype"
library(DESeq2)
## Loading required package: S4Vectors
## Loading required package: stats4
##
## Attaching package: 'S4Vectors'
## The following object is masked from 'package:base':
##
## expand.grid
## Loading required package: IRanges
## Loading required package: GenomicRanges
## Loading required package: GenomeInfoDb
## Loading required package: SummarizedExperiment
## Loading required package: DelayedArray
## Loading required package: matrixStats
##
## Attaching package: 'matrixStats'
## The following objects are masked from 'package:hpgltools':
##
## anyMissing, rowMedians
## The following objects are masked from 'package:Biobase':
##
## anyMissing, rowMedians
##
## Attaching package: 'DelayedArray'
## The following objects are masked from 'package:matrixStats':
##
## colMaxs, colMins, colRanges, rowMaxs, rowMins, rowRanges
## The following objects are masked from 'package:base':
##
## aperm, apply, rowsum
"~ celltype + celltype:genotype"
alt_model <- pData(mm_filt)
metadata <-"location"]] <- factor(metadata[["location"]], levels=c("retina", "scn", "dlgn"))
metadata[["genotype"]] <- factor(metadata[["genotype"]], levels=c("wt", "ko", "het"))
metadata[[ DESeqDataSetFromMatrix(countData=exprs(mm_filt),
dds <-colData=metadata,
design=~0+location+genotype+genotype:location)
## converting counts to integer mode
DESeq(dds) dds_run <-
## estimating size factors
## estimating dispersions
## gene-wise dispersion estimates
## mean-dispersion relationship
## final dispersion estimates
## fitting model and testing
resultsNames(dds_run)
## [1] "locationretina" "locationscn"
## [3] "locationdlgn" "genotypeko"
## [5] "genotypehet" "locationscn.genotypeko"
## [7] "locationdlgn.genotypeko" "locationscn.genotypehet"
## [9] "locationdlgn.genotypehet"
as.data.frame(results(dds_run, contrast=list("locationdlgn.genotypeko",
dlgn_vs_scn <-"locationscn.genotypeko")))
is.na(dlgn_vs_scn[["padj"]])
check <-"padj"] <- 1
dlgn_vs_scn[check, dlgn_vs_scn[["padj"]] <= 0.1 check <-
mm_limma$all_tables
tables <-for (i in 1:length(tables)) {
names(tables)[i]
table_name <- tables[[i]]
table <-print(table_name)
dlgn_vs_scn[, c("log2FoldChange", "pvalue")]
t1 <- table[, c("logFC", "P.Value")]
t2 <- merge(t1, t2, by="row.names")
t3 <-print(cor.test(t3[["log2FoldChange"]], t3[["logFC"]]))
}
## [1] "het_retina_vs_het_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -55, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4296 -0.4028
## sample estimates:
## cor
## -0.4163
##
## [1] "het_scn_vs_het_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 16, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.1162 0.1481
## sample estimates:
## cor
## 0.1322
##
## [1] "ko_dlgn_vs_het_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 17, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.1265 0.1583
## sample estimates:
## cor
## 0.1425
##
## [1] "ko_retina_vs_het_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -53, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4122 -0.3849
## sample estimates:
## cor
## -0.3986
##
## [1] "ko_scn_vs_het_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -5.9, df = 14599, p-value = 3e-09
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.06508 -0.03271
## sample estimates:
## cor
## -0.04891
##
## [1] "wt_dlgn_vs_het_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -63, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4773 -0.4519
## sample estimates:
## cor
## -0.4647
##
## [1] "wt_retina_vs_het_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -50, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.3955 -0.3678
## sample estimates:
## cor
## -0.3818
##
## [1] "wt_scn_vs_het_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 34, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.2570 0.2871
## sample estimates:
## cor
## 0.2721
##
## [1] "het_scn_vs_het_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 61, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.4381 0.4640
## sample estimates:
## cor
## 0.4511
##
## [1] "ko_dlgn_vs_het_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 60, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.433 0.459
## sample estimates:
## cor
## 0.4461
##
## [1] "ko_retina_vs_het_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -4.2, df = 14599, p-value = 3e-05
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.05094 -0.01854
## sample estimates:
## cor
## -0.03475
##
## [1] "ko_scn_vs_het_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 49, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.3595 0.3874
## sample estimates:
## cor
## 0.3735
##
## [1] "wt_dlgn_vs_het_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 34, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.2529 0.2830
## sample estimates:
## cor
## 0.268
##
## [1] "wt_retina_vs_het_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 8, df = 14599, p-value = 2e-15
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.04975 0.08205
## sample estimates:
## cor
## 0.06591
##
## [1] "wt_scn_vs_het_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 73, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.5043 0.5281
## sample estimates:
## cor
## 0.5163
##
## [1] "ko_dlgn_vs_het_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -9.4, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.09356 -0.06131
## sample estimates:
## cor
## -0.07746
##
## [1] "ko_retina_vs_het_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -59, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4490 -0.4228
## sample estimates:
## cor
## -0.436
##
## [1] "ko_scn_vs_het_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -40, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.3283 -0.2990
## sample estimates:
## cor
## -0.3137
##
## [1] "wt_dlgn_vs_het_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -46, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.3689 -0.3405
## sample estimates:
## cor
## -0.3548
##
## [1] "wt_retina_vs_het_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -56, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4345 -0.4078
## sample estimates:
## cor
## -0.4212
##
## [1] "wt_scn_vs_het_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 34, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.2574 0.2874
## sample estimates:
## cor
## 0.2725
##
## [1] "ko_retina_vs_ko_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -57, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4422 -0.4157
## sample estimates:
## cor
## -0.429
##
## [1] "ko_scn_vs_ko_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -13, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.1269 -0.0949
## sample estimates:
## cor
## -0.1109
##
## [1] "wt_dlgn_vs_ko_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -84, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.5813 -0.5594
## sample estimates:
## cor
## -0.5704
##
## [1] "wt_retina_vs_ko_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -55, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4275 -0.4006
## sample estimates:
## cor
## -0.4141
##
## [1] "wt_scn_vs_ko_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 28, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.2074 0.2382
## sample estimates:
## cor
## 0.2228
##
## [1] "ko_scn_vs_ko_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 47, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.3506 0.3787
## sample estimates:
## cor
## 0.3647
##
## [1] "wt_dlgn_vs_ko_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 32, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.2420 0.2723
## sample estimates:
## cor
## 0.2572
##
## [1] "wt_retina_vs_ko_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 9.4, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.06157 0.09382
## sample estimates:
## cor
## 0.07772
##
## [1] "wt_scn_vs_ko_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 69, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.4815 0.5060
## sample estimates:
## cor
## 0.4938
##
## [1] "wt_dlgn_vs_ko_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -26, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.2239 -0.1929
## sample estimates:
## cor
## -0.2085
##
## [1] "wt_retina_vs_ko_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -44, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.3546 -0.3259
## sample estimates:
## cor
## -0.3404
##
## [1] "wt_scn_vs_ko_scn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 78, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.5298 0.5527
## sample estimates:
## cor
## 0.5414
##
## [1] "wt_retina_vs_wt_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -30, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.2578 -0.2273
## sample estimates:
## cor
## -0.2426
##
## [1] "wt_scn_vs_wt_dlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 69, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.4834 0.5079
## sample estimates:
## cor
## 0.4957
##
## [1] "wt_scn_vs_wt_retina"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 69, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.4842 0.5086
## sample estimates:
## cor
## 0.4965
##
## [1] "normdlgn_vs_normret"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 61, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.4375 0.4633
## sample estimates:
## cor
## 0.4505
##
## [1] "normscn_vs_normret"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -22, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.1939 -0.1625
## sample estimates:
## cor
## -0.1782
##
## [1] "kodlgn_vs_koret"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 67, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.4705 0.4954
## sample estimates:
## cor
## 0.483
##
## [1] "koscn_vs_koret"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -55, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4253 -0.3984
## sample estimates:
## cor
## -0.4119
##
## [1] "koscn_vs_kodlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -211, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.8716 -0.8635
## sample estimates:
## cor
## -0.8676
##
## [1] "normko_vs_retdlgn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = -14, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.1335 -0.1015
## sample estimates:
## cor
## -0.1175
##
## [1] "normko_vs_retscn"
##
## Pearson's product-moment correlation
##
## data: t3[["log2FoldChange"]] and t3[["logFC"]]
## t = 28, df = 14599, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.2084 0.2392
## sample estimates:
## cor
## 0.2239
Until we get full replicates, I will do simple subtractions.
In an attempt to keep some clarity in the terms used, I want to define them now. There are three contexts in which we will consider the data:
The individual sample type. When considering individual samples, I will use three terms in this and only this context: wild-type (wt), het, and mut.
The individual translatome. These are defines as something / baseline. I will exclusively call the wt samples ‘baseline’ when speaking in this context. I will exclusively state ‘normal’ when referring to het / wt samples, and I will state ‘ko’ when referring to mut / wt samples in the translatome context.
Translatome vs. translatome. Whenever comparing translatomes, I will use the names as in #2 and always put the numerator first when writing the name of a comparison.
The most complex example of the above nomenclature is:
“normko_retdlgn is defined as normret_vs_normdlgn - koret_vs_kodlgn”
This states we are examining at the translatome context: (norm(retina translatome) - norm(dlgn translatome)) - (ko(retina translatome) - ko(dlgn translatome))
Which in turn is synonymous to the following at the sample context: ((rethet - retwt) - (dlgnhet - dlgnwt)) - ((retko - retwt) - (dlgnko - dlgnwt))
Now let us associate the various variable names with the appropriate samples:
"iprgc_01"
dlgnwt <- "iprgc_02"
retwt <- "iprgc_03"
scnwt <-
"iprgc_04"
dlgnhet <- "iprgc_05"
rethet <- NULL ## Does not yet exist.
scnhet <-
"iprgc_06"
dlgnmut <- "iprgc_07"
retmut <- "iprgc_08" scnmut <-
Give these variable names, now lets associate columns of the expression data with them. These are at the sample context, so the appropriate names are: ‘wt’, ‘het’, and ‘mut’. In each case I will prefix the genotype with the tissue type: ‘ret’, ‘dlgn’, and ‘scn’. Thus ‘retwt’ refers to the sample used to calculate the translatome retina baseline; in contrast ‘dlgnmut’ is the sample which provides the dlgn knockout.
## Sample context
mm_norm_sa mm38_norm <-
## Error in eval(expr, envir, enclos): object 'mm_norm_sa' not found
exprs(mm38_norm)[, dlgnwt]
dlgnwt <- exprs(mm38_norm)[, retwt]
retwt <- exprs(mm38_norm)[, scnwt]
scnwt <- exprs(mm38_norm)[, dlgnhet]
dlgnhet <- exprs(mm38_norm)[, rethet]
rethet <- exprs(mm38_norm)[, dlgnmut]
dlgnmut <- exprs(mm38_norm)[, retmut]
retmut <- exprs(mm38_norm)[, scnmut] scnmut <-
Each of the above 8 variables provides 1 column of information. We have 3 baseline comparisons available to us. In each of these we compare one wt sample to another.
## Baseline comparisons
dlgnwt - retwt
wt_dlgnret <- scnwt - retwt
wt_scnret <- dlgnwt - scnwt wt_dlgnscn <-
Simultaneously, we have 5 available translatomes. This are provided by comparing each het or mut to the associated wt. These will therefore receive names: ‘norm’ and ‘ko’ instead of ‘het’ and ‘mut’.
## Translatome context
rethet - retwt
normret <- retmut - retwt
koret <- scnmut - scnwt
koscn <- dlgnhet - dlgnwt
normdlgn <- dlgnmut - dlgnwt kodlgn <-
Given these translatomes, there are a few contrasts of likely interest. These are performed by comparing the relevant translatomes.
Will will split these into 4 separate categories: het vs het, ko vs ko, ko vs het, and ratio vs ratio.
Finally, note that we are being explicitly redundant in these definitions. I am making variable names for both the a/b ratio and the b/a ratio. Thus we have some redundantly redundant (haha) flexibility when deciding on what we want to plot.
## norm vs norm
normdlgn - normret
normdlgn_vs_normret <- normret - normdlgn normret_vs_normdlgn <-
## ko vs ko
koret - kodlgn
koret_vs_kodlgn <- kodlgn - koret
kodlgn_vs_koret <-
koret - koscn
koret_vs_koscn <- koscn - koret
koscn_vs_koret <-
kodlgn - koscn
kodlgn_vs_koscn <- koscn - kodlgn koscn_vs_kodlgn <-
On the other hand, I am assuming we always want the normals as denominators and kos as numerators.
## ko vs norm
koret - normret
koret_vs_normret <-
kodlgn - normdlgn kodlgn_vs_normdlgn <-
Finally, here is the ratio of ratios example I printed above:
I named it ‘normko_retdlgn’ in an attempt to make clear that it is actually: (normret/normdlgn)/(koret/kodlgn)
or stated differently: “norm divided by ko for ret divided by dlgn.”
## ratio of ratios
normret_vs_normdlgn - koret_vs_kodlgn normko_retdlgn <-
My matrix of data will now contain 1 column for each of the above 27 samples/comparisons.
cbind(
pair_mtrx <-## Individual samples
dlgnwt, retwt, scnwt, dlgnhet, rethet, dlgnmut, retmut, scnmut,## Baseline comparisons
wt_dlgnret, wt_scnret, wt_dlgnscn,## Baseline subtractions
normdlgn, normret, kodlgn, koret, koscn,## het_vs_het, of which there is only 1 because we do not have hetscn
normdlgn_vs_normret, normret_vs_normdlgn,## ko_vs_ko, of which we have 3
koret_vs_kodlgn, kodlgn_vs_koret,
koret_vs_koscn, koscn_vs_koret,
kodlgn_vs_koscn, koscn_vs_kodlgn,## ko_vs_het, 3 including one getting around missing hetscn
koret_vs_normret, kodlgn_vs_normdlgn,## ratio of ratios
normko_retdlgn)
I am not sure if we will use these indexes, but I am writing these out as subsets of genes to look at. These indexes are stating that, given a cutoff (0), we want to look at only the genes which have higher x / baseline values than the cutoff.
## Queries about gene subsets.
## These are all in the context of translatomes.
0
cutoff <- normret > cutoff & koret > cutoff
ret_kept_idx <- koscn > cutoff
scn_kept_idx <- normdlgn > cutoff & kodlgn > cutoff
dlgn_kept_idx <- ret_kept_idx & dlgn_kept_idx
ret_dlgn_kept_idx <- ret_kept_idx & scn_kept_idx
ret_scn_kept_idx <- dlgn_kept_idx & scn_kept_idx
dlgn_scn_kept_idx <-
##normdlgn_vs_normret[!ret_dlgn_kept_idx] <- NA
##normret_vs_normdlgn[!ret_dlgn_kept_idx] <- NA
##koret_vs_kodlgn[!ret_dlgn_kept_idx] <- NA
##kodlgn_vs_koret[!ret_dlgn_kept_idx] <- NA
##koret_vs_koscn[!ret_scn_kept_idx] <- NA
##koscn_vs_koret[!ret_scn_kept_idx] <- NA
##kodlgn_vs_koscn[!dlgn_scn_kept_idx] <- NA
##koscn_vs_kodlgn[!dlgn_scn_kept_idx] <- NA
##koret_vs_normret[!ret_kept_idx] <- NA
##kodlgn_vs_normdlgn[!dlgn_kept_idx] <- NA
##normko_retdlgn <- normko_retdlgn[!ret_dlgn_kept_idx] <- NA
I will use my function combine_de_tables() to add this information to my existing annotation data along with the results from the statistically valid comparison of the three tissue types.
sm(combine_de_tables(
mm_tables <-extra_annot=pair_mtrx,
mm_de_sa, excel=glue::glue("excel/{rundate}mm_salmon_tables-v{ver}.xlsx")))
## Put retina baseline on y axis as black, retina het on x axis as black.
## Then recolor a subset of these as red, the reds are when normret > 0
library(ggplot2)
as.data.frame(pair_mtrx[, c("rethet", "retwt")])
plotted <- normret > 0
red_idx <-"color"] <- ifelse(red_idx, "red", "black")
plotted[, "label"]] <- rownames(plotted)
plotted[[ ggplot(
ret_hetwt <-
plotted,aes_string(x="rethet", y="retwt", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
ret_hetwt ggplotly_url(
ret_hetwt_clicky <-"ret_hetwt.html", title="Retina expression, het vs. wt.",
ret_hetwt, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("retmut", "retwt")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
ret_mutwt <-
plotted,aes_string(x="retmut", y="retwt", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
ret_mutwt ggplotly_url(
ret_mutwt_clicky <-"ret_mutwt.html", title="Retina expression, mutant vs. wt.",
ret_mutwt, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("dlgnhet", "dlgnwt")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
dlgn_hetwt <-
plotted,aes_string(x="dlgnhet", y="dlgnwt", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
dlgn_hetwt ggplotly_url(
dlgn_hetwt_clicky <-"dlgn_hetwt.html", title="dlgn expression, het vs. wt.",
dlgn_hetwt, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("dlgnmut", "dlgnwt")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
dlgn_mutwt <-
plotted,aes_string(x="dlgnmut", y="dlgnwt", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
dlgn_mutwt ggplotly_url(
dlgn_mutwt_clicky <-"dlgn_mutwt.html", title="dlgn expression, mut vs. wt.",
dlgn_mutwt, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("scnmut", "scnwt")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
scn_mutwt <-
plotted,aes_string(x="scnmut", y="scnwt", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
scn_mutwt ggplotly_url(
scn_mutwt_clicky <-"scn_mutwt.html", title="scn expression, mut vs. wt.",
scn_mutwt, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
## x-axis: normdlgn_vs_normret or normret_vs_normdlgn,
## ^^^^
## y-axis: dlgnwt-retwt (baseline dlgn - baseline retina)
as.data.frame(pair_mtrx[, c("normdlgn_vs_normret", "wt_dlgnret")])
plotted <- normret > 0
red_idx <-## Note that this order is opposite of above.
"color"] <- ifelse(red_idx, "black", "red")
plotted[, "label"]] <- rownames(plotted)
plotted[[ ggplot(
axon_trans_ret_target <-
plotted,aes_string(x="normdlgn_vs_normret", y="wt_dlgnret", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
axon_trans_ret_target ggplotly_url(
axon_trans_ret_target_clicky <-"axon_trans_ret_target.html", title="Axon translatome, retina target.",
axon_trans_ret_target, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("normret", "normdlgn")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
normret_normdlgn <-
plotted,aes_string(x="normret", y="normdlgn", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
normret_normdlgn ggplotly_url(
normret_normdlgn_clicky <-"normret_normdlgn.html", title="Normal retina translatome vs normal dlgn translatome.",
normret_normdlgn, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("koret", "kodlgn")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
koret_kodlgn <-
plotted,aes_string(x="koret", y="kodlgn", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
koret_kodlgn ggplotly_url(
koret_kodlgn_clicky <-"koret_kodlgn.html", title="KO retina translatome vs KO dlgn translatome.",
koret_kodlgn, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("koret", "koscn")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
koret_koscn <-
plotted,aes_string(x="koret", y="koscn", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
koret_koscn ggplotly_url(
koret_koscn_clicky <-"koret_koscn.html", title="KO retina translatome vs KO scn translatome.",
koret_koscn, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("normdlgn", "kodlgn")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
normdlgn_kodlgn <-
plotted,aes_string(x="normdlgn", y="kodlgn", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
normdlgn_kodlgn ggplotly_url(
normdlgn_kodlgn_clicky <-"normdlgn_kodlgn.html", title="Normal dlgn translatome vs KO dlgn translatome.",
normdlgn_kodlgn, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("normret", "koret")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
normret_koret <-
plotted,aes_string(x="normret", y="koret", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
normret_koret ggplotly_url(
normret_koret_clicky <-"normret_koret.html", title="Normal retina translatome vs KO retina translatome.",
normret_koret, url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
as.data.frame(pair_mtrx[, c("normret_vs_normdlgn", "koret_vs_kodlgn")])
plotted <-"label"]] <- rownames(plotted)
plotted[["color"]] <- "black"
plotted[[ ggplot(
normal_ko_axon_translatome <-
plotted,aes_string(x="normret_vs_normdlgn", y="koret_vs_kodlgn", label="label", color="color")) +
geom_point(alpha=0.5) +
scale_color_manual(values=c("black", "red"))
normal_ko_axon_translatome ggplotly_url(
normal_ko_axon_translatome_clicky <-"normal_ko_axon_translatome.html",
normal_ko_axon_translatome, title="Normal retina ko axon translatome.",
url_data="http://useast.ensembl.org/Mus_musculus/Gene/Summary?g={ids}")
function(pair_mtrx, x_axis="koret", y_axis="koscn", lfc=NULL, fc=NULL, linewidth=1.5,
translatome_plotter <-up_color="red", down_color="#098534", line_color="#fcba03", alpha=0.5,
x_limit=c(-7, 7), y_limit=c(-7, 7)) {
if (is.null(fc) & is.null(lfc)) {
message("No fc/lfc was provided, defaulting to 10 fold.")
log2(10)
lfc <-else if (is.null(lfc)) {
} log2(fc)
lfc <-
} as.data.frame(pair_mtrx[, c(x_axis, y_axis)])
plotted <- is.na(plotted)
na_idx <- 0
plotted[na_idx] <- plotted[, y_axis] - plotted[, x_axis] >= lfc
up_idx <- plotted[, y_axis] - plotted[, x_axis] <= (lfc * -1)
down_idx <- rownames(plotted)[up_idx]
up_genes <- rownames(plotted)[down_idx]
down_genes <-## Note that this order is opposite of above.
"color"]] <- "black"
plotted[["color"] <- up_color
plotted[up_idx, "color"] <- down_color
plotted[down_idx, "color"]] <- as.factor(plotted[["color"]])
plotted[[levels(plotted[["color"]]) <- c("black", up_color, down_color)
ggplot2::ggplot(plotted, aes_string(x=x_axis,
plt <-y=y_axis,
color="color")) +
geom_abline(size=1.1, slope=1, intercept=lfc, color="orange") +
geom_abline(size=1.1, slope=1, intercept=(-1 * lfc), color="orange") +
scale_x_continuous(limits=c(x_limit)) +
scale_y_continuous(limits=c(y_limit)) +
geom_point(alpha=alpha) +
scale_color_manual(values=c(down_color, "black", up_color))
list(
retlist <-"mtrx" = plotted,
"ups" = up_genes,
"downs" = down_genes,
"plot" = plt)
return(retlist)
}
First plot: KO scn translatome on y axis vs. KO retina translatome on x axis.
translatome_plotter(pair_mtrx)
scnko_wrt_retko_translatome <-$plot
scnko_wrt_retko_translatome simple_gprofiler(sig_genes=scnko_wrt_retko_translatome$ups,
scnko_wrt_retko_up_go <-species="mmusculus")
simple_gprofiler(sig_genes=scnko_wrt_retko_translatome$downs,
scnko_wrt_retko_down_go <-species="mmusculus")
$go
scnko_wrt_retko_down_go$kegg
scnko_wrt_retko_down_go$corum scnko_wrt_retko_down_go
translatome_plotter(pair_mtrx,
dlgnnorm_wrt_retnorm_translatome <-x_axis="normret", y_axis="normdlgn")
$plot
dlgnnorm_wrt_retnorm_translatome simple_gprofiler(sig_genes=dlgnnorm_wrt_retnorm_translatome$ups,
dlgnnorm_wrt_retnorm_up_go <-species="mmusculus")
simple_gprofiler(sig_genes=dlgnnorm_wrt_retnorm_translatome$downs,
dlgnnorm_wrt_retnorm_down_go <-species="mmusculus")
$go
dlgnnorm_wrt_retnorm_down_go$pvalue_plots$bpp_plot_over
dlgnnorm_wrt_retnorm_down_go$pvalue_plots$ccp_plot_over dlgnnorm_wrt_retnorm_down_go
translatome_plotter(pair_mtrx,
dlgnko_wrt_retko_translatome <-x_axis="koret", y_axis="kodlgn")
$plot
dlgnko_wrt_retko_translatome simple_gprofiler(sig_genes=dlgnko_wrt_retko_translatome$ups,
dlgnko_wrt_retko_up_go <-species="mmusculus")
$go
dlgnko_wrt_retko_up_go$kegg
dlgnko_wrt_retko_up_go$reac
dlgnko_wrt_retko_up_go simple_gprofiler(sig_genes=dlgnko_wrt_retko_translatome$downs,
dlgnko_wrt_retko_down_go <-species="mmusculus")
$go
dlgnko_wrt_retko_down_go$kegg
dlgnko_wrt_retko_down_go$reac dlgnko_wrt_retko_down_go
translatome_plotter(pair_mtrx,
dlgnnorm_wrt_dlgnko_translatome <-x_axis="normdlgn",
y_axis="kodlgn")
$plot
dlgnnorm_wrt_dlgnko_translatome simple_gprofiler(sig_genes=dlgnnorm_wrt_dlgnko_translatome$ups,
dlgnnorm_wrt_dlgnko_up_go <-species="mmusculus")
simple_gprofiler(sig_genes=dlgnnorm_wrt_dlgnko_translatome$downs,
dlgnnorm_wrt_dlgnko_down_go <-species="mmusculus")
$go dlgnnorm_wrt_dlgnko_down_go
translatome_plotter(pair_mtrx,
dlgnko_wrt_scnko_translatome <-x_axis="kodlgn",
y_axis="koscn")
$plot
dlgnko_wrt_scnko_translatome simple_gprofiler(sig_genes=dlgnko_wrt_scnko_translatome$ups,
dlgnko_wrt_scnko_up_go <-species="mmusculus")
$go
dlgnko_wrt_scnko_up_go$pvalue_plots$bpp_plot_over
dlgnko_wrt_scnko_up_go simple_gprofiler(sig_genes=dlgnko_wrt_scnko_translatome$downs,
dlgnko_wrt_scnko_down_go <-species="mmusculus")
$go dlgnko_wrt_scnko_down_go
As I understand it, there is some interest in an ontology search using the ratio of ratios.
normko_retdlgn
ror <- ror >= 1
up_idx <- ror <= -1
down_idx <- ror[up_idx]
ror_up <-length(ror_up)
ror[down_idx]
ror_down <-length(ror_down)
simple_gprofiler(
ror_gprofiler_up <-sig_genes=ror_up, species="mmusculus",
excel=glue::glue("excel/{rundate}mm_ror_gpfoiler_up-v{ver}.xlsx"))
$pvalue_plots$mfp_plot_over
ror_gprofiler_up$pvalue_plots$bpp_plot_over
ror_gprofiler_up$pvalue_plots$ccp_plot_over
ror_gprofiler_up$pvalue_plots$tf_plot_over
ror_gprofiler_up$pvalue_plots$hp_plot_over
ror_gprofiler_up
simple_gprofiler(
ror_gprofiler_down <-sig_genes=ror_down, species="mmusculus",
excel=glue::glue("excel/{rundate}mm_ror_gpfoiler_down-v{ver}.xlsx"))
$pvalue_plots$mfp_plot_over
ror_gprofiler_down$pvalue_plots$bpp_plot_over
ror_gprofiler_down$pvalue_plots$reactome_plot_over
ror_gprofiler_down$pvalue_plots$ccp_plot_over
ror_gprofiler_down$pvalue_plots$tf_plot_over ror_gprofiler_down
::pander(sessionInfo())
pandermessage(paste0("This is hpgltools commit: ", get_git_commit()))
paste0(gsub(pattern="\\.Rmd", replace="", x=rmd_file), "-v", ver, ".rda.xz")
this_save <-message(paste0("Saving to ", this_save))
sm(saveme(filename=this_save)) tmp <-
loadme(filename=this_save)