open Phylogenetics

module Site_independent_mutsel : sig

  type simulation = {
    sequences : Dna.t list ;
    fitness_profiles : Amino_acid.vector array ;
    tree : Convergence_tree.t ;
    h0_profiles : int array ;
    ha_profiles : (int * int) array ;
    h0_params : (Mutsel.param * Mutsel.param) array ;
    ha_params : (Mutsel.param * Mutsel.param) array ;
    base_param : Mutsel.param ;
    gBGC : float * float ;
    ne_s : float * float ;
  }

  val simulation :
    ?branch_scale:float ->
    ?seed:int ->
    n_h0:int ->
    n_ha:int ->
    ne_s:float * float ->
    gBGC:float * float ->
    tree:Convergence_tree.t ->
    fitness_profiles:Amino_acid.vector array ->
    unit ->
    simulation
end