Testing for group effect in a 2 x k heteroscedastic ANOVA model

It is natural to want to relax the assumption of homoscedasticity and Gaussian error in ANOVA models. For a two-way ANOVA model with 2 x k cells, one can derive tests of main effect for the factor with two levels (referred to as group) without assuming homoscedasticity or Gaussian error. Empirical likelihood can be used to derive testing procedures. An approximate empirical likelihood ratio test (AELRT) is derived for the test of group main effect. To approximate the distributions of the test statistics under the null hypothesis, simulation from the approximate empirical maximum likelihood estimate (AEMLE) restricted by the null hypothesis is used. The homoscedastic ANOVA F -test and a Box-type approximation to the distribution of the heteroscedastic ANOVA F -test are compared to the AELRT in level and power. The AELRT procedure is shown by simulation to have appropriate type I error control (although possibly conservative) when the distribution of the test statistics are approximated by simulation from the constrained AEMLE. The methodology is motivated and illustrated by an analysis of folate levels in the blood among two alcohol intake groups while accounting for gender.     

A general approach for sample size and power calculations based on the Haseman-Elston method

Risch and Zhang (1995; Science 268: 1584-9) reported a simple sample size and power calculation approach for the Haseman-Elston method and based their computations on the null hypothesis of no genetic effect. We argue that the more reasonable null hypothesis is that of no recombination. For this null hypothesis, we provide a general approach for sample size and power calculations within the Haseman-Elston framework. We demonstrate the validity of our approach in a Monte-Carlo simulation study and illustrate the differences using data from published segregation analyses on body weight and heritability estimates on carotid artery artherosclerotic lesions.     

EMG median frequency changes in the neck–shoulder stabilizers of symptomatic office workers when challenged by different physical stressors

The problem of work-related neck and upper limb disorders among computer users has been reported extensively in the literature, and commonly cited risk factors include static posture, speed and force of keyboard operation. The present study examined changes in median frequency (MF) of the neck–shoulder muscles in symptomatic and asymptomatic office workers when they were exposed to these three physical stressors.

A quasi-experimental Case–Control design was used to examine MF changes in two groups of female office workers when they were subjected to controlled doses of computer work involving prolonged static posture, increased typing speed and increased typing force. The MF of four major neck–shoulder muscles were examined bilaterally and compared between groups.

The MF changes over time-at-task did not clearly illustrate any muscle fatigue mechanism. However, Case Group consistently showed trends for higher MF than the Control Group, and this pattern was observed in response to all three physical stressors. The consistent group differences in MF suggest different muscle recruitment strategies between symptomatic and asymptomatic office workers. These results implied that symptomatic individuals had altered motor control, which may have important implications in understanding the etiology of work-related musculoskeletal disorders.     

基于非线性方法的DNA序列分析

以人类DNA序列为原始数据,对其进行数字编码。在利用功率谱分析编码DNA序列的基础上,利用Hurst指数进一步分析序列的自相似性。从分析结果看出,DNA序列中的确存在长程相关现象,而且这种长程相关现象与DNA的组成结构有关,表现为它的结构基团的Hurst指数大于其功能基团的Hurst指数。同时,从内含子含量的角度分析序列的长程相关程度,结果表明,同一序列中,内含子含量与不同的化学基团具有不同的关系。这表明不同的化学结构对DNA序列特征具有不同的贡献。这些结论从非线性方法的角度对DNA序列的分析提供了新的思路。     

长江胭脂鱼的洄游问题及水利工程对其资源的影响

胭脂鱼为我国特有鱼类,自然分布地仅在长江和闽江。长江上游为胭脂鱼的主要分布区,孵化出的部分幼鱼可能会随水流漂至长江中下游及其附属水体生长发育,至性成熟前返回上游参加繁殖,其后留在上游生活。目前长江上游和葛洲坝以下江段已存在自然繁殖的不同的胭脂鱼种群;葛洲坝和三峡水利枢纽的建设不会直接威胁胭脂鱼物种的生存和繁衍。     

Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits

Understanding the genetic architecture of quantitative traits can greatly assist the design of strategies for their manipulation in plant-breeding programs. For a number of traits, genetic variation can be the result of segregation of a few major genes and many polygenes (minor genes). The joint segregation analysis (JSA) is a maximum-likelihood approach for fitting segregation models through the simultaneous use of phenotypic information from multiple generations. Our objective in this paper was to use computer simulation to quantify the power of the JSA method for testing the mixed-inheritance model for quantitative traits when it was applied to the six basic generations: both parents (P₁ and P₂), F₁, F₂, and both backcross generations (B₁ and B₂) derived from crossing the F₁ to each parent. A total of 1968 genetic model-experiment scenarios were considered in the simulation study to quantify the power of the method. Factors that interacted to influence the power of the JSA method to correctly detect genetic models were: (1) whether there were one or two major genes in combination with polygenes, (2) the heritability of the major genes and polygenes, (3) the level of dispersion of the major genes and polygenes between the two parents, and (4) the number of individuals examined in each generation (population size). The greatest levels of power were observed for the genetic models defined with simple inheritance; e.g., the power was greater than 90% for the one major gene model, regardless of the population size and major-gene heritability. Lower levels of power were observed for the genetic models with complex inheritance (major genes and polygenes), low heritability, small population sizes and a large dispersion of favourable genes among the two parents; e.g., the power was less than 5% for the two major-gene model with a heritability value of 0.3 and population sizes of 100 individuals. The JSA methodology was then applied to a previously studied sorghum data-set to investigate the genetic control of the putative drought resistance-trait osmotic adjustment in three crosses. The previous study concluded that there were two major genes segregating for osmotic adjustment in the three crosses. Application of the JSA method resulted in a change in the proposed genetic model. The presence of the two major genes was confirmed with the addition of an unspecified number of polygenes. Received: 18 August 2000 / Accepted: 9 March 2001     

Comparative chlorine and temperature tolerance of the oyster crassostrea madrasensis: Implications for cooling system fouling

Crassostrea madrasensis is an important fouling oyster in tropical industrial cooling water systems. C. madrasensis individuals attach to surfaces by cementing one of their two valves to the substratum. Therefore, oyster fouling creates more problems than mussel fouling in the cooling conduits of power stations, because unlike the latter, the shell of the former remains attached to the substratum even after the death of the animal. However, there are no published reports on the tolerance of this species to chlorination and heat treatment. The mortality pattern and physiological behaviour (oxygen consumption and filtration rate) of three size groups (13 mm, 44 mm and 64 mm mean shell length) of C. madrasensis were studied at different residual chlorine concentrations (0.25, 0.5, 0.75, 1, 2, 3 to 5 mg l^{m 1}) and temperatures (30°C to 45°C). The effect of shell size (=age) on C. madrasensis mortality in the presence of chlorine and taking into account temperature was significant, with the largest size group oysters showing highest resistance. At 1 mg l^{m 1} residual chlorine, the 13 mm and 64 mm size group oysters took 504 h (21 d) and 744 h (31 d), respectively to reach 100% mortality. At 39°C, the 13 mm size group oysters took 218 min to reach 100% mortality, whereas the 64 mm size group oysters took 325 min. The oxygen consumption and filtration rate of C. madrasensis showed progressive reduction with increasing residual chlorine concentrations. However, the filtration rate and oxygen consumption responses of C. madrasensis were not significantly different between 30°C (control) and 37.5°C. There was a sharp decrease in the filtration rate and oxygen consumption at 38.5°C. A comparison of the present mortality data with previous reports on other bivalves suggests that the chlorine tolerance of C. madrasensis lies in between that of Perna viridis and Perna perna, while its temperature tolerance is significantly higher than that of the other two bivalve species. However, in power station heat exchangers, where simultaneous chlorine and thermal stresses are existent, C. madrasensis may have an edge over other common foulants, because of its high temperature tolerance.     

Residual force enhancement contributes to increased performance during stretch-shortening cycles of human plantar flexor muscles in vivo

It is well known that muscular force production is history-dependent, which results in enhanced (RFE) and depressed (RFD) steady-state forces after stretching and shortening, respectively. However, it remains unclear if force-enhancing mechanisms can contribute to increased performance during in vivo stretch-shortening cycles (SSCs) of human locomotor muscles. The purpose of this study was to investigate whether RFE-related mechanisms contribute to enhanced force and power output during SSCs of the human plantar flexor muscles. Net ankle torques of fourteen participants were measured during and after pure isometric, pure stretch, pure shortening, and SSC contractions when the triceps surae muscles were electrically stimulated at a submaximal level that resulted in 30% of their maximum isometric torque. Dynamic contractions were performed over an amplitude of 15°, from 5° plantar flexion to 10° dorsiflexion, at a speed of 120° s⁻¹. External ankle work during shortening was 11.6% greater during SSCs compared to pure shortening contractions (p = .003). Additionally, RFD after SSCs (8.6%) was reduced compared to RFD after pure shortening contractions (12.0%; p < .05). It is therefore concluded that RFE-related mechanisms contribute to increased performance following SSCs of human locomotor muscles. Since RFD after SSCs decreased although work during shortening was increased, we speculate that the relevant mechanism lies outside actin-myosin interaction. Finally, our data suggests that RFE might be relevant and beneficial for human locomotion whenever a muscle is stretched, but this needs to be confirmed.