コード例 #1
0
        
        if max(row) > 0:
            psf_matias_single = temp_psf
            del indices_single, new_data, temp_psf, row
            break
        
    del best_cells   
    ###################     
     
    #for chi2
    expected = sample_size*phi_tilde_pdf*.01
    
    ################

    
    z_dist = df.max_finder(psf_base_single,phi_tilde_point)
    
    base_point_dist.append(z_dist)
    
    ################
    
    z_dist, c_val = df.integral_calc(initial_guess,psf_base_t,phi_tilde_pdf)
    z_temp = z_dist
    
    for j in xrange(max_iters):
        #expected = sample_size*phi_tilde_pdf*.01
        observed = sample_size*c_val*.01
        chisq_stat = chisquare(observed,expected)[0]
        
        if chisq_stat < 233:
            base_pdf_dist.append(z_temp)