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.

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.

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.

Multipartite entanglement and quantum state exchange

We investigate multipartite entanglement in relation to the process of quantum state exchange. In particular, we consider such entanglement for a certain pure state involving two groups of N trapped atoms. The state, which can be produced via quantum state exchange, is analogous to the steady-state intracavity state of the subthreshold optical nondegenerate parametric amplifier. We show that, first, it possesses some 2N-way entanglement. Second, we place a lower bound on the amount of such entanglement in the state using a measure called the entanglement of minimum bipartite entropy.

Mapping hippocampal and ventricular change in Alzheimer disease

We developed an anatomical mapping technique to detect hippocampal and ventricular changes in Alzheimer disease (AD). The resulting maps are sensitive to longitudinal changes in brain structure as the disease progresses. An anatomical surface modeling approach was combined with surface-based statistics to visualize the region and rate of atrophy in serial MRI scans and isolate where these changes link with cognitive decline. Fifty-two high-resolution MRI scans were acquired from 12 AD patients (age: 68.4 +/- 1.9 years) and 14 matched controls (age: 71.4 +/- 0.9 years), each scanned twice (2.1 +/- 0.4 years apart). 3D parametric mesh models of the hippocampus and temporal horns were created in sequential scans and averaged across subjects to identify systematic patterns of atrophy. As an index of radial atrophy, 3D distance fields were generated relating each anatomical surface point to a medial curve threading down the medial axis of each structure. Hippocampal atrophic rates and ventricular expansion were assessed statistically using surface-based permutation testing and were faster in AD than in controls. Using color-coded maps and video sequences, these changes were visualized as they progressed anatomically over time. Additional maps localized regions where atrophic changes linked with cognitive decline. Temporal horn expansion maps were more sensitive to AD progression than maps of hippocampal atrophy, but both maps correlated with clinical deterioration. These quantitative, dynamic visualizations of hippocampal atrophy and ventricular expansion rates in aging and AD may provide a promising measure to track AD progression in drug trials. (C) 2004 Elsevier Inc. All rights reserved.

Quantum nondemolition measurement of Fock states of mesoscopic mechanical oscillators

We investigate a scheme that makes a quantum nondemolition (QND) measurement of the excitation level of a mesoscopic mechanical oscillator by utilizing the anharmonic coupling between two beam bending modes. The nonlinear coupling between the two modes shifts the resonant frequency of the readout oscillator in proportion to the excitation level of the system oscillator. This frequency shift may be detected as a phase shift of the readout oscillation when driven on resonance. We derive an equation for the reduced density matrix of the system oscillator, and use this to study the conditions under which discrete jumps in the excitation level occur. The appearance of jumps in the actual quantity measured is also studied using the method of quantum trajectories. We consider the feasibility of the scheme for experimentally accessible parameters.

Spatial and temporal variability in somatic growth of green sea turtles (Chelonia mydas) resident in the Hawaiian Archipelago

The somatic growth dynamics of green turtles ( Chelonia mydas) resident in five separate foraging grounds within the Hawaiian Archipelago were assessed using a robust non-parametric regression modelling approach. The foraging grounds range from coral reef habitats at the north-western end of the archipelago, to coastal habitats around the main islands at the southeastern end of the archipelago. Pelagic juveniles recruit to these neritic foraging grounds from ca. 35 cm SCL or 5 kg ( similar to 6 years of age), but grow at foraging-ground-specific rates, which results in quite different size- and age-specific growth rate functions. Growth rates were estimated for the five populations as change in straight carapace length ( cm SCL year) 1) and, for two of the populations, also as change in body mass ( kg year) 1). Expected growth rates varied from ca. 0 - 2.5 cm SCL year) 1, depending on the foraging-ground population, which is indicative of slow growth and decades to sexual maturity, since expected size of first-time nesters is greater than or equal to 80 cm SCL. The expected size- specific growth rate functions for four populations sampled in the southeastern archipelago displayed a non-monotonic function, with an immature growth spurt at ca. 50 - 53 cm SCL ( similar to 18 - 23 kg) or ca. 13 - 19 years of age. The growth spurt for the Midway atoll population in the northwestern archipelago occurs at a much larger size ( ca. 65 cm SCL or 36 kg), because of slower immature growth rates that might be due to a limited food stock and cooler sea surface temperature. Expected age-at-maturity was estimated to be ca. 35 - 40 years for the four populations sampled at the south-eastern end of the archipelago, but it might well be > 50 years for the Midway population. The Hawaiian stock comprises mainly the same mtDNA haplotype, with no differences in mtDNA stock composition between foraging-ground populations, so that the geographic variability in somatic growth rates within the archipelago is more likely due to local environmental factors rather than genetic factors. Significant temporal variability was also evident, with expected growth rates declining over the last 10 - 20 years, while green turtle abundance within the archipelago has increased significantly since the mid-1970s. This inverse relationship between somatic growth rates and population abundance suggests a density-dependent effect on somatic growth dynamics that has also been reported recently for a Caribbean green turtle stock. The Hawaiian green turtle stock is characterised by slow growth rates displaying significant spatial and temporal variation and an immature growth spurt. This is consistent with similar findings for a Great Barrier Reef green turtle stock that also comprises many foraging-ground populations spanning a wide geographic range.

Ion trap transducers for quantum electromechanical oscillators

An enduring challenge for contemporary physics is to experimentally observe and control quantum behavior in macroscopic systems. We show that a single trapped atomic ion could be used to probe the quantum nature of a mesoscopic mechanical oscillator precooled to 4 K, and furthermore, to cool the oscillator with high efficiency to its quantum ground state. The proposed experiment could be performed using currently available technology.

The impact of changes to heroin supply on blood-borne virus notifications and injecting related harms in New South Wales, Australia

Background: In early 2001 Australia experienced a sudden and unexpected disruption to heroin availability, know as the 'heroin shortage'. This 'shortage has been linked to a decrease in needle and syringe output and therefore possibly a reduction in injecting drug use. We aimed to examine changes, if any, in blood-borne viral infections and presentations for injecting related problems related to injecting drug use following the reduction heroin availability in Australia, in the context of widespread harm reduction measures. Methods: Time series analysis of State level databases on HIV, hepatitis B, hepatitis C notifications and hospital and emergency department data. Examination of changes in HIV, hepatitis B, hepatitis C notifications and hospital and emergency department admissions for injection-related problems following the onset of the heroin shortage; non-parametric curve-fitting of number of hepatitis C notifications among those aged 15 - 19 years. Results: There were no changes observed in hospital visits for injection-related problems. There was no change related to the onset heroin shortage in the number of hepatitis C notifications among persons aged 15 - 19 years, but HCV notifications have subsequently decreased in this group. No change occurred in HIV and hepatitis B notifications. Conclusion: A marked reduction in heroin supply resulted in no increase in injection-related harm at the community level. However, a delayed decrease in HCV notifications among young people may be related. These changes occurred in a setting with widespread, publicly funded harm reduction initiatives.

Carotid artery intimal medial thickness, brachial artery flow-mediated vasodilation and cardiovascular risk factors in diabetic and non-diabetic indigenous Australians

Background: Indigenous Australians are at high risk for cardiovascular disease and type 2 diabetes. Carotid artery intimal medial thickness (CIMT) and brachial artery flow-mediated vasodilation (FMD) are ultrasound imaging based surrogate markers of cardiovascular risk. This study examines the relative contributions of traditional cardiovascular risk factors on CIMT and FMD in adult Indigenous Australians with and without type 2 diabetes mellitus. Method: One hundred and nineteen Indigenous Australians were recruited. Physical and biochemical markers of cardiovascular risk, together with CIMT and FMD were meausred for all subjects. Results: Fifty-three Indigenous Australians subjects (45%) had type 2 diabetes mellitus. There was a significantly greater mean CIMT in diabetic versus non-diabetic subjects (p = 0.049). In the non-diabetic group with non-parametric analyses, there were significant correlations between CIMT and: age (r = 0.64, p < 0.001), systolic blood pressure (r = 0.47, p < 0.001) and non-smokers (r = -0.30, p = 0.018). In the diabetic group, non-parametric analysis showed correlations between CIMT, age (r = 0.36, p = 0.009) and duration of diabetes (r = 0.30, p = 0.035) only. Adjusting forage, sex, smoking and history of cardiovascular disease, Hb(A1c) became the sole significant correlate of CIMT (r = 0.35,p = 0.01) in the diabetic group. In non-parametric analysis, age was the sole significant correlate of FMD (r = -0.31,p = 0.013), and only in non-diabetic subjects. Linear regression analysis showed significant associations between CIMT and age (t = 4.6,p < 0.001), systolic blood pressure (t = 2.6, p = 0.010) and Hb(A1c) (t = 2.6, p = 0.012), smoking (t = 2.1, p = 0.04) and fasting LDL-cholesterol (t = 2.1, p = 0.04). There were no significant associations between FMD and examined cardiovascular risk factors with linear regression analysis Conclusions: CIMT appears to be a useful surrogate marker of cardiovascular risk in this sample of Indigenous Australian subjects, correlating better than FMD with established cardiovascular risk factors. A lifestyle intervention programme may alleviate the burden of cardiovascular disease in Indigenous Australians by reducing central obesity, lowering blood pressure, correcting dyslipidaemia and improving glycaemic control. CIMT may prove to be a useful tool to assess efficacy of such an intervention programme. (c) 2004 Elsevier Ireland Ltd. All rights reserved.

Entanglement and the Einstein-Podolsky-Rosen paradox with coupled intracavity optical down-converters

We show that two evanescently coupled χ((2)) parametric down-converters inside a Fabry-Perot cavity provide a tunable source of quadrature squeezed light, Einstein-Podolsky-Rosen (EPR) correlations and quantum entanglement. Analyzing the operation in the below threshold regime, we show how these properties can be controlled by adjusting the coupling strengths and the cavity detunings. As this can be implemented with integrated optics, it provides a possible route to rugged and stable EPR sources.

Heavy tails, importance sampling and cross-entropy

We consider the problem of estimating P(Yi + (...) + Y-n > x) by importance sampling when the Yi are i.i.d. and heavy-tailed. The idea is to exploit the cross-entropy method as a toot for choosing good parameters in the importance sampling distribution; in doing so, we use the asymptotic description that given P(Y-1 + (...) + Y-n > x), n - 1 of the Yi have distribution F and one the conditional distribution of Y given Y > x. We show in some specific parametric examples (Pareto and Weibull) how this leads to precise answers which, as demonstrated numerically, are close to being variance minimal within the parametric class under consideration. Related problems for M/G/l and GI/G/l queues are also discussed.

Comment on "Evidence for quantized displacement in macroscopic nanomechanical oscillators"

A Comment on the Letter by Alexei Gaidarzhy et al., Phys. Rev. Lett. 94, 030402 (2005). The authors of the Letter offer a Reply.

Universality of quantum critical dynamics in a planar optical parametric oscillator

We analyze the critical quantum fluctuations in a coherently driven planar optical parametric oscillator. We show that the presence of transverse modes combined with quantum fluctuations changes the behavior of the quantum image critical point. This zero-temperature nonequilibrium quantum system has the same universality class as a finite-temperature magnetic Lifshitz transition.