Fitting models for the joint action of two drugs using SAS(R)

Combinations of drugs are increasingly being used for a wide variety of diseases and conditions. A pre-clinical study may allow the investigation of the response at a large number of dose combinations. In determining the response to a drug combination, interest may lie in seeking evidence of synergism, in which the joint action is greater than the actions of the individual drugs, or of antagonism, in which it is less. Two well-known response surface models representing no interaction are Loewe additivity and Bliss independence, and Loewe or Bliss synergism or antagonism is defined relative to these. We illustrate an approach to fitting these models for the case in which the marginal single drug dose-response relationships are represented by four-parameter logistic curves with common upper and lower limits, and where the response variable is normally distributed with a common variance about the dose-response curve. When the dose-response curves are not parallel, the relative potency of the two drugs varies according to the magnitude of the desired effect and the models for Loewe additivity and synergism/antagonism cannot be explicitly expressed. We present an iterative approach to fitting these models without the assumption of parallel dose-response curves. A goodness-of-fit test based on residuals is also described. Implementation using the SAS NLIN procedure is illustrated using data from a pre-clinical study. Copyright © 2007 John Wiley & Sons, Ltd.

A hands-on approach for fitting long-term survival models under the GAMLSS framework

In many data sets from clinical studies there are patients insusceptible to the occurrence of the event of interest. Survival models which ignore this fact are generally inadequate. The main goal of this paper is to describe an application of the generalized additive models for location, scale, and shape (GAMLSS) framework to the fitting of long-term survival models. in this work the number of competing causes of the event of interest follows the negative binomial distribution. In this way, some well known models found in the literature are characterized as particular cases of our proposal. The model is conveniently parameterized in terms of the cured fraction, which is then linked to covariates. We explore the use of the gamlss package in R as a powerful tool for inference in long-term survival models. The procedure is illustrated with a numerical example. (C) 2009 Elsevier Ireland Ltd. All rights reserved.

Transformed generalized linear models

The estimation of data transformation is very useful to yield response variables satisfying closely a normal linear model, Generalized linear models enable the fitting of models to a wide range of data types. These models are based on exponential dispersion models. We propose a new class of transformed generalized linear models to extend the Box and Cox models and the generalized linear models. We use the generalized linear model framework to fit these models and discuss maximum likelihood estimation and inference. We give a simple formula to estimate the parameter that index the transformation of the response variable for a subclass of models. We also give a simple formula to estimate the rth moment of the original dependent variable. We explore the possibility of using these models to time series data to extend the generalized autoregressive moving average models discussed by Benjamin er al. [Generalized autoregressive moving average models. J. Amer. Statist. Assoc. 98, 214-223]. The usefulness of these models is illustrated in a Simulation study and in applications to three real data sets. (C) 2009 Elsevier B.V. All rights reserved.

Egg production curve fitting using nonlinear models for selected and nonselected lines of White Leghorn hens

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

Análise dos processos de decomposição e sucessão ecológica em carcaças de suínos (Sus scrofa L.) mortos por disparo de arma de fogo e overdose de cocaína e protocolo de procedimento diante de corpo de delito

Pós-graduação em Biologia Geral e Aplicada - IBB

Genetic parameter estimation for pre- and post-weaning traits in Brahman cattle in Brazil

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

Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

NONLINEAR TUBE-FITTING FOR THE ANALYSIS OF ANATOMICAL AND FUNCTIONAL STRUCTURES

We are concerned with the estimation of the exterior surface of tube-shaped anatomical structures. This interest is motivated by two distinct scientific goals, one dealing with the distribution of HIV microbicide in the colon and the other with measuring degradation in white-matter tracts in the brain. Our problem is posed as the estimation of the support of a distribution in three dimensions from a sample from that distribution, possibly measured with error. We propose a novel tube-fitting algorithm to construct such estimators. Further, we conduct a simulation study to aid in the choice of a key parameter of the algorithm, and we test our algorithm with validation study tailored to the motivating data sets. Finally, we apply the tube-fitting algorithm to a colon image produced by single photon emission computed tomography (SPECT)and to a white-matter tract image produced using diffusion tensor `imaging (DTI).

Image registration for triggered and non-triggered DTI of the human kidney: Reduced variability of diffusion parameter estimation

BACKGROUND: To investigate if non-rigid image-registration reduces motion artifacts in triggered and non-triggered diffusion tensor imaging (DTI) of native kidneys. A secondary aim was to determine, if improvements through registration allow for omitting respiratory-triggering. METHODS: Twenty volunteers underwent coronal DTI of the kidneys with nine b-values (10-700 s/mm2 ) at 3 Tesla. Image-registration was performed using a multimodal nonrigid registration algorithm. Data processing yielded the apparent diffusion coefficient (ADC), the contribution of perfusion (FP ), and the fractional anisotropy (FA). For comparison of the data stability, the root mean square error (RMSE) of the fitting and the standard deviations within the regions of interest (SDROI ) were evaluated. RESULTS: RMSEs decreased significantly after registration for triggered and also for non-triggered scans (P < 0.05). SDROI for ADC, FA, and FP were significantly lower after registration in both medulla and cortex of triggered scans (P < 0.01). Similarly the SDROI of FA and FP decreased significantly in non-triggered scans after registration (P < 0.05). RMSEs were significantly lower in triggered than in non-triggered scans, both with and without registration (P < 0.05). CONCLUSION: Respiratory motion correction by registration of individual echo-planar images leads to clearly reduced signal variations in renal DTI for both triggered and particularly non-triggered scans. Secondarily, the results suggest that respiratory-triggering still seems advantageous.J. Magn. Reson. Imaging 2014. (c) 2014 Wiley Periodicals, Inc.

Deciduous tree reconstruction algorithm based on cylinder fitting from mobile terrestrial laser scanned point clouds

Vector reconstruction of objects from an unstructured point cloud obtained with a LiDAR-based system (light detection and ranging) is one of the most promising methods to build three dimensional models of orchards. The cylinder fitting method for woody structure reconstruction of leafless trees from point clouds obtained with a mobile terrestrial laser scanner (MTLS) has been analysed. The advantage of this method is that it performs reconstruction in a single step. The most time consuming part of the algorithm is generation of the cylinder direction, which must be recalculated at the inclusion of each point in the cylinder. The tree skeleton is obtained at the same time as the cluster of cylinders is formed. The method does not guarantee a unique convergence and the reconstruction parameter values must be carefully chosen. A balanced processing of clusters has also been defined which has proven to be very efficient in terms of processing time by following the hierarchy of branches, predecessors and successors. The algorithm was applied to simulated MTLS of virtual orchard models and to MTLS data of real orchards. The constraints applied in the method have been reviewed to ensure better convergence and simpler use of parameters. The results obtained show a correct reconstruction of the woody structure of the trees and the algorithm runs in linear logarithmic time

Microfiber Fabry-Perot interferometer for dual-parameter sensing

We propose and demonstrate a microfiber Fabry-Perot interferometer (MFPI) fabricated by taper-drawing microfiber at the center of a uniform fiber Bragg grating (FBG). The MFPI employing the two separated sections of FBG as reflectors and a length of microfiber as its cavity is derived. Theoretic study shows that the reflection spectrum of such MFPI is consisted of two parts-interference fringes induced by multi-beam interference and reflection spectrum envelope induced by FBGs. Temperature affects both interference fringes and reflection wavelength of FBGs while ambient refractive index (RI) only influences the interference fringes, i.e., MFPI has different response to temperature and RI. Therefore, MFPI for simultaneous sensing of RI and temperature is experimentally demonstrated by tracking a reflection peak of interference fringes and the Bragg wavelength of the FBGs, which are respectively assisted by frequency domain processing and Gaussian fitting of the optical spectrum. Consequently, wavelength measurement resolution of 0.5 pm is realized. © 1983-2012 IEEE.

Terapêutica interdisciplinar para fratura cominutiva de côndilo por projétil de arma de fogo: enfoque miofuncional

TEMA: ferimentos causados por projéteis de arma de fogo apresentam alta incidência na região da cabeça e face. A articulação temporomandibular pode estar envolvida, além de estruturas anatômicas importantes como o nervo facial, necessitando de equipe multidisciplinar para efetuar tratamento adequado. PROCEDIMENTOS: apresentação de caso clínico de fratura condilar cominutiva causada por projétil de arma de fogo tratado de forma não-cirúrgica associado à terapia miofuncional orofacial. Paciente encaminhado para avaliação e procedimentos fonoaudiológicos após conduta da equipe de cirurgia bucomaxilofacial, sem remoção do projétil, alojado superficialmente, próximo da origem do músculo esternocleidomastóideo à direita, com fratura condilar cominutiva e lesão do nervo facial. Foram aspectos observados em avaliação: mordida aberta anterior, importante redução da amplitude dos movimentos mandibulares com desvios para o lado acometido, ausência de lateralidade contralateral, dor muscular, paralisia e parestesia em terço médio e superior da hemiface direita. Realizadas sessões de terapia miofuncional seguindo protocolo específico para traumas de face constando de: drenagem de edema; manipulações na musculatura levantadora da mandíbula ipsilateral; ampliação e correção dos movimentos mandibulares; procedimentos específicos referentes à paralisia facial e reorganização funcional direcionada. RESULTADOS: após oito semanas de terapia os resultados obtidos mostram restabelecimento de amplitude e da simetria dos movimentos mandibulares, reorganização da mastigação, adequação da deglutição e fala, remissão da sintomatologia dolorosa e remissão da paralisia do terço médio. CONCLUSÃO: o tratamento conservador da fratura por meio da terapia miofuncional orofacial resultou na reabilitação funcional da mandíbula e face dirigindo os movimentos e estimulando a adequação das funções estomatognáticas.

Dynamic and steady: shear rheological properties of xanthan and guar gums dispersed in yellow passion fruit pulp (Passiflora edulis f. flavicarpa)

Yellow passion fruit pulp is unstable, presenting phase separation that can be avoided by the addition of hydrocolloids. For this purpose, xanthan and guar gum [0.3, 0.7 and 1.0% (w/w)] were added to yellow passion fruit pulp and the changes in the dynamic and steady - shear rheological behavior evaluated. Xanthan dispersions showed a more pronounced pseudoplasticity and the presence of yield stress, which was not observed in the guar gum dispersions. Cross model fitting to flow curves showed that the xanthan suspensions also had higher zero shear viscosity than the guar suspensions, and, for both gums, an increase in temperature led to lower values for this parameter. The gums showed different behavior as a function of temperature in the range of 5 - 35ºC. The activation energy of the apparent viscosity was dependent on the shear rate and gum concentration for guar, whereas for xanthan these values only varied with the concentration. The mechanical spectra were well described by the generalized Maxwell model and the xanthan dispersions showed a more elastic character than the guar dispersions, with higher values for the relaxation time. Xanthan was characterized as a weak gel, while guar presented a concentrated solution behavior. The simultaneous evaluation of temperature and concentration showed a stronger influence of the polysaccharide concentration on the apparent viscosity and the G' and G" moduli than the variation in temperature.

Caracterização das vítimas de ferimentos por arma de fogo, atendidas pelo Serviço de Atendimento Móvel de Urgência em Campo Grande-MS

A violência, em todas as suas apresentações, tem crescido de forma avassaladora no país, elevando os indicadores de morbidade e mortalidade por causas externas. O uso das armas de fogo, muitas vezes, faz vítimas fatais ou que podem permanecer sequeladas. Tal fato tem aumentado o ônus do estado com internações hospitalares e acréscimo dos anos de vida perdidos da população jovem, que constitui a grande maioria dessas vítimas. Nesse sentido, o presente estudo objetivou fazer um levantamento das vítimas de ferimentos por arma de fogo, atendidas pelo Serviço de Atendimento Móvel de Urgência no município de Campo Grande-MS, no período de abril de 2005 a abril de 2007, nos dois primeiros anos de funcionamento, desde a implantação desse serviço na capital do Estado de Mato Grosso do Sul. Realizou-se estudo descritivo, baseado em análise documental do sistema de informação do SAMU do município de Campo Grande-MS. Foram descritos 233 atendimentos. Os resultados evidenciaram 213 vítimas do sexo masculino, a faixa etária mais acometida foi dos 20 aos 24 anos de idade, a cabeça e o pescoço foram as partes do corpo mais atingidas e a região Sul do município de Campo Grande foi a que concentrou maior número de atendimentos. Conclui-se que a violência por arma de fogo em Campo Grande-MS atinge a camada economicamente ativa da população e provém de regiões de bolsões de pobreza e desigualdade social, justificando a implantação de um serviço como o SAMU.