Comparison of two population pharmacokinetic programs, NONMEM and P-PHARM, for tacrolimus

Objectives: To compare the population modelling programs NONMEM and P-PHARM during investigation of the pharmacokinetics of tacrolimus in paediatric liver-transplant recipients. Methods: Population pharmacokinetic analysis was performed using NONMEM and P-PHARM on retrospective data from 35 paediatric liver-transplant patients receiving tacrolimus therapy. The same data were presented to both programs. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F). Covariates screened for influence on these parameters were weight, age, gender, post-operative day, days of tacrolimus therapy, transplant type, biliary reconstructive procedure, liver function tests, creatinine clearance, haematocrit, corticosteroid dose, and potential interacting drugs. Results: A satisfactory model was developed in both programs with a single categorical covariate - transplant type - providing stable parameter estimates and small, normally distributed (weighted) residuals. In NONMEM, the continuous covariates - age and liver function tests - improved modelling further. Mean parameter estimates were CL/F (whole liver) = 16.3 1/h, CL/F (cut-down liver) = 8.5 1/h and V/F = 565 1 in NONMEM, and CL/F = 8.3 1/h and V/F = 155 1 in P-PHARM. Individual Bayesian parameter estimates were CL/F (whole liver) = 17.9 +/- 8.8 1/h, CL/F (cutdown liver) = 11.6 +/- 18.8 1/h and V/F = 712 792 1 in NONMEM, and CL/F (whole liver) = 12.8 +/- 3.5 1/h, CL/F (cut-down liver) = 8.2 +/- 3.4 1/h and V/F = 221 1641 in P-PHARM. Marked interindividual kinetic variability (38-108%) and residual random error (approximately 3 ng/ml) were observed. P-PHARM was more user friendly and readily provided informative graphical presentation of results. NONMEM allowed a wider choice of errors for statistical modelling and coped better with complex covariate data sets. Conclusion: Results from parametric modelling programs can vary due to different algorithms employed to estimate parameters, alternative methods of covariate analysis and variations and limitations in the software itself.

Exploring complex fitness surfaces: Multiple ornamentation and polymorphism in male guppies

The sexual ornamentation used by male guppies to attract females comprises many components, each of which varies considerably among males. Although natural and sexual selection have been shown to contribute to divergence among populations in male sexual ornaments, the role of sexual selection in maintaining polymorphism within populations is less clear. We used both parametric quadratic regression and nonparametric projection pursuit regression techniques to reveal the major axes of non-linear sexual selection on male ornaments. We visualized the fitness surfaces defined by these axes using thin-plate splines to allow a direct comparison of the two methodologies. Identification of the major axes of selection and their visualization was critical in determining the form and strength of nonlinear selection. Both types of analysis revealed fitness surfaces comprising three peaks, suggesting that there is more than one way to make an attractive guppy. Disruptive selection may be an important process underlying the presence of multiple sexual ornaments and may contribute to the maintenance of the high levels of polymorphism in male sexual ornaments found in guppy populations.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

A mechanistic model for carbon dioxide corrosion of mild steel in the presence of protective iron carbonate films - Part 3: Film growth model

A model of iron carbonate (FeCO3) film growth is proposed, which is an extension of the recent mechanistic model of carbon dioxide (CO2) corrosion by Nesic, et al. In the present model, the film growth occurs by precipitation of iron carbonate once saturation is exceeded. The kinetics of precipitation is dependent on temperature and local species concentrations that are calculated by solving the coupled species transport equations. Precipitation tends to build up a layer of FeCO3 on the surface of the steel and reduce the corrosion rate. On the other hand, the corrosion process induces voids under the precipitated film, thus increasing the porosity and leading to a higher corrosion rate. Depending on the environmental parameters such as temperature, pH, CO2 partial pressure, velocity, etc., the balance of the two processes can lead to a variety of outcomes. Very protective films and low corrosion rates are predicted at high pH, temperature, CO2 partial pressure, and Fe2+ ion concentration due to formation of dense protective films as expected. The model has been successfully calibrated against limited experimental data. Parametric testing of the model has been done to gain insight into the effect of various environmental parameters on iron carbonate film formation. The trends shown in the predictions agreed well with the general understanding of the CO2 corrosion process in the presence of iron carbonate films. The present model confirms that the concept of scaling tendency is a good tool for predicting the likelihood of protective iron carbonate film formation.

Localization of a novel melanoma susceptibility locus to 1p22

Over the past 20 years, the incidence of cutaneous malignant melanoma (CMM) has increased dramatically worldwide. A positive family history of the disease is among the most established risk factors for CMM; it is estimated that 10% of CMM cases result from an inherited predisposition. Although mutations in two genes, CDKN2A and CDK4, have been shown to confer an increased risk of CMM, they account for only 20%-25% of families with multiple cases of CMM. Therefore, to localize additional loci involved in melanoma susceptibility, we have performed a genomewide scan for linkage in 49 Australian pedigrees containing at least three CMM cases, in which CDKN2A and CDK4 involvement has been excluded. The highest two-point parametric LOD score (1.82; recombination fraction [theta] 0.2) was obtained at D1S2726, which maps to the short arm of chromosome 1 (1p22). A parametric LOD score of 4.65 (theta = 0) and a nonparametric LOD score of 4.19 were found at D1S2779 in nine families selected for early age at onset. Additional typing yielded seven adjacent markers with LOD scores 13 in this subset, with the highest parametric LOD score, 4.95 (theta = 0) ( nonparametric LOD score 5.37), at D1S2776. Analysis of 33 additional multiplex families with CMM from several continents provided further evidence for linkage to the 1p22 region, again strongest in families with the earliest mean age at diagnosis. A nonparametric ordered sequential analysis was used, based on the average age at diagnosis in each family. The highest LOD score, 6.43, was obtained at D1S2779 and occurred when the 15 families with the earliest ages at onset were included. These data provide significant evidence of a novel susceptibility gene for CMM located within chromosome band 1p22.

Rock fracture mean trace length estimation and confidence interval calculation using maximum likelihood methods

There has been a resurgence of interest in the mean trace length estimator of Pahl for window sampling of traces. The estimator has been dealt with by Mauldon and Zhang and Einstein in recent publications. The estimator is a very useful one in that it is non-parametric. However, despite some discussion regarding the statistical distribution of the estimator, none of the recent works or the original work by Pahl provide a rigorous basis for the determination a confidence interval for the estimator or a confidence region for the estimator and the corresponding estimator of trace spatial intensity in the sampling window. This paper shows, by consideration of a simplified version of the problem but without loss of generality, that the estimator is in fact the maximum likelihood estimator (MLE) and that it can be considered essentially unbiased. As the MLE, it possesses the least variance of all estimators and confidence intervals or regions should therefore be available through application of classical ML theory. It is shown that valid confidence intervals can in fact be determined. The results of the work and the calculations of the confidence intervals are illustrated by example. (C) 2003 Elsevier Science Ltd. All rights reserved.

A model to investigate hepatic extraction of oxygen during anaesthesia in the dog

Measurement of hepatic oxygen extraction was performed on six healthy Greyhound dogs over a two hour period. The Greyhounds were anaesthetised and a right subcostal surgical incision performed. Ultrasonic flow transducers were used to measure flow rate in the hepatic artery and the portal vein. The blood oxygen tensions in arterial blood and in the portal and hepatic veins were also measured. Hepatic oxygen extraction remained stable throughout the study, despite a steady decline in arterial blood pressure. The methodology described in this study provides a direct measure of oxygen uptake by the liver in the dog and could readily be used to investigate hepatic uptake of drugs. (C) 2003 Elsevier Ltd. All rights reserved.

Postbuckling of piezoelectric FGM plates subject to thermo-electro-mechanical loading

In this paper, we examine the postbuckling behavior of functionally graded material FGM rectangular plates that are integrated with surface-bonded piezoelectric actuators and are subjected to the combined action of uniform temperature change, in-plane forces, and constant applied actuator voltage. A Galerkin-differential quadrature iteration algorithm is proposed for solution of the non-linear partial differential governing equations. To account for the transverse shear strains, the Reddy higher-order shear deformation plate theory is employed. The bifurcation-type thermo-mechanical buckling of fully clamped plates, and the postbuckling behavior of plates with more general boundary conditions subject to various thermo-electro-mechanical loads, are discussed in detail. Parametric studies are also undertaken, and show the effects of applied actuator voltage, in-plane forces, volume fraction exponents, temperature change, and the character of boundary conditions on the buckling and postbuckling characteristics of the plates. (C) 2003 Elsevier Science Ltd. All rights reserved.

Dynamic stability of laminated FGM plates based on higher-order shear deformation theory

This paper conducts a dynamic stability analysis of symmetrically laminated FGM rectangular plates with general out-of-plane supporting conditions, subjected to a uniaxial periodic in-plane load and undergoing uniform temperature change. Theoretical formulations are based on Reddy's third-order shear deformation plate theory, and account for the temperature dependence of material properties. A semi-analytical Galerkin-differential quadrature approach is employed to convert the governing equations into a linear system of Mathieu-Hill equations from which the boundary points on the unstable regions are determined by Bolotin's method. Free vibration and bifurcation buckling are also discussed as subset problems. Numerical results are presented in both dimensionless tabular and graphical forms for laminated plates with FGM layers made of silicon nitride and stainless steel. The influences of various parameters such as material composition, layer thickness ratio, temperature change, static load level, boundary constraints on the dynamic stability, buckling and vibration frequencies are examined in detail through parametric studies.

Stokes-operator-squeezed continuous-variable polarization states

We investigate nonclassical Stokes-operator variances in continuous-wave polarization-squeezed laser light generated from one and two optical parametric amplifiers. A general expression of how Stokes-operator variances decompose into two-mode quadrature operator variances is given. Stokes parameter variance spectra for four different polarization-squeezed states have been measured and compared with a coherent state. Our measurement results are visualized by three-dimensional Stokes-operator noise volumes mapped on the quantum Poincare sphere. We quantitatively compare the channel capacity of the different continuous-variable polarization states for communication protocols. It is shown that squeezed polarization states provide 33% higher channel capacities than the optimum coherent beam protocol.

Optical experiments beyond the quantum limit: Squeezing entanglement and teleportation

Results of experiments recently performed are reported, in which two optical parametric amplifiers were set up to generate two independently quadrature squeezed continuous wave laser beams. The transformation of quadrature squeezed states into polarization squeezed states and into states with spatial quantum correlations is demonstrated. By utilizing two squeezed laser beams, a polarization squeezed state exhibiting three simultaneously squeezed Stokes operator variances was generated. Continuous variable polarization entanglement was generated and the Einstein-Podolsky-Rosen paradox was observed. A pair of Stokes operators satisfied both the inseparability criterion and the conditional variance criterion. Values of 0.49 and 0.77, respectively, were observed, with entanglement requiring values below unity. The inseparability measure of the observed quadrature entanglement was 0.44. This value is sufficient for a demonstration of quantum teleportation, which is the next experimental goal of the authors.

Quantum theory of the far-off-resonance continuous-wave Raman laser: Heisenberg-Langevin approach

We present the quantum theory of the far-off-resonance continuous-wave Raman laser using the Heisenberg-Langevin approach. We show that the simplified quantum Langevin equations for this system are mathematically identical to those of the nondegenerate optical parametric oscillator in the time domain with the following associations: pump pump, Stokes signal, and Raman coherence idler. We derive analytical results for both the steady-state behavior and the time-dependent noise spectra, using standard linearization procedures. In the semiclassical limit, these results match with previous purely semiclassical treatments, which yield excellent agreement with experimental observations. The analytical time-dependent results predict perfect photon statistics conversion from the pump to the Stokes and nonclassical behavior under certain operational conditions.

Anomalous water absorption in porous materials

The absorption of fluid by unsaturated, rigid porous materials may be characterized by the sorptivity. This is a simple parameter to determine and is increasingly being used as a measure of a material's resistance to exposure to fluids (especially moisture and reactive solutes) in aggressive environments. The complete isothermal absorption process is described by a nonlinear diffusion equation, with the hydraulic diffusivity being a strongly nonlinear function of the degree of saturation of the material. This diffusivity can be estimated from the sorptivity test. In a typical test the cumulative absorption is proportional to the square root of time. However, a number of researchers have observed deviation from this behaviour when the infiltrating fluid is water and there is some potential for chemo-mechanical interaction with the material. In that case the current interpretation of the test and estimation of the hydraulic diffusivity is no longer appropriate. Kuntz and Lavallee (2001) discuss the anomalous behaviour and propose a non-Darcian model as a more appropriate physical description. We present an alternative Darcian explanation and theory that retrieves the earlier advantages of the simple sorptivity test in providing parametric information about the material's hydraulic properties and allowing simple predictive formulae for the wetting profile to be generated.

Análise da fibrose do tecido adiposo de mulheres submetidas à abdominoplastia e cirurgia bariátrica

A fibrose é um acúmulo demasiado de matriz extracelular, resultante de um desequilíbrio entre a síntese e a degradação dos seus componentes. É associada às alterações metabólicas do tecido adiposo, contudo sua ocorrência nos diferentes depósitos e repercussões clínicas ainda não são totalmente compreendidas. O objetivo deste estudo foi analisar a fibrose no tecido adiposo em relação à presença de obesidade, localização do depósito [tecido adiposo subcutâneo abdominal (TASA) e visceral (TAV)] e sua associação a variáveis clínicas. Amostras de gordura do TASA e TAV foram obtidas de 21 mulheres submetidas à cirurgia bariátrica (IMC>40Kg/m2) e 25 amostras de TASA das submetidas à abdominoplastia (IMC<30Kg/m2). As amostras foram processadas para histologia convencional. O corante picrosirius foi utilizado para avaliação das fibras colágenas totais. As imagens obtidas foram analisadas no ADIPOSOFT®. O percentual de fibrose no TASA e no TAV foi analisado com testes estatísticos não paramétricos, adotando-se um valor de p<0,05. A fibrose no TASA foi maior em mulheres com obesidade (p<0.0006). A fibrose entre os depósitos de TASA e de TAV foi observada apenas em mulheres pardas e negras com obesidade (p<0,012). A fibrose no TASA não foi correlacionada com as variáveis clínicas nas mulheres sem obesidade. No entanto, nas submetidas à cirurgia bariátrica, foram observadas correlações da fibrose no TASA com Índice de Massa Corpórea (IMC), hemoglobina glicada (A1c), LDL e triglicerídeos; e no TAV com porcentagem de perda gordura pré-operatório, % de perda de gordura total, % de massa magra pré, Taxa Metabólica Basal (TBM) e Gasto Energético Basal (GEB). Os parâmetros metabólicos e de perfil antropométrico antes da cirurgia bariátrica foram associados à fibrose no TASA, enquanto os parâmetros após a cirurgia foram associados à fibrose no TAV.

Eventos de vida em mulheres com diagnóstico de câncer de mama

Introdução: A relação entre eventos de vida e o surgimento e progressão do câncer de mama tem sido investigada por alguns estudos. Objetivo: Identificar os eventos de vida ocorridos em mulheres com diagnóstico de câncer de mama, examinar o tempo transcorrido entre o evento e o diagnóstico, examinar a associação entre a sobrecarga ocasionada pelo evento no momento da ocorrência e após o diagnóstico e examinar a associação entre metástase e eventos de vida pós diagnóstico de câncer de mama. Metodologia: Estudo transversal realizado no Hospital Santa Rita de Cássia, Vitória – ES. Compõe-se a amostra por 300 mulheres. Coletaram-se os dados no período de setembro a dezembro de 2014. Utilizou-se o instrumento Life Events Units – LEU/VAS que se baseia na Escala de Avaliação de Reajustamento Social de Holmes e Rahe, que no Brasil foi adaptada por Vasconcellos. Utilizou-se o Pacote Estatístico para Ciências Sociais (SPSS), versão 20.0, para calcular a frequência, média, mediana, desvio padrão e aplicar os testes não paramétrico de Wilcoxon e qui-quadrado. Resultados: A média de idade foi de 53 anos. Predominou-se mulheres de raça/cor não branca (65%), com menos de 8 anos de estudo (64%) e casadas (54%). Identificou-se que a maioria da amostra relatou pelo menos um evento de vida (99,3%). O principal evento de vida relatado foi morte de alguém na família. As medianas do tempo transcorrido entre os eventos de vida mais relatados e o diagnóstico de câncer de mama variaram de 5 a 15 anos. Observou-se diferença significante (p< 0,05) entre a sobrecarga ocasionada pelos eventos de vida nos dois momentos examinados. Em relação à metástase, 31,7% das que tiveram recidivas relataram um ou mais eventos de vida (p= 0,001). Quando considerado o tempo transcorrido entre o diagnóstico e o surgimento da metástase nas 46 mulheres, observou-se uma mediana de 18,0 meses. Conclusão: Os resultados deste estudo são potencialmente importantes, pois dão suporte a uma possível interação entre eventos de vida pós diagnóstico de câncer de mama e metástase. Estudos futuros são necessários para melhor compreensão desta relação.