Python mult_likelihood Examples, module.protocol_test.mult_likelihood Python Examples

Example #1

0

Show file

File: evaluate.py Project: terbed/parameter-inference

    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)

Example #2

0

Show file

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"]

Example #3

0

Show file

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"]

Example #4

0

Show file

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,

Example #5

0

Show file

    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",

Example #6

0

Show file

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"]