code: def coda_output(pymc_object, name=None, chain=-1):
      
  
      print_()
      print_("Generating CODA output")
      print_('=' * 50)
  
      if name is None:
          name = pymc_object.__name__
  
      # Open trace file
      trace_file = open(name + '_coda.out', 'w')
  
      # Open index file
      index_file = open(name + '_coda.ind', 'w')
  
      variables = [pymc_object]
      if hasattr(pymc_object, 'stochastics'):
          variables = pymc_object.stochastics
  
      # Initialize index
      index = 1
  
      # Loop over all parameters
      for v in variables:
  
          vname = v.__name__
          print_("Processing", vname)
          try:
              index = _process_trace(
                  trace_file,
                  index_file,
                  v.trace(chain=chain),
                  vname,
                  index)
          except TypeError:
              pass
  
      # Close files
      trace_file.close()
      index_file.close()
docs: Generate output files that are compatible with CODA
  
      :Arguments:
        pymc_object : Model or Node
            A PyMC object containing MCMC output.