L. major differential expression.

index.html preprocessing.html

1 Differential Expression version: 20170706

2 L. major differential expression

lm_filt <- normalize_expt(lm_expt, filter=TRUE)

## This function will replace the expt$expressionset slot with:

## cbcb(data)

## It backs up the current data into a slot named:
##  expt$backup_expressionset. It will also save copies of each step along the way
##  in expt$normalized with the corresponding libsizes. Keep the libsizes in mind
##  when invoking limma.  The appropriate libsize is the 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 66 low-count genes (8232 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.

lm_pairwise <- all_pairwise(input=lm_filt, model_batch="svaseq")

## The be method chose 6 surrogate variable(s).

## This ignores the surrogates parameter and uses the be method to estimate surrogates.

## Using svaseq to test before/after batch correction.

## Finished running DE analyses, collecting outputs.

## Comparing analyses 1/15: pmn_inf_12hr_vs_amastigote

## Comparing analyses 2/15: pmn_inf_14d_vs_amastigote

## Comparing analyses 3/15: pmn_uninf_12hr_vs_amastigote

## Comparing analyses 4/15: pmn_uninf_14d_vs_amastigote

## Comparing analyses 5/15: promastigote_vs_amastigote

## Comparing analyses 6/15: pmn_inf_14d_vs_pmn_inf_12hr

## Comparing analyses 7/15: pmn_uninf_12hr_vs_pmn_inf_12hr

## Comparing analyses 8/15: pmn_uninf_14d_vs_pmn_inf_12hr

## Comparing analyses 9/15: promastigote_vs_pmn_inf_12hr

## Comparing analyses 10/15: pmn_uninf_12hr_vs_pmn_inf_14d

## Comparing analyses 11/15: pmn_uninf_14d_vs_pmn_inf_14d

## Comparing analyses 12/15: promastigote_vs_pmn_inf_14d

## Comparing analyses 13/15: pmn_uninf_14d_vs_pmn_uninf_12hr

## Comparing analyses 14/15: promastigote_vs_pmn_uninf_12hr

## Comparing analyses 15/15: promastigote_vs_pmn_uninf_14d

keepers <- list(
    "12hr_infuninf" = c("pmn_inf_12hr", "pmn_uninf_12hr"),
    "14d_infuninf" = c("pmn_inf_14d", "pmn_uninf_14d"),
    "ama_pro" = c("amastigote", "promastigote"),
    "inf_14d12hr" = c("pmn_inf_14d", "pmn_inf_12hr"),
    "uninf_14d12hr" = c("pmn_uninf_14d", "pmn_uninf_12hr")
)
lm_tables <- combine_de_tables(lm_pairwise, keepers=keepers,
                               excel=paste0("excel/lm_pairwise-v", ver, ".xlsx"))

## Deleting the file excel/lm_pairwise-v20170706.xlsx before writing the tables.

## Writing a legend of columns.

## Printing a pca plot before/after surrogates/batch estimation.

## Working on 1/5: 12hr_infuninf

## Found inverse table with pmn_uninf_12hr_vs_pmn_inf_12hr

## Working on 2/5: 14d_infuninf

## Found inverse table with pmn_uninf_14d_vs_pmn_inf_14d

## Working on 3/5: ama_pro

## Found inverse table with promastigote_vs_amastigote

## Working on 4/5: inf_14d12hr

## Found table with pmn_inf_14d_vs_pmn_inf_12hr

## Working on 5/5: uninf_14d12hr

## Found table with pmn_uninf_14d_vs_pmn_uninf_12hr

## Adding venn plots for 12hr_infuninf.

## Limma expression coefficients for 12hr_infuninf; R^2: 0.437; equation: y = 0.721x + 1.56

## Edger expression coefficients for 12hr_infuninf; R^2: 0.048; equation: y = -0.0495x + 63.4

## DESeq2 expression coefficients for 12hr_infuninf; R^2: 0.485; equation: y = 0.743x + 1.44

## Adding venn plots for 14d_infuninf.

## Limma expression coefficients for 14d_infuninf; R^2: 0.772; equation: y = 0.891x + 0.654

## Edger expression coefficients for 14d_infuninf; R^2: 0.476; equation: y = 0.612x + 2.77

## DESeq2 expression coefficients for 14d_infuninf; R^2: 0.854; equation: y = 0.864x + 0.649

## Adding venn plots for ama_pro.

## Limma expression coefficients for ama_pro; R^2: 0.628; equation: y = 0.812x + 1.03

## Edger expression coefficients for ama_pro; R^2: 0.788; equation: y = 0.788x + 1.57

## DESeq2 expression coefficients for ama_pro; R^2: 0.796; equation: y = 0.926x + 0.354

## Adding venn plots for inf_14d12hr.

## Limma expression coefficients for inf_14d12hr; R^2: 0.807; equation: y = 0.893x + 0.636

## Edger expression coefficients for inf_14d12hr; R^2: 0.916; equation: y = 0.962x + 0.113

## DESeq2 expression coefficients for inf_14d12hr; R^2: 0.809; equation: y = 0.99x - 0.272

## Adding venn plots for uninf_14d12hr.

## Limma expression coefficients for uninf_14d12hr; R^2: 0.448; equation: y = 0.597x + 2.16

## Edger expression coefficients for uninf_14d12hr; R^2: 0.0063; equation: y = 0.195x + 45.5

## DESeq2 expression coefficients for uninf_14d12hr; R^2: 0.497; equation: y = 0.65x + 0.503

## Writing summary information.

## The sheet: pairwise_summary is in legend, 12hr_infuninf, 14d_infuninf, ama_pro, inf_14d12hr, uninf_14d12hr, pairwise_summary.

## Attempting to add the comparison plot to pairwise_summary at row: 23 and column: 1

## Performing save of the workbook.

lm_sig <- extract_significant_genes(combined=lm_tables,
                                    excel=paste0("excel/lm_sig-v", ver, ".xlsx"))

## Writing a legend of columns.

## Writing excel data sheet 0/5: 12hr_infuninf

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 659 genes.

## After (adj)p filter, the down genes table has 115 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 659 genes.

## After fold change filter, the down genes table has 115 genes.

## Writing excel data sheet 1/5: 14d_infuninf

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 1 genes.

## After (adj)p filter, the down genes table has 2 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 1 genes.

## After fold change filter, the down genes table has 2 genes.

## Writing excel data sheet 2/5: ama_pro

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 636 genes.

## After (adj)p filter, the down genes table has 552 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 621 genes.

## After fold change filter, the down genes table has 532 genes.

## Writing excel data sheet 3/5: inf_14d12hr

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 2 genes.

## After (adj)p filter, the down genes table has 1 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 2 genes.

## After fold change filter, the down genes table has 1 genes.

## Writing excel data sheet 4/5: uninf_14d12hr

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 356 genes.

## After (adj)p filter, the down genes table has 66 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 356 genes.

## After fold change filter, the down genes table has 66 genes.

## Printing significant genes to the file: excel/lm_sig-v20170706.xlsx

## 1/5: Writing excel data sheet up_1limma_12hr_infuninf

## 1/5: Writing excel data sheet down_1limma_12hr_infuninf

## 2/5: Writing excel data sheet up_2limma_14d_infuninf

## 2/5: Writing excel data sheet down_2limma_14d_infuninf

## 3/5: Writing excel data sheet up_3limma_ama_pro

## 3/5: Writing excel data sheet down_3limma_ama_pro

## 4/5: Writing excel data sheet up_4limma_inf_14d12hr

## 4/5: Writing excel data sheet down_4limma_inf_14d12hr

## 5/5: Writing excel data sheet up_5limma_uninf_14d12hr

## 5/5: Writing excel data sheet down_5limma_uninf_14d12hr

## Writing excel data sheet 5/5: 12hr_infuninf

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 2721 genes.

## After (adj)p filter, the down genes table has 126 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 2721 genes.

## After fold change filter, the down genes table has 125 genes.

## Writing excel data sheet 6/5: 14d_infuninf

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 235 genes.

## After (adj)p filter, the down genes table has 238 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 234 genes.

## After fold change filter, the down genes table has 238 genes.

## Writing excel data sheet 7/5: ama_pro

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 128 genes.

## After (adj)p filter, the down genes table has 159 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 121 genes.

## After fold change filter, the down genes table has 152 genes.

## Writing excel data sheet 8/5: inf_14d12hr

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 40 genes.

## After (adj)p filter, the down genes table has 0 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 40 genes.

## After fold change filter, the down genes table has 0 genes.

## Writing excel data sheet 9/5: uninf_14d12hr

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 2746 genes.

## After (adj)p filter, the down genes table has 119 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 2746 genes.

## After fold change filter, the down genes table has 119 genes.

## Printing significant genes to the file: excel/lm_sig-v20170706.xlsx

## The sheet: number_changed is in legend, number_changed, up_1limma_12hr_infuninf, down_1limma_12hr_infuninf, up_2limma_14d_infuninf, down_2limma_14d_infuninf, up_3limma_ama_pro, down_3limma_ama_pro, up_4limma_inf_14d12hr, down_4limma_inf_14d12hr, up_5limma_uninf_14d12hr, down_5limma_uninf_14d12hr.

## 1/5: Writing excel data sheet up_1edger_12hr_infuninf

## 1/5: Writing excel data sheet down_1edger_12hr_infuninf

## 2/5: Writing excel data sheet up_2edger_14d_infuninf

## 2/5: Writing excel data sheet down_2edger_14d_infuninf

## 3/5: Writing excel data sheet up_3edger_ama_pro

## 3/5: Writing excel data sheet down_3edger_ama_pro

## 4/5: Writing excel data sheet up_4edger_inf_14d12hr

## 4/5: Writing excel data sheet down_4edger_inf_14d12hr

## 5/5: Writing excel data sheet up_5edger_uninf_14d12hr

## 5/5: Writing excel data sheet down_5edger_uninf_14d12hr

## Writing excel data sheet 10/5: 12hr_infuninf

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 992 genes.

## After (adj)p filter, the down genes table has 141 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 991 genes.

## After fold change filter, the down genes table has 141 genes.

## Writing excel data sheet 11/5: 14d_infuninf

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 3 genes.

## After (adj)p filter, the down genes table has 0 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 3 genes.

## After fold change filter, the down genes table has 0 genes.

## Writing excel data sheet 12/5: ama_pro

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 96 genes.

## After (adj)p filter, the down genes table has 64 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 96 genes.

## After fold change filter, the down genes table has 64 genes.

## Writing excel data sheet 13/5: inf_14d12hr

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 1 genes.

## After (adj)p filter, the down genes table has 0 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 1 genes.

## After fold change filter, the down genes table has 0 genes.

## Writing excel data sheet 14/5: uninf_14d12hr

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 1163 genes.

## After (adj)p filter, the down genes table has 102 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 1162 genes.

## After fold change filter, the down genes table has 102 genes.

## Printing significant genes to the file: excel/lm_sig-v20170706.xlsx

## The sheet: number_changed is in legend, number_changed, up_1limma_12hr_infuninf, down_1limma_12hr_infuninf, up_2limma_14d_infuninf, down_2limma_14d_infuninf, up_3limma_ama_pro, down_3limma_ama_pro, up_4limma_inf_14d12hr, down_4limma_inf_14d12hr, up_5limma_uninf_14d12hr, down_5limma_uninf_14d12hr, up_1edger_12hr_infuninf, down_1edger_12hr_infuninf, up_2edger_14d_infuninf, down_2edger_14d_infuninf, up_3edger_ama_pro, down_3edger_ama_pro, up_4edger_inf_14d12hr, down_4edger_inf_14d12hr, up_5edger_uninf_14d12hr, down_5edger_uninf_14d12hr.

## 1/5: Writing excel data sheet up_1deseq_12hr_infuninf

## 1/5: Writing excel data sheet down_1deseq_12hr_infuninf

## 2/5: Writing excel data sheet up_2deseq_14d_infuninf

## 2/5: Writing excel data sheet down_2deseq_14d_infuninf

## 3/5: Writing excel data sheet up_3deseq_ama_pro

## 3/5: Writing excel data sheet down_3deseq_ama_pro

## 4/5: Writing excel data sheet up_4deseq_inf_14d12hr

## 4/5: Writing excel data sheet down_4deseq_inf_14d12hr

## 5/5: Writing excel data sheet up_5deseq_uninf_14d12hr

## 5/5: Writing excel data sheet down_5deseq_uninf_14d12hr

## Writing excel data sheet 15/5: 12hr_infuninf

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 0 genes.

## After (adj)p filter, the down genes table has 0 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 0 genes.

## After fold change filter, the down genes table has 0 genes.

## Writing excel data sheet 16/5: 14d_infuninf

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 1 genes.

## After (adj)p filter, the down genes table has 9 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 0 genes.

## After fold change filter, the down genes table has 9 genes.

## Writing excel data sheet 17/5: ama_pro

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 215 genes.

## After (adj)p filter, the down genes table has 226 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 176 genes.

## After fold change filter, the down genes table has 191 genes.

## Writing excel data sheet 18/5: inf_14d12hr

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 0 genes.

## After (adj)p filter, the down genes table has 0 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 0 genes.

## After fold change filter, the down genes table has 0 genes.

## Writing excel data sheet 19/5: uninf_14d12hr

## Assuming the fold changes are on the log scale and so taking >< 0

## After (adj)p filter, the up genes table has 89 genes.

## After (adj)p filter, the down genes table has 214 genes.

## Assuming the fold changes are on the log scale and so taking -1.0 * fc

## After fold change filter, the up genes table has 89 genes.

## After fold change filter, the down genes table has 214 genes.

## Printing significant genes to the file: excel/lm_sig-v20170706.xlsx

## The sheet: number_changed is in legend, number_changed, up_1limma_12hr_infuninf, down_1limma_12hr_infuninf, up_2limma_14d_infuninf, down_2limma_14d_infuninf, up_3limma_ama_pro, down_3limma_ama_pro, up_4limma_inf_14d12hr, down_4limma_inf_14d12hr, up_5limma_uninf_14d12hr, down_5limma_uninf_14d12hr, up_1edger_12hr_infuninf, down_1edger_12hr_infuninf, up_2edger_14d_infuninf, down_2edger_14d_infuninf, up_3edger_ama_pro, down_3edger_ama_pro, up_4edger_inf_14d12hr, down_4edger_inf_14d12hr, up_5edger_uninf_14d12hr, down_5edger_uninf_14d12hr, up_1deseq_12hr_infuninf, down_1deseq_12hr_infuninf, up_2deseq_14d_infuninf, down_2deseq_14d_infuninf, up_3deseq_ama_pro, down_3deseq_ama_pro, up_4deseq_inf_14d12hr, down_4deseq_inf_14d12hr, up_5deseq_uninf_14d12hr, down_5deseq_uninf_14d12hr.

## 1/5: Writing excel data sheet up_1basic_12hr_infuninf

## 1/5: Writing excel data sheet down_1basic_12hr_infuninf

## 2/5: Writing excel data sheet up_2basic_14d_infuninf

## 2/5: Writing excel data sheet down_2basic_14d_infuninf

## 3/5: Writing excel data sheet up_3basic_ama_pro

## 3/5: Writing excel data sheet down_3basic_ama_pro

## 4/5: Writing excel data sheet up_4basic_inf_14d12hr

## 4/5: Writing excel data sheet down_4basic_inf_14d12hr

## 5/5: Writing excel data sheet up_5basic_uninf_14d12hr

## 5/5: Writing excel data sheet down_5basic_uninf_14d12hr

## Adding significance bar plots.

## The sheet: number_changed is in legend, number_changed, up_1limma_12hr_infuninf, down_1limma_12hr_infuninf, up_2limma_14d_infuninf, down_2limma_14d_infuninf, up_3limma_ama_pro, down_3limma_ama_pro, up_4limma_inf_14d12hr, down_4limma_inf_14d12hr, up_5limma_uninf_14d12hr, down_5limma_uninf_14d12hr, up_1edger_12hr_infuninf, down_1edger_12hr_infuninf, up_2edger_14d_infuninf, down_2edger_14d_infuninf, up_3edger_ama_pro, down_3edger_ama_pro, up_4edger_inf_14d12hr, down_4edger_inf_14d12hr, up_5edger_uninf_14d12hr, down_5edger_uninf_14d12hr, up_1deseq_12hr_infuninf, down_1deseq_12hr_infuninf, up_2deseq_14d_infuninf, down_2deseq_14d_infuninf, up_3deseq_ama_pro, down_3deseq_ama_pro, up_4deseq_inf_14d12hr, down_4deseq_inf_14d12hr, up_5deseq_uninf_14d12hr, down_5deseq_uninf_14d12hr, up_1basic_12hr_infuninf, down_1basic_12hr_infuninf, up_2basic_14d_infuninf, down_2basic_14d_infuninf, up_3basic_ama_pro, down_3basic_ama_pro, up_4basic_inf_14d12hr, down_4basic_inf_14d12hr, up_5basic_uninf_14d12hr, down_5basic_uninf_14d12hr.
## The sheet: number_changed is in legend, number_changed, up_1limma_12hr_infuninf, down_1limma_12hr_infuninf, up_2limma_14d_infuninf, down_2limma_14d_infuninf, up_3limma_ama_pro, down_3limma_ama_pro, up_4limma_inf_14d12hr, down_4limma_inf_14d12hr, up_5limma_uninf_14d12hr, down_5limma_uninf_14d12hr, up_1edger_12hr_infuninf, down_1edger_12hr_infuninf, up_2edger_14d_infuninf, down_2edger_14d_infuninf, up_3edger_ama_pro, down_3edger_ama_pro, up_4edger_inf_14d12hr, down_4edger_inf_14d12hr, up_5edger_uninf_14d12hr, down_5edger_uninf_14d12hr, up_1deseq_12hr_infuninf, down_1deseq_12hr_infuninf, up_2deseq_14d_infuninf, down_2deseq_14d_infuninf, up_3deseq_ama_pro, down_3deseq_ama_pro, up_4deseq_inf_14d12hr, down_4deseq_inf_14d12hr, up_5deseq_uninf_14d12hr, down_5deseq_uninf_14d12hr, up_1basic_12hr_infuninf, down_1basic_12hr_infuninf, up_2basic_14d_infuninf, down_2basic_14d_infuninf, up_3basic_ama_pro, down_3basic_ama_pro, up_4basic_inf_14d12hr, down_4basic_inf_14d12hr, up_5basic_uninf_14d12hr, down_5basic_uninf_14d12hr.
## The sheet: number_changed is in legend, number_changed, up_1limma_12hr_infuninf, down_1limma_12hr_infuninf, up_2limma_14d_infuninf, down_2limma_14d_infuninf, up_3limma_ama_pro, down_3limma_ama_pro, up_4limma_inf_14d12hr, down_4limma_inf_14d12hr, up_5limma_uninf_14d12hr, down_5limma_uninf_14d12hr, up_1edger_12hr_infuninf, down_1edger_12hr_infuninf, up_2edger_14d_infuninf, down_2edger_14d_infuninf, up_3edger_ama_pro, down_3edger_ama_pro, up_4edger_inf_14d12hr, down_4edger_inf_14d12hr, up_5edger_uninf_14d12hr, down_5edger_uninf_14d12hr, up_1deseq_12hr_infuninf, down_1deseq_12hr_infuninf, up_2deseq_14d_infuninf, down_2deseq_14d_infuninf, up_3deseq_ama_pro, down_3deseq_ama_pro, up_4deseq_inf_14d12hr, down_4deseq_inf_14d12hr, up_5deseq_uninf_14d12hr, down_5deseq_uninf_14d12hr, up_1basic_12hr_infuninf, down_1basic_12hr_infuninf, up_2basic_14d_infuninf, down_2basic_14d_infuninf, up_3basic_ama_pro, down_3basic_ama_pro, up_4basic_inf_14d12hr, down_4basic_inf_14d12hr, up_5basic_uninf_14d12hr, down_5basic_uninf_14d12hr.

pander::pander(sessionInfo())

R version 3.3.3 (2017-03-06)

**Platform:** x86_64-pc-linux-gnu (64-bit)

locale: LC_CTYPE=en_US.utf8, LC_NUMERIC=C, LC_TIME=en_US.utf8, LC_COLLATE=en_US.utf8, LC_MONETARY=en_US.utf8, LC_MESSAGES=en_US.utf8, LC_PAPER=en_US.utf8, LC_NAME=C, LC_ADDRESS=C, LC_TELEPHONE=C, LC_MEASUREMENT=en_US.utf8 and LC_IDENTIFICATION=C

attached base packages: stats, graphics, grDevices, utils, datasets, methods and base

other attached packages: Vennerable(v.3.1.0.9000), ruv(v.0.9.6) and hpgltools(v.2017.01)

loaded via a namespace (and not attached): nlme(v.3.1-131), bitops(v.1.0-6), devtools(v.1.13.2), bit64(v.0.9-7), doParallel(v.1.0.10), RColorBrewer(v.1.1-2), rprojroot(v.1.2), GenomeInfoDb(v.1.10.3), tools(v.3.3.3), backports(v.1.1.0), R6(v.2.2.2), rpart(v.4.1-11), Hmisc(v.4.0-3), DBI(v.0.7), lazyeval(v.0.2.0), BiocGenerics(v.0.20.0), mgcv(v.1.8-17), colorspace(v.1.3-2), nnet(v.7.3-12), withr(v.1.0.2), gridExtra(v.2.2.1), DESeq2(v.1.14.1), bit(v.1.1-12), compiler(v.3.3.3), preprocessCore(v.1.36.0), graph(v.1.52.0), Biobase(v.2.34.0), htmlTable(v.1.9), xml2(v.1.1.1), labeling(v.0.3), scales(v.0.4.1), checkmate(v.1.8.3), DEoptimR(v.1.0-8), robustbase(v.0.92-7), genefilter(v.1.56.0), RBGL(v.1.50.0), commonmark(v.1.2), stringr(v.1.2.0), digest(v.0.6.12), foreign(v.0.8-69), rmarkdown(v.1.6), XVector(v.0.14.1), base64enc(v.0.1-3), htmltools(v.0.3.6), limma(v.3.30.13), htmlwidgets(v.0.8), rlang(v.0.1.1), RSQLite(v.2.0), BiocParallel(v.1.8.2), gtools(v.3.5.0), acepack(v.1.4.1), RCurl(v.1.95-4.8), magrittr(v.1.5), Formula(v.1.2-1), Matrix(v.1.2-10), Rcpp(v.0.12.11), munsell(v.0.4.3), S4Vectors(v.0.12.2), stringi(v.1.1.5), yaml(v.2.1.14), edgeR(v.3.16.5), SummarizedExperiment(v.1.4.0), zlibbioc(v.1.20.0), plyr(v.1.8.4), grid(v.3.3.3), blob(v.1.1.0), parallel(v.3.3.3), ggrepel(v.0.6.5), crayon(v.1.3.2), lattice(v.0.20-35), splines(v.3.3.3), pander(v.0.6.0), annotate(v.1.52.1), locfit(v.1.5-9.1), knitr(v.1.16), GenomicRanges(v.1.26.4), corpcor(v.1.6.9), reshape2(v.1.4.2), geneplotter(v.1.52.0), codetools(v.0.2-15), stats4(v.3.3.3), XML(v.3.98-1.9), evaluate(v.0.10.1), latticeExtra(v.0.6-28), data.table(v.1.10.4), foreach(v.1.4.3), testthat(v.1.0.2), gtable(v.0.2.0), ggplot2(v.2.2.1), openxlsx(v.4.0.17), xtable(v.1.8-2), roxygen2(v.6.0.1), survival(v.2.41-3), tibble(v.1.3.3), iterators(v.1.0.8), AnnotationDbi(v.1.36.2), memoise(v.1.1.0), IRanges(v.2.8.2), cluster(v.2.0.6) and sva(v.3.22.0)

this_save <- paste0(gsub(pattern="\\.Rmd", replace="", x=rmd_file), "-v", ver, ".rda.xz")
message(paste0("Saving to ", this_save))

## Saving to 03_differential_expression_lmajor-v20170706.rda.xz

tt <- sm(saveme(filename=this_save))