# A stochastic differential equation analysis of cerebrospinal fluid dynamics

- Kalyan Raman
^{1}Email author

**8**:9

**DOI: **10.1186/2045-8118-8-9

© Raman; licensee BioMed Central Ltd. 2011

**Received: **1 September 2010

**Accepted: **18 January 2011

**Published: **18 January 2011

## Abstract

### Background

Clinical measurements of intracranial pressure (ICP) over time show fluctuations around the deterministic time path predicted by a classic mathematical model in hydrocephalus research. Thus an important issue in mathematical research on hydrocephalus remains unaddressed--modeling the effect of noise on CSF dynamics. Our objective is to mathematically model the noise in the data.

### Methods

The classic model relating the temporal evolution of ICP in pressure-volume studies to infusions is a nonlinear differential equation based on natural physical analogies between CSF dynamics and an electrical circuit. Brownian motion was incorporated into the differential equation describing CSF dynamics to obtain a nonlinear stochastic differential equation (SDE) that accommodates the fluctuations in ICP.

### Results

The SDE is explicitly solved and the dynamic probabilities of exceeding critical levels of ICP under different clinical conditions are computed. A key finding is that the probabilities display strong threshold effects with respect to noise. Above the noise threshold, the probabilities are significantly influenced by the resistance to CSF outflow and the intensity of the noise.

### Conclusions

Fluctuations in the CSF formation rate increase fluctuations in the ICP and they should be minimized to lower the patient's risk. The nonlinear SDE provides a scientific methodology for dynamic risk management of patients. The dynamic output of the SDE matches the noisy ICP data generated by the actual intracranial dynamics of patients better than the classic model used in prior research.

## Background

Intracranial dynamics play a central role in healthy brain function because disturbances in the internal fluid environment of the skull can lead to multiple complications such as, among other things, hydrocephalus [1]. Intracranial dynamics, driven by the circulation of CSF, are important because CSF protects the brain from injury, contains nutrients enabling normal functioning of the brain and, transports waste products away from the surrounding tissues. Much more is involved in hydrocephalus than a simple disorder of CSF circulation [2]; it is considered a complex spectrum of diseases, primarily defined by perturbation of the cranial contents--operationalized as CSF volume--and the intracranial pressure [3]. Given the complex nature of hydrocephalus, we define hydrocephalus as a disease associated with disturbances in the CSF dynamics, as in [1].

_{ss}):

_{ss}is considered a constant parameter, determined by central venous pressure. The coefficient R is the resistance to CSF reabsorbtion or outflow, measured in units of mmHg mL

^{-1}min. Storage of CSF is proportional to the cerebrospinal compliance C, measured in units of mL mmHg

^{-1}.

_{0}[8, 11], and is considered the most important law of cerebrospinal dynamic compensation [2]:

^{-1}[12]. Next, by exploiting an analogy between an electrical model of CSF compensation, as described in [5], and adapted in [2], the deterministic description of the dynamics of CSF flow are given by:

_{b}is a baseline pressure. The circuit diagram, reproduced from [2], is shown in Figure 1. An electrical circuit analogy is also used in [13] and [14] to study the dynamics of ICP in the ventricular compartments. The reference pressure parameter p

_{0}is sometimes taken to be zero, as for example, in [5] because, as noted in [2], the significance of p

_{0}is unclear. Consequently, we assume p

_{0}= 0, which results in the equation:

### The importance of modeling noise in CSF dynamics

## Methods

### Generalizing the CSF dynamics to incorporate noisy flow

Visual examination of the time-series of ICP recordings shows that the fluctuations are smooth (unlike electrons in a wire which generate shot noise, characterized by jumps [18]), and therefore continuous state space Markov processes are appropriate to capture the noisy dynamics of CSF flow. A large class of Markov processes can be represented by SDEs, and here a methodological choice must be made--noisy dynamic processes can be represented by stochastic differential equations of the Ito type or the Stratonovich type which correspond to two different ways of introducing noise into a dynamic system. A central difference between the two is that the Stratonovich SDE uses the usual deterministic calculus whereas the Ito SDE requires a completely new stochastic calculus. Extensive conceptual, empirical and philosophical discussions of this issue exist in the literature on mathematical models of electrical, biological and physical phenomena [19–21]. The overwhelming majority of these discussions conclude that Ito processes, generated by stochastic differential equations of the Ito type, are superior to Stratonovich processes, generated by stochastic differential equations of the Stratonovich type [22, 23]. Ito [24] extended standard deterministic calculus to a "stochastic calculus" applicable to functions of a wide class of continuous-time random processes, known as Ito processes. Given the SDE for the process under consideration, a result called Ito's Lemma yields the SDE driving the dynamics of a general transformation of the original process [24]. This utilitarian result allows deducing the stochastic properties of considerably complex models driven by Ito processes [23]. An essential property of Ito processes is that nonlinear functions of Ito processes remain Ito processes--a property called closure under nonlinear transformations, indispensable for practical reasons. From an empirical standpoint, a compelling advantage of Ito processes is that they often yield very precise statistical specifications for estimation [23]. An attractive property of Ito processes--on theoretical, mathematical, practical and computational grounds--is that they are Markov processes. Finally, the Ito calculus has been extended to embrace general martingale processes [25]--a development that permits joint consideration of both smooth noise and noise that occurs in jumps. Thus our modeling framework can accommodate neurological phenomena requiring noise that encompasses *both* smooth and jumpy variations in the state of the system, such as the firing of neurons [26].

### Modeling the CSF dynamics as an Ito process to incorporate noisy flow

*uncorrelated*, which substantively means that an error at one point in time does not influence the error at another point in time. This too is a standard assumption in the dynamic modeling literature.

Note that in order for equation (9) to be dimensionally consistent, the unit of σ is mL/min. Because the 'input' I(t) is the infusion rate which is under direct experimental control, therefore, in the language of control theory, I(t) is a 'control' variable. In the infusion studies conducted at Addenbrookes's Hospital in Cambridge, UK, I(t) is maintained at a constant rate of 1.5 mL/min. However, factors not within the experimenter's control also influence the input flow rate. In addition to the infusion rate of the experimenter which influences CSF formation, CSF is produced inside the brain, but much about its production remains unknown at the present time. Currently, there are no direct methods to measure the CSF production rate over short periods of time. Globally, the average secretion rate--used as a proxy for the production--is 0.35 mL/min with a 95% confidence range of 0.27 mL/min to 0.45 mL/min [2]. The lack of precise knowledge about the CSF production rate and the unmeasured factors that influence it are sources of noise in the total CSF formation rate. Consequently the stochastic Marmarou model may be conceptualized as the classic Marmarou model with a noisy input flow rate that reflects uncertainty about CSF formation.

### The clinical significance of the stochastic Marmarou model

By building the fluctuations right into the dynamics of the model structure, the stochastic model makes full use of the information in the variations of the ICP waveform. From this additional information, the time-varying probability distributions of the ICP waveform can be extracted, and it is these latter quantities that enable computation of the probabilities of clinically relevant events. It is the knowledge of these probabilities of clinically relevant events that facilitate dynamic risk management of the patient. Conceptually, the average value of p(t) at any given time 't' is the average ICP at that time in an ensemble of patients with a similar CSF flow profile, as reflected in the values of the CSF flow parameters.

### Analyzing the stochastic Marmarou model

In the Results section, we will display the exact analytical solution to the stochastic Marmarou model and derive insights from the solution into the influence of noise on the ICP at each point in time, and on average. Under the normal conditions described in [2], biological processes will settle down to a steady state after the transients have died out. In the deterministic Marmarou model [5], the steady state (equilibrium) is found by setting the time rate of change of the ICP equal to zero. What is the corresponding steady-state concept for a stochastic process? The stochastic counterpart to the time-independent steady-state level of the ICP is the time-independent probability distribution of the ICP, and the equilibrium probability distribution is to the stochastic environment as the stable equilibrium point is to the deterministic one [30]. We derive the equilibrium probability distribution for the ICP, and from it, draw conclusions for the influence of CSF flow parameters and noise intensity upon the average steady-state ICP level. We compute a measure relevant to the treatment and control of hydrocephalus: given the current value of the patient's ICP, what is the probability that it will exceed a critical high level? And how is that probability influenced by neurological characteristics of the patient such as their resistance to CSF flow and the noise intensity of the fluctuations in CSF formation rate which in turn drives the fluctuations in their ICP?

### Computing probabilities of clinically relevant events

The mathematical formulation of the problem posed in the previous paragraph is: given that a patient's ICP is currently x mmHg, where x is an arbitrary value, what is the probability that the ICP will exceed a critical threshold 'b' (mmHg) at a future time? Mathematically stated: given that p(s) = x, find the following transition probability-- P[p(t) > b | p(s) = x], t > s. Simple though the question seems, finding the answer requires computing the conditional probability distributions of the CSF process. Since the conditional probability distributions follow the Fokker-Planck partial differential equation, the problem is non-trivial, but Karlin and Taylor [31] circumvent the difficulty by solving a boundary-value problem associated with this dynamically changing probability. They show that the required probability satisfies a nonlinear ordinary differential equation which must be solved subject to two conditions on the probability that are natural consequences of the current ICP level when it is at one of the two extreme points of the range of ICP values under consideration. It is these conditions that give rise to the term 'boundary value problem.'

## Results and Discussion

We state and discuss the significance of the mathematical results, deferring their proofs to the Appendices (Additional file 1, Additional file 2 and Additional file 3) in the interest of maintaining clarity of exposition. Our first result is the exact analytical solution to the stochastic Marmarou model.

### Solution to stochastic Marmarou model with constant rate infusion

_{0}" (mmHg) at time t = 0, the future ICP value at any time "t" is given by:

The proof is provided in additional file 1: Solving the stochastic Marmarou model. Note that the solution to the stochastic Marmarou model is found through an "integrating factor" which involves an integration constant, the evaluation of which necessitates a unit of 1/min unit for the 2 inside the exponent of the exponential function. The noise intensity parameter σ and the Brownian Motion process W(t) in the solution show the explicit influence of noise on the evolution of the ICP, underscoring the importance of modeling the noise in the clinical ICP data. In addition to the practical utility of offering a closed-form analytical solution, this result has value for another reason: it shows explicitly that noise cannot be averaged away when the process is nonlinear. If the Brownian motion process W(t) entered the solution for p(t) in an additive linear way, its effect would disappear on average. But the Brownian motion process enters the solution in a highly nonlinear fashion, making it impossible to average out its effect to zero. Finally, the solution depends upon the noise intensity parameter σ in a mathematically continuous way, a fact that is meaningful because the result shows that the solution to the deterministic Marmarou model [5] emerges as the special case corresponding to σ = 0 mL/min, and so, it is natural to ask if the simpler deterministic model would suffice when the noise intensity is small. Should the influence of noise be negligible in a particular case, the value of σ will be very small, and, because of the mathematical continuity in its dependence upon σ, the stochastic solution will be very close to the deterministic solution in such a case, and we may use the simpler deterministic model with confidence. However, the stochastic model is preferable in general for two reasons: it captures the dynamics of the ICP data better than the deterministic model when the noise intensity is larger, and furthermore, the stochastic model characterizes the risk profile of the patient probabilistically. Almost tautologically, the deterministic model cannot evaluate the risks due to the errors that are an inseparable part of medical data because deterministic modeling philosophy sees the future as completely predictable from the present situation. These considerations suggest that, from a conservative modeling perspective, incorporating the influence of noise into the dynamics is conceptually more defensible.

In principle, the solution contains all the transient probability distributions of the ICP process that characterize it on its way to equilibrium. In practice, mathematical difficulties may make these transient distributions hard to extract from the solution. But we can still compute the probability of the critical events by using a methodology that does not depend on that knowledge. And we can still draw useful information about the nature of the process at steady-state. Next, we find the steady-state probability distribution of the ICP process.

### Steady-state probability distribution of ICP

The steady-state probability distribution of the ICP is gamma with the parameters shown p.149 in [29], and will exist provided that the noise intensity parameter σ satisfies the condition: ${\sigma}^{2}<\frac{2\left(\text{RI}+{\text{p}}_{\text{b}}\right)}{\text{RE}}$.

*jointly*determine its location and shape, because the mean and variance of this distribution are functions of

*both*the parameters. In practice, biological phenomena will converge to a steady-state, but nonetheless it is important to check that the technical condition stated for the existence of a steady-state distribution is satisfied by realistic values of the Marmarou model parameters. We obtained typical values of R, E, p

_{0}, p

_{b}and I from a combination of [2] and private communication with Dr. M. Czosnyka, yielding R = 7 mmHgmL

^{-1}min (reported values range between 6 and 10 mmHgmL

^{-1}min), p

_{0}= 0 mmHg and p

_{b}= 8 mmHg. Elevated elasticity is reported to be E > 0.18 ml

^{-1}and the rate of infusion is I = 1.5 mLmin

^{-1}[2]. The value of E was taken to be E = 0.15 mL

^{-1}, based on private communication with Dr. Czosnyka, and this value came from data gathered in infusion studies conducted at Addenbrooke's Hospital. p

_{b}is a baseline pressure which is different for each individual patient. Based on the ICP plots in [2], we set p

_{b}= 8 mmHg. This value is close to the average p

_{b}across all infusion studies conducted at Addenbrooke's Hospital which Dr. Czosnyka, in private communication, reported to be 6 mmHg. While the authors solved the deterministic Marmarou model for the general case of p

_{0}≠ 0 mmHg in [2], and found that the average value of p

_{0}in the infusion studies was p

_{0}= 4 (private communication with Dr. Czosnyka), the non-zero p

_{0}case is currently not analytically solvable for the stochastic Marmarou model. Our p

_{0}= 0 mmHg assumption is consistent with [5], in which the authors ignore the reference pressure. However, we acknowledge that the non-zero p

_{0}case is an important issue in mathematical modeling of hydrocephalus--and the author and his collaborators are working on an algorithm to solve the stochastic Marmarou model for non-zero p

_{0}. A typical value for σ is difficult to find since the input flow rate of CSF is not accessible to direct observation--only the fluctuations in the ICP are observable. We estimated σ roughly as, σ = 0.33*I = 0.33*1.5 mL/min = 0.5 mL/min, so that the flow fluctuations are 33% (1/3) of the flow rate. This is a rough estimate--the choice of a typical value for σ is unclear. Because of the uncertainty and approximations involved in the estimate, we did the computations in which σ was fixed, not just at σ = 0.5 mL/min, but also at values of σ lower as well as higher than 0.5 mL/min in order to check the robustness of the conclusions. We have reported the results for σ = 0.5 mL/min and for σ = 0.8 mL/min. The results of the computations of the probability as a function of R are robust across a wide range of σ, so that, even though the estimate of σ is only approximate, we can be reasonably confident about the conclusions of the analysis. To examine the influence of σ itself on the risk probability, we computed the probability across a wide range of σ as shown in Figure 4--again, in an effort to reduce the impact of our imprecise knowledge of σ upon the findings.. While our estimate of σ is only approximate, we note that there is imprecision and uncertainty about

*all*the parameter estimates, especially that of the CSF flow resistance 'R.' Estimation methods specifically developed for dynamic models are needed. In this paper, the primary objectives were to introduce SDE methodology to CSF research, demonstrate its analytical power, and show its clinical usefulness in dynamic risk management. Consequently, we used existing typical estimates of the model parameters even though some of them are imprecise and approximate. With σ = 0.5 mL/min, the condition for the existence of a steady-state probability distribution is met with ease.

Our next three results are motivated by the following considerations. A larger cerebrospinal fluid resistance R tends to increase ICP by increasing the pressure due to the circulatory CSF component. This is a direct consequence of Davson's equation [6]: ICP_{CSF} = (resistance to CSF outflow) × (CSF formation) + (pressure in sagittal sinus). This naturally leads to the following questions. How will the intensity of the fluctuations influence the relationship between resistance and ICP? The same relationship may hold on average, but, as anticipated in the solution to the stochastic Marmarou model, it may be moderated by the noise intensity parameter because of the nonlinearity of the ICP process. How will the intensity of fluctuations affect the average steady-state ICP--is the average steady-state ICP smaller or larger when the intensity of fluctuations increases? Finally, will the intensity of fluctuations attenuate or amplify the effect of resistance to CSF flow on the average steady-state ICP?

### Relationship between average steady-state ICP and cerebrospinal fluid resistance

The average steady-state ICP, denoted by μ, increases with the cerebrospinal fluid resistance R--thus the relationship between R and ICP holds on average.

The steady-state probability distribution of ICP is gamma with the parameters shown in the previous subsection. From well-known properties of the gamma distribution, it follows that the steady-state mean ICP level μ is given by: $\mu =\left(\text{RI}+{\text{p}}_{\text{b}}\right)-\frac{\text{RE}{\sigma}^{2}}{2}$. Therefore, $\frac{\partial \mu}{\partial \text{R}}=\text{I}-\frac{{\sigma}^{2}\text{E}}{2}$. From the expression for $\frac{\partial \mu}{\partial \text{R}}$, it is clear that the average ICP level does indeed increase with R, provided that $\frac{{\sigma}^{2}\text{E}}{2}<\text{I}$. This condition is satisfied, using the values of the parameters in the previous subsection. Thus, the increasing relationship between the actual ICP level and the cerebrospinal fluid resistance, predicted by Davson's equation when ICP is conceptualized as a deterministic process, also holds *on average* at steady-state when ICP is modeled as a stochastic process.

### Relationship between average steady-state ICP and noise intensity

The average steady-state ICP level, decreases with the intensity of fluctuations, measured by the infinitesimal variance parameter σ^{2}.

From the relationship derived in the previous subsection, $\mu =\left(\text{RI}+{\text{p}}_{\text{b}}\right)-\frac{\text{R}{\sigma}^{2}\text{E}}{2}$, it is clear that μ decreases as σ^{2} increases. A larger noise intensity corresponds to greater variation in the CSF input flow rate which translates into greater variation in ICP, and these larger fluctuations could cause the average ICP level to increase, decrease or remain unaffected. The nonlinear influence of the parameters of CSF flow dynamics on ICP level turns out to reduce the average ICP value when the fluctuations in ICP are greater. This is an outcome that one would expect to find when steady-state has been achieved--when the transition probabilities have settled down to constant levels so that the probability distribution of ICP is no longer changing over time. This mathematical finding could be tested by separating a random sample of patients into two groups, such that one group has more variability in its ICP levels (due to higher variability in its CSF input flow rate) than the other group, and then conducting a statistical test of significance--such as a t-test--on the difference in mean ICP levels in these two groups at steady-state.

### Effect of noise intensity on the relationship between average steady-state ICP and cerebrospinal fluid resistance

The resistance increases the ICP on average by a smaller amount when the intensity of fluctuations is higher.

From $\mu =\left(\text{RI}+{\text{p}}_{\text{b}}\right)-\frac{\text{R}{\sigma}^{2}\text{E}}{2}$, it is clear that a higher σ^{2} will dampen the effect of the cerebrospinal resistance on the average steady-state ICP level. This is an outcome that one would expect to find at steady-state. The mathematical finding could be tested by separating a random sample of patients into two groups, such that one group has more variability in its ICP levels than the other group (due to higher variability in its CSF input flow rate), and then correlating the mean ICP level with the cerebrospinal resistance in each group at steady-state. According to the mathematics, the correlation should be smaller in the group with more variable ICP. Given the linear relationship between the steady-state mean and the cerebrospinal resistance, a simple correlation coefficient such as the Pearson product moment should suffice.

Next we turn our attention to dynamic management of the patient's risk. Risk may be quantified in terms of the probability of the onset of some critical event, say the ICP exceeding a dangerously high level. Given the current value of the patient's ICP, what is the probability that it will exceed a high level? Such a probability is intrinsically dynamic because it depends upon the patient's current condition (their current ICP), the dynamics of the patient's CSF flow and the noise intensity σ^{2}. We want to understand how the probability is influenced by important clinical characteristics of the patient such as their resistance to CSF flow, and by the noise intensity.

### Computing clinically relevant dynamic probabilities

The integrals defining s(x) and S(x) are indefinite at the lower end because the final answer is unaffected by its choice. For our clinical applications, it is natural to take the lower end point to be zero.

*quasi-analytical*solution because the integral that defines s(x) cannot be obtained in closed-form. However, that is no limitation because we can integrate it numerically after substituting the empirically established values of the parameters. We used the parameter values shown in the subsection "Steady-state Probability Distribution of ICP." We took the critical level 'b' to be 40 mmHg, based on the clinical finding reported in [2] that patients were able to tolerate increases in ICP up to 40-50 mmHg. Our rationale was that, from a clinical perspective, a conservative approach to patient management would be consistent with assessing the probability of reaching the lower end point of the 40-50 mmHg range that patients are able to tolerate. Thus our critical event is defined as "reaching an ICP of 40 mmHg." In order to understand how the probability is influenced by the noise intensity parameter σ, we computed the probability over a range of σ = 0.4 mL/min to σ = 1.3 mL/min. For each value of σ in this range, we solved the boundary-value problem to find the probability of reaching 40 mmHg. Furthermore, in order to understand the influence of the patient's initial condition upon the probability of the critical event, we repeated this set of computations for three different starting levels of ICP; the curve shown in Figure 4 is for a starting level of ICP of 35 mmHg. In order to understand how the probability is influenced by the resistance to CSF outflow R, we computed the probability over a range of R = 4 mmHgmL

^{-1}min to R = 12 mmHgmL

^{-1}min. For each value of R in this range, we solved the boundary-value problem to find the probability of reaching 40 mmHg. Again, we repeated this set of computations for three different starting levels of ICP; the curve shown in Figure 5 is for a starting level of ICP of 35 mmHg. Across the three initial levels of ICP, the curves have a similar shape and are simply translated vertically. Figures 4 and 5 show that the probabilities increase at an increasing rate (convex functions). Furthermore, they tell an interesting neurological story--namely, that the probabilities of the critical events exhibit strong

*threshold*effects. In Figure 4, below a critical level of noise intensity, the probabilities are very low--almost zero--but beyond a threshold value of σ = 1.1 mL/min in Figure 4, they rise steeply. In Figure 5, as R (mmHgmL

^{-1}min) varies from 4 to below 10, the probabilities are almost zero, but beyond R = 10 mmHgmL

^{-1}min, they rise dramatically. Furthermore, at low levels of noise intensity, the probabilities remain close to zero

*throughout*the range of R (mmHgmL

^{-1}min) from 4 to 12. But as σ increases to 0.8 mL/min--the value assumed for it in Figure 5--R has a strong effect on the probability beyond the critical threshold of 10 mmHgmL

^{-1}min. The clinical significance of these findings is that erratic fluctuations in ICP (caused by a larger input flow rate noise intensity σ) will significantly increase the patient's risk, as measured by the probability of the critical event. Because the risk increases rapidly beyond the threshold value of σ, these results suggest that an essential component of risk management is to carefully minimize erratic fluctuations in the patient's CSF input flow rate at all times. Finally, Figure 6 shows the probability of the critical event as a function of both R and σ in a three-dimensional plot. The two-dimensional surface shows the value of the probability for each combination of values of R and σ. To facilitate interpretation of the surface, we used a mesh in which the dark lines are the probability plots as a function of R and the red lines are the probability plots as a function of σ. Figure 6 shows very clearly that threshold effects are sensitive to both R and σ, and beyond the threshold, the probabilities asymptotically approach one.

## Conclusions

The stochastic generalization of the Marmarou model offers a tractable analytical description of the noisy ICP dynamics and yields insights into the impact of noise. The SDE offers a rigorous analytical framework to study issues of clinical interest and neurological significance such as the patient's risk. A key clinical implication is that fluctuations in the CSF formation rate--which increase the fluctuations in ICP-- should be minimized to lower the patient's risk. Future work could extend the framework developed in this research to accommodate the non-zero reference pressure case. Finally, the stochastic differential equation framework, in conjunction with nonlinear control theory, can be used to develop a nonlinear automatic controller to regulate shunts to facilitate continuous CSF drainage.

## Declarations

### Acknowledgements

The author is grateful to Dr. Marek Czosnyka, Neuroscience Group, University of Cambridge, UK, for introducing him to the key issues in CSF dynamics and hydrocephalus, for making data available from infusion studies conducted at Addenbrooke's Hospital, UK, and to the late Professor Anthony Marmarou for suggesting the probabilistic computations for critical events.

## Authors’ Affiliations

## References

- Eklund A, Smielewski P, Chambers I, Alperin N, Malm J, Czosnyka M, Marmarou A: Assessment of cerebrospinal fluid outflow resistance. Med Bio Eng Comput. 2007, 45: 719-735. 10.1007/s11517-007-0199-5.View ArticleGoogle Scholar
- Czosnyka M, Czosnyka Z, Momjian S, Pickard JD: Cerebrospinal Fluid Dynamics. Physiol Meas. 2004, 25: R51-R76. 10.1088/0967-3334/25/5/R01.View ArticlePubMedGoogle Scholar
- Cohen B, Voorhees A, Vedel S, Wei T: Development of a theoretical framework for analyzing cerebrospinal fluid dynamics. Cerebrospinal Fluid Res. 2009, 6: 12-10.1186/1743-8454-6-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Sivaloganathan S, Tenti G, Drake JM: Mathematical pressure volume models of the cerebrospinal fluid. Appl Math Comp. 1998, 94: 243-266. 10.1016/S0096-3003(97)10093-5.View ArticleGoogle Scholar
- Marmarou A, Shulman K, Rosende RM: A nonlinear analysis of the cerebrospinal fluid system and intracranial pressure dynamics. J Neurosurg. 1978, 48: 332-344. 10.3171/jns.1978.48.3.0332.View ArticlePubMedGoogle Scholar
- Czosnyka M, Pickard JD: Monitoring and Interpretation of Intracranial Pressure. J Neurol Neurosurg Psychiatry. 2004, 75: 813-821. 10.1136/jnnp.2003.033126.PubMed CentralView ArticlePubMedGoogle Scholar
- Marmarou A: A theoretical model and experimental evaluation of the cerebrospinal fluid system. Thesis. 1973, Drexel University, Philadelphia, PAGoogle Scholar
- Avezaat CJJ, Eijndhoven JHM: Cerebrospinal fluid pulse pressure and craniospinal dynamics: a theoretical, clinical and experimental study. Thesis. Edited by: Jongbloed A. 1984, The HagueGoogle Scholar
- Frieden H, Ekstedt J: Instrumentation for cerebrospinal fluid hydrodynamic studies in man. Med Biol Eng Comput. 1982, 20: 167-80. 10.1007/BF02441352.View ArticleGoogle Scholar
- Feynman R: The Character of Physical Law. 1965, Cambridge, MA: MIT PressGoogle Scholar
- Raabe A, Czosnyka M, Piper I, Seifert V: Monitoring of intracranial compliance: correction for a change in body position. Acta Neurochir (Wien). 1999, 141: 31-36. 10.1007/s007010050263.View ArticleGoogle Scholar
- Czosnyka M, Batorski L, Roszkowski M, Tomaszewski J, Wocjan J, Walencik A, Zabolotny W: Cerebrospinal compensation in hydrocephalic children. Child's Nerv Syst. 1993, 9: 17-22.View ArticleGoogle Scholar
- Rekate HL: Circuit diagrams of the circulation of cerebrospinal fluid. Pediatr Neurosurg. 1994, 21: 248-253. 10.1159/000120845.View ArticlePubMedGoogle Scholar
- Rekate HL, Brodkey JA, Chizeck HJ, el Sakka WE, Ko WH: Ventricular volume regulation: a mathematical model and computer simulation. Pediatr Neurosci. 1988, 14: 77-84. 10.1159/000120367.View ArticlePubMedGoogle Scholar
- Clarke MJ, Meyer FB: The history of mathematical modeling in hydrocephalus. Neurosurg Focus. 2007, 22 (Suppl 4): 1-5. 10.3171/foc.2007.22.4.4.View ArticleGoogle Scholar
- Lavinio A, Czosnyka Z, Czosnyka M: Cerebrospinal Fluid Dynamics: Disturbances and Diagnostics. Eur J Anaesthesiol. 2008, 25 (Suppl 42): 137-141. 10.1017/S0265021507003298.View ArticleGoogle Scholar
- Petrella G, Czosnyka M, Keong N, Pickard JD, Czosnyka Z: How does CSF dynamics change after shunting. Acta Neurol Scand. 2008, 118: 182-188. 10.1111/j.1600-0404.2008.01041.x.View ArticlePubMedGoogle Scholar
- Gershenfeld N: The Physics of Information Technology. 2000, Cambridge: Cambridge University PressGoogle Scholar
- Van Kampen NG: Ito versus Stratanovich. J Statist Phys. 1981, 24: 175-187. 10.1007/BF01007642.View ArticleGoogle Scholar
- Turelli M: Random Environments and Stochastic Calculus. Theor Pop Biol. 1977, 12: 140-178. 10.1016/0040-5809(77)90040-5.View ArticleGoogle Scholar
- Turelli M: A Reexamination of stability in randomly varying versus deterministic environments with comments on the stochastic theory of limiting similarity. Theor Pop Biol. 1978, 13: 244-267. 10.1016/0040-5809(78)90045-X.View ArticleGoogle Scholar
- Kailath T, Poor HV: Detection of stochastic processes. IEEE Transactions on Information Theory. 1998, 44: 2230-2231. 10.1109/18.720538.View ArticleGoogle Scholar
- Lo AW: Maximum likelihood estimation of generalized Ito processes with discretely sampled data. Econometric Theory. 1988, 4: 231-247. 10.1017/S0266466600012044.View ArticleGoogle Scholar
- Ito K, Varadhan SRS, Stroock DW: On stochastic differential equations. Memoirs of the AMS. 1951, 4: 1-51.Google Scholar
- Skorohod AV: Studies in the Theory of Random Processes. 1965, New York: Dover PublicationsGoogle Scholar
- Harrison LM, David O, Friston KJ: Stochastic models of neuronal dynamics. Philos Trans R Soc Lond B Biol Sci. 2005, 360: 1075-1091. 10.1098/rstb.2005.1648.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsoularis A, Wallace J: Analysis of logistic models. Math Biosci. 2002, 179: 21-55. 10.1016/S0025-5564(02)00096-2.View ArticlePubMedGoogle Scholar
- Øksendal B: Stochastic Differential Equations: An Introduction with Applications. 2003, Berlin: SpringerView ArticleGoogle Scholar
- Gard TC: Introduction to Stochastic Differential Equations. 1988, New York: John Wiley & SonsGoogle Scholar
- May RM: Stability in randomly fluctuating versus deterministic environments. Am Nat. 1973, 107: 621-650. 10.1086/282863.View ArticleGoogle Scholar
- Karlin S, Taylor HM: A Second Course in Stochastic Processes. 1981, New York: Academic PressGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.