Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Adding Quickmodel script the option to use informed F_PL priors #138

Open
MehtapIsik opened this issue Feb 13, 2019 · 1 comment
Open

Adding Quickmodel script the option to use informed F_PL priors #138

MehtapIsik opened this issue Feb 13, 2019 · 1 comment

Comments

@MehtapIsik
Copy link
Contributor

We would like to be able to provide an informed F_PL prior based on previous experiments for Quickmodel analysis. Currently, it isn't possible to modify F_PL prior in Quickmodel analysis. We should make F_PL an accessible parameter that we can modify through input.py file for Quickmodel analysis.

@MehtapIsik
Copy link
Contributor Author

F_PL from an previous experiments Quickmodel run results can be extracted from the pickle files.
We should create a separate utility script for extracting F_PL from pickle files and place it under assaytools/scripts/ directory. Since pickle files are very large and difficult to store it would be more flexible to have this as a separate script.

Example script:

import cPickle
data_file = 'Src-Bosutinib-AB_mcmc-2017-06-18 08:42.pickle'
with open(r'%s'%data_file,'rb') as my_file:
     data = cPickle.load(my_file)
F_PL_mean = data['F_PL'].mean()

However, that this script above takes the mean from the entire trace, not just the equillibrated trace. We think that equilibration time (t_equil) should be extracted from json file, which was determined based on equilibriation of DeltaG value. We should throw away the same amount of equilibration time from the trace of fluorescence and calculate the mean F_PL.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant