def METHOD_NAME(file_name, obs_type="amplitudes", parent=None,
    n_bins=40, out=sys.stdout):
  from iotbx import file_reader
  hkl_in = file_reader.any_file(file_name, force_type="hkl")
  hkl_in.check_file_type("hkl")
  anom_data = []
  for array in hkl_in.file_server.miller_arrays :
    if (not array.anomalous_flag()) : continue
    if (array.is_xray_amplitude_array()) or (array.is_xray_intensity_array()):
      anom_data.append(array)
  if (len(anom_data) == 0):
    raise Sorry("No anomalous data arrays found.")
  x_label = "D_anom(F) / <F>"
  if (obs_type == "intensities"):
    x_label = "D_anom(I) / <I>"
  for array in anom_data :
    array_name = os.path.basename(file_name)+":"+array.info().label_string()
    frame = AnomHistPlotFrame(
      parent=parent,
      title="Bijvoet ratios for %s" % array_name,
      array=array,
      array_name=array_name)
    frame.Show()