Estimation of dynamic latent variable models using simulated nonparametric moments

Abstract. Given a model that can be simulated, conditional moments at a trial parameter value can be calculated with high accuracy by applying kernel smoothing methods to a long simulation. With such conditional moments in hand, standard method of moments techniques can be used to estimate the parameter. Because conditional moments are calculated using kernel smoothing rather than simple averaging, it is not necessary that the model be simulable subject to the conditioning information that is used to define the moment conditions. For this reason, the proposed estimator is applicable to general dynamic latent variable models. It is shown that as the number of simulations diverges, the estimator is consistent and a higher-order expansion reveals the stochastic difference between the infeasible GMM estimator based on the same moment conditions and the simulated version. In particular, we show how to adjust standard errors to account for the simulations. Monte Carlo results show how the estimator may be applied to a range of dynamic latent variable (DLV) models, and that it performs well in comparison to several other estimators that have been proposed for DLV models.

On estimation and identifiability issues of sex-linked inheritance with a case study of pigmentation in Swiss barn owl (Tyto alba)

Genetic evaluation using animal models or pedigree-based models generally assume only autosomal inheritance. Bayesian animal models provide a flexible framework for genetic evaluation, and we show how the model readily can accommodate situations where the trait of interest is influenced by both autosomal and sex-linked inheritance. This allows for simultaneous calculation of autosomal and sex-chromosomal additive genetic effects. Inferences were performed using integrated nested Laplace approximations (INLA), a nonsampling-based Bayesian inference methodology. We provide a detailed description of how to calculate the inverse of the X- or Z-chromosomal additive genetic relationship matrix, needed for inference. The case study of eumelanic spot diameter in a Swiss barn owl (Tyto alba) population shows that this trait is substantially influenced by variation in genes on the Z-chromosome (sigma(2)(z) = 0.2719 and sigma(2)(a) = 0.4405). Further, a simulation study for this study system shows that the animal model accounting for both autosomal and sex-chromosome-linked inheritance is identifiable, that is, the two effects can be distinguished, and provides accurate inference on the variance components.

Estimation of the radiation risks in diagnostic procedures to populations and patients

Since 1895, when X-rays were discovered, ionizing radiation became part of our life. Its use in medicine has brought significant health benefits to the population globally. The benefit of any diagnostic procedure is to reduce the uncertainty about the patient's health. However, there are potential detrimental effects of radiation exposure. Therefore, radiation protection authorities have become strict regarding the control of radiation risks.¦There are various situations where the radiation risk needs to be evaluated. International authority bodies point to the increasing number of radiologic procedures and recommend population surveys. These surveys provide valuable data to public health authorities which helps them to prioritize and focus on patient groups in the population that are most highly exposed. On the other hand, physicians need to be aware of radiation risks from diagnostic procedures in order to justify and optimize the procedure and inform the patient.¦The aim of this work was to examine the different aspects of radiation protection and investigate a new method to estimate patient radiation risks.¦The first part of this work concerned radiation risk assessment from the regulatory authority point of view. To do so, a population dose survey was performed to evaluate the annual population exposure. This survey determined the contribution of different imaging modalities to the total collective dose as well as the annual effective dose per caput. It was revealed that although interventional procedures are not so frequent, they significantly contribute to the collective dose. Among the main results of this work, it was shown that interventional cardiological procedures are dose-intensive and therefore more attention should be paid to optimize the exposure.¦The second part of the project was related to the patient and physician oriented risk assessment. In this part, interventional cardiology procedures were studied by means of Monte Carlo simulations. The organ radiation doses as well as effective doses were estimated. Cancer incidence risks for different organs were calculated for different sex and age-at-exposure using the lifetime attributable risks provided by the Biological Effects of Ionizing Radiations Report VII. Advantages and disadvantages of the latter results were examined as an alternative method to estimate radiation risks. The results show that this method is the most accurate, currently available, to estimate radiation risks. The conclusions of this work may guide future studies in the field of radiation protection in medicine.¦-¦Depuis la découverte des rayons X en 1895, ce type de rayonnement a joué un rôle important dans de nombreux domaines. Son utilisation en médecine a bénéficié à la population mondiale puisque l'avantage d'un examen diagnostique est de réduire les incertitudes sur l'état de santé du patient. Cependant, leur utilisation peut conduire à l'apparition de cancers radio-induits. Par conséquent, les autorités sanitaires sont strictes quant au contrôle du risque radiologique.¦Le risque lié aux radiations doit être estimé dans différentes situations pratiques, dont l'utilisation médicale des rayons X. Les autorités internationales de radioprotection indiquent que le nombre d'examens et de procédures radiologiques augmente et elles recommandent des enquêtes visant à déterminer les doses de radiation délivrées à la population. Ces enquêtes assurent que les groupes de patients les plus à risque soient prioritaires. D'un autre côté, les médecins ont également besoin de connaître le risque lié aux radiations afin de justifier et optimiser les procédures et informer les patients.¦Le présent travail a pour objectif d'examiner les différents aspects de la radioprotection et de proposer une manière efficace pour estimer le risque radiologique au patient.¦Premièrement, le risque a été évalué du point de vue des autorités sanitaires. Une enquête nationale a été réalisée pour déterminer la contribution des différentes modalités radiologiques et des divers types d'examens à la dose efficace collective due à l'application médicale des rayons X. Bien que les procédures interventionnelles soient rares, elles contribuent de façon significative à la dose délivrée à la population. Parmi les principaux résultats de ce travail, il a été montré que les procédures de cardiologie interventionnelle délivrent des doses élevées et devraient donc être optimisées en priorité.¦La seconde approche concerne l'évaluation du risque du point de vue du patient et du médecin. Dans cette partie, des procédures interventionnelles cardiaques ont été étudiées au moyen de simulations Monte Carlo. La dose délivrée aux organes ainsi que la dose efficace ont été estimées. Les risques de développer des cancers dans plusieurs organes ont été calculés en fonction du sexe et de l'âge en utilisant la méthode établie dans Biological Effects of Ionizing Radiations Report VII. Les avantages et inconvénients de cette nouvelle technique ont été examinés et comparés à ceux de la dose efficace. Les résultats ont montré que cette méthode est la plus précise actuellement disponible pour estimer le risque lié aux radiations. Les conclusions de ce travail pourront guider de futures études dans le domaine de la radioprotection en médicine.

Estimation of sex and age of "virtual skeletons"--a feasibility study.

This article presents a feasibility study with the objective of investigating the potential of multi-detector computed tomography (MDCT) to estimate the bone age and sex of deceased persons. To obtain virtual skeletons, the bodies of 22 deceased persons with known age at death were scanned by MDCT using a special protocol that consisted of high-resolution imaging of the skull, shoulder girdle (including the upper half of the humeri), the symphysis pubis and the upper halves of the femora. Bone and soft-tissue reconstructions were performed in two and three dimensions. The resulting data were investigated by three anthropologists with different professional experience. Sex was determined by investigating three-dimensional models of the skull and pelvis. As a basic orientation for the age estimation, the complex method according to Nemeskéri and co-workers was applied. The final estimation was effected using additional parameters like the state of dentition, degeneration of the spine, etc., which where chosen individually by the three observers according to their experience. The results of the study show that the estimation of sex and age is possible by the use of MDCT. Virtual skeletons present an ideal collection for anthropological studies, because they are obtained in a non-invasive way and can be investigated ad infinitum.

Quantitative real-time polymerase chain reaction for detection of adenovirus after T cell-replete hematopoietic cell transplantation: viral load as a marker for invasive disease.

BACKGROUND: The value of adenovirus plasma DNA detection as an indicator for adenovirus disease is unknown in the context of T cell-replete hematopoietic cell transplantation, of which adenovirus disease is an uncommon but serious complication. METHODS: Three groups of 62 T cell-replete hematopoietic cell transplant recipients were selected and tested for adenovirus in plasma by polymerase chain reaction. RESULTS: Adenovirus was detected in 21 (87.5%) of 24 patients with proven adenovirus disease (group 1), in 4 (21%) of 19 patients who shed adenovirus (group 2), and in 1 (10.5%) of 19 uninfected control patients. The maximum viral load was significantly higher in group 1 (median maximum viral load, 6.3x10(6) copies/mL; range, 0 to 1.0x10(9) copies/mL) than in group 2 (median maximum viral load, 0 copies/mL; range, 0 to 1.7x10(8) copies/mL; P<.001) and in group 3 (median maximum viral load, 0 copies/mL; range 0-40 copies/mL; P<.001). All patients in group 2 who developed adenoviremia had symptoms compatible with adenovirus disease (i.e., possible disease). A minimal plasma viral load of 10(3) copies/mL was detected in all patients with proven or possible disease. Adenoviremia was detectable at a median of 19.5 days (range, 8-48 days) and 24 days (range, 9-41 days) before death for patients with proven and possible adenovirus disease, respectively. CONCLUSION: Sustained or high-level adenoviremia appears to be a specific and sensitive indicator of adenovirus disease after T cell-replete hematopoietic cell transplantation. In the context of low prevalence of adenovirus disease, the use of polymerase chain reaction of plasma specimens to detect virus might be a valuable tool to identify and treat patients at risk for viral invasive disease.

Abrupt decrease in serum testosterone levels after an oral glucose load in men: implications for screening for hypogonadism.

OBJECTIVE: This study examines the physiological impact of a glucose load on serum testosterone (T) levels in men with varying glucose tolerance (GT). DESIGN: Cross-sectional study. PATIENTS AND METHODS: 74 men (19-74 years, mean 51·4 ± 1·4 years) underwent a standard 75-g oral glucose tolerance test with blood sampling at 0, 30, 60, 90 and 120 min. Fasting serum glucose, insulin, total T (and calculated free T), LH, SHBG, leptin and cortisol were measured. RESULTS: 57% of the men had normal GT, 30% had impaired GT and 13% had newly diagnosed type 2 diabetes. Glucose ingestion was associated with a 25% decrease in mean T levels (delta = -4·2 ± 0·3 nm, P < 0·0001). T levels remained suppressed at 120 min compared with baseline (13·7 ± 0·6 vs 16·5 ± 0·7 nm, P < 0·0001) and did not differ across GT or BMI. Of the 66 men with normal T levels at baseline, 10 (15%) had levels that decreased to the hypogonadal range (<9·7 nm) at one or more time points. SHBG, LH and cortisol levels were unchanged. Leptin levels decreased from baseline at all time points (P < 0·0001). CONCLUSIONS: Glucose ingestion induces a significant reduction in total and free T levels in men, which is similar across the spectrum of glucose tolerance. This decrease in T appears to be because of a direct testicular defect, but the absence of compensatory changes in LH suggests an additional central component. Men found to have low nonfasting T levels should be re-evaluated in the fasting state.

Estimation of parametric and nonparametric models for univariate claim severity distributions - an approach using R

This paper presents an analysis of motor vehicle insurance claims relating to vehicle damage and to associated medical expenses. We use univariate severity distributions estimated with parametric and non-parametric methods. The methods are implemented using the statistical package R. Parametric analysis is limited to estimation of normal and lognormal distributions for each of the two claim types. The nonparametric analysis presented involves kernel density estimation. We illustrate the benefits of applying transformations to data prior to employing kernel based methods. We use a log-transformation and an optimal transformation amongst a class of transformations that produces symmetry in the data. The central aim of this paper is to provide educators with material that can be used in the classroom to teach statistical estimation methods, goodness of fit analysis and importantly statistical computing in the context of insurance and risk management. To this end, we have included in the Appendix of this paper all the R code that has been used in the analysis so that readers, both students and educators, can fully explore the techniques described

Estimation of Genetic Divergence and Gene Flow between Culex pipiens and Culex quinquefasciatus (Diptera: Culicidae) in Argentina

Allele frequencies at seven polymorphic loci controlling the synthesis of enzymes were analyzed in six populations of Culex pipiens L. and Cx. quinquefasciatus Say. Sampling sites were situated along a north-south line of about 2,000 km in Argentina. The predominant alleles at Mdh, Idh, Gpdh and Gpi loci presented similar frequencies in all the samples. Frequencies at the Pgm locus were similar for populations pairs sharing the same geographic area. The loci Cat and Hk-1 presented significant geographic variation. The latter showed a marked latitudinal cline, with a frequency for allele b ranging from 0.99 in the northernmost point to 0.04 in the southernmost one, a pattern that may be explained by natural selection (FST = 0.46; p < 0.0001) on heat sensitive alleles. The average value of FST (0.088) and Nm (61.12) indicated a high gene flow between adjacent populations. A high correlation was found between genetic and geographic distance (r = 0.83; p < 0.001). The highest genetic identity (IN = 0.988) corresponded to the geographically closest samples from the central area. In one of these localities Cx. quinquefasciatus was predominant and hybrid individuals were detected, while in the other, almost all the specimens were identified as Cx. pipiens. To verify the fertility between Cx. pipiens and Cx. quinquefasciatus from the northern- and southernmost populations, experimental crosses were performed. Viable egg rafts were obtained from both reciprocal crosses. Hatching ranged from 76.5 to 100%. The hybrid progenies were fertile through two subsequent generations

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Quantitation of HIV-1 RNA Viral Load Using NASBA Methodology and Comparison with other Surrogate Markers for Disease Progression

In this study, HIV-1 viral load quantitation determined by Nucleic Acid Sequence Based Amplification (NASBA) was compared with other surrogate disease progression markers (antigen p24, CD4/CD8 cell counts and b-2 microglobulin) in 540 patients followed up at São Paulo, SP, Brazil. HIV-1 RNA detection was statistically associated with the presence of antigen p24, but the viral RNA was also detected in 68% of the antigen p24 negative samples, confirming that NASBA is much more sensitive than the determination of antigen p24. Regarding other surrogate markers, no statistically significant association with the detection of viral RNA was found. The reproducibility of this viral load assay was assessed by 14 runs of the same sample, using different reagents batches. Viral load values in this sample ranged from 5.83 to 6.27 log (CV = 36 %), less than the range (0.5 log) established to the determination of significant viral load changes.

A correlation sensitivity analysis of non-life underwriting risk in solvency capital requirement estimation

This paper analyses the impact of using different correlation assumptions between lines of business when estimating the risk-based capital reserve, the Solvency Capital Requirement (SCR), under Solvency II regulations. A case study is presented and the SCR is calculated according to the Standard Model approach. Alternatively, the requirement is then calculated using an Internal Model based on a Monte Carlo simulation of the net underwriting result at a one-year horizon, with copulas being used to model the dependence between lines of business. To address the impact of these model assumptions on the SCR we conduct a sensitivity analysis. We examine changes in the correlation matrix between lines of business and address the choice of copulas. Drawing on aggregate historical data from the Spanish non-life insurance market between 2000 and 2009, we conclude that modifications of the correlation and dependence assumptions have a significant impact on SCR estimation.

Age-related worm load and worm fecundity patterns in human populations, as indicated by schistosome circulating antigens

Recently, our group determined the relationship between serum CAA levels and fecal egg counts in two foci with very intense Schistosoma mansoni transmission: Maniema (Zaire), an area endemic for S. mansoni since several decades, and Ndombo (Senegal), where transmission has only been established since a few years. The objective was to study and compare age-related worm load and worm fecundity patterns in these two different endemic settings. Here, we will summarize the most important findings and conclusions of this study.

Effect of testosterone on immunocompetence, parasite load, and metabolism in the common wall lizard (Podarcis muralis)

Testosterone can benefit individual fitness by increasing ornament colour, aggressiveness, and sperm quality, but it can also impose both metabolic and immunological costs. However, evidence that testosterone causes immuno suppression in freely living populations is scant. We studied the effects of testosterone on one component of the immune system (i.e., the cell-mediated response to phytohaemagglutinin), parasite load, and metabolic rate in the common wall lizard, Podarcis muralis (Laurenti, 1768). For analyses of immunocompetence and parasitism, male lizards were implanted at the end of the breeding season with either empty or testosterone implants and were returned to their site of capture for 5-6 weeks before recapture. For analyses of the effects of testosterone on metabolic rate, male lizards were captured and implanted before hibernation and were held in the laboratory for 1 week prior to calorimetry. Experimental treatment with testosterone decreased the cell-mediated response to the T-cell mitogen phytohemagglutinin and increased mean metabolic rate. No effects of testosterone on the number of ectoparasites, hemoparasites, and resting metabolic rate could be detected. These results are discussed in the framework of the immunocompetence handicap hypothesis and the immuno-redistribution process hypothesis. [Authors]