Sparse kernel density estimation technique based on zero-norm constraint

A sparse kernel density estimator is derived based on the zero-norm constraint, in which the zero-norm of the kernel weights is incorporated to enhance model sparsity. The classical Parzen window estimate is adopted as the desired response for density estimation, and an approximate function of the zero-norm is used for achieving mathemtical tractability and algorithmic efficiency. Under the mild condition of the positive definite design matrix, the kernel weights of the proposed density estimator based on the zero-norm approximation can be obtained using the multiplicative nonnegative quadratic programming algorithm. Using the -optimality based selection algorithm as the preprocessing to select a small significant subset design matrix, the proposed zero-norm based approach offers an effective means for constructing very sparse kernel density estimates with excellent generalisation performance.

Quasi-optimal assimilation of satellite retrievals for incremental 4D-Var

An interface between satellite retrievals and the incremental version of the four-dimensional variational assimilation scheme is developed, making full use of the information content of satellite measurements. In this paper, expressions for the function that calculates simulated observations from model states (called “observation operator”), together with its tangent linear version and adjoint, are derived. Results from our work can be used for implementing a quasi-optimal assimilation of satellite retrievals (e.g., of atmospheric trace gases) in operational meteorological centres.

Inter-relationships between small, dense low-density lipoprotein (LDL), plasma triacylglycerol and LDL apoprotein B in an atherogenic lipoprotein phenotype in free-living subjects

A predominance of small, dense low-density lipoprotein (LDL) is a major component of an atherogenic lipoprotein phenotype, and a common, but modifiable, source of increased risk for coronary heart disease in the free-living population. While much of the atherogenicity of small, dense LDL is known to arise from its structural properties, the extent to which an increase in the number of small, dense LDL particles (hyper-apoprotein B) contributes to this risk of coronary heart disease is currently unknown. This study reports a method for the recruitment of free-living individuals with an atherogenic lipoprotein phenotype for a fish-oil intervention trial, and critically evaluates the relationship between LDL particle number and the predominance of small, dense LDL. In this group, volunteers were selected through local general practices on the basis of a moderately raised plasma triacylglycerol (triglyceride) level (>1.5 mmol/l) and a low concentration of high-density-lipoprotein cholesterol (<1.1 mmol/l). The screening of LDL subclasses revealed a predominance of small, dense LDL (LDL subclass pattern B) in 62% of the cohort. As expected, subjects with LDL subclass pattern B were characterized by higher plasma triacylglycerol and lower high-density lipoprotein cholesterol (<1.1 mmol/l) levels and, less predictably, by lower LDL cholesterol and apoprotein B levels (P<0.05; LDL subclass A compared with subclass B). While hyper-apoprotein B was detected in only five subjects, the relative percentage of small, dense LDL-III in subjects with subclass B showed an inverse relationship with LDL apoprotein B (r=-0.57; P<0.001), identifying a subset of individuals with plasma triacylglycerol above 2.5 mmol/l and a low concentration of LDL almost exclusively in a small and dense form. These findings indicate that a predominance of small, dense LDL and hyper-apoprotein B do not always co-exist in free-living groups. Moreover, if coronary risk increases with increasing LDL particle number, these results imply that the risk arising from a predominance of small, dense LDL may actually be reduced in certain cases when plasma triacylglycerol exceeds 2.5 mmol/l.

Variation in the FFAR1 gene modifies BMI, body composition and plasma lipids but not plasma insulin or Beta-cell function in overweight subjects

Background FFAR1 receptor is a long chain fatty acid G-protein coupled receptor which is expressed widely, but found in high density in the pancreas and central nervous system. It has been suggested that FFAR1 may play a role in insulin sensitivity, lipotoxicity and is associated with type 2 diabetes. Here we investigate the effect of three common SNPs of FFAR1 (rs2301151; rs16970264; rs1573611) on pancreatic function, BMI, body composition and plasma lipids. Methodology/Principal Findings For this enquiry we used the baseline RISCK data, which provides a cohort of overweight subjects at increased cardiometabolic risk with detailed phenotyping. The key findings were SNPs of the FFAR1 gene region were associated with differences in body composition and lipids, and the effects of the 3 SNPs combined were cumulative on BMI, body composition and total cholesterol. The effects on BMI and body fat were predominantly mediated by rs1573611 (1.06 kg/m2 higher (P = 0.009) BMI and 1.53% higher (P = 0.002) body fat per C allele). Differences in plasma lipids were also associated with the BMI-increasing allele of rs2301151 including higher total cholesterol (0.2 mmol/L per G allele, P = 0.01) and with the variant A allele of rs16970264 associated with lower total (0.3 mmol/L, P = 0.02) and LDL (0.2 mmol/L, P<0.05) cholesterol, but also with lower HDL-cholesterol (0.09 mmol/L, P<0.05) although the difference was not apparent when controlling for multiple testing. There were no statistically significant effects of the three SNPs on insulin sensitivity or beta cell function. However accumulated risk allele showed a lower beta cell function on increasing plasma fatty acids with a carbon chain greater than six. Conclusions/Significance Differences in body composition and lipids associated with common SNPs in the FFAR1 gene were apparently not mediated by changes in insulin sensitivity or beta-cell function.

Photochemical cross-linking of plastically compressed collagen gel produces an optimal scaffold for corneal tissue engineering

The experiments were designed to use photochemically cross-linked plastically compressed collagen (PCPCC) gel to support corneal epithelial cells. A plastically compressed collagen (PCC) scaffold was photo cross-linked by UVA in the presence of riboflavin to form a biomaterial with optimal mechanical properties. The breaking force, rheology, surgical suture strength, transparency, ultrastructure, and cell-based biocompatibility were compared between PCPCC and PCC gels. The breaking force increased proportionally with an increased concentration of riboflavin. The stress required to reach breaking point of the PCPCC scaffolds was over two times higher compared to the stress necessary to break PCC scaffolds in the presence of 0.1% riboflavin. Rheology results indicated that the structural properties of PCC remain unaltered after UVA cross-linking. The PCC gels were more easily broken than PCPCC gels when sutured on to bovine corneas. The optical density values of PCPCC and PCC showed no significant differences (p > 0.05). SEM analyses showed that the collagen fibres within the PCPCC gels were similar in morphology to PCC gels. No difference in cell-based biocompatibility was seen between the PCPCC and PCC scaffolds in terms of their ability to support the ex vivo expansion of corneal epithelial cells or their subsequent differentiation evidenced by similar levels of cytokeratin 14. In conclusion, PCPCC scaffold is an optimal biomaterial for use in therapeutic tissue engineering of the cornea.

A tutorial on pseudospectral methods for computational optimal control

[English] This paper is a tutorial introduction to pseudospectral optimal control. With pseudospectral methods, a function is approximated as a linear combination of smooth basis functions, which are often chosen to be Legendre or Chebyshev polynomials. Collocation of the differential-algebraic equations is performed at orthogonal collocation points, which are selected to yield interpolation of high accuracy. Pseudospectral methods directly discretize the original optimal control problem to recast it into a nonlinear programming format. A numerical optimizer is then employed to find approximate local optimal solutions. The paper also briefly describes the functionality and implementation of PSOPT, an open source software package written in C++ that employs pseudospectral discretization methods to solve multi-phase optimal control problems. The software implements the Legendre and Chebyshev pseudospectral methods, and it has useful features such as automatic differentiation, sparsity detection, and automatic scaling. The use of pseudospectral methods is illustrated in two problems taken from the literature on computational optimal control. [Portuguese] Este artigo e um tutorial introdutorio sobre controle otimo pseudo-espectral. Em metodos pseudo-espectrais, uma funcao e aproximada como uma combinacao linear de funcoes de base suaves, tipicamente escolhidas como polinomios de Legendre ou Chebyshev. A colocacao de equacoes algebrico-diferenciais e realizada em pontos de colocacao ortogonal, que sao selecionados de modo a minimizar o erro de interpolacao. Metodos pseudoespectrais discretizam o problema de controle otimo original de modo a converte-lo em um problema de programa cao nao-linear. Um otimizador numerico e entao empregado para obter solucoes localmente otimas. Este artigo tambem descreve sucintamente a funcionalidade e a implementacao de um pacote computacional de codigo aberto escrito em C++ chamado PSOPT. Tal pacote emprega metodos de discretizacao pseudo-spectrais para resolver problemas de controle otimo com multiplas fase. O PSOPT permite a utilizacao de metodos de Legendre ou Chebyshev, e possui caractersticas uteis tais como diferenciacao automatica, deteccao de esparsidade e escalonamento automatico. O uso de metodos pseudo-espectrais e ilustrado em dois problemas retirados da literatura de controle otimo computacional.

Retrieval of mesospheric electron densities using an optimal estimation inverse method

We present a new method to determine mesospheric electron densities from partially reflected medium frequency radar pulses. The technique uses an optimal estimation inverse method and retrieves both an electron density profile and a gradient electron density profile. As well as accounting for the absorption of the two magnetoionic modes formed by ionospheric birefringence of each radar pulse, the forward model of the retrieval parameterises possible Fresnel scatter of each mode by fine electronic structure, phase changes of each mode due to Faraday rotation and the dependence of the amplitudes of the backscattered modes upon pulse width. Validation results indicate that known profiles can be retrieved and that χ2 tests upon retrieval parameters satisfy validity criteria. Application to measurements shows that retrieved electron density profiles are consistent with accepted ideas about seasonal variability of electron densities and their dependence upon nitric oxide production and transport.

Oxidation of low density lipoprotein by iron at lysosomal pH: implications for atherosclerosis

Low density lipoprotein (LDL) has recently been shown to be oxidised by iron within the lysosomes of macrophages and this is a novel potential mechanism for LDL oxidation in atherosclerosis. Our aim was to characterise the chemical and physical changes induced in LDL by iron at lysosomal pH and to investigate the effects of iron chelators and α-tocopherol on this process. LDL was oxidised by iron at pH 4.5 and 37°C and its oxidation monitored by spectrophotometry and HPLC. LDL was oxidised effectively by FeSO4 (5-50 µM) and became highly aggregated at pH 4.5, but not at pH 7.4. Cholesteryl esters decreased and after a pronounced lag 7-ketocholesterol increased greatly. Total hydroperoxides (measured by tri-iodide assay) increased up to 24 h and then decreased only slowly. The lipid composition after 12 h at pH 4.5 and 37°C was similar to that of LDL oxidised by copper at pH 7.4 and 4°C, i.e. rich in hydroperoxides but low in oxysterols. Previously oxidised LDL aggregated rapidly and spontaneously at pH 4.5, but not at pH 7.4. Ferrous was much more effective than ferric iron at oxidising LDL when added after the oxidation was already underway. The iron chelators diethylenetriaminepentaacetic acid and, to a lesser extent, desferrioxamine inhibited LDL oxidation when added during its initial stages, but were unable to prevent LDL aggregating after it had been partially oxidised. Surprisingly, desferrioxamine increased the rate of LDL modification when added late in the oxidation process. α-Tocopherol enrichment of LDL initially increased the oxidation of LDL, but inhibited it later. The presence of oxidised and highly aggregated lipid within lysosomes has the potential to perturb the function of these organelles and to promote atherosclerosis.

Effects of oxidised low density lipoprotein on dendritic cells: a possible immunoregulatory component of the atherogenic micro-environment?

Objective: The objective of this study was to explore the relationship between low density lipoprotein (LDL) and dendritic cell (DC) activation, based upon the hypothesis that reactive oxygen species (ROS)-mediated modification of proteins that may be present in local DC microenvironments could be important as mediators of this activation. Although LDL are known to be oxidised in vivo, and taken up by macrophages during atherogenesis; their effect on DC has not been explored previously. Methods: Human DCs were prepared from peripheral blood monocytes using GM-CSF and IL-4. Plasma LDLs were isolated by sequential gradient centrifugation, oxidised in CuSO4, and oxidation arrested to yield mild, moderate and highly oxidised LDL forms. DCs exposed to these LDLs were investigated using combined phenotypic, functional (autologous T cell activation), morphological and viability assays. Results: Highly-oxidised LDL increased DC HLA-DR, CD40 and CD86 expression, corroborated by increased DC-induced T cell proliferation. Both native and oxidised LDL induced prominent DC clustering. However, high concentrations of highly-oxidised LDL inhibited DC function, due to increased DC apoptosis. Conclusions: This study supports the hypothesis that oxidised LDL are capable of triggering the transition from sentinel to messenger DC. Furthermore, the DC clustering–activation–apoptosis sequence in the presence of different LDL forms is consistent with a regulatory DC role in immunopathogenesis of atheroma. A sequence of initial accumulation of DC, increasing LDL oxidation, and DC-induced T cell activation, may explain why local breach of tolerance can occur. Above a threshold level, however, supervening DC apoptosis limits this, contributing instead to the central plaque core.

Ethnic differences in beta-cell function, dietary intake and expression of the metabolic syndrome among UK adults of south Asian, black African-Caribbean and white-European origin at high risk of metabolic syndrome

A cross-sectional analysis of ethnic differences in dietary intake, insulin sensitivity and beta-cell function, using the intravenous glucose tolerance test (IVGTT), was conducted on 497 healthy adult participants of the ‘Reading, Imperial, Surrey, Cambridge, and Kings’ (RISCK) study. Insulin sensitivity (Si) was significantly lower in African-Caribbean (AC) and South Asian (SA) participants [IVGTT-Si; AC: 2.13 vs SA: 2.25 vs white-European (WE): 2.84 (×10−4 mL µU min)2, p < 0.001]. AC participants had a higher prevalence of anti-hypertensive therapy (AC: 19.7% vs SA: 7.5%), the most cardioprotective lipid profile [total:high-density lipoprotein (HDL); AC: 3.52 vs SA: 4.08 vs WE: 3.83, p = 0.03] and more pronounced hyperinsulinaemia [IVGTT–acute insulin response (AIR)] [AC: 575 vs SA: 428 vs WE: 344 mL/µU/min)2, p = 0.002], specifically in female participants. Intake of saturated fat and carbohydrate was lower and higher in AC (10.9% and 50.4%) and SA (11.1% and 52.3%), respectively, compared to WE (13.6% and 43.8%, p < 0.001). Insulin resistance in ACs is characterised by ‘normal’ lipid profiles but high rates of hypertension and pronounced hyperinsulinaemia.

Estimation of optimal gravity wave parameters for climate models using data assimilation

There is a current need to constrain the parameters of gravity wave drag (GWD) schemes in climate models using observational information instead of tuning them subjectively. In this work, an inverse technique is developed using data assimilation principles to estimate gravity wave parameters. Because mostGWDschemes assume instantaneous vertical propagation of gravity waves within a column, observations in a single column can be used to formulate a one-dimensional assimilation problem to estimate the unknown parameters. We define a cost function that measures the differences between the unresolved drag inferred from observations (referred to here as the ‘observed’ GWD) and the GWD calculated with a parametrisation scheme. The geometry of the cost function presents some difficulties, including multiple minima and ill-conditioning because of the non-independence of the gravity wave parameters. To overcome these difficulties we propose a genetic algorithm to minimize the cost function, which provides a robust parameter estimation over a broad range of prescribed ‘true’ parameters. When real experiments using an independent estimate of the ‘observed’ GWD are performed, physically unrealistic values of the parameters can result due to the non-independence of the parameters. However, by constraining one of the parameters to lie within a physically realistic range, this degeneracy is broken and the other parameters are also found to lie within physically realistic ranges. This argues for the essential physical self-consistency of the gravity wave scheme. A much better fit to the observed GWD at high latitudes is obtained when the parameters are allowed to vary with latitude. However, a close fit can be obtained either in the upper or the lower part of the profiles, but not in both at the same time. This result is a consequence of assuming an isotropic launch spectrum. The changes of sign in theGWDfound in the tropical lower stratosphere, which are associated with part of the quasi-biennial oscillation forcing, cannot be captured by the parametrisation with optimal parameters.

Optimal frame length for keeping normalized goodput with lowest requirement on BER

This paper presents an adaptive frame length mechanism based on a cross-layer analysis of intrinsic relations between the MAC frame length, bit error rate (BER) of the wireless link and normalized goodput. The proposed mechanism selects the optimal frame length that keeps the service normalized goodput at required levels while satisfying the lowest requirement on the BER, thus increasing the transmission reliability. Numerical results are provided and show that an optimal frame length satisfying the lowest BER requirement does indeed exist. The performance of BER requirement as a function of the MAC frame length is evaluated and compared for transmission scenarios with and without automatic repeat request (ARQ). Furthermore, issues related to the MAC overhead length are also discussed to illuminate the functionality and performance of the proposed mechanism.

Surface water temperature observations of large lakes by optimal estimation

Optimal estimation (OE) and probabilistic cloud screening were developed to provide lake surface water temperature (LSWT) estimates from the series of (advanced) along-track scanning radiometers (ATSRs). Variations in physical properties such as elevation, salinity, and atmospheric conditions are accounted for through the forward modelling of observed radiances. Therefore, the OE retrieval scheme developed is generic (i.e., applicable to all lakes). LSWTs were obtained for 258 of Earth's largest lakes from ATSR-2 and AATSR imagery from 1995 to 2009. Comparison to in situ observations from several lakes yields satellite in situ differences of −0.2 ± 0.7 K for daytime and −0.1 ± 0.5 K for nighttime observations (mean ± standard deviation). This compares with −0.05 ± 0.8 K for daytime and −0.1 ± 0.9 K for nighttime observations for previous methods based on operational sea surface temperature algorithms. The new approach also increases coverage (reducing misclassification of clear sky as cloud) and exhibits greater consistency between retrievals using different channel–view combinations. Empirical orthogonal function (EOF) techniques were applied to the LSWT retrievals (which contain gaps due to cloud cover) to reconstruct spatially and temporally complete time series of LSWT. The new LSWT observations and the EOF-based reconstructions offer benefits to numerical weather prediction, lake model validation, and improve our knowledge of the climatology of lakes globally. Both observations and reconstructions are publically available from http://hdl.handle.net/10283/88.

Density functional theory study of rutile VO2 surfaces

We present the results of a density functional theory (DFT) investigation of the surfaces of rutile-like vanadium dioxide, VO2(R). We calculate the surface energies of low Miller index planes, and find that the most stable surface orientation is the (110). The equilibrium morphology of a VO2(R) particle has an acicular shape, laterally confined by (110) planes and topped by (011) planes. The redox properties of the (110) surface are investigated by calculating the relative surface free energies of the non-stoichiometric compositions as a function of oxygen chemical potential. It is found that the VO2(110) surface is oxidized with respect to the stoichiometric composition, not only at ambient conditions but also at the more reducing conditions under which bulk VO2 is stable in comparison with bulk V2O5. The adsorbed oxygen forms surface vanadyl species much more favorably than surface peroxo species.

On the application of optimal wavelet filter banks for ECG signal classification

This paper discusses ECG signal classification after parametrizing the ECG waveforms in the wavelet domain. Signal decomposition using perfect reconstruction quadrature mirror filter banks can provide a very parsimonious representation of ECG signals. In the current work, the filter parameters are adjusted by a numerical optimization algorithm in order to minimize a cost function associated to the filter cut-off sharpness. The goal consists of achieving a better compromise between frequency selectivity and time resolution at each decomposition level than standard orthogonal filter banks such as those of the Daubechies and Coiflet families. Our aim is to optimally decompose the signals in the wavelet domain so that they can be subsequently used as inputs for training to a neural network classifier.