Block-Conditional Missing at Random Models for Missing Data

Two major ideas in the analysis of missing data are (a) the EM algorithm [Dempster, Laird and Rubin, J. Roy. Statist. Soc. Ser. B 39 (1977) 1--38] for maximum likelihood (ML) estimation, and (b) the formulation of models for the joint distribution of the data ${Z}$ and missing data indicators ${M}$, and associated "missing at random"; (MAR) condition under which a model for ${M}$ is unnecessary [Rubin, Biometrika 63 (1976) 581--592]. Most previous work has treated ${Z}$ and ${M}$ as single blocks, yielding selection or pattern-mixture models depending on how their joint distribution is factorized. This paper explores "block-sequential"; models that interleave subsets of the variables and their missing data indicators, and then make parameter restrictions based on assumptions in each block. These include models that are not MAR. We examine a subclass of block-sequential models we call block-conditional MAR (BCMAR) models, and an associated block-monotone reduced likelihood strategy that typically yields consistent estimates by selectively discarding some data. Alternatively, full ML estimation can often be achieved via the EM algorithm. We examine in some detail BCMAR models for the case of two multinomially distributed categorical variables, and a two block structure where the first block is categorical and the second block arises from a (possibly multivariate) exponential family distribution.Comment: Published in at http://dx.doi.org/10.1214/10-STS344 the Statistical Science (http://www.imstat.org/sts/) by the Institute of Mathematical Statistics (http://www.imstat.org

Non-Coding RNA Sequencing of Equine Endometrium During Maternal Recognition of Pregnancy.

Maternal recognition of pregnancy (MRP) in the mare is not well defined. In a non-pregnant mare, prostaglandin F2α (PGF) is released on day 14 post-ovulation (PO) to cause luteal regression, resulting in loss of progesterone production. Equine MRP occurs prior to day 14 to halt PGF production. Studies have failed to identify a gene candidate for MRP, so attention has turned to small, non-coding RNAs. The objective of this study was to evaluate small RNA (<200 nucleotides) content in endometrium during MRP. Mares were used in a cross-over design with each having a pregnant and non-mated cycle. Each mare was randomly assigned to collection day 11 or 13 PO (n = 3/day) and endometrial biopsies were obtained. Total RNA was isolated and sequencing libraries were prepared using a small RNA library preparation kit and sequenced on a HiSeq 2000. EquCab3 was used as the reference genome and DESeq2 was used for statistical analysis. On day 11, 419 ncRNAs, representing miRNA, snRNA, snoRNA, scaRNA, and vaultRNA, were different between pregnancy statuses, but none on day 13. Equine endometrial ncRNAs with unknown structure and function were also identified. This study is the first to describe ncRNA transcriptome in equine endometrium. Identifying targets of these ncRNAs could lead to determining MRP

Inconsistency of the MLE for the joint distribution of interval censored survival times and continuous marks

This paper considers the nonparametric maximum likelihood estimator (MLE) for the joint distribution function of an interval censored survival time and a continuous mark variable. We provide a new explicit formula for the MLE in this problem. We use this formula and the mark specific cumulative hazard function of Huang and Louis (1998) to obtain the almost sure limit of the MLE. This result leads to necessary and sufficient conditions for consistency of the MLE which imply that the MLE is inconsistent in general. We show that the inconsistency can be repaired by discretizing the marks. Our theoretical results are supported by simulations.Comment: 27 pages, 4 figure

Improved Experimental Limits on the Production of Magnetic Monopoles

We present new limits on low mass accelerator-produced point-like Dirac magnetic monopoles trapped and bound in matter surrounding the D\O collision region of the Tevatron at Fermilab (experiment E-882). In the context of a Drell-Yan mechanism, we obtain cross section limits for the production of monopoles with magnetic charge values of 1, 2, 3, and 6 times the minimum Dirac charge of the order of picobarns, some hundred times smaller than found in similar previous Fermilab searches. Mass limits inferred from these cross section limits are presented.Comment: 5 pages, 4 eps figures, REVTe

Revision rates after primary hip and knee replacement in England between 2003 and 2006

<b>Background</b>: Hip and knee replacement are some of the most frequently performed surgical procedures in the world. Resurfacing of the hip and unicondylar knee replacement are increasingly being used. There is relatively little evidence on their performance. To study performance of joint replacement in England, we investigated revision rates in the first 3 y after hip or knee replacement according to prosthesis type. <b>Methods and Findings</b>: We linked records of the National Joint Registry for England and Wales and the Hospital Episode Statistics for patients with a primary hip or knee replacement in the National Health Service in England between April 2003 and September 2006. Hospital Episode Statistics records of succeeding admissions were used to identify revisions for any reason. 76,576 patients with a primary hip replacement and 80,697 with a primary knee replacement were included (51% of all primary hip and knee replacements done in the English National Health Service). In hip patients, 3-y revision rates were 0.9% (95% confidence interval [CI] 0.8%–1.1%) with cemented, 2.0% (1.7%–2.3%) with cementless, 1.5% (1.1%–2.0% CI) with “hybrid” prostheses, and 2.6% (2.1%–3.1%) with hip resurfacing (p < 0.0001). Revision rates after hip resurfacing were increased especially in women. In knee patients, 3-y revision rates were 1.4% (1.2%–1.5% CI) with cemented, 1.5% (1.1%–2.1% CI) with cementless, and 2.8% (1.8%–4.5% CI) with unicondylar prostheses (p < 0.0001). Revision rates after knee replacement strongly decreased with age. <b>Interpretation</b>: Overall, about one in 75 patients needed a revision of their prosthesis within 3 y. On the basis of our data, consideration should be given to using hip resurfacing only in male patients and unicondylar knee replacement only in elderly patients

A Perspective on the Scientific Registry of Transplant Recipients' Migration to Bayesian Methods

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/112268/1/ajt13354.pd

Observation of a New Charmed Strange Meson

Using the CLEO-II detector, we have obtained evidence for a new meson decaying to $D^0 K^+$. Its mass is $2573.2^{+1.7}_{-1.6}\pm 0.8\pm 0.5$ {}~MeV/$c^2$ and its width is $16^{+5}_{-4}\pm 3$~MeV/$c^2$. Although we do not establish its spin and parity, the new meson is consistent with predictions for an $L=1$, $S=1$, $J_P=2^+$ charmed strange state.Comment: 9 pages uuencoded compressed postscript (process with uudecode then gunzip). hardcopies with figures can be obtained by sending mail to: [email protected]

Revealing three-dimensional structure of individual colloidal crystal grain by coherent x-ray diffractive imaging

We present results of a coherent x-ray diffractive imaging experiment performed on a single colloidal crystal grain. The full three-dimensional (3D) reciprocal space map measured by an azimuthal rotational scan contained several orders of Bragg reflections together with the coherent interference signal between them. Applying the iterative phase retrieval approach, the 3D structure of the crystal grain was reconstructed and positions of individual colloidal particles were resolved. As a result, an exact stacking sequence of hexagonal close-packed layers including planar and linear defects were identified.Comment: 8 pages, 5 figure

Precision Measurement of the Ds∗+−Ds+D_s^{*+}- D_s^+ Mass Difference

We have measured the vector-pseudoscalar mass splitting $M(D_s^{*+})-M(D_s^+) = 144.22\pm 0.47\pm 0.37 MeV$, significantly more precise than the previous world average. We minimize the systematic errors by also measuring the vector-pseudoscalar mass difference $M(D^{*0})-M(D^0)$ using the radiative decay $D^{*0}\rightarrow D^0\gamma$, obtaining $[M(D_s^{*+})-M(D_s^+)]-[M(D^{*0})-M(D^0)] = 2.09\pm 0.47\pm 0.37 MeV$. This is then combined with our previous high-precision measurement of $M(D^{*0})-M(D^0)$, which used the decay $D^{*0}\rightarrow D^0\pi^0$. We also measure the mass difference $M(D_s^+)-M(D^+)=99.5\pm 0.6\pm 0.3$ MeV, using the $\phi\pi^+$ decay modes of the $D_s^+$ and $D^+$ mesons.Comment: 18 pages uuencoded compressed postscript (process with uudecode then gunzip). hardcopies with figures can be obtained by sending mail to: [email protected]

Monitoring international migration flows in Europe. Towards a statistical data base combining data from different sources

The paper reviews techniques developed in demography, geography and statistics that are useful for bridging the gap between available data on international migration flows and the information required for policy making and research. The basic idea of the paper is as follows: to establish a coherent and consistent data base that contains sufficiently detailed, up-to-date and accurate information, data from several sources should be combined. That raises issues of definition and measurement, and of how to combine data from different origins properly. The issues may be tackled more easily if the statistics that are being compiled are viewed as different outcomes or manifestations of underlying stochastic processes governing migration. The link between the processes and their outcomes is described by models, the parameters of which must be estimated from the available data. That may be done within the context of socio-demographic accounting. The paper discusses the experience of the U.S. Bureau of the Census in combining migration data from several sources. It also summarizes the many efforts in Europe to establish a coherent and consistent data base on international migration. The paper was written at IIASA. It is part of the Migration Estimation Study, which is a collaborative IIASA-University of Groningen project, funded by the Netherlands Organization for Scientific Research (NWO). The project aims at developing techniques to obtain improved estimates of international migration flows by country of origin and country of destination