A new method for the estimation of variance matrix with prescribed zeros in nonlinear mixed effects models

We propose a new method for the Maximum Likelihood Estimator (MLE) of nonlinear mixed effects models when the variance matrix of Gaussian random effects has a prescribed pattern of zeros (PPZ). The method consists in coupling the recently developed Iterative Conditional Fitting (ICF) algorithm with the Expectation Maximization (EM) algorithm. It provides positive definite estimates for any sample size, and does not rely on any structural assumption on the PPZ. It can be easily adapted to many versions of EM.Comment: Accepted for publication in Statistics and Computin

Environmental effects on progesterone profile measures of dairy cow fertility

Environmental effects on fertility measures early in lactation, such as the interval from calving to first luteal activity (CLA), proportion of samples with luteal activity during the first 60 days after calving (PLA) and interval to first ovulatory oestrus (OOE) were studied. In addition, traditional measurements of fertility, such as pregnancy to first insemination, number of inseminations per service period and interval from first to last insemination were studied as well as associations between the early and late measurements. Data were collected from an experimental herd during 15 years and included 1106 post-partum periods from 191 Swedish Holsteins and 325 Swedish Red and White dairy cows. Individual milk progesterone samples were taken twice a week until cyclicity and thereafter less frequently. First parity cows had 14.8 and 18.1 days longer CLA (LS-means difference) than second parity cows and older cows, respectively. Moreover, CLA was 10.5 days longer for cows that calved during the winter season compared with the summer season and 7.5 days longer for cows in tie-stalls than cows in loose-housing system. Cows treated for mastitis and lameness had 8.4 and 18.0 days longer CLA, respectively, compared with healthy cows. OOE was affected in the same way as CLA by the different environmental factors. PLA was a good indicator of CLA, and there was a high correlation (−0.69) between these two measurements. Treatment for lameness had a significant influence on all late fertility measurements, whereas housing was significant only for pregnancy to first insemination. All fertility traits were unfavourably associated with increased milk production. Regression of late fertility measurements on early fertility measurements had only a minor association with conception at first AI and interval from first to last AI for cows with conventional calving intervals, i.e. a 22 days later, CLA increased the interval from first to last insemination by 3.4 days. Early measurements had repeatabilities of 0.14–0.16, indicating a higher influence by the cow itself compared with late measurements, which had repeatabilities of 0.09–0.10. Our study shows that early fertility measurements have a possibility to be used in breeding for better fertility. To improve the early fertility of the cow, there are a number of important factors that have to be taken into account

A stitch in time: Efficient computation of genomic DNA melting bubbles

Background: It is of biological interest to make genome-wide predictions of the locations of DNA melting bubbles using statistical mechanics models. Computationally, this poses the challenge that a generic search through all combinations of bubble starts and ends is quadratic. Results: An efficient algorithm is described, which shows that the time complexity of the task is O(NlogN) rather than quadratic. The algorithm exploits that bubble lengths may be limited, but without a prior assumption of a maximal bubble length. No approximations, such as windowing, have been introduced to reduce the time complexity. More than just finding the bubbles, the algorithm produces a stitch profile, which is a probabilistic graphical model of bubbles and helical regions. The algorithm applies a probability peak finding method based on a hierarchical analysis of the energy barriers in the Poland-Scheraga model. Conclusions: Exact and fast computation of genomic stitch profiles is thus feasible. Sequences of several megabases have been computed, only limited by computer memory. Possible applications are the genome-wide comparisons of bubbles with promotors, TSS, viral integration sites, and other melting-related regions.Comment: 16 pages, 10 figure

Clustering and meso-level variables in cross-sectional surveys: an example of food aid during the Bosnian crisis

<p>Abstract</p> <p>Background</p> <p>Focus groups, rapid assessment procedures, key informant interviews and institutional reviews of local health services provide valuable insights on health service resources and performance. A long-standing challenge of health planning is to combine this sort of qualitative evidence in a unified analysis with quantitative evidence from household surveys. A particular challenge in this regard is to take account of the neighbourhood or clustering effects, recognising that these can be informative or incidental.</p> <p>Methods</p> <p>An example of food aid and food sufficiency from the Bosnian emergency (1995-96) illustrates two Lamothe cluster-adjustments of the Mantel Haenszel (MH) procedure, one assuming a fixed odds ratio and the other allowing for informative clustering by not assuming a fixed odds ratio. We compared these with conventional generalised estimating equations and a generalised linear mixed (GLMM) model, using a Laplace adjustment.</p> <p>Results</p> <p>The MH adjustment assuming incidental clustering generated a final model very similar to GEE. The adjustment that does not assume a fixed odds ratio produced a final multivariate model and effect sizes very similar to GLMM.</p> <p>Discussion</p> <p>In medium or large data sets with stratified last stage random sampling, the cluster adjusted MH is substantially more conservative than the naïve MH computation. In the example of food aid in the Bosnian crisis, the cluster adjusted MH that does not assume a fixed odds ratio produced similar results to the GLMM, which identified informative clustering.</p

Surface Covering of Downed Logs: Drivers of a Neglected Process in Dead Wood Ecology

Many species use coarse woody debris (CWD) and are disadvantaged by the forestry-induced loss of this resource. A neglected process affecting CWD is the covering of the surfaces of downed logs caused by sinking into the ground (increasing soil contact, mostly covering the underside of the log), and dense overgrowth by ground vegetation. Such cover is likely to profoundly influence the quality and accessibility of CWD for wood-inhabiting organisms, but the factors affecting covering are largely unknown. In a five-year experiment we determined predictors of covering rate of fresh logs in boreal forests and clear-cuts. Logs with branches were little covered because they had low longitudinal ground contact. For branchless logs, longitudinal ground contact was most strongly related to estimated peat depth (positive relation). The strongest predictor for total cover of branchless logs was longitudinal ground contact. To evaluate the effect on cover of factors other than longitudinal ground contact, we separately analyzed data from only those log sections that were in contact with the ground. Four factors were prominent predictors of percentage cover of such log sections: estimated peat depth, canopy shade (both increasing cover), potential solar radiation calculated from slope and slope aspect, and diameter of the log (both reducing cover). Peat increased cover directly through its low resistance, which allowed logs to sink and soil contact to increase. High moisture and low temperatures in pole-ward facing slopes and under a canopy favor peat formation through lowered decomposition and enhanced growth of peat-forming mosses, which also proved to rapidly overgrow logs. We found that in some boreal forests, peat and fast-growing mosses can rapidly cover logs lying on the ground. When actively introducing CWD for conservation purposes, we recommend that such rapid covering is avoided, thereby most likely improving the CWD's longevity as habitat for many species

Maternal Influences on the Transmission of Leukocyte Gene Expression Profiles in Population Samples from Brisbane, Australia

Two gene expression profiling studies designed to identify maternal influences on development of the neonate immune system and to address the population structure of the leukocyte transcriptome were carried out in Brisbane, Australia. In the first study, a comparison of 19 leukocyte samples obtained from mothers in the last three weeks of pregnancy with 37 umbilical cord blood samples documented differential expression of 7,382 probes at a false discovery rate of 1%, representing approximately half of the expressed transcriptome. An even larger component of the variation involving 8,432 probes, notably enriched for Vitamin E and methotrexate-responsive genes, distinguished two sets of individuals, with perfect transmission of the two profile types between each of 16 mother-child pairs in the study. A minor profile of variation was found to distinguish the gene expression profiles of obese mothers and children of gestational diabetic mothers from those of children born to obese mothers. The second study was of adult leukocyte profiles from a cross-section of Red Cross blood donors sampled throughout Brisbane. The first two axes in this study are related to the third and fourth axes of variation in the first study and also reflect variation in the abundance of CD4 and CD8 transcripts. One of the profiles associated with the third axis is largely excluded from samples from the central portion of the city. Despite enrichment of insulin signaling and aspects of central metabolism among the differentially expressed genes, there was little correlation between leukocyte expression profiles and body mass index overall. Our data is consistent with the notion that maternal health and cytokine milieu directly impact gene expression in fetal tissues, but that there is likely to be a complex interplay between cultural, genetic, and other environmental factors in the programming of gene expression in leukocytes of newborn children

Long term productivity and collaboration in information science

This is an accepted manuscript of an article published by Springer in Scientometrics on 02/07/2016, available online: https://doi.org/10.1007/s11192-016-2061-8 The accepted version of the publication may differ from the final published version.Funding bodies have tended to encourage collaborative research because it is generally more highly cited than sole author research. But higher mean citation for collaborative articles does not imply collaborative researchers are in general more research productive. This article assesses the extent to which research productivity varies with the number of collaborative partners for long term researchers within three Web of Science subject areas: Information Science & Library Science, Communication and Medical Informatics. When using the whole number counting system, researchers who worked in groups of 2 or 3 were generally the most productive, in terms of producing the most papers and citations. However, when using fractional counting, researchers who worked in groups of 1 or 2 were generally the most productive. The findings need to be interpreted cautiously, however, because authors that produce few academic articles within a field may publish in other fields or leave academia and contribute to society in other ways