code: def simulate_genotypes(
      minor_allele_frequency, n_samples, sample_list=sample_id_list
  ):
      """
      given an array that contains the minor_allele_frequency as the first element,
      return a genotypes struct of length=n_samples that conforms to the Glow variant schema,
      with genotypes that are in Hardy Weinberg Equilibrium
      """
      n_samples = int(n_samples)
      frequencies = hardy_weinberg_principle(minor_allele_frequency[0])
      calls = [[0, 0], [0, 1], [1, 1]]
      if len(minor_allele_frequency) > 1:
          genotype_list = generate_multiallelic_frequencies(
              minor_allele_frequency, n_samples
          )
      else:
          genotype_list = random.choices(
              calls, k=n_samples, weights=frequencies
          )
      new_lst = [list(x) for x in zip(sample_id_list, genotype_list)]
      genotypes = [{"sampleId": x, "calls": y} for x, y in new_lst]
      return genotypes
  
label_desc: Use image node_id 6 for calling a global function with example usage: simulate_genotypes(minor_allele_frequency, n_samples, sample_list) and returns: genotypes