def METHOD_NAME(maf, missingRate):
    bb = homAlt(maf)
    aa = homRef(maf)
    gps = []
    for j in range(nSamples):
        gt = random.random()
        missing = random.random()
        if missing < missingRate:
            gps += [0, 0, 0]
        else:
            d1 = 1.0 - random.uniform(0, 0.01)
            d2 = random.uniform(0, 1.0 - d1)
            d3 = 1.0 - d1 - d2
            if gt < aa:
                gps += [d1, d2, d3]
            elif gt >= aa and gt <= 1.0 - bb:
                gps += [d2, d1, d3]
            else:
                gps += [d3, d2, d1]
    return gps