QCYCLE SIMULATION DOCUMENTS


FORTRAN 77 program files

Mainline: read input files; set up simulation on Beowulf cluster; get chi^2 for each Latin hypercube sample; write results.  Download: sample15a.f.

Function: return chi^2 value.  Download: funk15a.f.

Subroutine: calculate new rate set based on old rate set and vector of multipliers.  Download: newrates15a.f.

Subroutine: carry out Monte Carlo simulation.  Download: qcycle15as.f.


Input files used by simulation whose results are shown in Fig. 13

Input parameter file.  Download: sampleIn15a.txt.

Base rates file.  Download qcycleRates15a.txt.

Latin hypercube samples file (8192 rows with 22 random numbers per row).  Download sample15a.txt.

Data from Kramer and Crofts (1993); Fig. 3A, cyt b data.  Download cut3ab3.txt.

Data from Kramer and Crofts (1993); Fig. 3B, cyt b data.  Download cut3bb3.txt.

Data from Kramer and Crofts (1993); Fig. 3A, cyt f data.  Download qcycle/cut3af3.txt.

Data from Kramer and Crofts (1993); Fig. 3B, cyt f data.  Download cut3bf3.txt.


Results:

chi^2 results corresponding to Latin hypercube samples. Download chi2.txt.

New rates file corresponding to smallest overall chi^2. Download qcycleRates15aNew.txt.