|Home | About | Journals | Submit | Contact Us | Français|
Investigators of homeopathy have proposed that nonlinear dynamical systems (NDS) and complex systems science offer conceptual and analytic tools for evaluating homeopathic remedy effects. Previous animal studies demonstrate that homeopathic medicines alter delta electroencephalographic (EEG) slow wave sleep. The present study extended findings of remedy-related sleep stage alterations in human subjects by testing the feasibility of using two different NDS analytic approaches to assess remedy effects on human slow wave sleep EEG.
Subjects (N=54) were young adult male and female college students with a history of coffee-related insomnia who participated in a larger 4-week study of the polysomnographic effects of homeopathic medicines on home-based all-night sleep recordings. Subjects took one bedtime dose of a homeopathic remedy (Coffea cruda or Nux vomica 30c). We computed multiscale entropy (MSE) and the correlation dimension (Mekler-D2) for stage 3 and 4 slow wave sleep EEG sampled in artifact-free 2-minute segments during the first two rapid-eye-movement (REM) cycles for remedy and post-remedy nights, controlling for placebo and post-placebo night effects.
MSE results indicate significant, remedy-specific directional effects, especially later in the night (REM cycle 2) (CC: remedy night increases and post-remedy night decreases in MSE at multiple sites for both stages 3 and 4 in both REM cycles; NV: remedy night decreases and post-remedy night increases, mainly in stage 3 REM cycle 2 MSE). D2 analyses yielded more sporadic and inconsistent findings.
Homeopathic medicines Coffea cruda and Nux vomica in 30c potencies alter short-term nonlinear dynamic parameters of slow wave sleep EEG in healthy young adults. MSE may provide a more sensitive NDS analytic method than D2 for evaluating homeopathic remedy effects on human sleep EEG patterns.
Homeopathy is a system of complementary and alternative medicine (CAM) that has generated much controversy over potential mechanisms of its purported effects as well as the nature of the animal-, mineral-, or plant-derived medicines.1 Some researchers have even questioned the validity of commonly-held practice theory tenets such as the existence of homeopathic aggravations (transient worsening prior to clinical improvement).2 Finding analytic techniques to detect and monitor subtle fluctuations in the course of remedy effects would improve understanding of the circumstances under which remedy responses are state-dependent, bidirectional, or even oscillatory in living systems.3–6
Time is a major feature of the homeopathic remedy response. Once the “correct” remedy, or simillimum, is administered, clinical responses change over time. The direction of these changes are sometimes reportedly not only multi-system, but also biphasic, i.e., nonlinear. That is, homeopathic patients may exhibit transient worsening prior to improvement in some symptoms and/or gradual reductions in susceptibility to and severity of flares in chronic conditions as other symptoms come to the fore.7 The changes appear to involve multiple subsystems of the person as well as the global state of the person as a whole. Treatment itself is an iterative individualized process using one remedy at a given moment, guided over time by repeated reassessments of the changes in the total clinical presentation.
Over the minutes, hours, days, and weeks following one dose of the remedy in an ill person, the clinical course evolves across many subsystems of the individual, at a pace related in part to the chronicity or acuity of the presenting problems and in a spatial pattern suggestive of self organization in the person as a complex system.8 Hering's Law of Cure in homeopathy predicts an improvement pattern for the “center of gravity of the illness” to shift within the person over time after initiation of remedy treatment from inside out (more important to less important organs), above downward (head toward toes), and from most recent to oldest onset symptoms in time.7,9 In contrast to weekly psychotherapy, homeopathic follow up visits for chronic diseases follow a more medical model, typically spaced 1–2 months or more apart to allow the clinical course to unfold after administration of a single dose.7,10 However, the time course is much faster and compressed in acute conditions1 and perhaps in proving or pathogenetic trials.11,12
Various contemporary scientists who study the nature of the clinical claims in homeopathy have postulated that medicines interact with the patient's state13 by modulating the nonlinear dynamics of the individual as a complex living system.5,6,8,14–24 Nonlinear dynamical systems (NDS) is a branch of the emerging holistic science of complex systems, chaos theory, and network models for physiological and social systems.25–27 NDS appears more applicable to the nature of homeopathic clinical change patterns than conventional medical linear models.18,22,23 For instance, in parallel with the complex nature of homeopathic healing outlined above, experts characterize NDS theory as “the study of how complex processes unfold over time…”.25 A complex system is more than the sum of its parts because of its capacity to generate emergent properties from the interactions of the interconnected parts through their relationships to one another and the larger whole.27–31
It is possible to use complex systems methods to study any level of organizational scale, from biochemical to organ system to organism to societal behaviors.26,31–33 However, these methods usually require sufficiently long time series of data points to evaluate the recurring patterns of change over time in a given behavior. Electroencephalographic (EEG) activity is one type of biological signal to which NDS analytic methods can apply.34–40
For research in a controversial field such as homeopathy, EEG as an outcome measure has advantages.41,42 Neurophysiological assessments enable utilization of more objective and sensitive tools than crude behavioral observations or self-report scales for assessing the real-time dynamics of an organism at integrative levels of functioning such as the central nervous system.19,43 In animal research, for example, Ruiz-Vega et al44 showed that homeopathically prepared medicines such as Coffea cruda at 30c potencies can produce increases in the electroencephalographic delta band power of slow wave sleep (stages 3 and 4) versus placebo. Other placebo-controlled animal studies revealed that (a) homeopathic histamine 30c can mobilize a dynamic biphasic process affecting sleep EEG delta activity that evolves over the course of the night between two locally stationary states according to a power law45,46; and (b) homeopathic Coffea cruda 30c produces different patterns of changes in parietal EEG delta power later versus earlier in the sleep period in a subset of animals pre- or post-treated with caffeine.46 The homeopathic remedy Nux vomica 30c has also produced alterations of sleep in other animal studies.47,48
The purpose of the present study was to explore the feasibility of applying nonlinear dynamical analytic methods to human slow wave sleep EEG from a placebo-controlled study on the effects of two different homeopathic remedies, Coffea cruda 30c and Nux vomica 30c, in relatively healthy young adult human subjects with histories of coffee-induced insomnia (conventional sleep stage findings reported elsewhere49). For the analysis, we adapted previously-published NDS methods for computation of both multiscale entropy (MSE)50–52 and EEG correlation dimension (Mekler method D2).39,53–55
The primary hypothesis was that the verum would exert a pattern of short-term modulatory effects on nonlinear dynamics of the slow wave sleep EEG (stages 3 and 4, which are by definition predominantly, but not exclusively, delta and theta spectral EEG frequencies) different from those seen on baseline nights or with placebo.8,56 The secondary hypothesis, based on homeopathic theory20,42,57 and accumulating basic science data on nanoparticles and succussion-related unique alterations in the dynamic structure of the ethanol-water solvent networks of specific remedies,58–69 was that Coffea cruda and Nux vomica would differ from one another in the pattern of their nonlinear EEG effects. Skeptics would argue the null hypothesis, given the 30c doses (potencies) used of both medicines would have been diluted past Avogadro's number and would, from a diluted bulk form source molecule perspective, contain no active ingredients and hence, effects no different from placebo.70
Subjects were young adult college students of both sexes (N=54, 50% female, mean age 20 years) recruited from an introductory psychology class after completing a screening questionnaire packet. Homeopathic clinical texts suggest that people who benefit from Coffea cruda are sensitive to all types of stimuli, including anxiety; whereas individuals who respond to Nux vomica are intensely driven, competitive, impatient Type A personalities.57 Both remedy types are reportedly associated with insomnia as an adverse effect of coffee. Thus, to increase the possibility of enrolling remedy-responsive subjects71,72, inclusion criteria were not only moderate to good self-rated global health and a history of coffee-induced insomnia in the past, but also an above-mean elevation of scores (using the class member scores of the first 1,036 students screened) for either the Cook-Medley Hostility Scale (high hostility trait, CMHO) or the Anxiety Sensitivity Index (high anxiety sensitivity trait, ASI). To be included in the polysomnography study, subjects with an above-mean value for CMHO or ASI had to score below the sample mean for the other scale (see also49). In order to control for the personality variables in regression analyses, a single combined personality index (CMASI) was created in which higher scores indicate elevated Type A traits and lower scores indicate elevated anxiety sensitivity traits. Specific findings relating to the interaction of personality type with remedy received are also published elsewhere.13
Because of the belief by some homeopaths that beverage coffee (but not caffeine) sometimes antidotes homeopathic medicines,7 subjects had to be willing to eliminate drinking coffee for the full duration of the study, although they were encouraged to change to and stabilize prior to the first baseline sleep recording night for at least one week on non-coffee caffeinated beverages to maintain their customary daily caffeine intake without interruption, if necessary. Exclusion criteria were pregnancy or plan to become pregnant, major psychiatric or serious chronic medical conditions, a history of anaphylactic shock, and/or chronic use of medications other than contraceptive drugs. This study was reviewed and approved by the Institutional Review Board of the University of Arizona. Subjects were paid for participation.
All participants underwent a total of eight all-night sleep recordings in their home, spaced over 4 weeks in pairs of successive days (same days of each week). The recording nights were as follows: Week 1 Baseline=nights 1,2; Week 2=nights 8,9 (placebo on night 8); Week 3 Baseline=nights15,16; Week 4=nights 22,23 (remedy on night 22). In specific, Weeks 1 and 3 involved baseline recordings. On Week 2, subjects dissolved 3 #30 lactose-sucrose placebo pellets under the tongue at bedtime (night 8 only) in the presence of the sleep recording technician. On Week 4, subjects dissolved 3 #30 lactose-sucrose verum remedy pellets under the tongue at bedtime (night 22 only) in the presence of the sleep recording technician.
Subjects were blinded to the ingredients of the pellets throughout the study. Staff were aware that the Week 2 pellets were placebo but fully blinded as to the identity of the remedy given in a randomly assigned manner on Week 4 (verum remedy: either Coffea cruda 30c or Nux vomica 30c). Medicines and placebos were purchased from an FDA-regulated homeopathic pharmacy, Hahnemann Laboratories, San Rafael, CA USA. The lactose-sucrose pellets used as placebos and on which specific remedy solution was sprayed per usual remedy preparation procedures were tested and certified as pure and contaminant-free by an independent testing laboratory.
Sleep recordings (polysomnograms including bilateral electrooculograms (EOG) and mental/submental electromyography (EMG), with EEG electrode sites at C3, C4, Cz, Pz, O1 and O2, referenced to contralateral mastoids (A1/A2)) were performed using Cadwell Laboratories (Kennewick, WA, USA) Easy Ambulatory 2 equipment. Electrode impedance levels were kept below 5 kiloohms. Data acquisition and data processing for sleep stage scoring involved an EEG sampling rate of 200 Hz, with a low pass filter set at 35 Hz and a high pass filter set at 0.53 Hz (permitting sleep EEG frequency band acquisition in the physiological range between 0.53 and 30 Hz). A 60 cycle notch filter was used to eliminate ambient electrical noise. Sleep stage data were scored for sleep stages using standard criteria (relying on EEG from C3/A2 or C4/A1, EOG, and EMG criteria) from the Rechtschaffen-Kales manual73 by a trained sleep technician blinded as to day and condition. Bell et al (2011) report the details of the conventional sleep stage methodology and findings elsewhere.49
For purposes of the present study, the technician then identified artifact-free 2-minute segments of EEG data within the previously-scored slow wave sleep stages (by definition, Stages 3 and 4, in which delta sleep is more commonly observed) during the first two rapid-eye-movement (REM) cycles of each of the eight recording nights (each REM cycle averages 90 minutes in length and includes some REM sleep stage events) for NDS analysis. Although the sleep stage scoring itself was based on conventional criteria involving in part the EEG C4/A1 or C3/A2 placements together with EOG and EMG,49 the concurrent EEG data from all electrode sites within the selected data segments were used for the current nonlinear dynamical exploratory analyses.
To control for baseline sleep patterns, results of baseline nights 1, 2, 15, and 16 were included for use in the regression models. Because of the animal evidence that the study medicines might have biphasic effects over time on any given biomarker3,46 and findings suggesting changes in normal human sleep EEG complexity over course of a night,39 we examined data from the first two REM cycles of the sleep period for both the night of and the night after taking the remedy, controlling for both baseline and placebo nights.
The analysis focused on the findings from nights 22 (verum remedy night) and 23 (post-remedy night). In particular, we examined the findings for night 22 and the change score from night 22 to night 23. We used a random effect regression model for analysis, with person as the random effect. The regression equation used the subject's complexity score (MSE or D2) as the dependent variable, and included dummy coded effects for each week of the study (baseline, placebo, second baseline, remedy), the first and second observed nights of the study week, whether the subject was given pellets or not, and interactions with the personality variables and whether they received the remedy Nux vomica or Coffea cruda. For example, we tested whether a subject's remedy night MSE differed from their placebo night MSE, and whether there was significant change in MSE on the post-remedy night compared to the post-placebo night change. Each comparison adjusted for each subject's own baseline MSE values and personality score. STATA 9.1 or 10.0 were used for all data reduction and analyses. This is an overall multivariate regression-based analytic approach.
Multiscale entropy (MSE) analysis50,52,74 is a method for computing how a certain entropy measure varies as one successively coarsens a time series over multiple scales, by computing averages over non-overlapping intervals of fixed size (the scale). We applied this to our 2-minute, artifact-free EEG signals, using 20 time scales, following Costa's example analyses. Costa has put her program into the public domain (http://physionet.org), and we translated it into somewhat less ornate code (see Appendix for specifics) for the large number of analyses we needed to perform. Although it may have been preferred to validate the original Costa C-language code involved, funding limitations prevented our doing so. As an alternative for this exploratory study, we have provided the code in the Appendix, for future validation and confirmatory studies.
In addition to computing entropies at 20 scales for an actual series, we then permuted the series at random 5 times, each time re-computing all 20 entropies. The final outcome measure was the area between the actual entropy-scale curve and the average of the 5 entropy-scale curves for the permuted data. Higher values of this measurement are taken to indicate higher scale-free entropy, which in turn may be due either to complex dynamics having effects over long time-scales, or to noise in the data. This duality of interpretation has been remarked by Costa.
Correlation dimension (D2) is considered one of severa key measures in chaotic time series analysis, assessing the degree of similarity between one observation and others later in the time series (http://www.societyforchaostheory.org/tutorials/00006/GlossaryTerms.html, accessed 12/21/11). Prior research indicates that D2 analysis of human EEG can distinguish between normals and patients with neurodegenerative disorders.37 Some recent research indicates that, at least in waking EEG, changes in D2 may relate to change in the theta frequency band (also a feature of slow wave sleep, i.e., stages 3 and 4).38
For computation of the correlation dimension D2 (originally developed by Grassberger and Procaccia 75), we closely followed the EEG analytic approach of Mekler.53 He advocates embedding in 12-dimensional space with relatively large lags, computing the correlation integral for a sequence of small parameter values, and then computing the slope from the middle third of the resulting curve. Mekler acknowledges that the resultant values are not D2 in the strict sense, but still capture nonlinear variability of brain functioning (hereafter, D2 in this paper refers to the calculated Mekler-D2). The advantages of the Mekler technique are that (a) EEG data per se does not meet all of the necessary assumptions or requirements for applying the Grassberger-Procaccia method without requiring subjective evaluation procedures for selecting calculation parameters (leading to difficulties in replication between different laboratories,53 an issue already a known challenge in other types of homeopathic laboratory research76); (b) Mekler's method enables capacity to process large amounts of experimental data.
We confirmed in test cases that his method appeared to apply to our data, in the sense that the correlation integral curves were regular in form, and the slope was a reasonable representation of the middle of the curve. Like Mekler, we also observed high correlation dimensions, near to the dimension of the embedding space. The conventional interpretation of a correlation dimension is that it is the (possibly fractional) dimension of the attractor of the nonlinear dynamic system, in an appropriate representation space. The suggestion is thus that lower correlation dimension represents simpler dynamics or lower complexity.
As in the MSE analysis, higher correlation dimension can also represent the presence of noise. These limitations are inherent in the state-of-the-science of various nonlinear dynamical systems computational techniques at present. Nonetheless, our attempt was to replicate methods already in the literature for EEG. The correlation dimension is a well-defined quantity related to certain nonlinear dynamical features of a complex system. We used the Physionet program (http://physionet.org) for computing the correlation integral, and then used the regression slope of the middle third (Mekler used the slope based on only two points, which is less accurate).
Of the 54 participants, one had unusable EEG and was dropped from the analysis. Because the analysis, as described above, requires complete data (all 8 nights for the given REM cycle and sleep stage), missing outcomes were linearly interpolated when that was possible, and carried forward or backward when interpolation was not possible. Data completion rates were about 70% for REM cycle 1, 60% for REM cycle 2, stage 3, and 45% for REM cycle 2, stage 4. This was an accumulation of missing nights for some participants, and unavailability of quality signal segments for others. The effect of the imputation of interpolated values should be to suppress any effects on the remedy night 22 or post-remedy night 23. Carrying out analysis without imputation was regarded as infeasible, since then different nights would consist of different subsamples, potentially creating analysis artifacts.
Table 1 gives the demographic information on the participants. Overall there were no baseline differences in age, gender, CMHO, ASI, or CMASI between the Coffea cruda and Nux vomica groups. Twenty-six subjects were randomized to receive Coffea cruda, and 28 subjects received Nux vomica. Conventional sleep stage findings are reported elsewhere.49 Tables 2a,b and 3a,b summarize the findings for the MSE and D2 analyses at each electrode site for each of the two REM cycles of remedy night (night 22) and the change from night 22 to post-remedy night 23, after controlling for factors noted above (including slow wave sleep EEG from baseline nights, placebo night, and change from placebo to post-placebo night, as well as personality trait index).
MSE revealed more extensive patterns of significant changes over multiple electrode sites for both medicines than did D2 analyses. For both remedies, the findings emerged later in the sleep period (REM cycle 2 more than REM cycle 1). For Coffea cruda 30c, MSE in stages 3 and 4 increased on night 22 remedy administration and then changed from night 22 to post-remedy night 23 with a decrease, especially toward the posterior (Pz) and right side (C4, O2) areas of the head (Table 2a). In contrast, for Nux vomica 30c, MSE findings were significant primarily for stage 3 (still REM cycle 2), but not stage 4 sleep (Table 3a). In the Nux vomica MSE data, the directionality was opposite to that of the Coffea effects. That is, Nux vomica led to a decrease in MSE on remedy night 22 and an increase from night 22 to night 23. Figures 1 and and22 provide a graphical presentation of MSE unstandardized beta coefficient values at the central parietal electrode site (Pz) to highlight the unique findings for each remedy. REM cycle 1 to REM cycle 2 shows change over time within the night for each sleep stage.
On the other hand, the Mekler-D2 analyses detected fewer, isolated site-specific significant changes than did MSE. With Coffea cruda, for instance, D2 showed a set of significant changes (increases from remedy night 22 to post-remedy night 23) toward the posterior electrode sites in Stage 4 REM cycle 2 (at Pz, O1, and O2), with similar significant effects in Stage 3 REM cycle 1 at Pz and O1 (Table 2b). Nux vomica, however, resulted in more widely scattered, isolated significant changes from remedy night 22 to post-remedy night 23 (increased values) in Stage 3 (REM cycle 1, at C3 and O1) and Stage 4 (REM cycle 1 at Cz and Pz; REM cycle 2 only at Cz)(Table 3b).
The multiscale entropy data showed that two homeopathic remedies in 30c potencies can induce short-term alterations, with carryover effects the next night in the nonlinear dynamics of human slow wave sleep EEG. Coffea cruda had a more extensive effect than did Nux vomica over multiple electrode sites. The direction of the MSE effects was opposite to those seen with Nux vomica. The Mekler-D2 findings were more isolated over scattered electrode sites for both medicines. The direction of changes for MSE and D2 also exhibited opposite patterns from one another within a given remedy.
It might appear that the pattern of significant electrode sites and the direction of changes for MSE and Mekler-D2 suggest differing conclusions for the extent and direction of remedy effects in terms of “complexity ” increasing versus decreasing. However, higher values of MSE and D2 are each considered to reflect different aspects of EEG complexity, i.e., within NDS, MSE evaluates similarity across time scales, whereas Mekler-D2 assesses attractor dimension. Moreover, bifurcations, bipolarity and variability in the direction of change are common and expected features of nonlinear dynamical systems.27,77 In an experimental context, the net direction depends not only on the history, sensitivity to initial conditions, and dynamics of the system itself, but also on the nature of a given nonlinear parameter under consideration (www.SocietyForChaosTheory.org; www.NECSI.org; www.PlexusInstitute.org; www.SantaFe.edu).25
Thus, it is not surprising to find different results using different NDS variables, as in this exploratory study. It will be important in future homeopathy research to choose more than one way to quantify the characteristics of complexity in a given system in order to more fully describe the nature of the treatment effects. “Complexity” of a living system such as a human brain is not a single uniform process or property easily captured with a single variable.78–80 This study is admittedly only a first step toward evaluating the feasibility of using nonlinear dynamical systems approaches to evaluate remedy effects on human sleep EEG.
Because of the inherently changeable nature of nonlinear dynamical systems,27 we also anticipate that reproducibility could be an issue in subsequent investigations. That is, although EEG effects different from those of placebo should still occur in replication studies, the specific affected electrode sites and directionality of remedy-related effects could differ from study to study. Individual differences between study participants in traits and dynamic states, 19,23 as well as possible variations in remedy preparation methods between homeopathic manufacturers 69 and/or batches by the same manufacturer67 could all contribute to inter-experimental variability. More systematic biological and psychological characterization of study participants at rest and under stress, as well as of the medicines themselves could improve reliability of study outcomes. 23,24
Despite these reservations, the findings support the primary and secondary hypotheses that homeopathically-prepared remedies at 30c can affect nonlinear slow wave sleep EEG dynamics beyond effects seen at baseline or with placebo administration. Such data further support the value of examining slow wave sleep EEG parameters in homeopathy. Previous animal studies also observed remedy effects on slow wave sleep EEG patterns.44–48,81,82 The convergent data indicate that additional research in this area is warranted. Because of the controversies surrounding homeopathy, EEG offers a promising objective biomarker for use in preclinical and clinical studies that subjective measures fail to provide.
An underlying question that NDS and other complex systems analytic methods attempt to address is whether or not there are self-organized repeating patterns detectable within a given time series or set of datapoints. That is, if one zooms in and out of different levels of organizational scale,32 does the systemic behavioral pattern repeat? Therefore, NDS analytic methods such as MSE evaluate whether or not the data reflect uncorrelated random signals or noise as opposed to “meaningful structural richness”83 in temporal organization. Other complex systems analytic methods focus more on structural organization32 (http://physionet.org/physiotools/mse/, accessed 02/01/12).
Homeopathic theory7,84,85 suggests that the simillimum should evoke similar themes across global and local levels of scale in the patient.24,43 This feature of repeating patterns in nonlinear systems and networks29 could be particularly relevant to the way in which homeopathic clinicians evaluate their patients. In clinical analysis of a case, for example, the homeopath is often looking for repetitive patterns and similar themes across the patient as a whole (global behavior of the organism at the general level of scale) as well as the more local subsystems of the body (e.g., the descriptive qualities or modalities of the physical symptoms at the organ system level of scale).7,83,84
Given the small sample size and other limitations of the study, the conservative interpretatio of the findings is that it is feasible to use nonlinear dynamical analyses to examine the effects of homeopathic remedies on human sleep patterns. MSE appears more sensitive to revealing homeopathic remedy effects than does the Mekler-method D2 computation for sleep EEG analyses, although both approaches reveal remedy effects beyond placebo. Although it would be desirable to perform explicit between-subject comparisons of the subgroups receiving the two different remedies (Nux vomica and Coffea cruda), the small sample limited power for such a statistical approach. Therefore, we limited our computations to within-subject analyses. Subsequent studies with larger samples should directly evaluate the between-subject effects.
The robust effects of Coffea cruda on MSE, in addition to mood13 and conventional polysomnographic measures49 in this study merit additional research. Many young adult college students have repeated, extensive use of beverage coffee (Coffea as a “tincture”) prior to their participation in the study. In contrast, the single bedtime dose of Nux vomica would likely have been a first-ever exposure to the remedy in any form for the study participants. In the phenomenon of time-dependent sensitization,86 past exposure history can lead to significant amplification of behavioral and physiological responses in susceptible individuals to various agents given intermittently but repeatedly over time (these subjects had a pre-existing history of beverage coffee-induced insomnia).87–89 Future sensitizers and non-sensitizer controls do not necessarily show differences from one another in response patterns to an initial, single exposure to a given agent on the first dose, or at rest, prior to a subsequent exposure.89
Also, for this initial complexity analysis, we controlled for personality and did not attempt to look for NDS sleep EEG effects within subsets of participants whose personality traits might have enhanced or dampened the potential to respond to a given remedy. Future studies will need to look at the adequacy of matching the remedy to the individual, with the prediction that a good match would generate not only net clinical benefit, but also corresponding changes in nonlinear EEG dynamics5 that might help predict or distinguish responders from non-responders to treatment.43,71 Data from this laboratory suggest that individual difference traits and remedy can interact in generating differential mood13 or waking EEG effects.42
Various investigators5,14 previously suggested that the state-dependency of homeopathic remedy effects could relate to effects on the loss of connectivity across the person as a complex network in the diseased state and enhanced system susceptibility to the effects of the signal properties of the remedy. In a recent evidence-based review of biological model systems and homeopathy, Baumgartner76 concluded that (a) remedies have stronger effects on more complex systems such as whole organisms than on simple cell culture model systems; (b) the state dependency of remedy effects involves ceiling/floor limits in the initial condition of individual organisms (e.g., smaller plants at baseline grow better than larger plants with homeopathic treatment)(see also87); (c) homeopathic treatment may produce reduction in the inter-experimental variability of data distribution more than consistent effects on group means.
It is not clear what are the clinical implications, if any, of the present data. In our subjects, the medicines appeared initially to increase total sleep time, but disrupt the quality of the sleep by conventional scoring criteria.49 The self-report data from this study suggest that the medicines also affected mood as a function of the individual personality traits of the participants,13 and that the quality of the dreams recorded in morning diaries may also have changed (unpublished data). Given the history of inter-experimental variability in results for homeopathic basic science and preclinical research findings,76 however, it would be essential to replicate and extend the present study before drawing more speculative conclusions about the type or direction of remedy effects clinically.
Homeopathic remedy effects on complex living systems may be regulatory or adaptogenic overall, rather than unidirectional, in restoring homeostasis.3,6,8,23,24,76,90–92 Hahnemann's Organon (Aphorism 65) discusses the initial action of a stressor, stimulus, drug, or remedy as inducing changes in one direction, followed by the organism's compensatory counter-action or after-action “for the restoration of the normal state.”93 Classical homeopaths refer to a loss of vitality in a less resilient patient as a reduction in “vital force.”94 Using very different terminology to characterize clinical status, contemporary NDS researchers such as Goldberger et al conceptualize most diseases and aging as conditions that reduce overall systemic complexity.52,95,96 Other investigators have postulated that there is an optimum range of dynamical complexity in a living system that allows human well-being, resilience,97 and flourishing to emerge.22,23,98–100 However, EEG complexity assessed by multiscale entropy analysis can change over shorter time scales as well.101 Overall, then, transiently reduced biological complexity might indicate a comparatively adverse effect in the direction of disease, whereas increased complexity might suggest a comparatively beneficial effect in the direction of recovery or healing.22,23
Homeopaths commonly report a biphasic remedy effect phenomenon of initial, short-term clinical exacerbation of symptoms (sometimes termed an “aggravation” or “healing crisis”), prior to gradual resolution of symptoms.7,102 Although our current subjects had histories of coffee-induced insomnia, they were generally healthy. The homeopathic prediction for healthy individuals would be a worsening rather than improvement of sleep quality.7 Exploring these issues from an NDS perspective could contribute a great deal to resolving confusion about the bidirectionality often observed from one homeopathic study to the next.76
Finding ways to document the initial dynamical state of the individual homeopathic subjects' physiology immediately prior to remedy administration may help improve predictions about the nature of any subsequent changes observed over both short- and long-term follow-up. The evolution of change in the dynamics of a complex nonlinear system is highly sensitive to initial starting conditions.25,103 Such research may parallel the issues raised by clinical observations of nonlinear, bidirectional response patterns to isopathic asthma treatment with dust mite 30c5,104 and of sudden worsenings or improvements early in psychotherapy for depression.105,106
As complex organisms further from equilibrium (`sicker' in some sense) may respond more to individually-tailored homeopathic treatment,76 it would also be important to assess the current medicines on patients with primary insomnia evaluated by clinical interview to need one of these specific medicines by an experienced homeopathic provider, rather than the relatively healthy young adults screened by questionnaire and enrolled in this study. At the same time, this study design reflected a modification of a phase 1-like drug development trial (homeopathic pathogenetic trial) in homeopathy. Relatively healthy volunteers usually take the test substance and record subjective effects over a period of days to weeks,12,72,107 with the expectation of a short-term transient increase in some remedy/person-characteristic symptoms.
In conclusion, the homeopathic remedies Coffea cruda and Nux vomica in 30c potencies alter short-term nonlinear dynamic parameters of slow wave sleep EEG in relatively healthy young adults. As an NDS analytic technique, MSE50,52,74 is apparently more sensitive than Mekler-method D253 for detecting these effects. Additional research to explore individual differences in state-dependency, direction and longer term patterns of homeopathic remedy effects on objectively-recorded human sleep in people with and without insomnia is indicated.
This study was supported by National Center for Complementary and Alternative Medicine Grants R21 AT000388, K24 AT000057, and T32 AT01287 (PI: IRB). The authors would like to thank Keith Fridel, PhD RPSGT for training and supervising the sleep scoring, and Elizabeth Acker, Michael Biuso, Erica Morey, Mallory Taylor, and Alivia Wieseler for their technical assistance in data collection and quality assurance procedures.
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
Conflicts of Interest Drs. Bell and Brooks serve as consultants to Standard Homeopathic Co./Hyland's Inc., a homeopathic manufacturer whose products were not used in the current study. Standard Homeopathic Co./Hyland's Inc. did not provide any funding for the current study.