2. Fitting statistical models

When input data and input parameter file are ready, use the following command line call in the working directory containing the files. Please note that the run time for MCMC sampler is not short, likely tens of minutes for a proper number of iterations (e.g., N SAMPLE=1,000).

For Linux/Mac OS X

Run the MCMC component of PECA program in the directory holding the input files:

> [path_to_peca_[name]]/peca_[name]_bin [peca-input-parameter-file]

After completion of the MCMC sampling, run the Python script to get the summary tables and gene-specific plots

> python3 [path_to_peca_[name]]/peca_[name].py

Following this step, run the GSA analysis on the PECA output:

> [path_to_peca_gsa]/gsa_bin [gsa-input-parameter-file]

For 64-bit Windows

Run the MCMC component of PECA program:

> [path_to_win64bin]\peca_[name].exe [peca-input-parameter-file]

After completion of the MCMC sampling, run the Python script to get the summary tables and gene-specific plots. Optionally, run the GSA analysis on the PECA output:

> [path_to_win64bin]\gsa.exe [gsa-input-parameter-file]