Published on in Vol 6, No 2 (2021): Apr-Jun

Preprints (earlier versions) of this paper are available at, first published .
Virtual Reality–Guided Meditation for Chronic Pain in Patients With Cancer: Exploratory Analysis of Electroencephalograph Activity

Virtual Reality–Guided Meditation for Chronic Pain in Patients With Cancer: Exploratory Analysis of Electroencephalograph Activity

Virtual Reality–Guided Meditation for Chronic Pain in Patients With Cancer: Exploratory Analysis of Electroencephalograph Activity

Original Paper

1School of Engineering Science, Simon Fraser University, Burnaby, BC, Canada

2School of Nursing, University of British Columbia, Vancouver, BC, Canada

3School of Rehabilitation Sciences, University of British Columbia, Vancouver, BC, Canada

4School of Interactive Arts and Technology, Simon Fraser University, Surrey, BC, Canada

5Surrey Memorial Hospital, Fraser Health Authority, Surrey, BC, Canada

6Behavioural and Cognitive Neuroscience Institute, Simon Fraser University, Burnaby, BC, Canada

Corresponding Author:

Bernie Garrett, DPhil

School of Nursing

University of British Columbia

T201-2211 Wesbrook Mall

Vancouver, BC, V6T 2B5


Phone: 1 604 822 7443


Background: Mindfulness-based stress reduction has demonstrated some efficacy for chronic pain management. More recently, virtual reality (VR)–guided meditation has been used to assist mindfulness-based stress reduction. Although studies have also found electroencephalograph (EEG) changes in the brain during mindfulness meditation practices, such changes have not been demonstrated during VR-guided meditation.

Objective: This exploratory study is designed to explore the potential for recording and analyzing EEG during VR experiences in terms of the power of EEG waveforms, topographic mapping, and coherence. We examine how these measures changed during a VR-guided meditation experience in participants with cancer-related chronic pain.

Methods: A total of 10 adult patients with chronic cancer pain underwent a VR-guided meditation experience while EEG signals were recorded during the session using a BioSemi ActiveTwo system (64 channels, standard 10-20 configuration). The EEG recording session consisted of an 8-minute resting condition (pre), a 30-minute sequence of 3 VR-guided meditation conditions (med), and a final rest condition (post). Power spectral density (PSD) was compared between each condition using a cluster-based permutation test and across conditions using multivariate analysis of variance. A topographic analysis, including coherence exploration, was performed. In addition, an exploratory repeated measures correlation was used to examine possible associations between pain scores and EEG signal power.

Results: The predominant pattern was for increased β and γ bandwidth power in the meditation condition (P<.025), compared with both the baseline and postexperience conditions. Increased power in the  bandwidth was evident, although not statistically significant. The pre versus post comparison also showed changes in the θ and α bands (P=.02) located around the frontal, central, and parietal cortices. Across conditions, multivariate analysis of variance tests identified 4 clusters with significant (P<.05) PSD differences in the δ, θ, β, and γ bands located around the frontal, central, and parietal cortices. Topographically, 5 peak channels were identified: AF7, FP2, FC1, CP5, and P5, and verified the changes in power in the different brain regions. Coherence changes were observed primarily between the frontal, parietal, and occipital regions in the θ, α, and γ bands (P<.0025). No significant associations were observed between pain scores and EEG PSD.

Conclusions: This study demonstrates the feasibility of EEG recording in exploring neurophysiological changes in brain activity during VR-guided meditation and its effect on pain reduction. These findings suggest that distinct altered neurophysiological brain signals are detectable during VR-guided meditation. However, these changes were not necessarily associated with pain. These exploratory findings may guide further studies to investigate the highlighted regions and EEG bands with respect to VR-guided meditation.

Trial Registration: NCT00102401;

JMIR Biomed Eng 2021;6(2):e26332




Chronic pain (CP) is a common condition occurring in 1 in 5 Canadians [1], and it has limited effective treatment approaches. Mindfulness-based stress reduction (MBSR) has shown some evidence of efficacy in this area [2,3] and has also been used to treat other clinical conditions, such as migraine, depression, addiction, and substance misuse [2,4].

Mindfulness meditation encompasses a range of mental exercises that share a common focus on regulating attention and awareness to improve well-being [2-4]. It is described as the quality of being completely engaged in the present moment, free from distractions and judgments, and being aware of bodily sensations, thoughts, and feelings without getting caught up by them, and it is used as a therapeutic technique in MBSR [5,6]. These practices involve mental training that allows practitioners to develop their minds in specific ways to help them deal with stress and anxiety [7,8]. Although the clinical benefits remain somewhat controversial, it is generally viewed as a beneficial practice for mental well-being, stress reduction, and pain management [9-11].

A recent trend in MBSR practice has been the use of immersive virtual reality (VR) to help participants focus on meditative exercise [12-16]. However, to date, neurological studies have not been performed with VR-guided meditation practices. Therefore, this exploratory study sought to identify any identifiable neurological effects of VR-guided meditation practices using electroencephalographs (EEGs).

Neurological Mindfulness Studies

Subjective reports of the benefits of mindfulness meditation have prompted investigations into the potential corresponding neurophysiological states. Exploration of fluctuations in brain wave voltage amplitude (power) topography and coherence (associated areas of activity) using EEG variations in neural activity assessed with functional magnetic resonance imaging and cortical evoked responses to visual and auditory stimuli that reflect the impact of meditation on attention [4,17-26]. However, findings remain speculative. EEG studies have previously reported modulation in α, θ, and γ band frequencies, generally with increased power and coherence during mediation [4,7,17,27]. Some studies theoretically conjectured how meditative states relate to EEG bandwidths. For example, Travis and Shear [28] suggested that focused attention (sustained attention is focused on a given object) increases γ band power, open monitoring (nonreactive monitoring of an ongoing experience) increases θ, and an automatic-self-transcending stage (transcending the procedures of the meditation) increases α [28]. Nevertheless, a consensus on what we currently know about how EEG forms correlate with meditation or how this may map onto stages of mental development or specific cognitive skills is yet to be reached [7,8].

The α frequency band lies between 8 and 12 Hz and is predominantly located in the occipital cortex. α waves are present in deep relaxation and sleep, usually when the eyes are closed. θ waves are characterized by oscillations in the 4-8 Hz band found in both cortical and subcortical structures. Increases in θ have been described during a variety of learning and recognition tasks, light sleep (including the rapid-eye-movement dream state), and deep meditation. θ and α power changes have been reported to increase in a number of meditation studies [4,17,19,29-31]. Moreover, γ is a higher frequency range generally regarded as between 30 and 50 Hz, although the range reported has varied substantially between 20 and 200 Hz across different studies [30]. This initial research suggests that γ is associated with a number of sensory and cognitive high-level information processing functions, such as semantic insights, learning, and neural plasticity. Peak γ frequencies around 40 Hz in bilateral hemispheres have only been observed in highly practiced meditators [7,30]. In addition, a posterior increase in γ activity may be related to enhanced perceptual clarity reported in some open monitoring (focusing on awareness itself) meditative processes [7,17,30].

VR-Assisted Mindfulness

VR is rapidly developing as a new form of media and uses computer-generated audio-visual displays and hand controller user interfaces to produce a sense of immersion in a digital 3D environment. Instead of watching an image on a typical computer or video display, VR technologies provide an increased sense of presence by engaging the senses (sight, sound, and touch) in real-time stereoscopic audio-visual media where users can move around and explore the environment as if they were there.

VR in Pain Management

One of the most common health care applications of VR is its use in pain management. Several studies have explored the value of MBSR in CP management, although the reported effect sizes for this technique have been typically mild to moderate, and regular adherence with meditative practice is problematic [3,32,33]. It has been suggested that combining mindfulness meditation within a VR intervention may help support acceptance and adherence to practice while having a synergistic effect on pain reduction through immersive VR distraction [34-36]. This remains an active area of research, as adjunctive VR strategies have been used successfully in the treatment of acute pain [37-41] and more recently have also been explored with CP [12,35,42,43]. The theoretical rationale behind why VR may enhance mindfulness skills is that VR provides the user with cognitive displacement by actively engaging in a coping activity that provides a profound sense of through presence in another world. Cognitive distraction is a common strategy in pain control and relies on creating competition for cognitive resources, that is, attention to a novel spatial orientation, and engaging within it reduces the perception of pain [32,33,44]. Therefore, immersive VR interventions using stereoscopic head-mounted displays (HMDs) have been proposed as powerful tools for providing visual, audio, cognitive, and emotional engagement [45-47]. In MBSR, VR experiences are typically accomplished using computer-simulated environments, stereoscopic headsets, and motion tracking to support a more immersive meditative experience. This was the approach taken in this study, and an EEG analysis was selected as a practical technique to assess neurophysiological activity while experiencing VR.


Overall, this inductive exploratory study sought to assess how EEG power of waveforms, their topographic mapping, and coherence measures altered in 3 main states during a VR-guided meditation experience in patients with cancer-related CP: at baseline (pre), during VR-guided meditation (med), and after VR-guided meditation (post). Moreover, we explored whether their pain level was associated with waveform power measures. We were particularly interested in the power, topography, and coherence of α, γ, and θ wave activity, and possible synchrony with VR activity, as other researchers have reported changes in these 3 waveforms with MBSR activities. Therefore, we were interested in determining whether prior findings were consistent with those observed during a VR experience. Specifically, the questions we sought to address in this study were as follows:

  1. Were there any observable or significant changes between pre and during VR-guided meditation experiences and between pre and post VR-guided meditation experiences with the power, topographic changes, or EEG coherence in specific waveforms?
  2. Was there any evidence that changes in the pain experienced by participants during a VR-guided meditation activity were associated with observable EEG changes?


An exploratory, single-subject design study was undertaken to compare EEG activity and pain levels before, during, and after VR-based meditation practice. This study was reviewed and approved by the University of British Columbia Clinical Research Ethics Board.


A convenience sample of 10 participants was used and recruited from those in an existing randomized controlled trial (RCT; NCT 02995434) where patients with cancer were using VR as an adjunctive therapy to help manage their CP (Textbox 1). These participants were completing or had completed cancer treatment and experienced a range of cancer-related pain, including neuropathy, fibromyalgia, postsurgical pain, or an exacerbation of pre-existing pain.

Randomized controlled trial eligibility criteria.

Eligibility Criteria

  • Aged ≥16 years, with a past or current diagnosis of cancer
  • Prior or ongoing cancer treatment
  • A chronic pain diagnosis (ongoing daily pain for ≥3 months, with a Neuropathic Pain Rating Scale of ≥4)
  • Able to understand English (read and write)
  • Normal stereoscopic vision
  • Able to move their head up, down, left, and right and able to wear a virtual reality head-mounted display
  • Sufficient fine motor control in one hand to use a game controller
  • Have space at home for a computer and monitor
Textbox 1. Randomized controlled trial eligibility criteria.

Participants were purposively recruited, focusing on recruiting those from the RCT cohort who had previously responded well to a VR-based meditation experience, with a self-reported reduction of a Visual Analog Scale for Pain ≥1. They were invited to participate in a single 2-hour VR-guided meditation experience, with EEG recorded in their home or at the university, and were offered a Can $100 (US $83) honorarium and expenses for taking part. As an exploratory study designed to establish methods and feasibility using a limited convenience sample, with an unknown effect size, the power to inform the sample size was not calculated a priori. This is acceptable for this type of inductive study [48].


EEG signals were recorded during the session using a BioSemi ActiveTwo system (BioSemi) with 64 channels in a standard 10-20 configuration. This system uses a head cap system with pin active silver chloride electrodes. The EEG ground (labeled DRL [Driven Right Leg]) on the Biosemi system was placed between POz and PO4, and the EEG reference (labeled CMS [Command Mode Sense]) was placed between POz and PO3. The antialiasing filter was a fixed first-order analog filter with −3 dB at 3.6 kHz, and the low-pass filter was a fifth-order cascaded integrator-comb digital filter with −3 dB at one-fifth the sample rate. A powerline notch filter was not applied because the system used active electrodes, a battery power supply, and optic fiber, greatly reducing noise from the powerline. Each channel consisted of a 24-bit analog-to-digital converter. Recordings were made at 1024 Hz (although one was recorded at 2048 Hz) using ActiView software (Biosemi Instrumentation).

An HTC Vive VR system (HTC) with a Deluxe Audio Strap fitted over the top of the EEG cap was used. This system features a 2160×1200 resolution, 90 Hz refresh rate, and 110° field of view. Positional tracking from 2 infrared cameras enabled 5-degrees-of-freedom motion tracking of the headset and hand controllers. Integrated stereo headphones supported 3D audio immersions. During the initial pilot testing, this system was found to work effectively. Minor noise was evident in the EEG from the VR system on rare occasions, but we were able to remove most of this noise by careful repositioning of the equipment.

For meditation practice, a commercially available Guided Meditation VR application (Cubicle Ninjas) was used. As each meditation was only 10 minutes in this app, we selected a single sequence of 3 unique meditations, Zen 2-4, to form a 30-minute block of meditation. In the guided meditation experience, users were situated in the Lost Woods virtual environment with the Calm music selected. They were able to explore a calm, forest 3D environment with running water features with soft chirping bird and gentle wind sounds. Using a controller, participants could move to different positions in the forest to explore or find a particular viewpoint they liked and found most conducive to their meditation. A narrative provided audio guidance on the meditative practice. This environment was selected to maximize the similarity with the participants’ prior VR experiences in the RCT.

Both EEG recording and the VR system were run on a Dell G7 17 7790 gaming laptop (Intel Core i7-8750H, 16 GB DDR4, RTX 2060) placed in front of the participant. During recording, the laptop screen was arranged to display the ActiView software and a mirror window of the Guided Meditation VR application experience, showing the participant’s perspective in the VR HMD. In addition, a Windows 10 (Microsoft) camera app was used to record the video. This arrangement was recorded using Snagit (Techsmith) screen capture software to accommodate the time synchronization of EEG and recordings of participants’ physical movements (Figure 1).

Figure 1. Electroencephalograph recording during a virtual reality experience.
View this figure


Participants were seated in front of the laptop, and the laptop webcam was framed to capture the participants’ head, arms, and upper body. Before putting on the EEG head cap, participants were familiarized with the VR app, its controls, and how to navigate and select meditations. The Cz, inion, and left and right preauricular locations were marked using standard EEG landmarking methods. These locations were used to align the EEG head cap on the participant. The electrode paste was then applied, and the electrode contacts were adjusted until the electrode offsets were less than 50 µV. The VR HMD was then placed in position over the EEG head cap, and the experience commenced.

The overall structure of data collection used an 8-minute resting condition, followed by the 30-minute sequence of meditations and followed by another 8-minute rest condition (Figure 2). The 8-minute period was considered a balance between being shorter to reduce strain and being longer for more data; this period has been used in EEG studies in a variety of fields for baselines [49-51].

Figure 2. (A) Diagram of equipment setup. (B) Timeline of recording, including rest and meditation conditions. (C) Diagram of desktop view. EEG: electroencephalograph; NRS: numerical rating scale.
View this figure

During the rest conditions, the participants were asked to rest quietly and observe a small white crosshair displayed on a black background in the VR HMD. Participants were instructed to keep their eyes open while blinking naturally, to keep their eyes on the crosshair, and to stay still while not thinking about anything. During the meditation practice, participants were instructed to engage with the guided meditation and look around or move around the virtual environment as they liked, to find an enjoyable perspective. Each guided meditation experience lasted 10 minutes, and participants were instructed to begin the next one immediately after the intervening rest condition. Pain was assessed before and after the first rest condition, after each 10-minute guided meditation, and after the second rest condition. Participants were asked to verbally rate their pain from 0 to 10 on the simple numerical rating scale (NRS) [52].


Data Preprocessing

The recording sessions of the resting condition, the 3 guided meditation conditions, and the final resting condition were referred to as the Pre, Med1, Med2, Med3, and Post, respectively, for the analysis. These 5 conditions were extracted based on the time stamps acquired from the EEG video recordings.

Figure 3 illustrates the EEG data preprocessing steps. The raw EEG data in the Biosemi format were imported and preprocessed using FieldTrip software (Donders Institute for Brain, Cognition and Behaviour) [53]. As the bandwidth of interest was less than 50 Hz, the data were downsampled to 512 Hz to reduce computing time in later processes. To prepare for downsampling, an antialiasing sixth-order low-pass Butterworth two-pass filter set to 70 Hz was first applied to the data. The primary 64 EEG channels were then re-referenced to the averaged left and right mastoid electrodes (T7 and T8). This re-referencing process effectively eliminated the noise introduced from the original reference electrodes. After this process, the reference electrodes T7 and T8 were eliminated from the data, and only 62 channels were retained.

Figure 3. Electroencephalograph data preprocessing chart. EEG: electroencephalograph; ICA: independent component analysis.
View this figure

In addition, a sixth-order high-pass Butterworth two-pass filter with a cut-off at 1 Hz was used to remove slow drifts and to prepare the data for independent component analysis (ICA) [54]. The data were further notch-filtered to remove any powerline noise and harmonics, with cut-offs set at 60, 120, and 180 Hz. Bad channels were identified during the EEG recording process using an EEG data browser and rejected.

These data for each condition were further segmented into nonoverlapping epochs of 2 seconds in length to enable signal averaging in the frequency domains [55]. A data cleaning process was then performed using the FieldTrip automatic artifact removal feature, which is based on a Z-transformation and the setting of a threshold to reject bad epochs. This artifact removal result was found to be unsatisfactory, so a manual cleaning process was then performed. Segments with participant movement observed in the videos were identified, and any residual bad segments and channels were rejected. These steps effectively eliminated the artifacts caused by head and body movements.

Any remaining artifacts caused by eye blinks, eye movements, and external noise were eliminated from the data using ICA. These manual artifact removal and ICA processes were repeated 3 times during the cleaning process to ensure that all movement-based artifacts had been captured. Finally, bad channels were repaired [56] using the established practice of spline interpolation of the neighborhood channels on the bad channels [57].

Power Spectral Density Analysis

The power spectral density (PSD) [58] was computed for the conditions using a fast Fourier transform method [59]. Hanning taper smoothing was applied to reduce spectral leakage owing to the discontinuity of the signal at the start and end points of the epochs. To improve the signal-to-noise ratio, the average PSD of the epochs was used. The PSD was further normalized by dividing the average power across all the frequency bins in each channel. The frequency range in the analysis was set to 2-50 Hz, which effectively covers the waveforms commonly used in EEG analysis. The PSD results of the premeditation experiences, during meditation experiences, and postmeditation experiences were then plotted and shown as pre, med, and post conditions, respectively. For med, the average PSD of the meditation conditions was computed using the average of the 3 meditation conditions (Med1, Med2, and Med3). For a graphical overview, a band power box and whisker plot were created. The plot revealed the signal power over the typical EEG frequency bands: 2-4 Hz (δ), 4-8 Hz (θ), 8-12 Hz (α), 12-30 Hz (β), and 30-50 Hz (γ) [60]. The lower boundary of 2 Hz was set to avoid any distortion introduced by the high-pass filter at 1 Hz. The upper boundary of 50 Hz was selected to avoid the heavy power-level drop due to the use of a 60 Hz notch filter. In addition, 50 Hz was considered a reasonable cut-off between the high and low γ ranges.

Topographic and Coherence Analysis

To explore the spatial properties of the signals on each EEG channel, a topographic mapping of the PSD was plotted according to the traditional frequency ranges. The PSD was plotted using a heat map visualization technique to display the magnitude of the PSD with a 2D color representation. An interpolation based on the MATLAB 4 grid data method was used to smooth the topography, where a surface was fitted to the scattered data points using a biharmonic spline interpolation [61]. This PSD plot was used to graphically illustrate the spatial properties associated with the pre, med, and post experiences.

A coherence analysis [62] was performed to visualize the functional connectivity between the electrodes as an indication of the brain areas that may be functionally integrated. Coherence of the pre, med, and post conditions were plotted for the frequency bands. A value between 0 and 1 was displayed for each channel pair. A value of 1 indicates full synchronization between the channel pair, and a value of 0 reveals that the channel pair does not work in the synchronized condition. The significance of the coherence difference between the conditions was computed using the two-tailed Student t test and then plotted.


Pain experience was measured using the NRS. Pain scores were collected from the patients as an initial baseline at the start, and after each condition, to explore any correlation between pain levels and any PSD changes identified during the meditative experience.

Statistical Analysis

For the PSD and coherence analysis, the null hypothesis set for inferential statistical analysis was that the probability distribution of the condition-specific averages for PSD and coherence would be identical for all conditions.

As a cluster-based permutation test [63] can be used to effectively resolve multiple comparisons in EEG signal statistical analysis, this test was used in this study to examine the overall PSD differences between the pre-med, med-post, and pre-post conditions by identifying the clusters of electrodes with significant changes. The tests were conducted using the Monte Carlo method, with 128 permutations, two tails, and α=.025 (negative and positive tails together equal .05). In addition, for comparison of all 3 conditions together, a multivariate analysis of variance (MANOVA) was used, with 2048 permutations, one tail, and α level set to .05.

The effect size was calculated after the cluster-based permutation test. First, a bounded rectangular area spanning each cluster was identified. This rectangular area was bounded by the frequency window of the cluster and all the channels in that cluster. The PSD in this area was then averaged. The maximum Cohen d effect size was then computed using FieldTrip for each of the conditions. The effect size was not computed for the MANOVA because it requires the specification of 2 mean groups for comparison. Generally, a Cohen d effect size around 0.2 is considered small and around 0.8 and higher is considered large to very large [64,65]. Finally, peak channel tests were performed to verify the PSD changes, and post hoc power analyses using G*Power [66] were performed to indicate statistical power based on the effect sizes observed in the sample.

For the coherence analysis, the Student t test was used to examine significant coherence differences between the conditions. Parameters were set to two-tailed, paired samples, and α=.05, with a 10% false discovery rate (FDR) adjustment set [67].

To explore the associations between pain scores and PSD, repeated measures correlation [68] was conducted. Some peak PSD channels identified with a high significance of PSD change from the cluster topoplot were input for each condition as the repeated independent variable, and the NRS pain score was input as the dependent variable.


A total of 10 participants were recruited (Table 1). As we went through the data cleaning process, the data from 1 participant (S05) were found to be excessively noisy, and on video review, 2 other participants (S08 and S09) were found not to be following the procedures. Hence, their results were excluded, and a sample size of 7 was used in the final analysis.

Table 1. Participant demographics.
SubjectAge (years)SexCancer history
S0159FemaleAbdominal tumors
S0266MaleNon-Hodgkin lymphoma
S0337MaleNon-Hodgkin lymphoma
S0458FemaleBreast cancer
S0750FemaleColon cancer
S1064MaleProstate cancer

PSD Analysis

Figure 4 illustrates the average of the power spectrum of the 3 conditions to provide a general overview of the power spectrum in the data, whereas Figure 5 shows the box and whisker plot of band power of the pre, med, and post conditions in the 5 traditional frequency ranges. From these figures, the differences in power in the frequency bands between the conditions were identified. The post condition power level had the lowest median and mean values in the δ range. The med condition changed from the lowest median and mean values in θ and α to the highest median and mean values in β and γ.

Figure 4. Grand average power spectrum.
View this figure
Figure 5. Power level changes between the pre, med, and post conditions in different frequency ranges. Box plot shows the median and range of power level of the participants. Line plot shows the changes in mean power.
View this figure

The topographic distribution of the PSD is shown in Figure 6. Differences in pre, med, and post conditions are shown spatially in all bands and in different brain regions. Changes in the power levels during meditation were observed in all frequency bands. In δ, an increase in power level in the central occipital region was observed, and a drop of power in the frontal-central region was observed in the post condition. In θ, there was a drop in the power level in the frontal cortex. In α, a decrease in power level in the central parietal region was noted. In the β band, an increase in power level was found in the bilateral central and prefrontal regions during meditation. In γ, an increase in power level was noted in the left frontal (LF) and right frontal (RF) regions, and a slight increase in the central parietal region was observed.

Figure 6. Topography of power spectrum shown in different frequency ranges. Different color scales are used for the frequency ranges to reveal the details in the central areas.
View this figure

Figure 7 shows the topoplot of the two-tailed cluster-based permutation test results. The conditions were compared in the test, including two-tailed pre versus med, med versus post, and pre versus post and a one-tailed MANOVA for all 3 conditions together. The clusters of electrodes with significant differences found in PSD changes were marked with a circle marker at the electrodes. The color bars indicate the t values computed from the test. For the pre versus med comparison, a cluster with a significance of P=.02 was found in the range of 24.5-31 Hz (high β and low γ) in the frontal cortex. For the med versus post comparison, significance was found in the β and γ ranges. A cluster with a significance of P=.001 was found in the range of 37-50 Hz (γ), and the cluster covered most brain regions; a cluster with a P value of .008 was found in the range of 23-36 Hz (high β and low γ) in the frontal and central cortices. For the pre versus post condition comparison, a cluster with P=.02 was obtained in the frontal, central, and parietal regions in the range of 8-9.5 Hz (high θ and low α). The MANOVA cluster-based permutation test for all the pre, med, and post conditions returned 4 clusters with P≤.05, where the test was one-tailed. The first cluster found has P=.002 in the range of 37.5-50 Hz (γ) in the frontal, central, and parietal cortices. The second cluster with P=.03 was found in the range of 31-36 Hz (low γ) in the frontal and central cortices. The third cluster with P=.03 was found in the range of 2-5 Hz (δ and low θ) in the frontal, central, and parietal cortices. The last cluster with P=.04 was found in the range of 24-30 Hz (high β) in the frontal and central cortices. Table 2 shows the test results for all the clusters of significance. The very large effect size indicated that the permutation tests were effective in rejecting the null hypothesis.

Figure 7. Power spectrum analysis using the cluster-based permutation test. Clusters of electrodes found with significant changes in power are marked with a circle marker at the electrodes. The color bars show the permutation test t-value level. MANOVA: multivariate analysis of variance.
View this figure
Table 2. Cluster-based permutation test results. All clusters with P≤.025 are shown for the first 3 tests, and clusters with P≤.05 are shown for the multivariate analysis of variance test.
Comparison and clusterP valueEffect size









aMANOVA: multivariate analysis of variance.

bNot calculated as it required to specify which 2 conditions to compare.

By inspecting the topoplots (Figures 6 and 7), 5 peak channels with noticeable changes in the power levels between the conditions were selected for additional analysis and for better understanding the observed power level changes. The peak channels selected were AF7 and Fp2 in the prefrontal region, FC1 in the frontal region, CP5 in the left central (LC) region, and P5 in the left parietal (LP) region. Two-tailed, paired-sample t tests were performed to examine the overall PSD changes in different conditions. These tests were conducted on the pre and med, med and post, and pre and post conditions. The resulting P values were reported with an FDR adjusted to 10%. The t test results were plotted graphically for ease of interpretation. For each of the peak channels, a combined plot of the results of the 3 conditions was used.

Figure 8 shows the P values obtained in a graphical form for the comparison of the 3 conditions. The upper shaded areas represent P≤.025. The shaded areas above the dashed line represent the adjusted significance with P≤.0025.

Figure 8. Graphical depiction of P values of changes in power spectral density between conditions for the selected channels. Shading indicates significance found at the .025 level, and shading above the dashed line indicates the adjusted significance at the .0025 level (with the false discovery rate set to 10%).
View this figure

A statistical power analysis was also performed to verify the probability of detecting a true effect in the t tests. Table 3 shows the maximum effect size value and the respective frequency point found for each channel and for each comparison. Table 4 shows the statistical power based on the effect sizes listed in Table 3. It is shown that the t test results have a high power ≥0.7302 even when using the adjusted α value of .005.

Table 3. The maximum channel effect size and the frequency at which it was located.
Channel (N=7)Pre-medMed-postPre-post

Max effect sizeFrequency (Hz)Max effect sizeFrequency (Hz)Max effect sizeFrequency (Hz)

aItalicized values indicate where P≤.0025 was found for that frequency value.

Table 4. Power analysis for the effect size found.
ChannelPre-med, α=.005Med-post, α=.005Pre-post, α=.005

aItalicized values indicate where P≤.0025 was achieved.

Topographic and Coherence Analysis

Figure 9 demonstrates the coherence detected in the α band during meditation.

The functional connectivity level is represented by a red to blue color scale. The red color with a value close to 1 indicates that a channel pair is highly synchronized in the signal transfer. The blue color with a value close to 0 indicates that the channel pair is working independently.

Figure 10 shows the significant difference in coherence for the comparison of the med and post conditions in the α band. Channel pairs with P≤.025 are highlighted in green. Those with P≤.0025 (FDR adjusted) are highlighted in yellow. For clarity, the EEG channels were divided into regions according to their respective locations in the brain, namely, LF, LC, LP, left occipital, central parietal and occipital, RF, right central, right parietal and right occipital. Table 5 summarizes the results of a series of plots, as shown in Figure 10, for all the bands and all 3 comparisons. The table shows the significant coherence changes with P≤.0025 between the brain regions, where 1 denotes the pre and med comparison, 2 denotes the med and post comparison, and 3 denotes the pre and post comparison. The frequency band of the region pair at which a significance was found is shown with the Greek letter of the band. The number inside parentheses indicates the number of significant channel pairs in the frequency band.

Figure 9. Coherence between the electroencephalograph channels during the virtual reality–guided meditation and in the α frequency range.
View this figure
Figure 10. Channel pairs with significant coherence difference between the meditation and post meditation conditions in the α band. t test was used in the analysis, and channel pairs with a P≤.025 are highlighted in green, and channel pairs with a P≤.0025 (with the false discovery rate set to 10%) are highlighted in yellow. CPO: central parietal and occipital; LC: left central; LF: left frontal; LO: left occipital; LP: left parietal; RC: right central; RF: right frontal; RO: right occipital; RP: right parietal.
View this figure
Table 5. False discovery rate–adjusted significant coherence changes found in the region pairs for the 3 comparisons: 1: pre versus med, 2: med versus post, and 3: pre versus post.
RO3: δj3: θkl
RP2: δ(2)m2: δ2: δ2: α; 3: β1: θ
RC1: β3: α (✓)
RF1: γ(2); 2: θ(2), α(5); 3: δ(2)2: α
CPO3: δ α1: α; 2: α(2)
LP2: θ, α(2), β

aLF: left frontal.

bLC: left central.

cLP: left parietal.

dLO: left occipital.

eCPO: central parietal and occipital.

fRF: right frontal.

gRC: right central.

hRP: right parietal.

iRO: right occipital.

jGreek letter indicates the frequency band of the region pair with significant change.

kThe empty cells in the upper left of the check mark diagonal show no significant channel pair was found. The cells in the lower right of the check mark diagonal are not used, as they are just mirrored duplicates of the cells in the upper left of the check mark diagonal.

lThe checkmark indicates a region connects to the same region in the region pair.

mNumber within parentheses indicates the number of channel pairs with significant changes in that frequency band. If there is only 1 channel pair, the number is not shown.

As shown in Table 5, significant changes were mostly found between the frontal and parietal regions, namely, RF-LP (11 channel pairs) and LF-LP (4 channel pairs). In addition, the most active regions were LP (related to 17 channel pairs), RF (15 channel pairs), LF (7 channel pairs), right parietal (7 channel pairs), and central parietal and occipital (7 channel pairs). For the 3 comparisons overall, significant changes were mostly found between the frontal and parietal and occipital regions, particularly in the θ, α, and β bands in the med-post comparison. For the pre-med comparison, the frontal-parietal region-pair coherence changes were observed in the γ band. For the pre-post comparison, the frontal-parietal and occipital region-pair coherence changes were observed in the δ, α, and β bands. Table 6 shows the effect sizes of the 2 regions of interest, namely, RF-LP and LF-LP. Some large to very large effect sizes (≥0.8) were found in the δ, θ, α, and γ bands.

Table 6. Effect sizes of the 2 region pairs of interest: right frontal-left parietal and left frontal-left parietal.
BandRegion-pair RF-LPa average, effect sizeRegion-pair LF-LPb average, effect size


aRF-LP: right frontal-left parietal.

bLF-LP: left frontal-left parietal.

cItalicized values indicate large and very large effect sizes.

Table 7 shows the NRS pain scores collected immediately after the pre, Med1, Med2, Med3, and post conditions. A repeated measures correlation analysis was performed on some selected peak channels and clusters. In the test, no significant relationship between the collected pain scores and PSD was found.

Table 7. Numerical rating scale scores collected after each condition.
ParticipantNRSa score


S01422 21

aNRS: numerical rating scale.

Principal Findings


In answer to the primary research question, there were significant changes in EEG power and coherence among three conditions (pre, VR-guided meditation, and post). Therefore, the null hypothesis of no difference was rejected.

The visual inspection of the global normalized power spectrum analyses revealed various changes in all bandwidths. The predominant pattern was for increased δ, β, and γ bandwidth power in the meditation condition, compared with both the pre and post conditions. In the θ and α bandwidths, the changes in power were more varied within the 3 conditions.

Pre Medication Versus VR-Guided Meditation

Visual inspection of the topographic distribution showed 2 main patterns comparing meditation with the prior resting condition. The first was an increase in power of δ (mainly in the central and occipital areas), β (mainly in the bilateral prefrontal areas), and γ (mainly in the frontal and bilateral prefrontal areas) during VR-guided meditation.

The second pattern that emerged from the visual topography map was decreased low-frequency range power of the θ (mainly in frontal areas) and α (mainly in occipital and parietal areas) bandwidths in the med condition compared with the pre condition. However, these visually observed changes were not significant in the permutation test.

Among the significant changes identified, the permutation test showed that a cluster of increased signals occurring across the high β and low γ range (24.5-31 Hz) in frontal areas was significantly different in the pre condition than in the med condition. In addition, comparison of single selective channels between conditions showed a significant difference in this bandwidth in the frontal areas recorded from AF7 and FP2.

β waves generally replace α waves when participants open their eyes, and in the motor cortex, β waves are associated with muscle motor activity [58,60]. They are normally most prominent in the frontal and central head regions and attenuate posteriorly. This may have been the case for the increase observed here. β activity is also commonly associated with drowsiness, stage nonrapid eye movement 1 sleep, and subsequently decreases in deeper sleep, and β activity is not affected by eye opening [60]. Interestingly, sedative medications are also known to increase the amplitude and quantity of β activity [69]. This finding suggests that an increase in the power of the β range might potentially be useful as a neurophysiological correlate of VR-guided meditation. Nevertheless, findings regarding changes in the β band with meditation have been inconsistent. Several studies have reported no significant changes associated with meditation in the β range [4,31]. However, an increase in β and  band power was reported in one study after a longer period of 6 weeks of meditation compared with baseline [70]. It is possible that the changes in β power in the meditation condition compared with the premeditation condition here may be more specific to the use of VR-guided meditation and the visual activity involved, but further work is required to explore this.

In contrast, the activity of the γ band has been reported to be associated with activation of the default mode network. The default mode network is most frequently detected during the resting conditions and reflects the neural activity of different brain areas, such as the cingulate cortex, hippocampus, medial frontal lobes, inferior parietal lobes, and temporal lobes. It is thought to be involved in self-consciousness; self-processing; and introspective functions, including emotional awareness and processing [71-74]. An increase in γ band activity in the frontal and prefrontal areas during VR-guided meditation could reflect the activation of such introspective experiences through meditation. The increase in γ power may also have been due to the activation of attentional networks and visual processing of the meditative VR environment [75].

It is also noteworthy that others have reported an increase in γ band power related to meditation [17,75,76]. For example, one study reported an increase in the γ band (25-45 Hz) during meditation compared with the resting state in the temporal and parieto-occipital areas in mindfulness meditation practitioners [76]. Another study reported that the proficiency level of the meditator is associated with the increased level of the γ band (60-110 Hz) in the parieto-occipital region in meditative states relative to the mind-wandering state in experienced meditators compared with healthy controls. Although it is noteworthy, in this study, no significant difference was found between the states of meditation and mind wandering [75].

VR-Guided Meditation Versus Post Meditation

In contrast to the pattern of increase in high β and low γ activity in the meditation condition compared with the prior rest condition, here, we observed a pattern of reduced β and γ power in the post condition compared with the meditation condition in the topographic maps. This was followed by permutation test results in terms of 2 significant clusters of differences in the β cluster (23-36 Hz) and γ cluster (37-50 Hz) in the post versus med condition. These findings suggest a potential regression back to the baseline state activity and further suggest that changes in high β and low γ activity are associated with a VR meditative state. Analysis of single selective channels also supported widespread significant differences in the power of β and γ bandwidths in the frontal, central, and parietal channels (FP2, FC1, CP5, and P5).

As γ band oscillation has also been reported to be associated with attention toward pain and hypervigilance [77,78], the significant reduction in γ band activity following the VR meditation experience could potentially show that less attentional capacity is directed toward pain after using a meditative VR environment. However, this is conjectural and requires further verification.

Pre Meditation Versus Post VR-Guided Meditation

The comparison of the pre and post conditions could provide an indication of a VR-guided meditation effect in our study. These changes were mainly observed in the α frequency range in terms of an increase in α power in the frontal and central areas in the post condition compared with the pre condition. This was accompanied by a significant cluster-based permutation analysis finding over a cluster of channels in the frontal and central areas (8-9.5 Hz) and significant differences in α power in channels such as Fp2, FC1, CP5, and P5. An increase in the θ band in the central areas and a decrease in the θ band in the posterior occipital areas were also observed in the power spectrum analysis; however, these changes were not found to be significant in single-channel analysis. The posterior dominant α rhythm is characteristically present in normal conscious EEG recordings in the occipital region. It is a defining feature of the normal background rhythm of the adult EEG, best observed with the eyes closed and during mental relaxation and is attenuated by eye opening and mental effort. θ waveforms are characteristically observed more in drowsiness and in the early stages of sleep, such as light sleep (the nonrapid eye movement 1 and nonrapid eye movement 2 sleep phases) [60]. Increases in α-θ bandwidths have previously been reported to be associated with mindfulness meditation, and the α-θ border (7-8 Hz) has also been suggested as an optimal range for indicating visualization activity [4,79,80].

A final observation worthy of note is that there was a reduction in the δ range power in the post condition compared with the pre condition, which was significantly different in the LC and parietal channels (CP5 and P5). Moreover, although changes in the δ range were not significant in the pairwise pre-post comparison, a cluster of channels was identified to be significant in MANOVA. δ is seen more in deep, dreamless sleep, and meditation activities where awareness is more detached [4,60]. Such a general pattern of reduction in δ power in the post condition compared with the pre condition could possibly be related to the effect of VR-guided meditation on brain activity and would need further work to explore if this is a significant trend.


The significant coherence changes suggested that variations in brain connectivity occurred between the different test conditions. Coherence was found predominantly between the frontal and parietal and occipital cortices and in different wave bands, namely, in γ for the pre and med conditions; in θ, α, and β for the med and post conditions; and in δ, α, and β for the pre and post conditions. The γ coherence changes between the pre and med conditions were likely associated with activity during the VR-guided meditation. Between the med and post conditions, it shifted to slower frequencies, possibly suggesting a postmeditation effect. The pre-post comparison showed coherence in the δ, α, and β ranges. The reasons for this are unclear and could be related to the individual responses of the participants after the VR-guided meditation.

Pain and EEG Signals

In terms of the secondary focus of the study to explore pain changes correlated with EEG variations, we found no significant association between pain reduction and changes in electrophysiological signal. This could be due to the limited sample size of this exploratory study. More electrophysiological studies on a larger sample population could potentially identify EEG correlates associated with pain reduction after VR-guided meditation.


As an exploratory study, the sample numbers were small and not necessarily representative of the wider population of patients with CP, which limited the power to identify differences. Therefore, there is a need for neurophysiological studies with larger samples to validate these results and to better explore this phenomenon. In addition, this was a single cohort study with no comparison group, although resting states before and after meditation were used as a no-mindfulness within-subjects control. Finally, the study focused on short-term neurophysiological alterations in the electrophysiology of the brain, and the long-term effects of meditative VR environments are still unknown, which will require longitudinal studies.


These findings suggest that distinct altered neurophysiological brain signals are detectable during VR-guided meditation, predominantly in terms of an increase in the power of the β and γ bands. Changes in the α and θ bands were also identified, predominantly as a pattern in VR-guided meditation compared with the resting baseline, possibly reflecting the specific impact of visual activity during VR-guided meditation. Some changes in coherence were also observed between the frontal and parietal and occipital cortices during VR-guided meditation. No significant association between NRS pain scores and changes in EEG signals was observed. Although this is an exploratory study, the results of this work clearly demonstrate the feasibility of EEG recording and subsequent data processing and analysis during VR experiences in patients using modern VR HMDs. To our knowledge, this is the first exploration of EEG alterations in the brain’s electrophysiological signals associated with VR-guided meditation in patients with CP and should provide some valuable initial data to inform future work in this field.


The authors would like to acknowledge that the EEG system used in this study was purchased by Dr Urs Ribary using the Canadian Foundation of Innovation funding. The authors would also like to thank Dr Diane Gromala, Director of the Pain Studies Lab, School of Interactive Arts and Technology, Simon Fraser University, for supporting this study, and all participants for their time and effort in taking part in the study. This work was supported by the Lotte and John Hecht Memorial Foundation (grant 4110).

Conflicts of Interest

None declared.

  1. Chronic pain in Canada: laying a foundation for action. Government of Canada. 2019.   URL: [accessed 2021-05-16]
  2. Kabat-Zinn J. Mindfulness-based interventions in context: past, present, and future. Clin Psychol Sci Pract 2003;10(2):144-156. [CrossRef]
  3. Rosenzweig S, Greeson JM, Reibel DK, Green JS, Jasser SA, Beasley D. Mindfulness-based stress reduction for chronic pain conditions: variation in treatment outcomes and role of home meditation practice. J Psychosom Res 2010 Jan;68(1):29-36. [CrossRef] [Medline]
  4. Lomas T, Ivtzan I, Fu CH. A systematic review of the neurophysiology of mindfulness on EEG oscillations. Neurosci Biobehav Rev 2015 Oct;57:401-410. [CrossRef] [Medline]
  5. Hanson R, Mendius R. Buddha's Brain: The Practical Neuroscience of Happiness, Love & Wisdom. Oakland, California, United States: New Harbringer Publications; 2009:1-272.
  6. Nyklíček I, Dijksman SC, Lenders PJ, Fonteijn WA, Koolen JJ. A brief mindfulness based intervention for increase in emotional well-being and quality of life in percutaneous coronary intervention (PCI) patients: the MindfulHeart randomized controlled trial. J Behav Med 2014 Feb 23;37(1):135-144. [CrossRef] [Medline]
  7. Lee DJ, Kulubya E, Goldin P, Goodarzi A, Girgis F. Review of the neural oscillations underlying meditation. Front Neurosci 2018;12:178 [FREE Full text] [CrossRef] [Medline]
  8. Schoenberg PL, Vago DR. Mapping meditative states and stages with electrophysiology: concepts, classifications, and methods. Curr Opin Psychol 2019 Aug;28:211-217. [CrossRef] [Medline]
  9. Farias M, Wikholm C. Has the science of mindfulness lost its mind? BJPsych Bull 2016 Dec;40(6):329-332 [FREE Full text] [CrossRef] [Medline]
  10. Keng S, Smoski MJ, Robins CJ. Effects of mindfulness on psychological health: a review of empirical studies. Clin Psychol Rev 2011 Aug;31(6):1041-1056 [FREE Full text] [CrossRef] [Medline]
  11. Seay A, Gromala D, Hodges L, Shaw C. The meditation chamber: a debriefing. In: Proceedings of the The 29th International Conference on Computer Graphics and Interactive Techniques. 2002 Presented at: SIGGRAPH02: The 29th International Conference on Computer Graphics and Interactive Techniques; July, 2002; San Antonio Texas p. 263. [CrossRef]
  12. Botella C, Garcia-Palacios A, Vizcaíno Y, Herrero R, Baños RM, Belmonte MA. Virtual reality in the treatment of fibromyalgia: a pilot study. Cyberpsychol Behav Soc Netw 2013 Mar;16(3):215-223. [CrossRef] [Medline]
  13. Erlacher D, Chapin H. Lucid dreaming: Neural virtual reality as a mechanism for performance enhancement: Commentary on “The neurobiology of consciousness: Lucid dreaming wakes up” by J. Allan Hobson. Int J Dream Res 2010;3(1):7-10 [FREE Full text] [CrossRef]
  14. Gromala D, Tong X, Choo A, Karamnejad M, Shaw C. The virtual meditative walk: virtual reality therapy for chronic pain management. In: Proceedings of the 33rd Annual ACM Conference on Human Factors in Computing Systems. 2015 Presented at: CHI '15: CHI Conference on Human Factors in Computing Systems; April, 2015; Seoul Republic of Korea p. 521-524. [CrossRef]
  15. Navarro-Haro MV, Modrego-Alarcón M, Hoffman HG, López-Montoyo A, Navarro-Gil M, Montero-Marin J, et al. Evaluation of a mindfulness-based intervention with and without virtual reality dialectical behavior therapy mindfulness skills training for the treatment of generalized anxiety disorder in primary care: a pilot study. Front Psychol 2019 Jan 28;10:55 [FREE Full text] [CrossRef] [Medline]
  16. Venuturupalli RS, Chu T, Vicari M, Kumar A, Fortune N, Spielberg B. Virtual reality-based biofeedback and guided meditation in rheumatology: a pilot study. ACR Open Rheumatol 2019 Dec 11;1(10):667-675 [FREE Full text] [CrossRef] [Medline]
  17. Cahn BR, Polich J. Meditation states and traits: EEG, ERP, and neuroimaging studies. Psychol Bull 2006 Mar;132(2):180-211. [CrossRef] [Medline]
  18. Davidson RJ, Lutz A. Buddha's brain: neuroplasticity and meditation. IEEE Signal Process Mag 2008 Jan 01;25(1):172-176 [FREE Full text] [CrossRef] [Medline]
  19. Ivanovski B, Malhi GS. The psychological and neurophysiological concomitants of mindfulness forms of meditation. Acta Neuropsychiatr 2007 Apr 24;19(2):76-91. [CrossRef] [Medline]
  20. Panjwani U, Selvamurthy W, Singh S, Gupta H, Mukhopadhyay S, Thakur L. Effect of Sahaja Yoga meditation on Auditory Evoked Potentials (AEP) and Visual Contrast Sensitivity (VCS) in epileptics. Appl Psychophysiol Biofeedback 2000;25(1):1-12. [CrossRef]
  21. Heeger DJ, Ress D. What does fMRI tell us about neuronal activity? Nat Rev Neurosci 2002 Feb;3(2):142-151. [CrossRef] [Medline]
  22. Chiesa A, Serretti A. A systematic review of neurobiological and clinical features of mindfulness meditations. Psychol Med 2010 Aug;40(8):1239-1252. [CrossRef] [Medline]
  23. Magalhaes AA, Oliveira L, Pereira MG, Menezes CB. Does meditation alter brain responses to negative stimuli? A systematic review. Front Hum Neurosci 2018;12:448 [FREE Full text] [CrossRef] [Medline]
  24. Mishra SK, Khosa S, Singh S, Moheb N, Trikamji B. Changes in functional magnetic resonance imaging with meditation: a pilot study. Ayu 2017;38(3-4):108-112 [FREE Full text] [CrossRef] [Medline]
  25. Froeliger BE, Garland EL, Modlin LA, McClernon FJ. Neurocognitive correlates of the effects of yoga meditation practice on emotion and cognition: a pilot study. Front Integr Neurosci 2012;6:48 [FREE Full text] [CrossRef] [Medline]
  26. Goldin P, Ziv M, Jazaieri H, Hahn K, Gross JJ. MBSR vs aerobic exercise in social anxiety: fMRI of emotion regulation of negative self-beliefs. Soc Cogn Affect Neurosci 2013 Jan;8(1):65-72 [FREE Full text] [CrossRef] [Medline]
  27. Delmonte MM. Electrocortical activity and related phenomena associated with meditation practice: a literature review. Int J Neurosci 1984 Nov 07;24(3-4):217-231. [CrossRef] [Medline]
  28. Travis F, Shear J. Focused attention, open monitoring and automatic self-transcending: categories to organize meditations from Vedic, Buddhist and Chinese traditions. Conscious Cogn 2010 Dec;19(4):1110-1118. [CrossRef] [Medline]
  29. Dunn B, Hartigan J, Mikulas W. Concentration and mindfulness meditations: unique forms of consciousness? Appl Psychophysiol Biofeedback 1999;24(3):147-165. [CrossRef]
  30. Fell J, Axmacher N, Haupt S. From alpha to gamma: electrophysiological correlates of meditation-related states of consciousness. Med Hypotheses 2010 Aug;75(2):218-224. [CrossRef] [Medline]
  31. Lagopoulos J, Xu J, Rasmussen I, Vik A, Malhi GS, Eliassen CF, et al. Increased theta and alpha EEG activity during nondirective meditation. J Altern Complement Med 2009 Nov;15(11):1187-1192. [CrossRef] [Medline]
  32. Eccleston C. Chronic pain and distraction: an experimental investigation into the role of sustained and shifting attention in the processing of chronic persistent pain. Behav Res Ther 1995 May;33(4):391-405. [CrossRef]
  33. Johnson MH. How does distraction work in the management of pain? Curr Pain Headache Rep 2005 Apr;9(2):90-95. [CrossRef] [Medline]
  34. Garrett B, Taverner T, Gromala D, Tao G, Cordingley E, Sun C. Virtual reality clinical research: promises and challenges. JMIR Serious Games 2018 Oct 17;6(4):e10839 [FREE Full text] [CrossRef] [Medline]
  35. Jin W, Choo A, Gromala D, Shaw C, Squire P. A virtual reality game for chronic pain management: a randomized, controlled clinical study. In: Westwood JD, Westwood SW, Fellander-Tsai L, Fidopiastis CM, Liu A, editors. Medicine Meets Virtual Reality 22: NextMed / MMVR22. Amsterdam, Netherlands: IOS Press; 2016:154-160.
  36. Shaw C, Gromala D, Seay A. The meditation chamber: enacting autonomic senses. In: Proceedings of the 4th International Conference on Enactive Interfaces. 2007 Presented at: ENACTIVE/07 4th International Conference on Enactive Interfaces; November 19-22, 2007; Grenoble, France p. 1-4   URL:
  37. Garrett B, Taverner T, Masinde W, Gromala D, Shaw C, Negraeff M. A rapid evidence assessment of immersive virtual reality as an adjunct therapy in acute pain management in clinical practice. Clin J Pain 2014 Dec;30(12):1089-1098. [CrossRef] [Medline]
  38. Maani CV, Hoffman HG, Morrow M, Maiers A, Gaylord K, McGhee LL, et al. Virtual reality pain control during burn wound debridement of combat-related burn injuries using robot-like arm mounted VR goggles. J Trauma 2011 Jul;71(1 Suppl):125-130 [FREE Full text] [CrossRef] [Medline]
  39. Scapin S, Echevarría-Guanilo ME, Fuculo Jr PR, Gonçalves N, Rocha PK, Coimbra R. Virtual Reality in the treatment of burn patients: a systematic review. Burns 2018 Sep;44(6):1403-1416. [CrossRef] [Medline]
  40. Schmitt YS, Hoffman HG, Blough DK, Patterson DR, Jensen MP, Soltani M, et al. A randomized, controlled trial of immersive virtual reality analgesia, during physical therapy for pediatric burns. Burns 2011 Feb;37(1):61-68 [FREE Full text] [CrossRef] [Medline]
  41. Wolitzky K, Fivush R, Zimand E, Hodges L, Rothbaum BO. Effectiveness of virtual reality distraction during a painful medical procedure in pediatric oncology patients. Psychol Health 2005 Dec;20(6):817-824. [CrossRef]
  42. Garrett B, Taverner T, McDade P. Virtual reality as an adjunct home therapy in chronic pain management: an exploratory study. JMIR Med Inform 2017 May 11;5(2):e11 [FREE Full text] [CrossRef] [Medline]
  43. Li S, Kay S, Hardicker N. Virtual reality: Towards a novel treatment environment for ankylosing spondylitis. In: Bos L, Blobel B, editors. Medical and Care Compunetics 4. Amsterdam, Netherlands: IOS Press; 2007:190-196.
  44. Moore DJ, Keogh E, Eccleston C. The interruptive effect of pain on attention. Q J Exp Psychol (Hove) 2012 Jan;65(3):565-586. [CrossRef] [Medline]
  45. Dascal J, Reid M, IsHak WW, Spiegel B, Recacho J, Rosen B, et al. Virtual reality and medical inpatients: a systematic review of randomized, controlled trials. Innov Clin Neurosci 2017;14(1-2):14-21 [FREE Full text] [Medline]
  46. Li A, Montaño Z, Chen VJ, Gold JI. Virtual reality and pain management: current trends and future directions. Pain Manag 2011 Mar;1(2):147-157 [FREE Full text] [CrossRef] [Medline]
  47. Malloy KM, Milling LS. The effectiveness of virtual reality distraction for pain reduction: a systematic review. Clin Psychol Rev 2010 Dec;30(8):1011-1018. [CrossRef] [Medline]
  48. Johanson GA, Brooks GP. Initial scale development: sample size for pilot studies. Educ Psychol Meas 2009 Dec 18;70(3):394-400. [CrossRef]
  49. Boonstra TW, Nikolin S, Meisener A, Martin DM, Loo CK. Change in mean frequency of resting-state electroencephalography after transcranial direct current stimulation. Front Hum Neurosci 2016;10:270 [FREE Full text] [CrossRef] [Medline]
  50. Konicar L, Radev S, Silvoni S, Bolinger E, Veit R, Strehl U, et al. Balancing the brain of offenders with psychopathy? Resting state EEG and electrodermal activity after a pilot study of brain self-regulation training. PLoS One 2021;16(1):e0242830 [FREE Full text] [CrossRef] [Medline]
  51. Sikka P, Revonsuo A, Noreika V, Valli K. EEG frontal alpha asymmetry and dream affect: alpha oscillations over the right frontal cortex during REM sleep and presleep wakefulness predict anger in REM sleep dreams. J Neurosci 2019 Jun 12;39(24):4775-4784 [FREE Full text] [CrossRef] [Medline]
  52. Williamson A, Hoggart B. Pain: a review of three commonly used pain rating scales. J Clin Nurs 2005 Aug;14(7):798-804. [CrossRef] [Medline]
  53. Oostenveld R, Fries P, Maris E, Schoffelen J. FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput Intell Neurosci 2011;2011:156869 [FREE Full text] [CrossRef] [Medline]
  54. Makeig S, Jung T, Bell A, Sejnowski T. Independent component analysis of electroencephalographic data. In: Proceedings of the 8th International Conference on Neural Information Processing Systems. 1995 Presented at: 8th International Conference on Neural Information Processing Systems; November 27-30, 1995; Denver, CO, USA p. 145-151. [CrossRef]
  55. Levy W. Effect of epoch length on power spectrum analysis of the EEG. Anesthesiology 1987 Apr;66(4):489-495 [FREE Full text] [CrossRef] [Medline]
  56. Bigdely-Shamlo N, Mullen T, Kothe C, Su K, Robbins KA. The PREP pipeline: standardized preprocessing for large-scale EEG analysis. Front Neuroinform 2015;9:16 [FREE Full text] [CrossRef] [Medline]
  57. Perrin F, Pernier J, Bertrand O, Echallier J. Spherical splines for scalp potential and current density mapping. Electroencephalogr Clin Neurophysiol 1989 Feb;72(2):184-187. [CrossRef]
  58. Johnson L, Lubin A, Naitoh P, Nute C, Austin M. Spectral analysis of the EEG of dominant and non-dominant alpha subjects during waking and sleeping. Electroencephalogr Clin Neurophysiol 1969 Apr;26(4):361-370. [CrossRef]
  59. Cooley JW, Tukey JW. An algorithm for the machine calculation of complex Fourier series. Math Comp 1965 May 01;19(90):297. [CrossRef]
  60. Nayak CS, Anilkumar AC. EEG Normal Waveforms. Treasure Island, FL, U.S: StatPearls Publishing; 2020.
  61. Sandwell DT. Biharmonic spline interpolation of GEOS-3 and SEASAT altimeter data. Geophys Res Lett 2012 Dec 07;14(2):139-142. [CrossRef]
  62. Classen J, Gerloff C, Honda M, Hallett M. Integrative visuomotor behavior is associated with interregionally coherent oscillations in the human brain. J Neurophysiol 1998 Mar 01;79(3):1567-1573 [FREE Full text] [CrossRef] [Medline]
  63. Maris E, Oostenveld R. Nonparametric statistical testing of EEG- and MEG-data. J Neurosci Methods 2007 Aug 15;164(1):177-190. [CrossRef] [Medline]
  64. Cohen J. Statistical Power Analysis for the Behavioural Sciences. New York: Routledge; 1988:1-567.
  65. Sawilowsky SS. New effect size rules of thumb. J Mod App Stat Meth 2009 Nov 01;8(2):597-599. [CrossRef]
  66. Faul F, Erdfelder E, Lang A, Buchner A. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Method 2007 May;39(2):175-191. [CrossRef]
  67. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B 2018 Dec 05;57(1):289-300. [CrossRef]
  68. Bakdash JZ, Marusich LR. Repeated measures correlation. Front Psychol 2017 Apr 7;8:456 [FREE Full text] [CrossRef] [Medline]
  69. Louis E, Frey L. Electroencephalography (EEG): An Introductory Text and Atlas of Normal and Abnormal Findings in Adults, Children, and Infants. Chicago: American Epilepsy Society; 2016.
  70. Ahani A, Wahbeh H, Nezamfar H, Miller M, Erdogmus D, Oken B. Quantitative change of EEG and respiration signals during mindfulness meditation. J Neuroeng Rehabil 2014 May 14;11(1):87 [FREE Full text] [CrossRef] [Medline]
  71. Imperatori C, Della Marca G, Brunetti R, Carbone GA, Massullo C, Valenti EM, et al. Default Mode Network alterations in alexithymia: an EEG power spectra and connectivity study. Sci Rep 2016 Nov 15;6:36653 [FREE Full text] [CrossRef] [Medline]
  72. Garrison KA, Zeffiro TA, Scheinost D, Constable RT, Brewer JA. Meditation leads to reduced default mode network activity beyond an active task. Cogn Affect Behav Neurosci 2015 Sep;15(3):712-720 [FREE Full text] [CrossRef] [Medline]
  73. Buckner R, Andrews-Hanna J, Schacter D. The brain's default network: anatomy, function, and relevance to disease. Ann N Y Acad Sci 2008 Mar;1124:1-38. [CrossRef] [Medline]
  74. Greicius MD, Krasnow B, Reiss AL, Menon V. Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc Natl Acad Sci U S A 2003 Jan 07;100(1):253-258 [FREE Full text] [CrossRef] [Medline]
  75. Braboszcz C, Cahn BR, Levy J, Fernandez M, Delorme A. Increased gamma brainwave amplitude compared to control in three different meditation traditions. PLoS One 2017;12(1):e0170647 [FREE Full text] [CrossRef] [Medline]
  76. Berkovich-Ohana A, Glicksohn J, Goldstein A. Mindfulness-induced changes in gamma band activity - implications for the default mode network, self-reference and attention. Clin Neurophysiol 2012 Apr;123(4):700-710. [CrossRef] [Medline]
  77. Tiemann L, Schulz E, Winkelmann A, Ronel J, Henningsen P, Ploner M. Behavioral and neuronal investigations of hypervigilance in patients with fibromyalgia syndrome. PLoS One 2012;7(4):e35068 [FREE Full text] [CrossRef] [Medline]
  78. Tiemann L, Schulz E, Gross J, Ploner M. Gamma oscillations as a neuronal correlate of the attentional effects of pain. Pain 2010 Aug;150(2):302-308. [CrossRef] [Medline]
  79. Belkofer CM, Van Hecke AV, Konopka LM. Effects of drawing on alpha activity: a quantitative EEG Study with implications for art therapy. Art Ther 2014 Jun 13;31(2):61-68. [CrossRef]
  80. Faber PL, Lehmann D, Gianotti LR, Milz P, Pascual-Marqui RD, Held M, et al. Zazen meditation and no-task resting EEG compared with LORETA intracortical source localization. Cogn Process 2015 Feb 5;16(1):87-96. [CrossRef] [Medline]

CP: chronic pain
CMS: Command Mode Sense
DRL: Driven Right Leg
EEG: electroencephalograph
FDR: false discovery rate
HMD: head-mounted display
ICA: independent component analysis
LC: left central
LF: left frontal
LP: left parietal
MBSR: mindfulness-based stress reduction
MANOVA: multivariate analysis of variance
NRS: numerical rating scale
PSD: power spectral density
RCT: randomized controlled trial
RF: right frontal
VR: virtual reality

Edited by G Eysenbach; submitted 07.12.20; peer-reviewed by C Fowler, G Lewis, E Kutafina, A Mastropietro; comments to author 03.02.21; revised version received 23.04.21; accepted 29.04.21; published 24.06.21


©Henry Fu, Bernie Garrett, Gordon Tao, Elliott Cordingley, Zahra Ofoghi, Tarnia Taverner, Crystal Sun, Teresa Cheung. Originally published in JMIR Biomedical Engineering (, 24.06.2021.

This is an open-access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Biomedical Engineering, is properly cited. The complete bibliographic information, a link to the original publication on, as well as this copyright and license information must be included.