A Monte Carlo Study to Evaluate the Performance of Selected Tests of Homogeneity in Sparse Contingency Tables

We address the problem of tests of homogeneity in two-way contingency tables in case-control studies when the case category is subdivided into k subcategories. In this situation, we have two cells with large frequencies and 2 X k cells with frequencies that become small as k increases. We propose two ad hoc statistics in which a statistic for the sparse cells is combined with a statistic for the cells with large frequencies. We will study these tests along with the Pearson test (using a chi-square approximation) in a Monte Carlo simulation study. Two sets of null hypothesis models and two sets of alternative hypothesis models are considered. The best test for the models considered is the usual Pearson test (using an approximate chi-square distribution) although the ad hoc models are more powerful under one alternative model considered.     

Statistical quality control of variety trial data

I propose detection criteria for identifying an abnormal or erroneous data vector provided by a single variety trial in a longer series of variety trials. The test criteria are based on the linear effects estimated separately for each studied trial using global variance components estimated from the whole series of variety trials. The criteria comprise three mutually independent test statistics. The first one is a quadratic form in the estimated fixed effects, the second one is a quadratic form in the estimated realized linear random effects not including the residual effects, and the third one is a quadratic form in the estimated realized residual effects. Under the null hypothesis defining a valid data vector, the three quadratics have independent chi2 distributions. Under natural alternative hypotheses, they have noncentral chi2 distributions. Decomposing the total variation of the data vector studied into quadratic forms due to different types of the realized linear effects intuitively justifies the resulting test criteria. The decomposition may also be used to show that the resulting tests are likelihood ratio tests. I further present computational procedures that allow us to dispense with the need for prior estimation of the linear effects.     

Spectra-based estimates of certain time-domain parameters of a bivariate stationary-point process.

We consider estimates of certain time-domain parameters of a bivariate stationary-point process based on modified periodogram statistics. These estimates are shown to be asymptotically normal under regularity conditions. In the computations of the estimates, the advantage of using the FFT algorithm is demonstrated. Three examples, obtained by analyzing two data sets from the field of neurophysiology, are illustrated.     

Numerical Investigation of a Branch-Shaped Filter Based on Metal-Insulator-Metal Waveguide

Using the finite difference time-domain method, we present a comprehensive numerical investigation of a branch-shaped filter based on the metal-insulator-metal (MIM) waveguide. The results show that several passbands and stopbands appear in the transmission spectra, which are resulted by the phase differences between the surface plasmon polaritons (SPPs) propagating along the straight waveguide and the SPPs resonating in the circuit formed by the branch and the straight waveguide. The effects of the structural parameters of the branch-shaped filters on their transmission properties are also studied. These results not only present an alternative plasmonic filter for the MIM waveguides but also help us to understand the transmission properties of the circuit-shaped structures.     

Features of anti‐inflammatory effects of modulated extremely high‐frequency electromagnetic radiation

Using a model of acute zymosan‐induced paw edema in NMRI mice, we test the hypothesis that anti‐inflammatory effects of extremely high‐frequency electromagnetic radiation (EHF EMR) can be essentially modified by application of pulse modulation with certain frequencies. It has been revealed that a single exposure of animals to continuous EHF EMR for 20 min reduced the exudative edema of inflamed paw on average by 19% at intensities of 0.1–0.7 mW/cm² and frequencies from the range of 42.2–42.6 GHz. At fixed effective carrier frequency of 42.2 GHz, the anti‐inflammatory effect of EHF EMR did not depend on modulation frequencies, that is, application of different modulation frequencies from the range of 0.03–100 Hz did not lead to considerable changes in the effect level. On the contrary, at “ineffective” carrier frequencies of 43.0 and 61.22 GHz, the use of modulation frequencies of 0.07–0.1 and 20–30 Hz has allowed us to restore the effect up to a maximal level. The results obtained show the critical dependence of anti‐inflammatory action of low‐intensity EHF EMR on carrier and modulation frequencies. Within the framework of this study, the possibility of changing the level of expected biological effect of modulated EMR by a special selection of combination of carrier and modulation frequencies is confirmed. Bioelectromagnetics 30:454–461, 2009. © 2009 Wiley‐Liss, Inc.     

Effect of low frequency modulated microwave exposure on human EEG: individual sensitivity

The aim of this study was to evaluate the effect of modulated microwave exposure on human EEG of individual subjects. The experiments were carried out on four different groups of healthy volunteers. The 450 MHz microwave radiation modulated at 7 Hz (first group, 19 subjects), 14 and 21 Hz (second group, 13 subjects), 40 and 70 Hz (third group, 15 subjects), 217 and 1000 Hz (fourth group, 19 subjects) frequencies was applied. The field power density at the scalp was 0.16 mW/cm(2). The calculated spatial peak SAR averaged over 1 g was 0.303 W/kg. Ten cycles of the exposure (1 min off and 1 min on) at fixed modulation frequencies were applied. All subjects completed the experimental protocols with exposure and sham. The exposed and sham-exposed subjects were randomly assigned. A computer also randomly assigned the succession of modulation frequencies. Our results showed that microwave exposure increased the EEG energy. Relative changes in the EEG beta1 power in P3-P4 channels were selected for evaluation of individual sensitivity. The rate of subjects significantly affected is similar in all groups except for the 1000 Hz group: in first group 3 subjects (16%) at 7 Hz modulation; in second group 4 subjects (31%) at 14 Hz modulation and 3 subjects (23%) at 21 Hz modulation; in third group 3 subjects (20%) at 40 Hz and 2 subjects (13%) at 70 Hz modulation; in fourth group 3 subjects (16%) at 217 Hz and 0 subjects at 1000 Hz modulation frequency.     

Effect of Ca ion concentration on cross-bridge kinetics in rabbit psoas fibers. Evidence for the presence of two Ca-activated states of thin filament

The effect of Ca ion concentration on cross-bridge kinetics in a small bundle (one to three fibers) of chemically skinned rabbit psoas muscle is studied. The length of the muscle is oscillated in small amplitude sine waves (0.2% L0 peak-to-peak) at varying frequencies (0.125 -- 167 Hz), and the resulting amplitude and phase shift in tension are measured. The frequency response function (complex stiffness) thus obtained can be divided into three parts, which we name process (A) (centered at 1 Hz), process (B) (3--17 Hz), and process (C) (50 Hz). Process (B), which represents oscillatory work, further splits into two processes (B' and B) at partial Ca activation (less than 50% P0), where the phase-frequency plot appears W-shaped. The slower of the two processes (B') disappears by full activation, at which time the plot appears V-shaped. The characteristic frequencies associated with the minima of the plot do not shift in a graded way with Ca concentration, indicating that there is no change in apparent rate constants. Apparent rate constants of processes (A) and (C) are minimally affected by Ca. The above results are not altered when ionic strength is changed between 128 and 265 mM. We propose that activated thin filaments can have two "on" states and that Ca concentration controls the distribution of these two states. This mechanism generally supports the "switch" hypothesis of Ca regulation.     

Comparison between low-level 50 Hz and 900 MHz electromagnetic stimulation on single channel ionic currents and on firing frequency in dorsal root ganglion isolated neurons

Alteration of membrane surface charges represents one of the most interesting effects of the electromagnetic exposure on biological structures. Some evidence exists in the case of extremely low frequency whereas the same effect in the radiofrequency range has not been detected. Changes in transmembrane voltages are probably responsible for the mobilization of intracellular calcium described in some previous studies but not confirmed in others. These controversial results may be due to the cell type under examination and/or to the permeability properties of the membranes. According to such a hypothesis, calcium oscillations would be a secondary effect due to the induced change in the membrane voltage and thus dependent on the characteristics of ionic channels present in a particular preparation. Calcium increases could suggest more than one mechanism to explain the biological effects of exposure due to the fact that all the cellular pathways using calcium ions as a second messenger could be, in theory, disturbed by the electromagnetic field exposure. In the present work, we investigate the early phase of the signal transmission in the peripheral nervous system. We present evidence that the firing rate of rat sensory neurons can be modified by 50/60 Hz magnetic field but not by low level 900 MHz fields. The action of the 50/60 Hz magnetic field is biphasic. At first, the number of action potentials increases in time. Following this early phase, the firing rate decreases more rapidly than in control conditions. The explanation can be found at the single-channel level. Dynamic action current recordings in dorsal root ganglion neurons acutely exposed to the electromagnetic field show increased functionality of calcium channels. In parallel, a calcium-activated potassium channel is able to increase its mean open time.     

Effect of electromagnetic field on cyclic adenosine monophosphate (cAMP) in a human mu-opioid receptor cell model

ABSTRACT

During the cell communication process, endogenous and exogenous signaling affect normal as well as pathological developmental conditions. Exogenous influences such as extra-low-frequency electromagnetic field (EMF) have been shown to effect pain and inflammation by modulating G-protein receptors, down-regulating cyclooxygenase-2 activity, and affecting the calcium/calmodulin/nitric oxide pathway. Investigators have reported changes in opioid receptors and second messengers, such as cyclic adenosine monophosphate (cAMP), in opiate tolerance and dependence by showing how repeated exposure to morphine decreases adenylate cyclase activity causing cAMP to return to control levels in the tolerant state, and increase above control levels during withdrawal. Resonance responses to biological systems using exogenous EMF signals suggest that frequency response characteristics of the target can determine the EMF biological response. In our past research we found significant down regulation of inflammatory markers tumor necrosis factor alpha (TNF-α) and nuclear factor kappa B (NFκB) using 5?Hz EMF frequency. In this study cAMP was stimulated in Chinese Hamster Ovary (CHO) cells transfected with human mu-opioid receptors, then exposed to 5?Hz EMF, and outcomes were compared with morphine treatment. Results showed a 23% greater inhibition of cAMP-treating cells with EMF than with morphine. In order to test our results for frequency specific effects, we ran identical experiments using 13?Hz EMF, which produced results similar to controls. This study suggests the use of EMF as a complementary or alternative treatment to morphine that could both reduce pain and enhance patient quality of life without the side-effects of opiates.     

Comparing Power Behaviours for Some Goodness-of-Fit Test Statistics in Sparse Multinomials

This paper is concerned with the power behaviour of four goodness-of-fit test statistics in sparse multinomials with k cells. Most previous work has been concerned only with both Pearson's X² and the likelihood ratio test statistics. We consider in this study, two additional test statistics, namely, the Cressie-Read test statistic – I(2/3) and the modified Freeman-Tukey test (FT) statistic. Because k ≥ 10 in this study, a Monte Carlo procedure based on 1000 simulated samples is used to estimate the powers for the four test statistics. Alternatives on various line segments are employed. Results suggest that none of the test statistics completely dominate the other and that the choice of which test to use depends on the nature of the alternative hypothesis. These results are consistent with those obtained by West and Kempthorne (1972), although, the Pearson's χ² test statistic may be preferred because of its closer approximation to the χ² distribution in terms of the attained α levels.     

Influence of Stimuli Spatial Proximity on a SSVEP-Based BCI Performance

Objective: Steady-State Visual Evoked Potentials based Brain-Computer Interfaces (SSVEP-based BCIs) systems have been shown as promising technology due to their short response time and ease of use. SSVEP-based BCIs use brain responses to a flickering visual stimulus as an input command to an external application or device, and it can be influenced by stimulus properties, signal recording, and signal processing. We aim to investigate the system performance varying the stimuli spatial proximity (a stimulus property).Material and methods: We performed a comparative analysis of two visual interface designs (named cross and square) for an SSVEP-based BCI. The power spectrum density (PSD) was used as feature extraction and the Support Machine Vector (SVM) as classification method. We also analyzed the effects of five flickering frequencies (6.67, 8.57, 10, 12 e 15 Hz) between and within interfaces.Results: We found higher accuracy rates for the flickering frequencies of 10, 12, and 15 Hz. The stimulus of 10 Hz presented the highest SSVEP amplitude response for both interfaces. The system presented the best performance (highest classification accuracy and information transfer rate) using the cross interface (lower visual angle).Conclusion: Our findings suggest that the system has the highest performance in the spatial proximity range from 4° to 13° (visual angle). In addition, we conclude that as the stimulus spatial proximity increases, the interference from other stimuli reduces, and the SSVEP amplitude response decreases, which reduces system accuracy. The inter-stimulus distance is a visual interface parameter that must be chosen carefully to increase the efficiency of an SSVEP-based BCI.     

Non-linear tests for identifying differentially expressed genes or genetic networks

MOTIVATION: One of the recently developed statistics for identifying differentially expressed genetic networks is Hotelling T2 statistic, which is a quadratic form of difference in linear functions of means of gene expressions between two types of tissue samples, and so their power is limited. RESULTS: To improve the power of test statistics, a general statistical framework for construction of non-linear tests is presented, and two specific non-linear test statistics that use non-linear transformations of means are developed. Asymptotical distributions of the non-linear test statistics under the null and alternative hypothesis are derived. It has been proved that under some conditions the power of the non-linear test statistics is higher than that of the T2 statistic. Besides theory, to evaluate in practice the performance of the non-linear test statistics, they are applied to two real datasets. The preliminary results demonstrate that the P-values of the non-linear statistics for testing differential expressions of the genetic networks are much smaller than those of the T2 statistic. And furthermore simulations show the Type I errors of the non-linear statistics agree with the threshold used and the statistics fit the chi2 distribution. SUPPLEMENTARY INFORMATION: Supplementary data are available on Bioinformatics online.     

Mean airway pressure and mean alveolar pressure during high-frequency jet ventilation in rabbits

Mean airway pressure underestimates mean alveolar pressure during high-frequency oscillatory ventilation. We hypothesized that high inspiratory flows characteristic of high-frequency jet ventilation may generate greater inspiratory than expiratory pressure losses in the airways, thereby causing mean airway pressure to overestimate, rather than underestimate, mean alveolar pressure. To test this hypothesis, we ventilated anesthetized paralyzed rabbits with a jet ventilator at frequencies of 5, 10, and 15 Hz, constant inspiratory-to-expiratory time ratio of 0.5 and mean airway pressures of 5 and 10 cmH2O. We measured mean total airway pressure in the trachea with a modified Pitot probe, and we estimated mean alveolar pressure as the mean pressure corresponding in the static pressure-volume relationship to the mean volume of the respiratory system measured with a jacket plethysmograph. We found that mean airway pressure was similar to mean alveolar pressure at frequencies of 5 and 10 Hz but overestimated it by 1.1 and 1.4 cmH2O at mean airway pressures of 5 and 10 cmH2O, respectively, when frequency was increased to 15 Hz. We attribute this finding primarily to the combined effect of nonlinear pressure frictional losses in the airways and higher inspiratory than expiratory flows. Despite the nonlinearity of the pressure-flow relationship, inspiratory and expiratory net pressure losses decreased with respect to mean inspiratory and expiratory flows at the higher rates, suggesting rate dependence of flow distribution. Redistribution of tidal volume to a shunt airway compliance is thought to occur at high frequencies.(ABSTRACT TRUNCATED AT 250 WORDS)     

Renaturation effects of cis and trans platinum II and IV compounds on calf thymus deoxyribonucleic acid.

The effects of cis dichlorodiammine platinum [cis Pt(II)], trans dichlorodiammine platinum (trans Pt(II)], cis tetrachlorodiammine platinum [cis Pt(IV)], trans tetrachlorodiammine platinum [trans Pt(IV)], and ethylenediaminedichloride platinum [Pt(II)en] on the absorption spectra, and thermal hyper- and hypochromicity of calf thymus DNA were investigated. Platinum-induced renaturation was studied as one parameter of interstrand cross-linking. Based on a DNA cross-linking hypothesis, the tumor-inhibitory platinum compounds cis Pt(II), cis Pt(IV) and Pt(II)en would be expected to induce renaturation following thermal denaturation, whereas the ineffective drugs, trans Pt(II) and trans Pt(IV) would not. All five bind to DNA in such a way as to induce renaturation. However, cis Pt(IV) requires at least a 3- to 4-fold longer incubation time than is required by the other compounds to form the coordination bonds necessary for renaturation. Maximum renaturation with all compounds was observed at a molar Pt/base ratio of 0.05 except cis Pt(IV), with which it was 0.25. The rate of the formation of the platinum-coordinated cross-links by fresh cis Pt(II) suggests two reactions or types of reactions occur. The first is rapid and destabilizes the DNA helix, whereas the second is slow and responsible for renaturation following thermal denaturation. These results suggest that cis Pt(IV) may be activated cellularly and that cross-linking is not the primary mechanism of action of the tumor-inhibitory platinum compounds.     

Influence of pulsed electromagnetic and pulsed vector magnetic potential field on the growth of tumor cells

Aims and Background: Tumor diseases cause 20% of deaths in Europe and they are the second most common cause of death and morbidity after cardiovascular diseases. Thus, tumor cells are target of many therapeutic strategies and tumor research is focused on searching more efficient and specific drugs as well as new therapeutic approaches. One of the areas of tumor research is an issue of external fields. In our work, we tested influence of a pulsed electromagnetic field (PEMF) and a hypothetic field of the pulsed vector magnetic potential (PVMP) on the growth of tumor cells; and further the possible growth inhibition effect of the PVMP. Methods: Both unipolar and bipolar PEMF fields of 5?mT and PVMP fields of 0?mT at frequencies of 15?Hz, 125?Hz and 625?Hz were tested on cancer cell lines derived from various types of tumors: CEM/C2 (acute lymphoblastic leukemia), SU-DHL-4 (B-cell lymphoma), COLO-320DM (colorectal adenocarcinoma), MDA-BM-468 (breast adenocarcinoma), and ZR-75-1 (ductal carcinoma). Cell morphology was observed, proliferation activity using WST assay was measured and simultaneous proportion of live, early apoptotic and dead cells was detected using flow cytometry. Results: A PEMF of 125?Hz and 625?Hz for 24?h–48?h increased proliferation activity in the 2 types of cancer cell lines used, i.e. COLO-320DM and ZR-75-1. In contrast, any of employed methods did not confirm a significant inhibitory effect of hypothetic PVMP field on tumor cells.     

The effects of motor unit synchronization on the power spectrum of the electromyogram

A realistic model for two synchronized motor unit action potential trains (MUAPT) is presented in which the variability of the time difference between corresponding action potentials (hereafter denoted by delay) is taken into account. Specifically, this delay is modeled as a continuous random variable that may assume both positive and negative values.Expressions are derived for the auto- and cross-power spectra of two such trains using their relations with the auto- and cross-correlation functions, respectively, with which they form Fourier transform pairs.The results show that the auto- and the cross-power spectra of two such synchronized MUAPTs differ from the auto- and the cross-spectra of two independent MUAPTs. The contribution of the statistics of the interpulse intervals to one of the autopower spectra is smaller and the cross-power spectra no longer reduce to a Dirac -function at the origin but are now determined by the other auto-power spectrum and by the Fourier transform of the density function associated with the time difference between corresponding action potentials. As a consequence of this change in the cross-power spectra synchronization leads to an absolute increase of power at low frequencies and to a relative decrease of power at high frequencies.The results are then generalized to electromyograms (EMG) composed of more than just two MUAPTs and illustrated with simulated power spectra with which the theory shows excellent agreement.     

Comparison of maximum statistics for hypothesis testing when a nuisance parameter is present only under the alternative

In many practical problems, a hypothesis testing involves a nuisance parameter which appears only under the alternative hypothesis. Davies (1977, Biometrika 64, 247-254) proposed the maximum of the score statistics over the whole range of the nuisance parameter as a test statistic for this type of hypothesis testing. Freidlin, Podgor, and Gastwirth (1999, Biometrics 55, 883-886) studied two other simpler maximum test statistics, the maximum of the score statistics at two extreme points of the nuisance parameter, and the maximum of the score statistics at three points of the nuisance parameter including the two extreme points. In this article, we compare the powers of these three maximum-type statistics in the context of three genetic problems.     

The F1-F2 vowel chart for Czech whispered vowels a, e, i, o, u

Aim: The aim of this contribution is to present the formant chart of the Czech vowels a, e, i, o, u and show that this can be achieved by means of digital methods of sound processing. Method: A group of 35 Czech students of the Pedagogical Faculty of Palacky University was tested and a record of whispered vowels was taken from each of them. The record was digitalized and processed by the Discrete Fourier Trasform. The result is the power spectrum of the individual vocals - the graphic output consists of a plot of the relative power of individual frequencies in the original sound. The values of the first two maxima which represent the first and the second formants were determined from the graph. The values were plotted on a formant chart. Results: Altogether, 175 spectral analyses of individual vowels were performed. In the resulting power spectrum, the first and the second formant frequencies were identified. The first formant was plotted against the second one and pure vocal formant regions were identified. Conclusion: Frequency bands for the Czech vowel "a" were circumscribed between 850 and 1150 Hz for first formant (F1) and between 1200 and 2000 Hz for second formant (F2). Similarly, borders of frequency band for vowel "e" they were 700 and 950 Hz for F1 and 1700 and 3000 Hz for F2. For vowel "i" 300 and 450 Hz for F1 and 2000 and 3600 Hz for F2, for vowel "o" 600 and 800 Hz for F1 and 600 and 1400 Hz for F2, for vowel "u" 100 and 400 Hz for F1 and 400 and 1200 Hz for F2. Discussion: At low frequencies it is feasible to invoke the source-filter model of voice production and associate vowel identity with frequencies of the first two formants in the voice spectrum. On the other hand, subject to intonation, singing or other forms of exposed voice (such as emotional speech, focused speech), the formant regions tend to spread. In spectral analysis other frequencies dominate, so specific formant frequency bands are not easily recognizable. Although the resulting formant map is not much different from the formant map of Peterson, it carries basic information about specific Czech vowels. The results may be used in further research and in education.     

Photoacoustic spectra of needles as an indicator of the activity of the photosynthetic apparatus of healthy and damaged conifers

Photoacoustic spectra were measured at chopping frequencies of the excitation light at 22 and 515 Hz in order to detect differences between needles from healthy and damaged conifers ( Abies alba Mill, and Picea abies Karst.) with respect to pigment distribution and possibly also to photosynthetic activity. Fully green needles of healthy trees exhibit photoacoustic spectra at 22 Hz with a maximum in the red absorption region of the chlorophyll. This maximum is lost with increasing damage to needles and chlorophyll breakdown. The photoacoustic spectra at 22 Hz of the damaged needles therefore are characterized by a higher signal in the blue-light region as compared to that in the red-light region. This can be quantified by forming the ratio of the photoacoustic signals at 675 and 475 nm (ratio PA 675/475). The needles of the damaged trees possess a lower photosynthetic activity, as seen from the CO₂ fixation rate and the variable fluorescence (Rfd-values). It is assumed that the changes in the PA-spectra of the needles from damaged trees are the result of this decline of the photosynthetic apparatus and its function. In contrast to the PA-spectrum at 22 Hz the PA-spectra at 515 Hz, where the PA-signals primarily emanate from the epidermal layer, exhibit a different shape with a maximum near 550 nm. In the needles from the damaged trees the PA-signals are higher, particularly in this maximum range of the spectrum, than those of fully green needles from healthy trees.     

Similar Spectral Power Densities Within the Schumann Resonance and a Large Population of Quantitative Electroencephalographic Profiles: Supportive Evidence for Koenig and Pobachenko

In 1954 and 1960 Koenig and his colleagues described the remarkable similarities of spectral power density profiles and patterns between the earth-ionosphere resonance and human brain activity which also share magnitudes for both electric field (mV/m) and magnetic field (pT) components. In 2006 Pobachenko and colleagues reported real time coherence between variations in the Schumann and brain activity spectra within the 6–16 Hz band for a small sample. We examined the ratios of the average potential differences (~3 μV) obtained by whole brain quantitative electroencephalography (QEEG) between rostral-caudal and left-right (hemispheric) comparisons of 238 measurements from 184 individuals over a 3.5 year period. Spectral densities for the rostral-caudal axis revealed a powerful peak at 10.25 Hz while the left-right peak was 1.95 Hz with beat-differences of ~7.5 to 8 Hz. When global cerebral measures were employed, the first (7–8 Hz), second (13–14 Hz) and third (19–20 Hz) harmonics of the Schumann resonances were discernable in averaged QEEG profiles in some but not all participants. The intensity of the endogenous Schumann resonance was related to the ‘best-of-fitness’ of the traditional 4-class microstate model. Additional measurements demonstrated real-time coherence for durations approximating microstates in spectral power density variations between Schumann frequencies measured in Sudbury, Canada and Cumiana, Italy with the QEEGs of local subjects. Our results confirm the measurements reported by earlier researchers that demonstrated unexpected similarities in the spectral patterns and strengths of electromagnetic fields generated by the human brain and the earth-ionospheric cavity.