Home-range size in juveniles of the temperate reef fish, the common triplefin (Forsterygion lapillum)

An organism’s home range dictates the spatial scale on which important processes occur (e.g. competition and predation) and directly affects the relationship between individual fitness and local habitat quality. Many reef fish species have very restricted home ranges after settlement and, here, we quantify home-range size in juveniles of a widespread and abundant reef fish in New Zealand, the common triplefin (Forsterygion lapillum). We conducted visual observations on 49 juveniles (mean size = 35-mm total length) within the Wellington harbour, New Zealand. Home ranges were extremely small, 0.053 m2 ± 0.029 (mean ± s.d.) and were unaffected by adult density, body size or substrate composition. A regression tree indicated that home-range size sharply decreased ~4.5 juveniles m–2 and a linear mixed model confirmed that home-range sizes in high-density areas (>4.5 juveniles m–2) were significantly smaller (34%) than those in low-density areas (after accounting for a significant effect of fish movement on our home-range estimates). Our results suggest that conspecific density may have negative and non-linear effects on home-range size, which could shape the spatial distribution of juveniles within a population, as well as influence individual fitness across local density gradients.

Outlier effects on robust joint modelling of longitudinal and survival data

Robust joint modelling is an emerging field of research. Through the advancements in electronic patient healthcare records, the popularly of joint modelling approaches has grown rapidly in recent years providing simultaneous analysis of longitudinal and survival data. This research advances previous work through the development of a novel robust joint modelling methodology for one of the most common types of standard joint models, that which links a linear mixed model with a Cox proportional hazards model. Through t-distributional assumptions, longitudinal outliers are accommodated with their detrimental impact being down weighed and thus providing more efficient and reliable estimates. The robust joint modelling technique and its major benefits are showcased through the analysis of Northern Irish end stage renal disease patients. With an ageing population and growing prevalence of chronic kidney disease within the United Kingdom, there is a pressing demand to investigate the detrimental relationship between the changing haemoglobin levels of haemodialysis patients and their survival. As outliers within the NI renal data were found to have significantly worse survival, identification of outlying individuals through robust joint modelling may aid nephrologists to improve patient's survival. A simulation study was also undertaken to explore the difference between robust and standard joint models in the presence of increasing proportions and extremity of longitudinal outliers. More efficient and reliable estimates were obtained by robust joint models with increasing contrast between the robust and standard joint models when a greater proportion of more extreme outliers are present. Through illustration of the gains in efficiency and reliability of parameters when outliers exist, the potential of robust joint modelling is evident. The research presented in this thesis highlights the benefits and stresses the need to utilise a more robust approach to joint modelling in the presence of longitudinal outliers.

Generalized compound transport of charged particles in turbulent magnetized plasmas

The transport of charged particles in partially turbulent magnetic systems is investigated from first principles. A generalized compound transport model is proposed, providing an explicit relation between the mean-square deviation of the particle parallel and perpendicular to a magnetic mean field, and the mean-square deviation which characterizes the stochastic field-line topology. The model is applied in various cases of study, and the relation to previous models is discussed.

Plant-wide predictive control for a thermal power plant based on a physical plant model

A constrained non-linear, physical model-based, predictive control (NPMPC) strategy is developed for improved plant-wide control of a thermal power plant. The strategy makes use of successive linearisation and recursive state estimation using extended Kalman filtering to obtain a linear state-space model. The linear model and a quadratic programming routine are used to design a constrained long-range predictive controller One special feature is the careful selection of a specific set of plant model parameters for online estimation, to account for time-varying system characteristics resulting from major system disturbances and ageing. These parameters act as nonstationary stochastic states and help to provide sufficient degrees-of-freedom to obtain unbiased estimates of controlled outputs. A 14th order non-linear plant model, simulating the dominant characteristics of a 200 MW oil-fired pou er plant has been used to test the NPMPC algorithm. The control strategy gives impressive simulation results, during large system disturbances and extremely high rate of load changes, right across the operating range. These results compare favourably to those obtained with the state-space GPC method designed under similar conditions.

Modelling the mechanical response of the reed-mouthpiece-lip system of a clarinet. Part II: A lumped model approximation

A non-linear lumped model of the reed-mouthpiece-lip system of a clarinet is formulated, in which the lumped parameters are derived from numerical experiments with a finite-difference simulation based on a distributed reed model. The effective stiffness per unit area is formulated as a function of the pressure signal driving the reed, in order to simulate the effects of the reed bending against the lay, and mass and damping terms are added as a first approximation to the dynamic behaviour of the reed. A discrete-time formulation is presented, and its response is compared to that of the distributed model. In addition, the lumped model is applied in the simulation of clarinet tones, enabling the analysis of the effects of using a pressure-dependent stiffness per unit area on sustained oscillations. The analysed effects and features are in qualitative agreement with players' experiences and experimental results obtained in prior studies.

Heritability and Genome-Wide Association Study to Assess Genetic Differences between Advanced Age-Related Macular Degeneration Subtypes

PURPOSE: To investigate whether the 2 subtypes of advanced age-related macular degeneration (AMD), choroidal neovascularization (CNV), and geographic atrophy (GA) segregate separately in families and to identify which genetic variants are associated with these 2 subtypes. DESIGN: Sibling correlation study and genome-wide association study (GWAS). PARTICIPANTS: For the sibling correlation study, 209 sibling pairs with advanced AMD were included. For the GWAS, 2594 participants with advanced AMD subtypes and 4134 controls were included. Replication cohorts included 5383 advanced AMD participants and 15 240 controls. METHODS: Participants had the AMD grade assigned based on fundus photography, examination, or both. To determine heritability of advanced AMD subtypes, a sibling correlation study was performed. For the GWAS, genome-wide genotyping was conducted and 6 036 699 single nucleotide polymorphisms (SNPs) were imputed. Then, the SNPs were analyzed with a generalized linear model controlling for genotyping platform and genetic ancestry. The most significant associations were evaluated in independent cohorts. MAIN OUTCOME MEASURES: Concordance of advanced AMD subtypes in sibling pairs and associations between SNPs with GA and CNV advanced AMD subtypes. RESULTS: The difference between the observed and expected proportion of siblings concordant for the same subtype of advanced AMD was different to a statistically significant degree (P = 4.2×10(-5)), meaning that in siblings of probands with CNV or GA, the same advanced subtype is more likely to develop. In the analysis comparing participants with CNV to those with GA, a statistically significant association was observed at the ARMS2/HTRA1 locus (rs10490924; odds ratio [OR], 1.47; P = 4.3×10(-9)), which was confirmed in the replication samples (OR, 1.38; P = 7.4×10(-14) for combined discovery and replication analysis). CONCLUSIONS: Whether CNV versus GA develops in a patient with AMD is determined in part by genetic variation. In this large GWAS meta-analysis and replication analysis, the ARMS2/HTRA1 locus confers increased risk for both advanced AMD subtypes, but imparts greater risk for CNV than for GA. This locus explains a small proportion of the excess sibling correlation for advanced AMD subtype. Other loci were detected with suggestive associations that differ for advanced AMD subtypes and deserve follow-up in additional studies. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

New generalized dispersion relation for low-frequency electromagnetic waves in Hall-magnetohydrodynamic dusty plasma

A generalized linear theory for electromagnetic waves in a homogeneous dusty magnetoplasma is presented. The waves described are characterized by a frequency which is much smaller (larger) than the electron gyrofrequency (dust plasma and dust gyrofrequencies), and a long wavelength (in comparison with the ion gyroradius and the electron skin depth). The generalized Hall- magnetohydrodynamic (GH-MHD) equations are derived by assuming massive charged dust macroparticles to be immobile, and Fourier transformed to obtain a general dispersion relation. The latter is analyzed to understand the influence of immobile charged dust grains on various electromagnetic wave modes in a magnetized dusty plasma.

Nonlinear Field Line Random Walk and Generalized Compound Diffusion of Charged Particles

A new nonlinear theory for the perpendicular transport of charged particles is presented. This approach is based on an improved nonlinear treatment of field line random walk in combination with a generalized compound diffusion model. The generalized compound diffusion model is much more systematic and reliable, in comparison to previous theories. Furthermore, the new theory shows remarkably good agreement with test-particle simulations and heliospheric observations.

Regression analysis of spatial data

Many of the most interesting questions ecologists ask lead to analyses of spatial data. Yet, perhaps confused by the large number of statistical models and fitting methods available, many ecologists seem to believe this is best left to specialists. Here, we describe the issues that need consideration when analysing spatial data and illustrate these using simulation studies. Our comparative analysis involves using methods including generalized least squares, spatial filters, wavelet revised models, conditional autoregressive models and generalized additive mixed models to estimate regression coefficients from synthetic but realistic data sets, including some which violate standard regression assumptions. We assess the performance of each method using two measures and using statistical error rates for model selection. Methods that performed well included generalized least squares family of models and a Bayesian implementation of the conditional auto-regressive model. Ordinary least squares also performed adequately in the absence of model selection, but had poorly controlled Type I error rates and so did not show the improvements in performance under model selection when using the above methods. Removing large-scale spatial trends in the response led to poor performance. These are empirical results; hence extrapolation of these findings to other situations should be performed cautiously. Nevertheless, our simulation-based approach provides much stronger evidence for comparative analysis than assessments based on single or small numbers of data sets, and should be considered a necessary foundation for statements of this type in future.

Role of Predalors, Habitat Attributes, and Spatial Autocorrelation on the Distribution of Eggs in the Northern Spectacled Salamander (Salamandrina perspicillata)

We predicted that the probability of egg occurrence of salamander Salamandrina perspicillata depended on stream features and predation by native crayfish Austropotamobius fulcisianus and the introduced trout Salmo trutta. We assessed the presence of S. perspicillata at 54 sites within a natural reserve of southern Tuscany, Italy. Generalized linear models with binomial errors were constructed using egg presence/absence and altitude, stream mean size and slope, electrical conductivity, water pH and temperature, and a predation factor, defined according to the presence/absence of crayfish and trout. Some competing models also included an autocovariate term, which estimated how much the response variable at any one sampling point reflected response values at surrounding points. The resulting models were compared using Akaike's information criterion. Model selection led to a subset of 14 models with Delta AIC(c) <7 (i.e., models ranging from substantial support to considerably less support), and all but one of these included an effect of predation. Models with the autocovariate term had considerably more support than those without the term. According to multimodel inference, the presence of trout and crayfish reduced the probability of egg occurrence from a mean level of 0.90 (SE limits: 0.98-0.55) to 0.12 (SE limits: 0.34-0.04). The presence of crayfish alone had no detectable effects (SE limits: 0.86-0.39). The results suggest that introduced trout have a detrimental effect on the reproductive output of S. perspicillata and confirm the fundamental importance of distinguishing the roles of endogenous and exogenous forces that act on population distribution.

Large-scale spatial patterns in the distribution of Collembola (Hexapoda) species in Antarctic terrestrial ecosystems

We tested whether the distribution of three common springtail species (Gressittacantha terranova, Gomphiocephalus hodgsoni and Friesea grisea) in Victoria Land (Antarctica) could be modelled as a function of latitude, longitude, altitude and distance from the sea.

Victoria Land, Ross Dependency, Antarctica.

Generalized linear models were constructed using species presence/absence data relative to geographical features (latitude, longitude, altitude, distance from sea) across the species' entire ranges. Model results were then integrated with the known phylogeography of each species and hypotheses were generated on the role of climate as a major driver of Antarctic springtail distribution.

Based on model selection using Akaike's information criterion, the species' distributions were: hump-shaped relative to longitude and monotonic with altitude for Gressittacantha terranova; hump-shaped relative to latitude and monotonic with altitude for Gomphiocephalus hodgsoni; and hump-shaped relative to longitude and monotonic with latitude, altitude and distance from the sea for Friesea grisea.

No single distributional pattern was shared by the three species. While distributions were partially a response to climatic spatial clines, the patterns observed strongly suggest that past geological events have influenced the observed distributions. Accordingly, present-day spatial patterns are likely to have arisen from the interaction of historical and environmental drivers. Future studies will need to integrate a range of spatial and temporal scales to further quantify their respective roles.

Population pharmacokinetic and pharmacogenetic analysis of tacrolimus in paediatric liver transplant patients

Aims: To build a population pharmacokinetic model that describes the apparent clearance of tacrolimus and the potential demographic, clinical and genetically controlled factors that could lead to inter-patient pharmacokinetic variability within children following liver transplantation.

Methods: The present study retrospectively examined tacrolimus whole blood pre-dose concentrations (n = 628) of 43 children during their first year post-liver transplantation. Population pharmacokinetic analysis was performed using the non-linear mixed effects modelling program (nonmem) to determine the population mean parameter estimate of clearance and influential covariates.

Results: The final model identified time post-transplantation and CYP3A5*1 allele as influential covariates on tacrolimus apparent clearance according to the following equation:

TVCL=12.9×(Weight /13.2)^0.75×EXP(-0.00158×TPT)×EXP(0.428×CYP3A5)

where TVCL is the typical value for apparent clearance, TPT is time post-transplantation in days and the CYP3A5 is 1 where*1 allele is present and 0 otherwise. The population estimate and inter-individual variability (%CV) of tacrolimus apparent clearance were found to be 0.977 l h kg (95% CI 0.958, 0.996) and 40.0%, respectively, while the residual variability between the observed and predicted concentrations was 35.4%.

Conclusion: Tacrolimus apparent clearance was influenced by time post-transplantation and CYP3A5 genotypes. The results of this study, once confirmed by a large scale prospective study, can be used in conjunction with therapeutic drug monitoring to recommend tacrolimus dose adjustments that take into account not only body weight but also genetic and time-related changes in tacrolimus clearance. © 2013 The British Pharmacological Society.

An Observational Study of Blood Concentrations and Kinetics of Methyl- and Propyl-Parabens in Neonates

Purpose: Systemic exposure to parabens in the neonatal population, in particular propyl-parabens (PPB), remains a concern. Blood concentrations and kinetics of methyl-parabens (MPB) and PPB were therefore determined in neonates receiving medicines containing these excipients.

Methods: A multi-centre, non-interventional, observational study of excipient-kinetics in neonates. ‘Dried Blood Spot’ samples were collected opportunistically at the same time as routine samples and the observations modelled using a non-linear mixed effects approach.

Results: A total of 841 blood MPB and PPB concentration data were available for evaluation from 181 pre- and term-neonates. Quantifiable blood concentrations of MPB and PPB were observed in 99% and 49% of patients, and 55% and 25% of all concentrations were above limit of detection (10 ng/ml), respectively. Only MPB data was amenable to modelling. Oral bioavailability was influenced by type of formulation and disposition was best described by a two compartment model with clearance (CL) influenced by post natal age (PNA); CL_{PNA<21 days} 0.57 versus CL_PNA>21days 0.88 L/h.

Conclusions: Daily repeated administration of parabens containing medicines can result in prolonged systemic exposure to the parent compound in neonates. Animal toxicology studies of PPB that specifically address the neonatal period are required before a permitted daily exposure for this age group can be established.

PREDICTING THE SPATIAL-DISTRIBUTION OF CLIMATE - TEMPERATURE IN GREAT-BRITAIN

1. The prediction and mapping of climate in areas between climate stations is of increasing importance in ecology.

2. Four categories of model, simple interpolation, thin plate splines, multiple linear regression and mixed spline-regression, were tested for their ability to predict the spatial distribution of temperature on the British mainland. The models were tested by external cross-verification.

3. The British distribution of mean daily temperature was predicted with the greatest accuracy by using a mixed model: a thin plate spline fitted to the surface of the country, after correction of the data by a selection from 16 independent topographical variables (such as altitude, distance from the sea, slope and topographic roughness), chosen by multiple regression from a digital terrain model (DTM) of the country.

4. The next most accurate method was a pure multiple regression model using the DTM. Both regression and thin plate spline models based on a few variables (latitude, longitude and altitude) only were comparatively unsatisfactory, but some rather simple methods of surface interpolation (such as bilinear interpolation after correction to sea level) gave moderately satisfactory results. Differences between the methods seemed to be dependent largely on their ability to model the effect of the sea on land temperatures.

5. Prediction of temperature by the best methods was greater than 95% accurate in all months of the year, as shown by the correlation between the predicted and actual values. The predicted temperatures were calculated at real altitudes, not subject to sea-level correction.

6. A minimum of just over 30 temperature recording stations would generate a satisfactory surface, provided the stations were well spaced.

7. Maps of mean daily temperature, using the best overall methods are provided; further important variables, such as continentality and length of growing season, were also mapped. Many of these are believed to be the first detailed representations at real altitude.

8. The interpolated monthly temperature surfaces are available on disk.

An integrated approach for real-time model-based state-of-charge estimation of lithium-ion batteries

Lithium-ion batteries have been widely adopted in electric vehicles (EVs), and accurate state of charge (SOC) estimation is of paramount importance for the EV battery management system. Though a number of methods have been proposed, the SOC estimation for Lithium-ion batteries, such as LiFePo4 battery, however, faces two key challenges: the flat open circuit voltage (OCV) vs SOC relationship for some SOC ranges and the hysteresis effect. To address these problems, an integrated approach for real-time model-based SOC estimation of Lithium-ion batteries is proposed in this paper. Firstly, an auto-regression model is adopted to reproduce the battery terminal behaviour, combined with a non-linear complementary model to capture the hysteresis effect. The model parameters, including linear parameters and non-linear parameters, are optimized off-line using a hybrid optimization method that combines a meta-heuristic method (i.e., the teaching learning based optimization method) and the least square method. Secondly, using the trained model, two real-time model-based SOC estimation methods are presented, one based on the real-time battery OCV regression model achieved through weighted recursive least square method, and the other based on the state estimation using the extended Kalman filter method (EKF). To tackle the problem caused by the flat OCV-vs-SOC segments when the OCV-based SOC estimation method is adopted, a method combining the coulombic counting and the OCV-based method is proposed. Finally, modelling results and SOC estimation results are presented and analysed using the data collected from LiFePo4 battery cell. The results confirmed the effectiveness of the proposed approach, in particular the joint-EKF method.