def compare_protocols(self, pathlist, xticks):
"""
Compare the results of different protocols
:param pathlist: list of protocol subdirectories to compare \ example: ["steps/3", "steps/100", "steps/comb"]
:param xticks: list of protocol names will show on the axe of the plot \ example: ["3ms", "100ms", "scomb"]
:return: save the plots of different protocol comparison statistics
"""
check_directory(self.rootdir+"Results")
pathlist = [(self.rootdir+item) for item in pathlist]
protocol_comparison(pathlist, self.nfp, self.nrep, self.pnames, self.rootdir+"Results", self.pinit, xticks)
# Create combinations
combine_likelihood(zaps_list, numfp=10, num_mult_single=10,
out_path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/comb")
combine_likelihood(steps_list, numfp=10, num_mult_single=10,
out_path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/comb")
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/100", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/10", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/1", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/comb", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/3", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/20", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/200", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/comb", 10, dbs=pinit)
path_list = ["/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/3",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/20",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/200",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/1",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/10",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/100",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/comb",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/comb"]
protocol_comparison(path_list, numfp=10, inferred_params=['Ra', 'cm', 'gpas'],
out_path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4", dbs=pinit)
pinit.close()
runningTime = (time.time()-startTime)/60
print("\n\nThe script was running for %f minutes" % runningTime)
combine_likelihood(sins_list, numfp=10, num_mult_single=10,
out_path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/comb")
combine_likelihood(steps_list, numfp=10, num_mult_single=10,
out_path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/comb")
plot_combined_results("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/100", 10, dbs=pinit)
plot_combined_results("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/10", 10, dbs=pinit)
plot_combined_results("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/1", 10, dbs=pinit)
plot_combined_results("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/comb", 10, dbs=pinit)
plot_combined_results("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/3", 10, dbs=pinit)
plot_combined_results("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/20", 10, dbs=pinit)
plot_combined_results("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/200", 10, dbs=pinit)
plot_combined_results("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/comb", 10, dbs=pinit)
path_list = ["/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/3",
"/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/20",
"/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/200",
"/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/1",
"/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/10",
"/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/100",
"/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/comb",
"/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/comb"]
protocol_comparison(path_list, numfp=10, inferred_params=p_names,
out_path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens", dbs=pinit)
pinit.close()
runningTime = (time.time()-startTime)/60
print("\n\nThe script was running for %f minutes" % runningTime)
num_mult_single=10,
out_path=".../sins/comb")
combine_likelihood(steps_list,
numfp=10,
num_mult_single=10,
out_path=".../steps/comb")
plot_combined_results(".../sins/100", 10, dbs=pinit)
plot_combined_results(".../sins/10", 10, dbs=pinit)
plot_combined_results(".../sins/1", 10, dbs=pinit)
plot_combined_results(".../sins/comb", 10, dbs=pinit)
plot_combined_results(".../steps/3", 10, dbs=pinit)
plot_combined_results(".../steps/20", 10, dbs=pinit)
plot_combined_results(".../steps/200", 10, dbs=pinit)
plot_combined_results(".../steps/comb", 10, dbs=pinit)
path_list = [
".../steps/3", ".../steps/20", ".../steps/200", ".../sins/1",
".../sins/10", ".../sins/100", ".../steps/comb", ".../sins/comb"
]
protocol_comparison(path_list,
numfp=10,
inferred_params=p_names,
out_path="...",
dbs=pinit)
pinit.close()
runningTime = (time.time() - startTime) / 60
print("\n\nThe script was running for %f minutes" % runningTime)
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/steps/200",
10,
dbs=pinit)
plot_combined_results(
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/steps/comb",
10,
dbs=pinit)
path_list = [
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/steps/3",
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/steps/20",
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/steps/200",
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/1",
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/10",
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/100",
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/steps/comb",
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/comb"
]
protocol_comparison(
path_list,
numfp=10,
inferred_params=p_names,
out_path=
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored",
dbs=pinit)
pinit.close()
runningTime = (time.time() - startTime) / 60
print("\n\nThe script was running for %f minutes" % runningTime)
plot_combined_results(
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/steps/200",
10,
dbs=pinit)
plot_combined_results(
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/steps/comb",
10,
dbs=pinit)
path_list = [
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/steps/3",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/steps/20",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/steps/200",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/sins/1",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/sins/10",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/sins/100",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/steps/comb",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug/sins/comb"
]
protocol_comparison(
path_list,
numfp=10,
inferred_params=p_names,
out_path=
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored/debug",
dbs=pinit)
pinit.close()
runningTime = (time.time() - startTime) / 60
print("\n\nThe script was running for %f minutes" % runningTime)
