Sorry. It's obsolete page. New address of page
http://mashinva.narod.ru/researches_en.html



The Projects and Studies.

20.01.2004 [ Russian Version ]

[ Home ] [ NPP safety ] [ Links ] [ About author ] [ Send mail ]

  • Program MABP.DBase-PFO v.2.01 (20.01.2004).

  • The process of the psychological support to activity of professionals (the selection, the current examinations) requires from psychologist the analysis of dynamics of test individual results, the compare test individual results to referent sample results, the building of the psychological conclusion on base of the analysis. For this purpose us was designed the progarm MABP.DBase-PFO v.2.01, in which are solved put above problems: entering, editing, analysis and making up the psychological conclusion on results of the examinations.
    Source information is kept in format Microsoft Access-2000. The Database is designed by Delphi-5. The data access is realized by means of component ADO (MDAC versions 2.5). In more detail gets acquainted with possibility presented program possible on: Short Help of Program MABP.DBase-PFO v.2.01. (In russian language)

  • Analysis of additional indexes of heart rate graph (24.09.2002)

  • Researching structure of heart rate graph, we came to conclusion about practicability of consideration of the following additional indexes: The Square of least distance from graph node before axis ó calculated on formula: QX=0.5*(RR[i]-RR[i+1])^2. Not difficult notice the generality Sum(QX) and RMSSD (Root Mean Sum Successful Deviation) formulas.
    The fact that medium of square (in which heart rate graph was inserted) not always comply with mean value of RR-intervals (RRav) was necessary take into account for calculation QY. In other words, there were events, when values of VR/2 were more or less RRav-RRmin. Herewith "center of gravity" of graph was displaced on axis ó or upwards, or downwards from the center of square. With provision for correction on RRav formula for calculation of square of least distance from graph node before axis Y has gained the following type: QY=0.5*(RR[i]+RR[i+1]-2*RRmin-2*RRav)^2.
    In table 1 are presented values of Pearson correlation (r) for different indexes of heart rate variability (the registration RR-intervals was conducted in rest state; more full description of procedure a data processing is presented in articles Analysis heart rate variability by means of graph method and Factorial analysis of heart rate graph indexes). N=10036, ln-natural logarithm, HF-power of spectrum of heart rate in high-frequency range, PW-general power of spectrum, SD-standard deviation.

    Table 1.
    r, p=0.00lnRMSSDHFlnSDlnPW
    Sum(QX)0.56800.96620.47960.4823
    lnSum(QX)0.99350.64940.88430.8869
    lnSum(QY)0.82870.48010.98900.9835

    According to Table 1, sum of squares of least distances from graph nodes before axis ó has very high correlations with tone of parasympathetic nervous system indexes - RMSSD and HF (correlation between lnRMSSD and lnHF is 0.9631). At the same time, sum of squares of least distances from graph nodes before axis Y has very high correlations with the general heart rate variability indexes - SD and PW (correlation between SD^2 and PW is 0.9938, between lnSum(QY) and ln(lnSD) is 0.9920).
    The high correlation of relation Sum(QX)/Sum(QY) with HFn (HFn=100*HF/PW) and LF/HF (LF-power of spectrum of heart rate in low-frequency range) indexes logistical was expect. The correlation analysis has confirmed this. The results are presented in Table 2 (PXY=Sum(QX)/Sum(QY)). The high values of correlation for DFA (Detrended Fluctuation Analysis) are presented in addition.

    Table 2.
    r, p=0.00HFnlnHFnln(LF/HF)DFA
    lnPXY0.85600.8607-0.7671-0.7344
    DFA-0.7369-0.75740.4881-

    On the following step we have calculated the values of PXY index for harmonic fluctuation since different period and amplitude. The analysis has shown presence of relationship PXY with frequency of harmonic fluctuations. The fluctuations since odd periods (for example, 5, 15 sec.) have formed exception that follows from nature of graph (asymmetry of graph comparatively axises X and Y are appeared).
    But real heart rate contains the ensemble of dynamic changing fluctuations of different frequency and amplitudes. Following hypothesis was formulated: PXY-index reflects leading (central) fluctuation frequency, having most influence upon heart rate. Spectrograms of RR-intervals were analyzed, frequencies of peaks with maximum spectral density were determined and were compared with PXY for ground of this hypothesis . Correlation between PXY and frequency a pica with maximum spectral density (fmax) though was enough high (r=0.7035, p=0.00), but has required the more detailed analysis of real rhythmograms. This has allowed selecting three variants of relations between PXY and fmax.
    1) Exists very high coincidence of PXY and fmax values (for instance, fmax=0.082, PXY=0.080). It is characterized of active rest state (readiness to execution of activity). In this state the high influence of autonomic nervous system (ANS) on heart rate regulation can expect.
    2) Growing PXY vastly overtakes the growing fmax (for instance, fmax=0.301, PXY=1.082). Given process is typical for condition of deep relaxation. This is accompanied the reduction role of ANS in heart rate regulation and growing of chaotic nature of process. PXY values strive to 1.0. Exactly such values were received at modeling "white noise".
    3) Reduction fmax is accompanied the growing PXY (for instance, fmax=0.016, SD=3257 msec2, PXY=0.328 - chronic overstrain, registration in rest state). Given process is typical for high psychic strain and different pathological processes. It is probably in this case normal regulation of heart rate through ANS becomes impossible, that it is also accompanied the growing of chaotic nature of process, but with very low amplitude (the denominated stabilization of rhythm).
    To consider the events, when spectrogram contains the significant peaks in different ranges of frequencies. In this case PXY index presents itself resulting frequencies two fluctuations. For instance, at presence of peaks of spectral density with frequency 0.102 Hz (LF-range, SD=337420 msec2) and 0.207 Hz (HF-range, SD=31345 msec2), PXY value was 0.132 Hz. Note also that PXY index is most sensitive to fluctuations in ranges LF and HF.
    For four groups of man: 1 - normal state, 2 - neurotic excitement, 3 - functional fatigue, 4 - psychic overstrain were calculated average values of analyzed graph indexes on each group (Table 3, Nv - number of samples, Sum(QX) and Sum(QY) values are given in msec2).

    Table 3.
    Group1234
    Sum(QX)345576105006217543626
    Sum(QY)85391618036299131925783
    PXY0.4080.0550.2700.1675
    Nv6311508452478

    Either as was expected, maximum Sum(QX) value (the indicator to activities of parasympathetic nervous system) was received for first group ("Norma"), and minimum - for fourth group ("Overstrain"). Maximum Sum(QY) value (the factor of general heart rate variability) is characteristic of neurotic excitement (the influence of corticolimbic brain structures systems, note the significant reduction PXY index of leading (resulting) of frequency of fluctuations of heart rate). High PXY value for group "Overstrain" reflects earlier considered dynamic of increasing PXY at reduction fmax for events of rhythm stabilizations.
    Use the SumQX, SumQY and PXY indexes in discriminate analysis have allowed greatly to raise accuracy of differential diagnostics presented above functional states.
    In conclusion we will afford to return to fact of difference medium of variational range and RRav-RRmin. There are three possible variants:
    1) VR/2-(RRav-RRmin)>0 - it is characteristic of for connecting process: fluent (long-time) change one functional state others, firm emotional excitement.
    2) VR/2-(RRav-RRmin)<0 - it is characteristic of for violation of heart rate regulation: sharp (short time) change one functional state others, sharp wave of emotional excitement.
    3) VR/2-(RRav-RRmin)=0 - it is characteristic of for stationary process.
    For validation this hypothesis we have compared own results with standard estimations of stationarity (through checking the constancy of mean and dispersion). The difference between VR/2 and (RRav-RRmin) have expressed in percent from (RRav-RRmin): Dvm=100*(VR/2-(RRav-RRmin))/(RRav-RRmin). For estimation of relations between two categorical variables on criterion Pearson Chi-square we used the registration of RR-intervals in rest state (N=10036). We have calculated stationarity on Dvm when difference did not exceed 1, 2 and 3 percent (Table 4).

    Table 4.
    Dvm1%2%3%
    Pearson Chi-square2.15588.935611.0789
    p0.142040.002800.00087

    According to results are presented in Table 4, Dvm index really can serve the estimation of stationarity time series. The analysis of real rhythmogram has shown that choice Dvm=1% brings about more exact diagnostics stationarity than in the event of standard estimation. Certainly, following studies must confirm or refuse our hypothesis.

  • Heart rate graph

  • According to studies heart rate graph by means of method a graph (in Russian),heart rate graph is characterized the number of nodes (ND) and ribs (NRib). The decimal logarithm ND is nor than other, as index of system complexity (HmND). Knowing probability of the nodes distribution, easy calculate the index of entropy - the measures of system uncertainty (ìND). The relative organization of system (RND) will is: 1 - HND/HmND. Enumerated indexes were calculated both for ND, and for NRib.
    For four groups: 1 - norm, 2 - neurotic state, 3 - asthenic state, 4 - tensity state the average values were calculated (the sliding 256 RR-intervals samples with a step of 10 RR-intervals were used). Average indexes on each groups are presented in table 1, as well as F-criterion values for analysis of variance (ANOVA, test for evaluation of significant differences between averages). (Nv - number of samples).

    Table 1.
    GroupHNDHmNDRNDHRibHmRibRRib Nv
    12.152.210.0292.372.380.004 1042
    22.052.140.0432.292.310.011 1610
    31.781.900.0652.192.240.022 441
    41.251.460.1551.621.750.080 322
    F380038152508335027523225 322

    According to table 1, complexity level of system (both on number nodes, and on number ribs) is most high in calm state (norm) and decreases at growing of tensity. The indexes of entropy (the level to uncertainties) and relative organization of system (both on number nodes, and on number ribs) are minimum for first group and maximum for tensity state.
    In addition the index ND/NRib was calculated. Correlation analysis has revealled the interesting regularity between this index and number of nodes. If for first three groups the values of r (Pearson coefficient of correlation) indicated to high positive coorelation between ND and ND/NRib (0.98, 0.94 and 0.87, accordingly), then for tensity state r=-0.75. This possible comment as follows: in usual conditions a growing or reduction of ND is accompanied the identical growing or reduction of NRib. In tensity state, reduction of ND is accompanied comparatively more denominated reduction of NRib.In the future is planned to continue the studies of heart rate graph indexes by means of discriminant analysis.

  • Detrended Fluctuation Analysis

  • One of the methods of nonlinear estimation serves Detrended Fluctuation Analysis (DFA is quantify the fractal correlation properties of time series, Peng et al. 1995). In article Peng C-K, Hausdorff JM, Goldberger AL. Fractal mechanisms in neural control: Human heartbeat and gait dynamics in health and disease. In: Walleczek J, ed. Nonlinear Dynamics, Self-Organization, and Biomedicine. Cambridge: Cambridge University Press, 1999. is kept enough full description DFA. But for calculations better to use the programs were kindly given by authors: the detrended fluctuation analysis (DFA) software package. (Except programs, zip-file contains the source texts in C-code.) Recent studies have shown that calculated by means of procedures DFA index alfa in norm is aproximately 1.00, reduction for patients with cardiovascular disorders, but increase for elderly people (Peng et al. 1995, Ho et al. 1997, Hausdorff et al. 1995, Iyengar et al. 1996). But note that these results were received for long-term records of RR-intervals (24 hours).
    Using data on next groups: 1 - norm, 2 - neurotic state, 3 - asthenic state, 4 - tensity state, alfa indexes were calculated on procedure DFA for short-term records (the sliding 256 RR-intervals samples with a step of 10 RR-intervals were used). Then correlations of alfa with variables of statistical and spectral analysis, as well as heart rate graph indexes were counted for each group. The significant positive correlation (p=0.000) alfa with normalized values of power spectrum in Þ very low-frequency (VLFn=100*VLF/TP, TP - overall power of heart rate spectrum, in msec2) was discovered. In table 2 the average alfa indexes, r - values of Pearson correlations between alfa and VLFn, as well as number of samples on each group (Nv) were brought.

    Table 1.
    Group alfa r Nv
    1 0.688 0.931042
    2 1.086 0.78 1610
    3 0.820 0.89 441
    4 0.921 0.52 322

    From tables 1 is seen that reduction of tensity is connected with a decrease of alfa (to 0.5). In tensity state alfa approaches to 1.00. In neurotic state alfa exceeds 1.00.
    The values of alfa and VLFn have a high positive correlation (the increase of VLFn-values are connected with growing of alfa) for all groups, except fourth. In tensity state a correlation between alfa and VLFn sharply falls (the dependency between variables becomes not obvious). Curiously enough note that alfa has negative correlation (but more denominated) with HFn (HFn = 100*HF/TP).
    According to ANOVA analysis (on four groups) for alfa F=2214, for NRib F=4468, for ND F=3814.
    Hereinafter the studies of possibility of using DFA on short samples for differential diagnostics of different functional conditions is planned to continue.

  • The Remarks to Project

  • At present time the group of specialists under the direction of R.ë.Baevsky is developping The Project of "Methodical recommendations on analysis heart rate variability use different ECG-systems (the analysis "short" record). Moscow, 2001. (In Russian)"
    On the grounds of studies have conducted in psychophysiological laboratory of Atom Energy Training Center, Novovoronezh, Russia, The Remarks to Project of "Methodical recommendations on analysis heart rate variability use different ECG-systems (the analysis "short" record). Moscow, 2001 (in Russian)" were prepared.



    [ Home ] [ NPP safety ] [ Links ] [ About author ] [ Send mail ]

    ©2001-2004 Vladimir Mashin