[Oscill8-devel] Re: Issues with Oscill8

SourceForge Headquarters 1320 Columbia Street Suite 310 San Diego, CA 92101 +1 (858) 422-6466

<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Vijay,<br>
<br>
Update to the latest Oscill8, v1.9.23. It has better default values set
for the Sample Period/Min/Max. Tersely, here are<br>
what the parameters mean:<br>
<br>
&nbsp;&nbsp;&nbsp; double period_min;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; min period to look for<br>
&nbsp;&nbsp;&nbsp; double period_max;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; max period to look for<br>
&nbsp;&nbsp;&nbsp; double h_min;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; min step for ODE solver (and use
this close to finding period for accuracy)<br>
&nbsp;&nbsp;&nbsp; double h_max;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; max step for ODE solver (when we
think we can make big jumps)<br>
&nbsp;&nbsp;&nbsp; double h_init;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; initial step for ODE solver<br>
&nbsp;&nbsp;&nbsp; double x_match_tol;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; (xi-xi_at_last_period) &lt;
match_tol ==&gt; match<br>
&nbsp;&nbsp;&nbsp; double period_match_tol; //!&lt; (peri-peri_at_last_match) &lt;
period_tol ==&gt; period is detected match<br>
&nbsp;&nbsp;&nbsp; int period_match_max;&nbsp;&nbsp;&nbsp; //!&lt; (peri-peri_at_last_match) &lt;
period_tol ==&gt; period is detected match<br>
&nbsp;&nbsp;&nbsp; int p_index;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; index of parameter to continue<br>
&nbsp;&nbsp;&nbsp; double dp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; sample step for parameter in
question<br>
&nbsp;&nbsp;&nbsp; int n_max;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; maximum number of valid (positive)
period calculations before quit<br>
&nbsp;&nbsp;&nbsp; int edge_max;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; number of bad (negative) results
before quit<br>
&nbsp;&nbsp;&nbsp; int x_index;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; //!&lt; variable on which to calculate
period<br>
<br>
In answer to your questions below:<br>
<br>
1) when you say "steady states", this refers to a 1 parameter
continuation. You would set the parameter range to .2 to 1.2 and run a
continuation (parameter ranges are set through the "Set Model Data"
button in the Run Config dialog).<br>
<br>
BUT, do you mean the max and min of the oscillation? If so, then you
want to start with a parameter value in the .2 to 1.2 range, and choose
n_max large and dp appropriately (to get a smooth number of samples).
The sampling of the period/min/max starts from your initial value and
goes as FAR as it can in each direction before it can't calculate the
period "edge_max" times in a row. This should take you to the
boundaries of your oscillatory domain.<br>
<br>
2) When the graph shows nothing, it wasn't able to calculate the period
at ALL.. this mean either 1) you didn't set the parameter for the
sampling well or 2) you forgot to put the initial parameter value into
the oscillatory domain.<br>
<br>
3) More information on a HB is in the works; but, this is not something
that AUTO provides, so I haven't yet added this feature. It's on my
list.<br>
<br>
Don't every worry about asking questions!!<br>
<br>
Emery<br>
<br>
<br>
<br>
Vijay Chickarmane wrote:
<blockquote
 cite="mid...@de..."
 type="cite">
  <meta http-equiv="Content-Type" content="text/html; ">
  <meta content="MSHTML 6.00.2800.1515" name="GENERATOR">
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">Hi Emory,</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">Sorry for the delay, I did set up the x_index as
2 and I was</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">able to obtain the limits. But I have no feel
for what the </font></span></div>
  <div><span class="182293416-31082005"></span><span
 class="182293416-31082005"><font color="#0000ff" face="Arial" size="2">parameters
mean. Could you tell me how I can do the foll.</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">1. If&nbsp; I want to&nbsp;get the&nbsp;steady states for a
system&nbsp;with the</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">&nbsp;parameter ranging from .2 to 1.2, how could I
set this up.</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">2. Sometimes the program cannot compute the
steady states,</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">the graph does not show anything, usually I
reset the parameter</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">and try again. But since there is no error
message, one does not</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">know what is hapenning.</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">3. Can we get more information on a HB, ex: sub
or supercritical?</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">Hope this deluge of questions will not take up
too much of your&nbsp;time.</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">thanks for the help,</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">best,</font></span></div>
  <div><span class="182293416-31082005"><font color="#0000ff"
 face="Arial" size="2">vijay</font></span></div>
  <div><span class="182293416-31082005"></span>&nbsp;</div>
  <div><span class="182293416-31082005"></span>&nbsp;</div>
  <blockquote style="margin-right: 0px;" dir="ltr">
    <div class="OutlookMessageHeader" align="left" dir="ltr"><font
 face="Tahoma" size="2">-----Original Message-----<br>
    <b>From:</b> Emery Conrad [<a class="moz-txt-link-freetext" href="mailto:ec...@vt...">mailto:ec...@vt...</a>]<br>
    <b>Sent:</b> Friday, August 26, 2005 11:19 AM<br>
    <b>To:</b> Vijay Chickarmane<br>
    <b>Subject:</b> Re: Issues with Oscill8<br>
    <br>
    </font></div>
Vijay,<br>
    <br>
Which parameters where you interested in the most? the sampling works
fine for me for some of the parameters, but you need to choose your
algorithmic pars correctly (for example, x_index should be set to 1 or
2 given that x0 is a very "stiff" variable).<br>
    <br>
Let me know which parameters you're interested in and I'll tell you how
to set the algorithm to run correctly.<br>
    <br>
Emery<br>
    <br>
    <br>
    <br>
Vijay Chickarmane wrote:
    <blockquote
 cite="mid...@de..."
 type="cite">
      <pre wrap="">Hi Emory,
Greetings, I am sending you a model
which has oscillatory behavior, but I was
unable to obtain the max/min of the oscillations.
It could very well be that I did not use the
tool properly. Could you let me know how
you are able to use this feature. 
best wishes,
vijay

-----Original Message-----
From: Emery Conrad [<a class="moz-txt-link-freetext"
 href="mailto:ec...@vt...">mailto:ec...@vt...</a>]
Sent: Wednesday, August 24, 2005 12:34 PM
To: Herbert Sauro
Cc: Vijay Chickarmane
Subject: Re: Issues with Oscill8

Vijay and Herb,

Sure thing... the period/min/max thingy is a bit exploratory right now 
and very sensitive to the algorithms parameters. I've updated things in 
the last few days (I need this feature for temperature compensation code 
which I'm working on), so it should be in better shape.

Vijay, can you send me the models you were dealing with so I can check 
things out?

Thanks,

Emery

Herbert Sauro wrote:

  </pre>
      <blockquote type="cite">
        <pre wrap="">Vijay tried to use Oscill8 to determine the max and min for an
oscillatory regime but was unable to for some reason, instead he did it
himself using Matlab. Could you email him to ask him for the models in
question and figure out what was wrong?

Thanks
Herbert

    </pre>
        <pre wrap=""><hr size="4" width="90%">
&lt;?xml version="1.0" encoding="UTF-8"?&gt;
&lt;!-- Created by XMLPrettyPrinter on 8/5/2005 from JDesigner 2.0.11 Beta --&gt;
&lt;sbml level = "2" version = "1" xmlns = <a
 class="moz-txt-link-rfc2396E" href="http://www.sbml.org/sbml/level2">"http://www.sbml.org/sbml/level2"</a> xmlns:jd = <a
 class="moz-txt-link-rfc2396E" href="http://www.sys-bio.org/sbml">"http://www.sys-bio.org/sbml"</a>&gt;
   &lt;!--                     --&gt;
   &lt;!--  Model Starts Here  --&gt;
   &lt;!--                     --&gt;
   &lt;model id = "untitled"&gt;
      &lt;annotation&gt;
         &lt;jd:header&gt;
            &lt;VersionHeader JDesignerVersion = "2.0"/&gt;
            &lt;ModelHeader Author = "Mr Unknown" ModelTitle = "No Title" ModelVersion = "0.0"/&gt;
         &lt;/jd:header&gt;
         &lt;jd:display&gt;
            &lt;SBMLGraphicsHeader BackGroundColor = "FFEFFFFF"/&gt;
         &lt;/jd:display&gt;
      &lt;/annotation&gt;
      &lt;listOfCompartments&gt;
         &lt;compartment id = "compartment" size = "1"&gt;
            &lt;annotation&gt;
               &lt;jd:display&gt;
                  &lt;boundingBox h = "0" w = "0" x = "0" y = "0"/&gt;
               &lt;/jd:display&gt;
            &lt;/annotation&gt;
         &lt;/compartment&gt;
      &lt;/listOfCompartments&gt;
      &lt;listOfSpecies&gt;
         &lt;species boundaryCondition = "false" compartment = "compartment" id = "x1" initialConcentration = "0"&gt;
            &lt;annotation&gt;
               &lt;jd:display&gt;
                  &lt;boundingBox h = "32" w = "48" x = "427" y = "62"/&gt;
                  &lt;pictogram boundarySpeciesStyle = "bsVertVice" complexName0 = "" complexName1 = "" complexName2 = ""
                             complexName3 = "" edgeColor = "FF969696" edgeThickness = "1" endColor = "FFFFFFFF"
                             gradientType = "gtLinear" namePosition = "npCenter" pictoComplexType = "cstSimple" pictoShape = "ptRoundSquare"
                             startColor = "FFFFFFCC"/&gt;
                  &lt;font fontColor = "0" fontName = "Arial" fontSize = "8" fontStyle = ""/&gt;
               &lt;/jd:display&gt;
            &lt;/annotation&gt;
         &lt;/species&gt;
         &lt;species boundaryCondition = "false" compartment = "compartment" id = "y1" initialConcentration = "1"&gt;
            &lt;annotation&gt;
               &lt;jd:display&gt;
                  &lt;boundingBox h = "32" w = "48" x = "442" y = "118"/&gt;
                  &lt;pictogram boundarySpeciesStyle = "bsVertVice" complexName0 = "" complexName1 = "" complexName2 = ""
                             complexName3 = "" edgeColor = "FF969696" edgeThickness = "1" endColor = "FFFFFFFF"
                             gradientType = "gtLinear" namePosition = "npCenter" pictoComplexType = "cstSimple" pictoShape = "ptRoundSquare"
                             startColor = "FFFFFFCC"/&gt;
                  &lt;font fontColor = "0" fontName = "Arial" fontSize = "8" fontStyle = ""/&gt;
               &lt;/jd:display&gt;
            &lt;/annotation&gt;
         &lt;/species&gt;
         &lt;species boundaryCondition = "false" compartment = "compartment" id = "b" initialConcentration = "0"&gt;
            &lt;annotation&gt;
               &lt;jd:display&gt;
                  &lt;boundingBox h = "32" w = "48" x = "513" y = "90"/&gt;
                  &lt;pictogram boundarySpeciesStyle = "bsVertVice" complexName0 = "" complexName1 = "" complexName2 = ""
                             complexName3 = "" edgeColor = "FF969696" edgeThickness = "1" endColor = "FFFFFFFF"
                             gradientType = "gtLinear" namePosition = "npCenter" pictoComplexType = "cstSimple" pictoShape = "ptRoundSquare"
                             startColor = "FFFFFFCC"/&gt;
                  &lt;font fontColor = "0" fontName = "Arial" fontSize = "8" fontStyle = ""/&gt;
               &lt;/jd:display&gt;
            &lt;/annotation&gt;
         &lt;/species&gt;
         &lt;species boundaryCondition = "true" compartment = "compartment" id = "Node3" initialConcentration = "0.1"&gt;
            &lt;annotation&gt;
               &lt;jd:display&gt;
                  &lt;boundingBox h = "32" w = "48" x = "291" y = "76"/&gt;
                  &lt;pictogram boundarySpeciesStyle = "bsVertVice" complexName0 = "" complexName1 = "" complexName2 = ""
                             complexName3 = "" edgeColor = "FF969696" edgeThickness = "1" endColor = "FFFFFFFF"
                             gradientType = "gtLinear" namePosition = "npCenter" pictoComplexType = "cstSimple" pictoShape = "ptRoundSquare"
                             startColor = "FFFFFFCC"/&gt;
                  &lt;font fontColor = "0" fontName = "Arial" fontSize = "8" fontStyle = ""/&gt;
               &lt;/jd:display&gt;
            &lt;/annotation&gt;
         &lt;/species&gt;
      &lt;/listOfSpecies&gt;
      &lt;listOfReactions&gt;
         &lt;reaction id = "J2" reversible = "false"&gt;
            &lt;annotation&gt;
               &lt;jd:display arrowFillColor = "FF648483" lineColor = "FF0099FF" lineThickness = "2"&gt;
                  &lt;lineType type = "ltBezier"&gt;
                     &lt;pt x = "0" y = "0"/&gt;
                     &lt;pt x = "382" y = "81"/&gt;
                     &lt;pt x = "382" y = "81"/&gt;
                  &lt;/lineType&gt;
               &lt;/jd:display&gt;
            &lt;/annotation&gt;
            &lt;listOfReactants&gt;
               &lt;speciesReference species = "Node3" stoichiometry = "1"/&gt;
            &lt;/listOfReactants&gt;
            &lt;listOfProducts&gt;
               &lt;speciesReference species = "x1" stoichiometry = "1"/&gt;
            &lt;/listOfProducts&gt;
            &lt;kineticLaw&gt;
               &lt;math xmlns = <a class="moz-txt-link-rfc2396E"
 href="http://www.w3.org/1998/Math/MathML">"http://www.w3.org/1998/Math/MathML"</a>&gt;
                  &lt;apply&gt;
                     &lt;minus/&gt;
                     &lt;apply&gt;
                        &lt;divide/&gt;
                        &lt;apply&gt;
                           &lt;times/&gt;
                           &lt;ci&gt;
                                  a1 
                           &lt;/ci&gt;
                           &lt;ci&gt;
                                  b 
                           &lt;/ci&gt;
                        &lt;/apply&gt;
                        &lt;apply&gt;
                           &lt;plus/&gt;
                           &lt;ci&gt;
                                  k1 
                           &lt;/ci&gt;
                           &lt;apply&gt;
                              &lt;power/&gt;
                              &lt;ci&gt;
                                     y1 
                              &lt;/ci&gt;
                              &lt;ci&gt;
                                     n 
                              &lt;/ci&gt;
                           &lt;/apply&gt;
                        &lt;/apply&gt;
                     &lt;/apply&gt;
                     &lt;apply&gt;
                        &lt;times/&gt;
                        &lt;ci&gt;
                               g1 
                        &lt;/ci&gt;
                        &lt;ci&gt;
                               x1 
                        &lt;/ci&gt;
                     &lt;/apply&gt;
                  &lt;/apply&gt;
               &lt;/math&gt;
               &lt;listOfParameters&gt;
                  &lt;parameter id = "a1" value = "0.2"/&gt;
                  &lt;parameter id = "k1" value = "1"/&gt;
                  &lt;parameter id = "n" value = "4"/&gt;
                  &lt;parameter id = "g1" value = "0.1"/&gt;
               &lt;/listOfParameters&gt;
            &lt;/kineticLaw&gt;
         &lt;/reaction&gt;
         &lt;reaction id = "J1" reversible = "false"&gt;
            &lt;annotation&gt;
               &lt;jd:display arrowFillColor = "FF648483" lineColor = "FF0099FF" lineThickness = "2"&gt;
                  &lt;lineType type = "ltBezier"&gt;
                     &lt;pt x = "0" y = "0"/&gt;
                     &lt;pt x = "390" y = "111"/&gt;
                     &lt;pt x = "390" y = "111"/&gt;
                  &lt;/lineType&gt;
               &lt;/jd:display&gt;
            &lt;/annotation&gt;
            &lt;listOfReactants&gt;
               &lt;speciesReference species = "Node3" stoichiometry = "1"/&gt;
            &lt;/listOfReactants&gt;
            &lt;listOfProducts&gt;
               &lt;speciesReference species = "y1" stoichiometry = "1"/&gt;
            &lt;/listOfProducts&gt;
            &lt;kineticLaw&gt;
               &lt;math xmlns = <a class="moz-txt-link-rfc2396E"
 href="http://www.w3.org/1998/Math/MathML">"http://www.w3.org/1998/Math/MathML"</a>&gt;
                  &lt;apply&gt;
                     &lt;minus/&gt;
                     &lt;apply&gt;
                        &lt;divide/&gt;
                        &lt;ci&gt;
                               a2 
                        &lt;/ci&gt;
                        &lt;apply&gt;
                           &lt;plus/&gt;
                           &lt;ci&gt;
                                  k2 
                           &lt;/ci&gt;
                           &lt;apply&gt;
                              &lt;power/&gt;
                              &lt;ci&gt;
                                     x1 
                              &lt;/ci&gt;
                              &lt;ci&gt;
                                     n 
                              &lt;/ci&gt;
                           &lt;/apply&gt;
                        &lt;/apply&gt;
                     &lt;/apply&gt;
                     &lt;apply&gt;
                        &lt;times/&gt;
                        &lt;ci&gt;
                               g2 
                        &lt;/ci&gt;
                        &lt;ci&gt;
                               y1 
                        &lt;/ci&gt;
                     &lt;/apply&gt;
                  &lt;/apply&gt;
               &lt;/math&gt;
               &lt;listOfParameters&gt;
                  &lt;parameter id = "a2" value = "0.15"/&gt;
                  &lt;parameter id = "k2" value = "1"/&gt;
                  &lt;parameter id = "n" value = "4"/&gt;
                  &lt;parameter id = "g2" value = "0.1"/&gt;
               &lt;/listOfParameters&gt;
            &lt;/kineticLaw&gt;
         &lt;/reaction&gt;
         &lt;reaction id = "J15" reversible = "false"&gt;
            &lt;annotation&gt;
               &lt;jd:display arrowFillColor = "FF648483" lineColor = "FF0099FF" lineThickness = "2"&gt;
                  &lt;lineType type = "ltBezier"&gt;
                     &lt;pt x = "0" y = "0"/&gt;
                     &lt;pt x = "426" y = "97"/&gt;
                     &lt;pt x = "426" y = "97"/&gt;
                  &lt;/lineType&gt;
               &lt;/jd:display&gt;
            &lt;/annotation&gt;
            &lt;listOfReactants&gt;
               &lt;speciesReference species = "Node3" stoichiometry = "1"/&gt;
            &lt;/listOfReactants&gt;
            &lt;listOfProducts&gt;
               &lt;speciesReference species = "b" stoichiometry = "1"/&gt;
            &lt;/listOfProducts&gt;
            &lt;kineticLaw&gt;
               &lt;math xmlns = <a class="moz-txt-link-rfc2396E"
 href="http://www.w3.org/1998/Math/MathML">"http://www.w3.org/1998/Math/MathML"</a>&gt;
                  &lt;apply&gt;
                     &lt;minus/&gt;
                     &lt;apply&gt;
                        &lt;times/&gt;
                        &lt;ci&gt;
                               c1 
                        &lt;/ci&gt;
                        &lt;ci&gt;
                               y1 
                        &lt;/ci&gt;
                     &lt;/apply&gt;
                     &lt;apply&gt;
                        &lt;times/&gt;
                        &lt;ci&gt;
                               g 
                        &lt;/ci&gt;
                        &lt;ci&gt;
                               b 
                        &lt;/ci&gt;
                     &lt;/apply&gt;
                  &lt;/apply&gt;
               &lt;/math&gt;
               &lt;listOfParameters&gt;
                  &lt;parameter id = "c1" value = "0.01"/&gt;
                  &lt;parameter id = "g" value = "0.01"/&gt;
               &lt;/listOfParameters&gt;
            &lt;/kineticLaw&gt;
         &lt;/reaction&gt;
      &lt;/listOfReactions&gt;
   &lt;/model&gt;
&lt;/sbml&gt;
    </pre>
      </blockquote>
    </blockquote>
  </blockquote>
</blockquote>
</body>
</html>