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

Could we just set a slightly more restrictive prior on F_PL? #108

Open
sonyahanson opened this issue Jun 8, 2017 · 1 comment
Open

Could we just set a slightly more restrictive prior on F_PL? #108

sonyahanson opened this issue Jun 8, 2017 · 1 comment

Comments

@sonyahanson
Copy link
Contributor

Tested the effect of setting F_PL using values taken from the mean after equilibration from a previous run of assaytools, and dF_PL=0.50, which was approximated from the standard deviation. Tested this for Src:Erlotinib and Src: Gefitinib.

F_PL_def = {'Src-Erlotinib-EF': 59601584903.0,
            'Src-Gefitinib-GH': 81835891154.8}

Similar to issue #106 we see an expected removal of long tails and better reproducibility compared to without setting the F_PL. Note here, we are still looking at a dPstated = 0.35 * inputs['Pstated'].

Without setting F_PL:

src-erlotinib-ef_binding_iter0

Setting F_PL and dF_PL:

src-erlotinib-ef_binding_iter0

Without setting F_PL:

comparing_src_geferl_3iter

Setting F_PL and dF_PL:

comparing_src_geferl_3iter

Also we no longer see any sort of correlation between F_PL and DeltaG:

Without setting F_PL:

delg_v_fpl_erl0

Setting F_PL and dF_PL:
delg_v_fpl_erl0

While this is nice, I'm worried this is maybe too restrictive? Maybe we can just have a slightly more informed prior instead of this (link):

F_PL_guess = (Fmax - Fmin) / min(Pstated.max(), Lstated.max())
model['F_PL'] = pymc.Uniform('F_PL', lower=0.0, upper=2*Fmax/min(Pstated.max(),Lstated.max()), value=F_PL_guess)
@jchodera
Copy link
Member

jchodera commented Jun 8, 2017

Wow!

Can you show this to Nick? We might ask his thoughts on using an estimate of F_PL and its uncertainty from our spectra assay and/or a separate cuvette assay.

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

2 participants