Estimation of dynamic latent variable models using simulated nonparametric moments

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. Since 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. Monte Carlo results show that the estimator performs well in comparison to other estimators that have been proposed for estimation of general DLV models.

BORN AT 27 WEEKS OF GESTATION WITH CLASSICAL PKU: CHALLENGES OF DIETETIC MANAGEMENT IN A VERY PRETERM INFANT

Only few cases of classical phenylketonuria (PKU) in premature infants have been reported. Treatment of these patients is challenging due to the lack of a phenylalanine-free amino acid solution for parenteral infusion. The boy was born at 27 weeks of gestation with a weight of 1000 g (P10). He received parenteral nutrition with a protein intake of 3 g/kg/day. On day 7 he was diagnosed with classical PKU (genotype IVS10-11G>A/IVS12+ 1G>A) due to highly elevated phenylalanine (Phe) level in newborn screening (2800 micromol/L). His maximum plasma Phe level reached 3696 micromol/L. Phe intake was stopped for 4 days. During this time the boy received intravenous glucose and lipids as well as little amounts of Phe-free formula by a nasogastric tube. Due to a deficit of essential amino acids and insufficient growth, a parenteral nutrition rich in branched-chain amino-acids and relatively poor in Phe was added, in order to promote protein synthesis without overloading in Phe. Under this regimen, Phe plasma levels normalized on day 19 when intake of natural protein was started. The boy has now a corrected age of 2 years. He shows normal growth parameters and psychomotor development. Despite a long period of highly elevated Phe levels in the postnatal period our patient shows good psychomotor development. The management of premature infants with PKU depends on the child's tolerance to enteral nutrition. It demands an intensive follow-up by an experienced team and dedicated dietician. Appropriate Phe-free parenteral nutrition would be necessary especially in case of gastro-intestinal complications of prematurity.

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.

Genomics of wine yeast: functional analysis of new and variable genes

Projecte de recerca elaborat a partir d’una estada a l’Institut National de la Recherche Agronomique, França, entre 2007 i 2009. Saccharomyces cerevisiae ha estat el llevat utilitzat durant mil.lenis en l'elaboració de vins. Tot i així, es té poc coneixement sobre les pressions de selecció que han actuat en la modelització del genoma dels llevats vínics. S’ha seqüenciat el genoma d'una soca vínica comercial, EC1118, obtenint 31 supercontigs que cobreixen el 97% del genoma de la soca de referència, S288c. S’ha trobat que el genoma de la soca vínica es diferencia bàsicament en la possessió de 3 regions úniques que contenen 34 gens implicats en funcions claus per al procés fermentatiu. A banda, s’han dut a terme estudis de filogènia i synteny (ordre dels gens) que mostren que una d'aquestes tres regions és pròxima a una espècie relacionada amb el gènere Saccharomyces, mentre que les altres dos regions tenen un origen no-Saccharomyces. S’ha identificat mitjançant PCR i seqüenciació a Zygosaccharomyces bailii, una espècie contaminant de les fermentacions víniques, com a espècie donadora d'una de les dues regions. Les hibridacions naturals entre soques de diferents espècies dins del grup Saccharomyces sensu stricto ja han estat descrites. El treball és el primer que presenta hibridacions entre espècies Saccharomyces i no-Saccharomyces (Z. bailii, en aquest cas). També s’assenyala que les noves regions es troben freqüent i diferencialment presents entre els clades de S. cerevisiae, trobant-se de manera gairebé exclusiva en el grup de les soques víniques, suggerint que es tracta d'una adquisició recent de transferència gènica. En general, les dades demostren que el genoma de les soques víniques pateix una constant remodelació mitjançant l'adquisició de gens exògens. Els resultats suggereixen que aquests processos estan afavorits per la proximitat ecològica i estan implicats en l'adaptació molecular de les soques víniques a les condicions d'elevada concentració en sucres, poc nitrogen i elevades concentracions en etanol.

Assessing the Relationship between the Baseline Value of a Continuous Variable and Subsequent Change over Time

Analyzing the relationship between the baseline value and subsequent change of a continuous variable is a frequent matter of inquiry in cohort studies. These analyses are surprisingly complex, particularly if only two waves of data are available. It is unclear for non-biostatisticians where the complexity of this analysis lies and which statistical method is adequate.With the help of simulated longitudinal data of body mass index in children,we review statistical methods for the analysis of the association between the baseline value and subsequent change, assuming linear growth with time. Key issues in such analyses are mathematical coupling, measurement error, variability of change between individuals, and regression to the mean. Ideally, it is better to rely on multiple repeated measurements at different times and a linear random effects model is a standard approach if more than two waves of data are available. If only two waves of data are available, our simulations show that Blomqvist's method - which consists in adjusting for measurement error variance the estimated regression coefficient of observed change on baseline value - provides accurate estimates. The adequacy of the methods to assess the relationship between the baseline value and subsequent change depends on the number of data waves, the availability of information on measurement error, and the variability of change between individuals.

Sobolev orthogonal polynomials on a simplex

Vegeu el resum a l'inici del document de l'arxiu adjunt

Optical counterpart of galactic plane variable radio sources

We present I-band deep CCD exposures of the fields of galactic plane radio variables. An optical counterpart, based on positional coincidence, has been found for 15 of the 27 observed program objects. The Johnson I magnitude of the sources identified is in the range 18-21.

On the use of spatially discrete data to compute energy and mass balance

In many practical applications the state of field soils is monitored by recording the evolution of temperature and soil moisture at discrete depths. We theoretically investigate the systematic errors that arise when mass and energy balances are computed directly from these measurements. We show that, even with no measurement or model errors, large residuals might result when finite difference approximations are used to compute fluxes and storage term. To calculate the limits set by the use of spatially discrete measurements on the accuracy of balance closure, we derive an analytical solution to estimate the residual on the basis of the two key parameters: the penetration depth and the distance between the measurements. When the thickness of the control layer for which the balance is computed is comparable to the penetration depth of the forcing (which depends on the thermal diffusivity and on the forcing period) large residuals arise. The residual is also very sensitive to the distance between the measurements, which requires accurately controlling the position of the sensors in field experiments. We also demonstrate that, for the same experimental setup, mass residuals are sensitively larger than the energy residuals due to the nonlinearity of the moisture transport equation. Our analysis suggests that a careful assessment of the systematic mass error introduced by the use of spatially discrete data is required before using fluxes and residuals computed directly from field measurements.

Interaction of a discrete breather with a lattice junction

We study the scattering of a moving discrete breather (DB) on a junction in a Fermi-Pasta-Ulam chain consisting of two segments with different masses of the particles. We consider four distinct cases: (i) a light-heavy (abrupt) junction in which the DB impinges on the junction from the segment with lighter mass, (ii) a heavy-light junction, (iii) an up mass ramp in which the mass in the heavier segment increases continuously as one moves away from the junction point, and (iv) a down mass ramp. Depending on the mass difference and DB characteristics (frequency and velocity), the DB can either reflect from, or transmit through, or get trapped at the junction or on the ramp. For the heavy-light junction, the DB can even split at the junction into a reflected and a transmitted DB. The latter is found to subsequently split into two or more DBs. For the down mass ramp the DB gets accelerated in several stages, with accompanying radiation (phonons). These results are rationalized by calculating the Peierls-Nabarro barrier for the various cases. We also point out implications of our results in realistic situations such as electron-phonon coupled chains.

Quantification of permanent and variable charges in reference soils of the State of Pernambuco, Brazil

The electrical charges in soil particles are divided into structural or permanent charges and variable charges. Permanent charges develop on the soil particle surface by isomorphic substitution. Variable charges arise from dissociation and association of protons (H+), protonation or deprotonation, and specific adsorption of cations and anions. The aim of this study was to quantify the permanent charges and variable charges of Reference Soils of the State of Pernambuco, Brazil. To do so, 24 subsurface profiles from different regions (nine in the Zona da Mata, eight in the Agreste, and seven in the Sertão) were sampled, representing approximately 80 % of the total area of the state. Measurements were performed using cesium chloride solution. Determination was made of the permanent charges and the charges in regard to the hydroxyl functional groups through selective ion exchange of Cs+ by Li+ and Cs+ by NH4+, respectively. All the soils analyzed exhibited variable cation exchange capacity, with proportions from 0.16 to 0.60 and an average of 0.40 when related to total cation exchange capacity.

Classical predictive electrodynamics of two charges with radiation: General framework. I.

Outgoing radiation is introduced in the framework of the classical predictive electrodynamics using LorentzDiracs equation as a subsidiary condition. In a perturbative scheme in the charges the first radiative self-terms of the accelerations, momentum and angular momentum of a two charge system without external field are calculated.

Classical predictive electrodynamics of two charges with radiation: Energy and 3. momentum balance and scattering crossections-II

We deal with a classical predictive mechanical system of two spinless charges where radiation is considered and there are no external fields. The terms (2,2)Paa of the expansion in the charges of the HamiltonJacobi momenta are calculated. Using these, together with known previous results, we can obtain the paa up to the fourth order. Then we have calculated the radiated energy and the 3-momentum in a scattering process as functions of the impact parameter and the incident energy for the former and 3-momentum for the latter. Scattering cross-sections are also calculated. Good agreement with well known results, including those of quantum electrodynamics, has been found.

Orthogonal protection of peptides and peptoids for cyclization by the thiol-ene reaction and conjugation

Cyclic peptides and peptoids were prepared using the thiolene Michael-type reaction. The linear precursors were provided with additional functional groups allowing for subsequent conjugation: an orthogonally protected thiol, a protected maleimide, or an alkyne. The functional group for conjugation was placed either within the cycle or in an external position. The click reactions employed for conjugation with suitably derivatized nucleoside or oligonucleotides were either cycloadditions (Diels-Alder, Cu(I)-catalyzed azide-alkyne) or the same Michael-type reaction as for cyclization.