CellML Discussion List

Text archives Help


[cellml-discussion] Time delays in CellML models


Chronological Thread 
  • From: matt at elyt.com (Matt)
  • Subject: [cellml-discussion] Time delays in CellML models
  • Date: Wed, 13 Feb 2008 13:17:37 +1300

My guess is to keep it simple:

1) don't prescribe how simulation software is to calculate this -
i.e. don't prescribe a particular mathematical formulation for
representing time delays
2) calculate the variable(s) representing the output of delayed
function separately - e.g. in a separate component or math element -
i.e. don't promote variables to being able to be functions also

I think it is necessarily going to require:

1) math elements that point to user defined code
e.g. in the simplest case the code will buffer of a particular
variable and perhaps create a local approximation curve around the
time required and return a result.

2) referencing math elements from other components
e.g. to calculate f(y-T) you need to at least have a reference to the
math of y in the case that someone has an analytical method to apply
instead of some algorithm over a history of values



On 13/02/2008, at 12:08 PM, Catherine Lloyd wrote:

> Dear All
>
> The issue of how a CellML model should handle time delays has been
> discussed previously by the Auckland CellML group, and at the time
> there
> was a tracker item filed on this subject:
>
> https://tracker.physiomeproject.org/show_bug.cgi?id=70
>
> The topic has been recently raised again (see the email discussion
> below), and it was discussed at today's CellML meeting (13th Feb,
> 2008).
>
> We'd like to open up this discussion to the CellML community to see if
> anyone has any thoughts, ideas, or preferences as to how time delays
> should be a) expressed in a CellML model and b) handled by the
> software
> tools.
>
> Thank you in advance for your feedback.
>
> Best wishes
> Catherine
>
>> Hi Rui,
>>
>> If you are already familiar with Virtual Cell, I'd stick to that. The
>> initial version of SBML was largely based on the internal VCML
>> used by
>> the Virtual Cell, and I suspect the VCell developers remain involved
>> in recent developments within the SBML community.
>>
>> In the past, VCell was capable of importing CellML models but I seem
>> to recall that there might be issues with that in recent versions. If
>> it is still possible, you should be able to import the CellML version
>> of the model and then touch it up from there (add in the time delay).
>> You should then be able to export the SBML version for use in other
>> tools if required.
>>
>> CellML2SBML is a collection of XSL stylesheets which are capable of
>> processing most models from the repository. As Catherine mentioned,
>> you'll lose the metadata from the CellML model in the resulting SBML
>> model, but the math should largely be intact. From my understanding,
>> the CellML2SBML stylesheets are mainly targeted at CellML models
>> which
>> use the now deprecated reaction element, so if the model you are
>> interested in uses reactions then you'll get the best SBML version.
>> Otherwise I suspect you'll end up with a generic SBML encoding of the
>> mathematics.
>>
>> Oh, and I'm pretty sure VCell will allow you to specify time delays
>> although I can't remember the specific syntax or interface options
>> and
>> you can see https://tracker.physiomeproject.org/show_bug.cgi?id=70
>> for
>> a bit of discussion on this issue within CellML - feel free to add to
>> that discussion!
>>
>> Not sure if I can offer much more help than that but let me know if
>> you have any more questions...
>>
>>
>> David.
>>
>>
>> Catherine Lloyd wrote:
>> Dear Rui
>>
>> You are welcome!
>>
>> Hmmm, much as I would hate to lose a potential CellML user to SBML, I
>> can see the advantage for you if indeed, SBML can handle time delays.
>> Although I'm not familiar with the structure of SBML, my
>> understanding
>> is that the cellml2sbml converter is not perfect - if nothing
>> else, it
>> strips all the metadata out of a CellML model. But you might be
>> able to
>> use it as a foundation which you can then edit by hand.
>>
>> I have cc'd David Nickerson in on this email - if anyone deserves the
>> title of "CellML expert" it is David!... David do you have any idea
>> which SBML tool might be best to use in this case?
>>
>> I've also cc'd Poul Nielsen and Peter Hunter who lead the CellML
>> project. I'm sure they will be interested in this case and they may
>> have an idea how long it might take to develop the CellML tools such
>> that they are capable of handling time delays.
>>
>> I'm sorry I can't help you more myself, but I imagine we can work
>> together to get this up and running for you.
>>
>> Best wishes
>> Catherine
>>
>>
>>
>> Rui Zhu wrote:
>>> Dear Catherine,
>>>
>>> Thank you so much for your reply! It is very helpful! Learning
>>> CellML
>>> modeling from an expert is a really delightful experience.
>>>
>>> I am having another thought. I searched online, and it says SBML
>>> level2 can deal with the time delay term, and there is a software
>>> (cellml2sbml) to convert CellML to SBML. Do you think it is a good
>>> idea to use SBML to solve the problem? And from SBML website, there
>>> are a lot of tools that can run SBML model, like virtual cell. I am
>>> wondering if you have any idea which tool can handle the time delay
>>> term. Thanks!
>>>
>>> Your suggestion will be appreaciated!
>>>
>>> Many thanks,
>>>
>>> Rui ZHU
>>>
>>>
>>> ----- Original Message -----
>>> From: Catherine Lloyd <c.lloyd at auckland.ac.nz>
>>> Date: Sunday, February 10, 2008 11:19 am
>>> Subject: Re: time delay term in CellML model
>>> To: Rui Zhu <rzhu at usc.edu>
>>> Cc: Poul Nielsen <p.nielsen at auckland.ac.nz>, Peter Hunter
>>> <p.hunter at auckland.ac.nz>
>>>
>>>
>>>> Dear Rui
>>>>
>>>> Thank you for your interest in CellML. I am familiar with the
>>>> model
>>>> you
>>>> are interested in, and I know the original model author (Michael
>>>> Zagar)helped me to translate this model into CellML. I'm sorry
>>>> if our
>>>> curation status confused you - it does indeed state that the CellML
>>>> model recreates the results of the published model - but it should
>>>> statethat the models almost match - with the exception of the time
>>>> delay(I'll change the text now).
>>>> Unfortunately this is a current limitation of CellML. While it is
>>>> possible to express a time delay in the CellML language, the tools
>>>> whichrun CellML models (PCEnv and COR) are unable to handle time
>>>> delays at
>>>> present. This is something which we expect to be a feature in
>>>> future
>>>> releases of the software.
>>>>
>>>> Please let me know if you have any further questions and I will
>>>> try to
>>>> answer them for you.
>>>>
>>>> Best wishes
>>>> Catherine
>>>>
>>>>
>>>>
>>>> Rui Zhu wrote:
>>>>
>>>>> Dear Catherine,
>>>>>
>>>>> This is Rui ZHU, a Biomedical Engineering Ph.D student in
>>>> University of Southern California. I am very interested in using
>>>> CellML to do modeling. Recently, I saw the CellML model you made
>>>> based on the paper "A Delayed Nonlinear PBPK Model for Genistein
>>>> Dosimetry in Rats" on the CellML website.
>>>> http://www.cellml.org/models/zager_schlosser_tran_2007_version03
>>>>> And I tried to use that model to do simulation using PCEnv to re-
>>>>>
>>>> generate the results in the original paper. However, I found the
>>>> result was different from the one the in original paper, and the
>>>> reason it that, in the original paper, there is a time delay
>>>> term in
>>>> the differential equation d( Acon_b )/d(time) = kbile * H_Ccon_l
>>>> * Ccon_l ( t - d(t, C_gen_l) ) * Vl, but in the CellML model the
>>>> differential equation is without the time delay term d( Acon_b
>>>> )/d(time) = kbile * H_Ccon_l * Ccon_l * Vl. I was trying to
>>>> add a time delay term in the CellML model, but I don't know how. If
>>>> you have any idea of how to do it, would you please let me know.
>>>> Thank you!
>>>>> Your kind help is really appreciated! Look forward to your reply!
>>>>>
>>>>> Rui ZHU
>>>>>
>>>>>
>>>>
>
> _______________________________________________
> cellml-discussion mailing list
> cellml-discussion at cellml.org
> http://www.cellml.org/mailman/listinfo/cellml-discussion





Archive powered by MHonArc 2.6.18.

Top of page