CMD: R --no-save --no-restore --no-site-file --no-init-file -q < __tmp_runTMM.R > library(edgeR) > > rnaseqMatrix = read.table("/home/jforment/biovice/internal_projects/200120_mblazquez_rnaseq/04-diffexpr.Ler/all.htseq-count.unstranded.tsv", header=T, row.names=1, com='') > rnaseqMatrix = round(rnaseqMatrix) > exp_study = DGEList(counts=rnaseqMatrix, group=factor(colnames(rnaseqMatrix))) > exp_study = calcNormFactors(exp_study) > exp_study$samples$eff.lib.size = exp_study$samples$lib.size * exp_study$samples$norm.factors > write.table(exp_study$samples, file="/home/jforment/biovice/internal_projects/200120_mblazquez_rnaseq/04-diffexpr.Ler/all.htseq-count.unstranded.tsv.TMM_info.txt", quote=F, sep="\t", row.names=F) >