Alexander/P-unit_model
2
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

update weights, and to new fits final_2

Browse Source
This commit is contained in:
a.ott 2020-09-08 18:32:07 +02:00
parent 7914d00978
commit 72b61cc5f7
3 changed files with 6 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
analysis.py
View File

@ -17,7 +17,7 @@ def main():
# parser.add_argument("dir", help="folder containing the cell folders with the fit results")
# args = parser.parse_args()
dir_path = "results/final_1/" # args.dir
dir_path = "results/final_2/" # args.dir
# if not os.path.isdir(dir_path):
# print("Argument dir is not a directory.")
@ -42,7 +42,7 @@ def main():
# labels, corr_values, corrected_p_values = parameter_correlations(fits_info)
# create_correlation_plot(labels, corr_values, corrected_p_values)
create_parameter_distributions(get_parameter_values(fits_info))
# create_parameter_distributions(get_parameter_values(fits_info))
cell_b, model_b = get_behaviour_values(fits_info)
create_behaviour_distributions(cell_b, model_b)
pass

1
glm_prediction.py
View File

@ -16,6 +16,7 @@ def main():
'f_inf_slope', 'f_zero_slope', 'serial_correlation', 'vector_strength']
behaviour, error = get_variables(folder, variable_order)
df_behaviour = pandas.DataFrame(behaviour, columns=variable_order)
# print(df)
gamma_glm = sm.GLM(error, df_behaviour, sm.families.Gamma())

9
run_Fitter.py
View File

@ -14,13 +14,10 @@ from helperFunctions import plot_errors
import multiprocessing as mp
# SAVE_DIRECTORY = "./results/invivo_results/"
SAVE_DIRECTORY = "./results/final_2/"
# SAVE_DIRECTORY_BEST = "./results/invivo_best/"
SAVE_DIRECTORY_BEST = "./results/final_2_best/"
# [bf, vs, sc, cv, isi_hist, bursty, f_inf, f_inf_slope, f_zero, f_zero_slope, f0_curve]
ERROR_WEIGHTS = (2, 2, 1, 1, 0, 1, 1, 1, 0, 1)
# [vs, sc, cv, isi_hist, bursty, f_inf, f_inf_slope, f_zero, f_zero_slope, f0_curve]
ERROR_WEIGHTS = (1, 1, 1, 1, 1, 1, 1, 1, 0, 1)
def main():
@ -34,7 +31,7 @@ def main():
quit()
# test_single_cell("data/invivo/2010-11-08-al-invivo-1/")
test_single_cell("data/invivo_bursty/2014-03-19-ae-invivo-1/")
test_single_cell("data/invivo_bursty/2013-04-09-ac-invivo-1/")
# start_parameters = [p for p in iget_start_parameters()]
# cell_data = CellData("data/invivo_bursty/2014-03-19-ae-invivo-1/")