np.array([np.array([run_single_regression(inputs) for  x in xrange(iterations)])]).transpose()

def Monte_Carlo_regressions(filename, iterations, do_plot = False):
    inputs = data_assignment_regression(filename)
    total_pow, total_energy = np.array([np.array([run_single_regression(inputs) for  x in xrange(iterations)])]).transpose()
    if do_plot:
        plot(total_pow, 'Total Power Capacity (GW)')
        plot(total_energy, 'Total Energy Storage Capacity (TWh)')
    return total_pow.mean(0), total_pow.std(0), total_energy.mean(0), total_energy.std(0)