experimentalFolder = "~/ArmModelPython/Data/ExperimentalHits/" experiment1 = np.loadtxt(experimentalFolder + "hits0.005")[:, 1] experiment2 = np.loadtxt(experimentalFolder + "hits0.01")[:, 1] experiment3 = np.loadtxt(experimentalFolder + "hits0.02")[:, 1] experiment4 = np.loadtxt(experimentalFolder + "hits0.04")[:, 1] print(experiment1.shape[0]) resultPath = "~/ArmModelPython/Data/CMAESK20/" resultName = "50/finalX/x.last" result1 = np.loadtxt(resultPath + "0.005/" + resultName) result2 = np.loadtxt(resultPath + "0.01/" + resultName) result3 = np.loadtxt(resultPath + "0.02/" + resultName) result4 = np.loadtxt(resultPath + "0.04/" + resultName) print("chi2, taille 0.005 : " + str(cp.chi2(result1, experiment1))) print("KS, taille 0.005 : " + str(cp.ks(result1, experiment1))) print("MWWRankSum, taille 0.005 : " + str(cp.MWWRankSum(result1, experiment1))) print("KL, taille 0.005 : " + str(cp.kl(result1, experiment1))) print cp.mw(np.random.normal(0, 0.005, 10000), np.random.normal(0, 0.005, 10000)) print cp.ks(np.random.normal(0, 0.005, 10000), np.random.normal(0, 0.005, 10000)) print( cp.kl(np.random.normal(0, 0.005, 10000), np.random.normal(0, 0.005, 10000)))
import Utils.Compare as cp import numpy as np import matplotlib.pyplot as plt experimentalFolder="~/ArmModelPython/Data/ExperimentalHits/" experiment1 = np.loadtxt(experimentalFolder+"hits0.005")[:,1] experiment2 = np.loadtxt(experimentalFolder+"hits0.01")[:,1] experiment3 = np.loadtxt(experimentalFolder+"hits0.02")[:,1] experiment4 = np.loadtxt(experimentalFolder+"hits0.04")[:,1] print(experiment1.shape[0]) resultPath="~/ArmModelPython/Data/CMAESK20/" resultName="50/finalX/x.last" result1= np.loadtxt(resultPath+"0.005/"+resultName) result2= np.loadtxt(resultPath+"0.01/"+resultName) result3= np.loadtxt(resultPath+"0.02/"+resultName) result4= np.loadtxt(resultPath+"0.04/"+resultName) print ("chi2, taille 0.005 : " + str(cp.chi2(result1, experiment1))) print ("KS, taille 0.005 : " + str(cp.ks(result1, experiment1))) print ("MWWRankSum, taille 0.005 : " + str(cp.MWWRankSum(result1, experiment1))) print ("KL, taille 0.005 : " + str(cp.kl(result1, experiment1))) print cp.mw(np.random.normal(0,0.005,10000),np.random.normal(0,0.005,10000)) print cp.ks(np.random.normal(0,0.005,10000),np.random.normal(0,0.005,10000)) print (cp.kl(np.random.normal(0,0.005,10000),np.random.normal(0,0.005,10000)))