示例#1
0
    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)
示例#2
0
    def combinations(self, comb_lists, output_dirs):
        """
        Multiplies likelihoods from different protocols and plots results.
        :param comb_lists: list of list of protocol names to combine \ struct: [ ["step3", "step100"], ["sin10", "sin100"] ]
        :param output_dirs: directory for the output of combinations (inside root dir) \ for examp: ["steps/comb", "sins/comb"]
        :return: Saves and plots the combinations of protocols to disk
        """

        # create lists for combine_likelihoods function
        comb_lists = [[(self.rootdir+elem) for elem in comb_lists[i]] for i in range(len(comb_lists))]
        out_dirs = [(self.rootdir+subdir) for subdir in output_dirs]

        print("Combination lists:")
        print(comb_lists)

        for idx, comblist in enumerate(comb_lists):
            combine_likelihood(comblist, numfp=self.nfp, num_mult_single=self.nrep/len(comblist), out_path=out_dirs[idx])
            plot_combined_results(out_dirs[idx], self.nfp, self.pinit)
示例#3
0
# 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"]

# 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",
示例#4
0
# 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"]

# Create combinations
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",
示例#5
0
# 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,
                   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,
示例#6
0
# 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"]

# Create combinations
combine_likelihood(sins_list, numfp=10, num_mult_single=10,
                   out_path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/comb")
combine_likelihood(steps_list, numfp=10, num_mult_single=10,
                   out_path="/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/comb")

plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/100", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/10", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/1", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/sins/comb", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/3", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/20", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/200", 10, dbs=pinit)
plot_combined_results("/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/comb", 10, dbs=pinit)

path_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",
             "/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",
             "/Users/Dani/TDK/parameter_estim/stim_protocol2/combining_colored_dend-soma/steps/comb",