The default output from LiBiNorm consists of a single text file that lists the expression values for each gene. There are a number of other optional text files that can be produced and which provide more detail about the bias normalisation using -u <fileroot>
Default gene expression output
By default the only output from LiBiNorm is a list of genes and their expression in text format which has the following format and is sent to stdout, mirroring the operation of htseq-count:
The first two columns are identical to those of htseq-count. There then follows the RNA length, the predicted bias relative to a linear model (= reads scale proportional to length, e.g. FPKM) and the expression in TPM after correction for the bias as predicted by the bias model that has been selected.
LiBiNorm can also output the results directly to a file using -c <filename>.
If LiBiNorm is run in the htseq-count compatible mode (-z), then only the first two columns are output and there is no header information.
The full mode option (-f) provides additional columns showing a range of different expression values, with and without bias normalisation.
When the -u <fileroot> option is used three additional files are created which provide more details on the results of the bias normalisation
This provides the best parameters and associated log likelihoods for each of the models, together with an indication of the variation seen within the MCMC chain used to derive the parameters. The endpoints of each of the chains are also listed. The following shows an example for one of the models. The columns show the results for each of the parameters in the model and the (negative) Log Likelihood that is obtained.
The first two rows give the log10 and absolute values for the model parameters. The "Initial" row identifies the values obtained by Nelder Mead that were used as a starting point for the MCMC runs. The next four rows give various measures of the dispersion (median absolute deviation) in the values as determined by the MCMC runs. Only a single measure of the (negative) Log Likelihood is shown. The following rows show the end points of each of the MCMC runs which allow checking agreement between these.
NM iterations indicate how many iterations of the Nelder Mead algorthm were run before the parameters had stabilised sufficiently for LiBiNorm to move to the second, MCMC based stage of parameter determination.
This outputs the predicted bias for RNA lengths between 200 bp and 20000 bp for each of the 6 models. Parameter estimates are given for each model, followed by two rows indicating the predicted bias at 100 bp intervals between 100 and 20000 bp. A length of 1000 bp is used as the reference.
This shows the bias in the underlying read distribution as the raw data for a heat map that can then be plotted using LiBiNormPlot.R
This file mirrors the counts data that is produced by LiBiNorm