Hglm: A package for fitting hierarchical generalized linear models

We present the hglm package for fitting hierarchical generalized linear models. It can be used for linear mixed models and generalized linear mixed models with random effects for a variety of links and a variety of distributions for both the outcomes and the random effects. Fixed effects can also be fitted in the dispersion part of the model.

Genetic heterogeneity of residual variance : estimation of variance components using double hierarchical generalized linear models

Background: The sensitivity to microenvironmental changes varies among animals and may be under genetic control. It is essential to take this element into account when aiming at breeding robust farm animals. Here, linear mixed models with genetic effects in the residual variance part of the model can be used. Such models have previously been fitted using EM and MCMC algorithms. Results: We propose the use of double hierarchical generalized linear models (DHGLM), where the squared residuals are assumed to be gamma distributed and the residual variance is fitted using a generalized linear model. The algorithm iterates between two sets of mixed model equations, one on the level of observations and one on the level of variances. The method was validated using simulations and also by re-analyzing a data set on pig litter size that was previously analyzed using a Bayesian approach. The pig litter size data contained 10,060 records from 4,149 sows. The DHGLM was implemented using the ASReml software and the algorithm converged within three minutes on a Linux server. The estimates were similar to those previously obtained using Bayesian methodology, especially the variance components in the residual variance part of the model. Conclusions: We have shown that variance components in the residual variance part of a linear mixed model can be estimated using a DHGLM approach. The method enables analyses of animal models with large numbers of observations. An important future development of the DHGLM methodology is to include the genetic correlation between the random effects in the mean and residual variance parts of the model as a parameter of the DHGLM.

Mixed-signal analog-digital circuits design on the pre-diffused digital array using trapezoidal association of transistors

The mixed-signal and analog design on a pre-diffused array is a challenging task, given that the digital array is a linear matrix arrangement of minimum-length transistors. To surmount this drawback a specific discipline for designing analog circuits over such array is required. An important novel technique proposed is the use of TAT (Trapezoidal Associations of Transistors) composite transistors on the semi-custom Sea-Of-Transistors (SOT) array. The analysis and advantages of TAT arrangement are extensively analyzed and demonstrated, with simulation and measurement comparisons to equivalent single transistors. Basic analog cells were also designed as well in full-custom and TAT versions in 1.0mm and 0.5mm digital CMOS technologies. Most of the circuits were prototyped in full-custom and TAT-based on pre-diffused SOT arrays. An innovative demonstration of the TAT technique is shown with the design and implementation of a mixed-signal analog system, i. e., a fully differential 2nd order Sigma-Delta Analog-to-Digital (A/D) modulator, fabricated in both full-custom and SOT array methodologies in 0.5mm CMOS technology from MOSIS foundry. Three test-chips were designed and fabricated in 0.5mm. Two of them are IC chips containing the full-custom and SOT array versions of a 2nd-Order Sigma-Delta A/D modulator. The third IC contains a transistors-structure (TAT and single) and analog cells placed side-by-side, block components (Comparator and Folded-cascode OTA) of the Sigma-Delta modulator.

Low cost BIST techniques for linear and non-linear analog circuits

With the ever increasing demands for high complexity consumer electronic products, market pressures demand faster product development and lower cost. SoCbased design can provide the required design flexibility and speed by allowing the use of IP cores. However, testing costs in the SoC environment can reach a substantial percent of the total production cost. Analog testing costs may dominate the total test cost, as testing of analog circuits usually require functional verification of the circuit and special testing procedures. For RF analog circuits commonly used in wireless applications, testing is further complicated because of the high frequencies involved. In summary, reducing analog test cost is of major importance in the electronic industry today. BIST techniques for analog circuits, though potentially able to solve the analog test cost problem, have some limitations. Some techniques are circuit dependent, requiring reconfiguration of the circuit being tested, and are generally not usable in RF circuits. In the SoC environment, as processing and memory resources are available, they could be used in the test. However, the overhead for adding additional AD and DA converters may be too costly for most systems, and analog routing of signals may not be feasible and may introduce signal distortion. In this work a simple and low cost digitizer is used instead of an ADC in order to enable analog testing strategies to be implemented in a SoC environment. Thanks to the low analog area overhead of the converter, multiple analog test points can be observed and specific analog test strategies can be enabled. As the digitizer is always connected to the analog test point, it is not necessary to include muxes and switches that would degrade the signal path. For RF analog circuits, this is specially useful, as the circuit impedance is fixed and the influence of the digitizer can be accounted for in the design phase. Thanks to the simplicity of the converter, it is able to reach higher frequencies, and enables the implementation of low cost RF test strategies. The digitizer has been applied successfully in the testing of both low frequency and RF analog circuits. Also, as testing is based on frequency-domain characteristics, nonlinear characteristics like intermodulation products can also be evaluated. Specifically, practical results were obtained for prototyped base band filters and a 100MHz mixer. The application of the converter for noise figure evaluation was also addressed, and experimental results for low frequency amplifiers using conventional opamps were obtained. The proposed method is able to enhance the testability of current mixed-signal designs, being suitable for the SoC environment used in many industrial products nowadays.

Sobre o sistema de amortização linear crescente

Aiming to reach a compromise solution to the issues of default risk and the payment capacity of takers of housing loans, Jorge Oscar de Mello Flôres submited to the Banco Nacional de Habitação, which was then in charge of the Brazilian System of Housing Financing, what he named as the Linearly Increasing System of Amortization. (LISA). Following a critical analysis of the LISA, it is proposed the alternative named as the Generalyzed System of Mixed Amortization (GSMA).

Bifurcation of limit cycles from an n-dimensional linear center inside a class of piecewise linear differential systems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A Weak Maximum Principle for Optimal Control Problems with Nonsmooth Mixed Constraints

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Miscibility in a degenerate fermionic mixture induced by linear coupling

We consider a one-dimensional mean-field-hydrodynamic model of a two-component degenerate Fermi gas in an external trap, each component representing a spin state of the same atom. We demonstrate that the interconversion between them (linear coupling), imposed by a resonant electromagnetic wave, transforms the immiscible binary gas into a miscible state, if the coupling constant, kappa, exceeds a critical value, kappa(cr). The effect is predicted in a variational approximation, and confirmed by numerical solutions. Unlike the recently studied model of a binary Bose-Einsten condensate with the linear coupling, the components in the immiscible phase of the binary fermion mixture never fill two separated domains with a wall between them, but rather form antilocked (pi-phase-shifted) density waves. Another difference from the bosonic mixture is spontaneous breaking of symmetry between the two components in terms of the numbers of atoms in them, N(1) and N(2). The latter effect is characterized by the parameter nu equivalent to(N(1)-N(2))/(N(1)+N(2)) (only N(1)+N(2) is a conserved quantity), the onset of miscibility at kappa >=kappa(cr) meaning a transition to nu equivalent to 0. At kappa linear coupling, and the developing spatial structure resembles the known Larkin-Ovchinnikov-Fulde-Ferrell states in the Fermi mixture. Dynamical states, when kappa is suddenly switched from zero to a value exceeding kappa(cr), are considered too. In the latter case, the system features oscillatory relaxation to the mixed state.

Valorização da reserva operativa em um ambiente competitivo para geradores hidrelétricos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Procedimento MIXED do SAS® para análise de modelos mistos

O modelo misto consiste numa importante classe de modelos que tem sido tradicionalmente analisada por meio de procedimentos da análise de variância. Nos modelos mistos, três aspectos são fundamentais: estimação e testes de hipóteses dos efeitos fixos, predição dos efeitos aleatórios e estimação dos componentes de variância. Na análise de modelos lineares mistos desbalanceados, a estimação dos componentes de variância é de fundamental importância e depende da estrutura de covariâncias e dos métodos de estimação utilizados. Nesse contexto, este artigo pretende apresentar os principais métodos de estimação e de análise utilizados no estudo de modelos lineares mistos com estruturas gerais de covariâncias nos efeitos aleatórios, disponíveis no procedimento MIXED, do SAS (Statistical Analysis System).

Bayesian longitudinal data analysis with mixed models and thick-tailed distributions using MCMC

Linear mixed effects models are frequently used to analyse longitudinal data, due to their flexibility in modelling the covariance structure between and within observations. Further, it is easy to deal with unbalanced data, either with respect to the number of observations per subject or per time period, and with varying time intervals between observations. In most applications of mixed models to biological sciences, a normal distribution is assumed both for the random effects and for the residuals. This, however, makes inferences vulnerable to the presence of outliers. Here, linear mixed models employing thick-tailed distributions for robust inferences in longitudinal data analysis are described. Specific distributions discussed include the Student-t, the slash and the contaminated normal. A Bayesian framework is adopted, and the Gibbs sampler and the Metropolis-Hastings algorithms are used to carry out the posterior analyses. An example with data on orthodontic distance growth in children is discussed to illustrate the methodology. Analyses based on either the Student-t distribution or on the usual Gaussian assumption are contrasted. The thick-tailed distributions provide an appealing robust alternative to the Gaussian process for modelling distributions of the random effects and of residuals in linear mixed models, and the MCMC implementation allows the computations to be performed in a flexible manner.

A linear optimization approach to the combined production planning model

Two fundamental processes usually arise in the production planning of many industries. The first one consists of deciding how many final products of each type have to be produced in each period of a planning horizon, the well-known lot sizing problem. The other process consists of cutting raw materials in stock in order to produce smaller parts used in the assembly of final products, the well-studied cutting stock problem. In this paper the decision variables of these two problems are dependent of each other in order to obtain a global optimum solution. Setups that are typically present in lot sizing problems are relaxed together with integer frequencies of cutting patterns in the cutting problem. Therefore, a large scale linear optimizations problem arises, which is exactly solved by a column generated technique. It is worth noting that this new combined problem still takes the trade-off between storage costs (for final products and the parts) and trim losses (in the cutting process). We present some sets of computational tests, analyzed over three different scenarios. These results show that, by combining the problems and using an exact method, it is possible to obtain significant gains when compared to the usual industrial practice, which solve them in sequence. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Selection of alfalfa cultivars adapted for tropical environments with repeated measures using PROC MIXED of SAS (R) System

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Effect of a heat-treatment on the linear dimensional change of a hard chairside reline resin

Statement of problem. Little data are available regarding the effect of heat-treatments on the dimensional stability of hard chairside reline resins. Purpose. The objective of this in vitro study was to evaluate whether a heat-treatment improves the dimensional stability of the reline resin Duraliner II and to compare the linear dimensional changes of this material with the heat-polymerized acrylic resin Lucitone 550. Material and methods. The materials were mixed according to the manufacturer's instructions and packed into a stainless steel split mold (50.0 mm diameter and 0.5 mm thickness) with reference points (A, B, C, and D). Duraliner II specimens were polymerized for 12 minutes in water at 37°C and bench cooled to room temperature before being removed from the mold. Twelve specimens were made and divided into 2 groups: group 1 specimens (n=6) were left untreated, and group 2 specimens (n=6) were submitted to a heat-treatment in a water bath at 55°C for 10 minutes and then bench cooled to room temperature. The 6 Lucitone specimens (control group) were polymerized in a water bath for 9 hours at 71°C. The specimens were removed after the mold reached the room temperature. A Nikon optical comparator was used to measure the distances between the reference points (AB and CD) on the stainless steel mold (baseline readings) and on the specimens to the nearest 0.001 mm. Measurements were made after processing and after the specimens had been stored in distilled water at 37°C for 8 different periods of time. Data were subjected to analysis of variance with repeated measures, followed by Tukey's multiple comparison test (P<.05). Results. All specimens exhibited shrinkage after processing (control, -0.41%; group 1, -0.26%; and group 2, -0.51%). Group 1 specimens showed greater shrinkage (-1.23%) than the control (-0.23%) and group 2 (-0.81%) specimens after 60 days of storage in water (P<.05). Conclusion. Within the limitations of this study, a significant improvement of the long-term dimensional stability of the Duraliner II reline resin was observed when the specimens were heat-treated. However, the shrinkage remained considerably higher than the denture base resin Lucitone 550. Copyright © 2002 by The Editorial Council of The Journal of Prosthetic Dentistry.

Stochastic Stability for Markovian Jump Linear Systems Subject to a Crucial Failure Event

This paper is concerned with the stability of discrete-time linear systems subject to random jumps in the parameters, described by an underlying finite-state Markov chain. In the model studied, a stopping time τ Δ is associated with the occurrence of a crucial failure after which the system is brought to a halt for maintenance. The usual stochastic stability concepts and associated results are not indicated, since they are tailored to pure infinite horizon problems. Using the concept named stochastic τ-stability, equivalent conditions to ensure the stochastic stability of the system until the occurrence of τ Δ is obtained. In addition, an intermediary and mixed case for which τ represents the minimum between the occurrence of a fix number N of failures and the occurrence of a crucial failure τ Δ is also considered. Necessary and sufficient conditions to ensure the stochastic τ-stability are provided in this setting that are auxiliary to the main result.