Statistical Analysis of Election Manipulation in Republican Primaries(themoneyparty.org)

    #!/usr/bin/python
    import random
    import matplotlib.pyplot as plt

    precinct_size = [1, 1, 1, 3, 5, 7, 8, 9, 10, 12, 18, 19, 21, 26, 28, 29, 30, 30, 34, 41,
    46, 50, 56, 57, 57, 58, 58, 61, 68, 69, 72, 78, 79, 85, 88, 94, 99, 100,
    103, 109, 109, 120, 126, 126, 129, 132, 133, 133, 136, 139, 139, 141, 147, 152,
    157, 158, 162, 162, 165, 166, 169, 172, 173, 175, 177, 177, 179, 181, 192,
    201, 224, 231, 235, 238, 246, 249, 249, 251, 258, 264, 268, 270, 272, 276,
    277, 281, 281, 293, 315, 322, 327, 333, 334, 337, 346, 348, 349, 350, 363, 366,
    367, 369, 370, 374, 374, 384, 385, 386, 387, 394, 397, 405, 407, 412, 413, 419, 420,
    421, 425, 429, 434, 438, 439, 449, 467, 474, 483, 483, 496, 502, 505, 507, 508, 513,
    518, 525, 526, 538, 545, 548, 555, 560, 571, 577, 581, 583, 584, 606, 606, 612, 620,
    620, 625, 635, 641, 646, 647, 650, 652, 662, 663, 666, 669, 673, 674, 710, 711, 721, 728,
    737, 744, 745, 747, 747, 780, 786, 787, 790, 791, 818, 824, 833, 834, 878, 899,
    903, 927, 928, 949, 961, 1002, 1031, 1059, 1070, 1133, 1587]

    vote_dist  = [.446, .314, .160, .080]    # romney, santorum, gingrich, paul
    vote_dist2 = [.296, .364, .210, .130]
    vote_dist3 = [.146, .384, .275, .195]

    def vote(d):
        r = random.random()
        s = 0
        for v in range(len(d)):
            s += d[v]
            if s >= r:
                return v


    def trial():
        total = [0, 0, 0, 0]
        x = []
        romney = []
        santorum = []
        gingrich = []
        paul = []
        i = 1
        vt = 0
        for s in precinct_size:
            for voter in range(s):
                if s < 50:
                    v = vote(vote_dist3)
                elif s < 150:
                    v = vote(vote_dist2)
                else:
                    v = vote(vote_dist)
                total[v] += 1
                vt += 1
            x.append(float(i))
            i += 1
            romney.append(float(total[0])/vt)
            santorum.append(float(total[1])/vt)
            gingrich.append(float(total[2])/vt)
            paul.append(float(total[3])/vt)

        plt.plot(x, paul, 'ro')
        plt.plot(x, gingrich, 'go')
        plt.plot(x, santorum, 'yo')
        plt.plot(x, romney, 'bo')
        plt.show()

    trial()