Groups keyboard shortcuts have been updated
Dismiss
See shortcuts

Parameter estimation for different reaction rates

36 views
Skip to first unread message

cel.c...@gmail.com

unread,
Jun 14, 2024, 7:17:33 AM6/14/24
to COPASI User Forum
Hello,
I want to estimate two parameters and I want to do that with timecourse data which changes based on one reaction parameter.
So e.g. we want to estimate the values for parameter k1 and k2.
And I am interested in doing that w.r.t the concentration curve of AC_activity , which is an assignment type (not sure if that is important).
Anyways as a rule of thumb, the concentration curve of AC_activity, in particular the peaks for AC_activity look different with different k_releases.
I generated my dataset like this:

Time AC_activity k_release
0              0                1e-10
..              ...                1e-10
100          0                1e-10
################# (empty line)
0               0.5              1e-2
...              ...                 1e-2
100           0.5              1e-2

But somehow the evaluation runs fail. Can someone help with what I am doing wrong ?
Best Celvic

Frank Bergmann

unread,
Jun 14, 2024, 7:53:39 AM6/14/24
to COPASI User Forum
Hello Celvic, 

I'm afraid we would need a little bit more information, especially how you mapped the data to your model. What are the dependent / independent model elements mapped to? If AC_activity is of type assignment, then they ought to be mapped as dependent element. K_release sounds like an independent value that you specify. 

I can offer that you send me the COPASI file + data and i can look over it to see why the evaluation fails. 

best
Frank

cel.c...@gmail.com

unread,
Jun 14, 2024, 8:07:30 AM6/14/24
to COPASI User Forum
Hello Frank,
yes I chose AC_activity to be dependent and k_release to be independent and  I would appreciate your offer
Best Celvic
Reply all
Reply to author
Forward
0 new messages