Provided by: bitseq_0.7.5+dfsg-3ubuntu1_amd64 bug


       estimateExpression  -  estimate  expression  given precomputed probabilities of (observed)
       reads alignments


       estimateExpression -o <outFilePrefix>  [OPTIONS] [prob file]


       Estimates expression given precomputed probabilities of (observed) reads' alignments.

              Uses MCMC sampling algorithm to produce relative abundance or RPKM.



              Show this help information.


              Length of sampler's burn in period. (default: 1000)


              Number of parallel chains used. At least two chains will be used. (default: 4)


              Alpha parameter for the Dirichlet distribution. (default: 1)


              Produce maximum number of samples  (samplesNmax)  in  second  iteration  and  quit.
              (default: Off)


              Initial number of samples produced. Doubles after every iteration. (default: 1000)


              Maximum  number  of  samples produced in one iteration. After producing samplesNmax
              samples sampler finishes. (default: 50000)


              Number of samples recorder in total. (default: 1000)


              Target scale reduction, sampler finishes after this value is met. (default: 1.2)

       -G ,   --gibbs

              Use Gibbs sampling instead of collapsed Gibbs sampling. (default: Off)

       -o <outFilePrefix> ,   --outPrefix=<outFilePrefix>

              Prefix for the output files.

       -O <outputType> ,   --outType=<outputType>

              Output type (theta, RPKM, counts, tau). (default: theta)

       -p <parFileName> ,   --parFile=<parFileName>

              File containing parameters for the sampler, which can  be  otherwise  specified  by
              --MCMC*  options. As the file is checked after every MCMC iteration, the parameters
              can be adjusted while running.

       -P <procN> ,   --procN=<procN>

              Limit the maximum number of threads to be used. (Default  is  the  number  of  MCMC


              Use  scale  reduction  as stopping criterion, instead of computing effective sample
              size. (default: Off)

       -s <seed> ,   --seed=<seed>

              Random initialization seed.


              File for logging noise parameter theta^{act}.

       -t <trInfoFileName> ,   --trInfoFile=<trInfoFileName>

              File containing transcript information. (Necessary for RPKM)

       -v ,   --verbose

              Verbose output. (default: Off)