def single_result_plot(self, which_rep=0):
"""
This function save the plots of likelihoods and posteriors for each fixed parameter and one chosen repetition
(because it would be too many to plot all of them...)
:param which_rep: at which repetition number to plot results (should be less than the number of repetitions...)
:return: save the plots of the single likelihood and posterior distributions at given repnum for all fix params
"""
for protocol in self.subdirs:
plot_single_results(path=(self.rootdir + protocol), numfp=self.nfp, which=which_rep, dbs=self.pinit)
from module.protocol_test import plot_single_results, plot_combined_results, mult_likelihood,\
combine_likelihood, protocol_comparison
import time
import tables as tb
startTime = time.time()
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",
from module.protocol_test import plot_single_results, plot_combined_results, mult_likelihood,\
combine_likelihood, protocol_comparison
import time
import tables as tb
startTime = time.time()
nfp = 10 # Number of fixed parameters
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",
from module.protocol_test import plot_single_results, plot_combined_results, mult_likelihood,\
combine_likelihood, protocol_comparison
import time
import tables as tb
startTime = time.time()
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"]
