path_git, path_datos_global = definir_path() os.chdir(path_git) sys.path.append('./Software/Funcionales/') from funciones_sampleo import MCMC_sampler from funciones_data import leer_data_pantheon, leer_data_cronometros, leer_data_BAO, leer_data_AGN from funciones_alternativos import params_to_chi2 #ORDEN DE PRESENTACION DE LOS PARAMETROS: Mabs,omega_m,b,H_0,n #%% ## Supernovas os.chdir(path_git + '/Software/Estadística/Datos/Datos_pantheon/') ds_SN = leer_data_pantheon('lcparam_full_long_zhel.txt') # Cronómetros os.chdir(path_git + '/Software/Estadística/Datos/') ds_CC = leer_data_cronometros('datos_cronometros.txt') # BAO os.chdir(path_git + '/Software/Estadística/Datos/BAO/') ds_BAO = [] archivos_BAO = [ 'datos_BAO_da.txt', 'datos_BAO_dh.txt', 'datos_BAO_dm.txt', 'datos_BAO_dv.txt', 'datos_BAO_H.txt' ] for i in range(5): aux = leer_data_BAO(archivos_BAO[i]) ds_BAO.append(aux) # AGN os.chdir(path_git + '/Software/Estadística/Datos/Datos_AGN') ds_AGN = leer_data_AGN('table3.dat')
M_true = -19.2 omega_m_true = 0.4 b_true = -0.5 H0_true = 73.48 #Unidades de (km/seg)/Mpc n = 1 params_fijos = [H0_true,n] #%% #Datos de SN os.chdir(path_git+'/Software/Estadística/Datos/Datos_pantheon/') zcmb,zhel, Cinv, mb = leer_data_pantheon('lcparam_full_long_zhel.txt') #Datos de crnómetros os.chdir(path_git+'/Software/Estadística/Datos/') z_data, H_data, dH = leer_data_cronometros('datos_cronometros.txt') #%% #Parametros a ajustar nll = lambda theta: params_to_chi2_taylor(theta,params_fijos,zcmb, zhel, Cinv, mb,z_data, H_data, dH,chi_riess=False) initial = np.array([M_true,omega_m_true,b_true]) soln = minimize(nll, initial, options = {'eps': 0.01}, bounds =((-20,-18),(0.1,0.5),(-1, 1))) M_ml, omega_m_ml, b_ml = soln.x print(M_ml,omega_m_ml,b_ml) os.chdir(path_git + '/Software/Estadística/Resultados_simulaciones') np.savez('valores_medios_HS_CC+SN_3params_taylor', sol=soln.x)