def likelihood_mult(self):
"""
Multiplies likelihoods for all repetition at each fixed parameters and plots some result
:return: saves and plots the multiplied likelihoods into disk
"""
for protocol in self.subdirs:
mult_likelihood(path=(self.rootdir+protocol), numfp=self.nfp, num_mult=self.nrep)
plot_combined_results(self.rootdir+protocol, self.nfp, dbs=self.pinit)
nfp = 10 # Number of fixed parameters
# Load parameter space initializator
pinit = tb.open_file("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/paramsetup.hdf5", mode="r")
# Plot some single results
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/3", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/20", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/200", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/1", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/10", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/100", numfp=nfp, which=0, dbs=pinit)
# Multiply likelihoods for each fixed parameter
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/100", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/10", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/zaps/1", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/3", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/20", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining4/steps/200", numfp=10, num_mult=30)
# Create combine path_lists:
steps_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"]
zaps_list = ["/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"]
p_names = ['gpas', 'k']
# Load parameter space initializator
pinit = tb.open_file("/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/paramsetup.hdf5", mode="r")
# Plot some single results
plot_single_results(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/3", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/20", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/200", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/1", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/10", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/100", numfp=nfp, which=0, dbs=pinit)
# Multiply likelihoods for each fixed parameter
mult_likelihood(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/100", numfp=10, num_mult=30)
mult_likelihood(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/10", numfp=10, num_mult=30)
mult_likelihood(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/sins/1", numfp=10, num_mult=30)
mult_likelihood(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/3", numfp=10, num_mult=30)
mult_likelihood(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/20", numfp=10, num_mult=30)
mult_likelihood(path="/home/szabolcs/parameter_inference/stim_protocol2_v24/comb_colored_dend_gpas-dens/steps/200", numfp=10, num_mult=30)
# Create combine path_lists:
steps_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"]
sins_list = ["/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"]
nfp = 10 # Number of fixed parameters
p_names = ['gpas', 'k']
# Load parameter space initializator
pinit = tb.open_file(".../paramsetup.hdf5", mode="r")
# Plot some single results
plot_single_results(path=".../steps/3", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path=".../steps/20", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path=".../steps/200", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path=".../sins/1", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path=".../sins/10", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path=".../sins/100", numfp=nfp, which=0, dbs=pinit)
# Multiply likelihoods for each fixed parameter
mult_likelihood(path=".../sins/100", numfp=10, num_mult=30)
mult_likelihood(path=".../sins/10", numfp=10, num_mult=30)
mult_likelihood(path=".../sins/1", numfp=10, num_mult=30)
mult_likelihood(path=".../steps/3", numfp=10, num_mult=30)
mult_likelihood(path=".../steps/20", numfp=10, num_mult=30)
mult_likelihood(path=".../steps/200", numfp=10, num_mult=30)
# Create combine path_lists:
steps_list = [".../steps/3", ".../steps/20", ".../steps/200"]
sins_list = [".../sins/1", ".../sins/10", ".../sins/100"]
# Create combinations
combine_likelihood(sins_list,
numfp=10,
num_mult_single=10,
path=
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/10",
numfp=nfp,
which=0,
dbs=pinit)
plot_single_results(
path=
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/100",
numfp=nfp,
which=0,
dbs=pinit)
# Multiply likelihoods for each fixed parameter
mult_likelihood(
path=
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/100",
numfp=10,
num_mult=30)
mult_likelihood(
path=
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/10",
numfp=10,
num_mult=30)
mult_likelihood(
path=
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/sins/1",
numfp=10,
num_mult=30)
mult_likelihood(
path=
"/home/szabolcs/parameter_inference/stim_protocol2_v21/combining_colored/steps/3",
p_names = ['Ra', 'gpas']
# Load parameter space initializator
pinit = tb.open_file("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/paramsetup.hdf5", mode="r")
# Plot some single results
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/3", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/20", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/200", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/1", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/10", numfp=nfp, which=0, dbs=pinit)
plot_single_results(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/100", numfp=nfp, which=0, dbs=pinit)
# Multiply likelihoods for each fixed parameter
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/100", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/10", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/1", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/3", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/20", numfp=10, num_mult=30)
mult_likelihood(path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/200", numfp=10, num_mult=30)
# Create combine path_lists:
steps_list = ["/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/3",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/20",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/200"]
sins_list = ["/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/1",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/10",
"/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/100"]