def simulate_sources(self, gen_conf_1, gen_conf_2): [simulated_gen_1, gen_conf_1, electrodes] = \ gen_simulation(gen_conf_1, num_sim=self.number_of_subjects, gen_magnitude_stddev=\ self.generator_magnitude_stddev) [simulated_gen_2, gen_conf_2, electrodes] = \ gen_simulation(gen_conf_2, num_sim=self.number_of_subjects, gen_magnitude_stddev=\ self.generator_magnitude_stddev) return [simulated_gen_1, gen_conf_1, simulated_gen_2, gen_conf_2, electrodes]
def simulate_noise(self, size): if self.normal_noise_amplitude: normal_noise_1 = normal(0, self.normal_noise_amplitude, size) normal_noise_2 = normal(0, self.normal_noise_amplitude, size) else: normal_noise_1 = None normal_noise_2 = None if self.constant_noise_amplitude: constant_noise = self.constant_noise_amplitude * ones(size) else: constant_noise = None if self.topographic_noise_amplitude: noise_conf_1 = \ random_generator_configuration( self.limits_number, self.limits_depth, self.limits_orientation, self.limits_magnitude, self.topographic_noise_amplitude)[0] noise_conf_2 = \ random_generator_configuration( self.limits_number, self.limits_depth, self.limits_orientation, self.limits_magnitude, self.topographic_noise_amplitude)[0] [topographic_noise_1, noise_conf_1, electrodes] = \ gen_simulation(noise_conf_1, num_sim=self.number_of_subjects) [topographic_noise_2, noise_conf_2, electrodes] = \ gen_simulation(noise_conf_2, num_sim=self.number_of_subjects) else: topographic_noise_1 = None topographic_noise_2 = None return [normal_noise_1, topographic_noise_1, normal_noise_2, topographic_noise_2, constant_noise]
plot_topographic_map_array(map_array, scale=scale_array, label_y=label_y_array, title=title_array) savefig(folder + "/subject_" + str(i + 1) + ".png") if parameters["todo_plot_grand_averages"]: print("Drawing grand averages for all imported data") figure() plot_ERP(ERP_avg["old"][0]) if parameters["todo_plot_simulated_noise"]: print("Running simulation to show topographic properties of noise") gen_conf = random_generator_configuration() [simulated_data, gen_conf, electrodes] = gen_simulation( gen_conf, num_sim=parameters["number_of_subjects"], gen_magnitude_stddev=parameters["generator_magnitude_stddev"], ) figure() plot_topographic_map(mean(simulated_data, 0)) # savefig('Results/' + timestamp + '/simulated_noise') savefig("Results/" + label + "/simulated_noise") # Here we will be running the scaling experiments # se <-- holds all simulation experiments, for all parameter combinations se = {} # All simulation experiment parameters we will be looping through experiments = [1, 2, 3, 4] configurations = ["random", "uandk"]