Estimating equations with nonignorably missing response data

Troxel, Lipsitz, and Brennan (1997, Biometrics 53, 857-869) considered parameter estimation from survey data with nonignorable nonresponse and proposed weighted estimating equations to remove the biases in the complete-case analysis that ignores missing observations. This paper suggests two alternative modifications for unbiased estimation of regression parameters when a binary outcome is potentially observed at successive time points. The weighting approach of Robins, Rotnitzky, and Zhao (1995, Journal of the American Statistical Association 90, 106-121) is also modified to obtain unbiased estimating functions. The suggested estimating functions are unbiased only when the missingness probability is correctly specified, and misspecification of the missingness model will result in biases in the estimates. Simulation studies are carried out to assess the performance of different methods when the covariate is binary or normal. For the simulation models used, the relative efficiency of the two new methods to the weighting methods is about 3.0 for the slope parameter and about 2.0 for the intercept parameter when the covariate is continuous and the missingness probability is correctly specified. All methods produce substantial biases in the estimates when the missingness model is misspecified or underspecified. Analysis of data from a medical survey illustrates the use and possible differences of these estimating functions.     

Photosynthetic Utilization of Radiant Energy by Temperate Lettuce Grown Under Natural Tropical Condition with Manipulation of Root-Zone Temperature

Photosynthetic utilization of radiant energy was studied by chlorophyll (Chl) fluorescence and maximum photosynthetic O₂ evolution (P _max) in temperate lettuce (Lactuca sativa L.) grown under natural tropical fluctuating ambient temperatures but with their roots exposed to two different root-zone temperatures (RZTs): a constant 20 °C-RZT (RZT₂₀) and a fluctuating ambient RZT (RZT_a) from 23 to 40 °C. On a sunny day, irrespective of RZT, F/_Fm [ratio of the variable to maximal fluorescence under irradiation (the maximal photosystem 2 quantum yield with actinic light)] decreased and non-photochemical quenching (NPQ) increased parallel to the increase of photosynthetic photon flux density (PPFD). However, RZT_a plants showed lower F/F_m and higher NPQ than RZT₂₀ plants. The electron transport rate (ETR) was much higher in RZT₂₀ plants than in RZT_a plants especially during moderately sunny days. There were no significant diurnal changes in P _max although these values of RZT₂₀ plants were much higher than those of RZT_a plants. On cloudy days, no significant diurnal changes in F/F_m and NPQ occurred, but F/F_m was higher and NPQ was lower in RZT₂₀ plants than in RZT_a plants. Diurnal changes in ETR were also observed in all plants while P _max values throughout the whole cloudy days in both RZT₂₀ and RZT_a plants were constant. Again, RZT₂₀ plants had much higher values of P _max than RZT_a plants. During RZT transfer period, all Chl fluorescence parameters measured at midday fluctuated with PPFD. Impact of RZT on these parameters was observed 2–3 d after RZT transfer. ETR and P _max measured with saturating PPFD in the laboratory did not vary with the fluctuating PPFD in the greenhouse but the effects of RZT on these two parameters were observed 3–4 d after RZT transfer. Thus RZT affects photosynthetic utilization of photon energy in temperate lettuce grown under natural tropical condition.     

Trends in butterfly species richness in response to the peninsular effect in South Korea

Aim The Korean peninsula is elongated in shape and is connected to the Asian continent on the north. The peninsular effect – a decline in species density or richness as a function of distance from the mainland base (towards the distal tip) of a peninsula – was evaluated for plants and animals in different peninsulas. The aims of the present study were to describe the pattern of butterfly species diversity and to determine what factors may be responsible for this pattern along the Korean peninsula. The distribution pattern of butterfly species in South Korea before and after the Korean War was also investigated. Location South Korea (34–38° N, 126–129° E). Methods Forty‐three quadrats, each 1/2° latitude by 1/2° longitude, and three data sets – butterfly distribution data from 1938 to 1950, butterfly distribution data from 1976 to 1999, and the combined data – were analysed. The influence of four variables – latitude, longitude, area and maximum altitude – on each quadrat was investigated using multiple regression analysis. Results and conclusion The analyses revealed a marked peninsular effect: there was a significant positive correlation between butterfly species richness and latitude. Additionally, habitat diversity, expressed as maximum altitude, was significantly correlated with butterfly species richness. I conclude that both the geographical orientation and habitat diversity contribute to butterfly species diversity across South Korea. Comparison of ranges between the older and recent data sets suggests that geographical distributions of several species are dramatically reduced in size. These species may be used for future conservation activities in South Korea.     

Quartet mapping and the extent of lateral transfer in bacterial genomes

Several recent analyses have used quartet-based methods to assess the congruence among phylogenies derived for large sets of genes from prokaryotic genomes. The principal conclusion from these studies is that lateral gene transfer (LGT) has blurred prokaryotic phylogenies to such a degree that the darwinian scheme of treelike evolution might be abandoned in favor of a net or web. Here, we focus on one of these methods, quartet mapping, and show that its application can lead to overestimation of the extent of inferred LGT in prokaryotes, particularly when applied to distantly related taxa.     

Protein-water interactions in ribonuclease A and angiogenin: a molecular dynamics study

It is known that water molecules play an important role in the biological functioning of proteins. The members of the ribonuclease A (RNase A) family of proteins, which are sequentially and structurally similar, are known to carry out the obligatory function of cleaving RNA and individually perform other diverse biological functions. Our focus is on elucidating whether the sequence and structural similarity lead to common hydration patterns, what the common hydration sites are and what the differences are. Extensive molecular dynamics simulations followed by a detailed analysis of protein-water interactions have been carried out on two members of the ribonuclease A superfamily-RNase A and angiogenin. The water residence times are analyzed and their relationship with the characteristic properties of the protein polar atoms, such as their accessible surface area and mean hydration, is studied. The capacity of the polar atoms to form hydrogen bonds with water molecules and participate in protein-water networks are investigated. The locations of such networks are identified for both proteins.     

Analysis of survival data with cross-effects of survival functions

We present a model and semiparametric estimation procedures for analysis of survival data with cross-effects (CE) of survival functions. Finite sample properties of the estimators are analyzed by simulation. A goodness-of-fit test for the proportional hazards model against the CE model is proposed. The well known data concerning effects of chemotherapy and radiotherapy on the survival times of gastric cancer patients is analyzed as an example.     

Evaluating the performance of likelihood methods for detecting population structure and migration

A plethora of statistical models have recently been developed to estimate components of population genetic history. Very few of these methods, however, have been adequately evaluated for their performance in accurately estimating population genetic parameters of interest. In this paper, we continue a research program of evaluation of population genetic methods through computer simulation. Specifically, we examine the software MIGRATEE-N 1.6.8 and test the accuracy of this software to estimate genetic diversity (Theta), migration rates, and confidence intervals. We simulated nucleotide sequence data under a neutral coalescent model with lengths of 500 bp and 1000 bp, and with three different per site Theta values of (0.00025, 0.0025, 0.025) crossed with four different migration rates (0.0000025, 0.025, 0.25, 2.5) to construct 1000 evolutionary trees per-combination per-sequence-length. We found that while MIGRATEE-N 1.6.8 performs reasonably well in estimating genetic diversity (Theta), it does poorly at estimating migration rates and the confidence intervals associated with them. We recommend researchers use this software with caution under conditions similar to those used in this evaluation.     

Phylogenetic and biogeographic diversification of Rhus (Anacardiaceae) in the Northern Hemisphere

Sequences of internal transcribed spacers (ITS) of nuclear ribosomal DNA, the chloroplast ndhF gene, and chloroplast trnL-F regions (trnL intron, and trnL [UAA] 3' exon-trnF [GAA] intergenic spacer) were used for phylogenetic analyses of Rhus, a genus disjunctly distributed in Asia, Europe, Hawaii, North America, and Northern Central America. Both ITS and cpDNA data sets support the monophyly of Rhus. The monophyly of subgenus Rhus was suggested by the combined cpDNA and ITS data, and largely supported in the cpDNA data except that Rhus microphylla of subgenus Lobadium was nested within it. The monophyly of subgenus Lobadium was strongly supported in the ITS data, whereas the cpDNA data revealed two main clades within the subgenus, which formed a trichotomy with the clade of subgenus Rhus plus R. microphylla. The ITS and cpDNA trees differ in the positions of Rhus michauxii, R. microphylla, and Rhus rubifolia, and hybridization may have caused this discordance. Fossil evidence indicates that Rhus dates back to the early Eocene. The penalized likelihood method was used to estimate divergence times, with fossils of Rhus subgenus Lobadium, Pistacia and Toxicodendron used for age constraints. Rhus diverged from its closest relative at 49.1+/-2.1 million years ago (Ma), the split of subgenus Lobadium and subgenus Rhus was at 38.1+/-3.0 Ma. Rhus most likely migrated from North America into Asia via the Bering Land Bridge during the Late Eocene (33.8+/-3.1 Ma). Rhus coriaria from southern Europe and western Asia diverged from its relatives in eastern Asia at 24.4+/-3.2 Ma. The Hawaiian Rhus sandwicensis diverged from the Asian Rhus chinensis at 13.5+/-3.0 Ma. Subgenus Lobadium was inferred to be of North American origin. Taxa of subgenus Lobadium then migrated southward to Central America. Furthermore, we herein make the following three nomenclatural combinations: (1) Searsia leptodictya (Diels) T. S. Yi, A. J. Miller and J. Wen, comb. nov., (2) Searsia pyroides (A. Rich.) T. S. Yi, A. J. Miller and J. Wen, comb. nov., and (3) Searsia undulata (Jacq.) T. S. Yi, A. J. Miller and J. Wen, because our analyses support the segregation of Searsia from Rhus.     

AFLP utility for population assignment studies: analytical investigation and empirical comparison with microsatellites

Individual-based population assignment tests have thus far mainly relied on the use of microsatellite loci. However, the logistic difficulty of screening large numbers of loci required to reach sufficient statistical power hampers the usefulness of microsatellites in situations of weak population structuring. Amplified fragment length polymorphisms (AFLP) represents an alternative for overcoming this logistical issue as the technique allows the user to characterize a much larger number of loci with a comparable analytical effort. In this study, an assignment test based on maximum likelihood for dominant markers was used to investigate the potential usefulness of AFLP for population assignment. We also compared assignment success achieved with AFLP with that obtained using microsatellites in a case study of low population differentiation involving whitefish (Coregonus clupeaformis) sympatric ecotypes. The analytical investigation showed that the minimum number of AFLP loci required to reach an assignment success of 95% stood within values that are easily achievable in many situations. This also showed how assignment success varied according to the number of AFLP loci used, their absolute frequency and their frequency differential and sampling errors, as well as the number of putative source populations. The case study showed that given a comparable analytical effort in the laboratory, AFLP were much more efficient than the microsatellite loci in discriminating the source of an individual among putative populations. AFLP resulted in higher assignment success at all levels of stringency and the log-likelihood differences between populations obtained with AFLP for each individual were much larger than those obtained with microsatellites. These results indicate that research involving individual-based population assignment methods should benefit importantly from the use of AFLP markers, especially in systems characterized by weak population structuring.     

Methods of parentage analysis in natural populations

The recent proliferation of hypervariable molecular markers has ushered in a surge of techniques for the analysis of parentage in natural and experimental populations. Consequently, the potential for meaningful studies of paternity and maternity is at an all-time high. However, the details and implementation of the multifarious techniques often differ in subtle ways that can influence the results of parentage analyses. Now is a good time to reflect on the available techniques and to consider their strengths and weaknesses. Here, we review the leading techniques in parentage analysis, with a particular emphasis on those that have been implemented in readily useable software packages. Our survey leads to some important insights with respect to the utility of the different approaches. This review should serve as a useful guide to anyone who wishes to embark on the study of parentage.