A modified PNN algorithm with optimal PD modeling using the orthogonal least squares method

In this paper a modified algorithm is suggested for developing polynomial neural network (PNN) models. Optimal partial description (PD) modeling is introduced at each layer of the PNN expansion, a task accomplished using the orthogonal least squares (OLS) method. Based on the initial PD models determined by the polynomial order and the number of PD inputs, OLS selects the most significant regressor terms reducing the output error variance. The method produces PNN models exhibiting a high level of accuracy and superior generalization capabilities. Additionally, parsimonious models are obtained comprising a considerably smaller number of parameters compared to the ones generated by means of the conventional PNN algorithm. Three benchmark examples are elaborated, including modeling of the gas furnace process as well as the iris and wine classification problems. Extensive simulation results and comparison with other methods in the literature, demonstrate the effectiveness of the suggested modeling approach.

An optimal inverse method using Doppler lidar measurements to estimate the surface sensible heat flux

Inverse methods are widely used in various fields of atmospheric science. However, such methods are not commonly used within the boundary-layer community, where robust observations of surface fluxes are a particular concern. We present a new technique for deriving surface sensible heat fluxes from boundary-layer turbulence observations using an inverse method. Doppler lidar observations of vertical velocity variance are combined with two well-known mixed-layer scaling forward models for a convective boundary layer (CBL). The inverse method is validated using large-eddy simulations of a CBL with increasing wind speed. The majority of the estimated heat fluxes agree within error with the proscribed heat flux, across all wind speeds tested. The method is then applied to Doppler lidar data from the Chilbolton Observatory, UK. Heat fluxes are compared with those from a mast-mounted sonic anemometer. Errors in estimated heat fluxes are on average 18 %, an improvement on previous techniques. However, a significant negative bias is observed (on average −63%) that is more pronounced in the morning. Results are improved for the fully-developed CBL later in the day, which suggests that the bias is largely related to the choice of forward model, which is kept deliberately simple for this study. Overall, the inverse method provided reasonable flux estimates for the simple case of a CBL. Results shown here demonstrate that this method has promise in utilizing ground-based remote sensing to derive surface fluxes. Extension of the method is relatively straight-forward, and could include more complex forward models, or other measurements.

A Variance-Minimizing Filter for Large-Scale Applications

A truly variance-minimizing filter is introduced and its per for mance is demonstrated with the Korteweg– DeV ries (KdV) equation and with a multilayer quasigeostrophic model of the ocean area around South Africa. It is recalled that Kalman-like filters are not variance minimizing for nonlinear model dynamics and that four - dimensional variational data assimilation (4DV AR)-like methods relying on per fect model dynamics have dif- ficulty with providing error estimates. The new method does not have these drawbacks. In fact, it combines advantages from both methods in that it does provide error estimates while automatically having balanced states after analysis, without extra computations. It is based on ensemble or Monte Carlo integrations to simulate the probability density of the model evolution. When obser vations are available, the so-called importance resampling algorithm is applied. From Bayes’ s theorem it follows that each ensemble member receives a new weight dependent on its ‘ ‘distance’ ’ t o the obser vations. Because the weights are strongly var ying, a resampling of the ensemble is necessar y. This resampling is done such that members with high weights are duplicated according to their weights, while low-weight members are largely ignored. In passing, it is noted that data assimilation is not an inverse problem by nature, although it can be for mulated that way . Also, it is shown that the posterior variance can be larger than the prior if the usual Gaussian framework is set aside. However , i n the examples presented here, the entropy of the probability densities is decreasing. The application to the ocean area around South Africa, gover ned by strongly nonlinear dynamics, shows that the method is working satisfactorily . The strong and weak points of the method are discussed and possible improvements are proposed.

Estimation of genetic variance for macro- and micro-environmental sensitivity using double hierarchical generalized linear models

Background: Genetic variation for environmental sensitivity indicates that animals are genetically different in their response to environmental factors. Environmental factors are either identifiable (e.g. temperature) and called macro-environmental or unknown and called micro-environmental. The objectives of this study were to develop a statistical method to estimate genetic parameters for macro- and micro-environmental sensitivities simultaneously, to investigate bias and precision of resulting estimates of genetic parameters and to develop and evaluate use of Akaike’s information criterion using h-likelihood to select the best fitting model. Methods: We assumed that genetic variation in macro- and micro-environmental sensitivities is expressed as genetic variance in the slope of a linear reaction norm and environmental variance, respectively. A reaction norm model to estimate genetic variance for macro-environmental sensitivity was combined with a structural model for residual variance to estimate genetic variance for micro-environmental sensitivity using a double hierarchical generalized linear model in ASReml. Akaike’s information criterion was constructed as model selection criterion using approximated h-likelihood. Populations of sires with large half-sib offspring groups were simulated to investigate bias and precision of estimated genetic parameters. Results: Designs with 100 sires, each with at least 100 offspring, are required to have standard deviations of estimated variances lower than 50% of the true value. When the number of offspring increased, standard deviations of estimates across replicates decreased substantially, especially for genetic variances of macro- and micro-environmental sensitivities. Standard deviations of estimated genetic correlations across replicates were quite large (between 0.1 and 0.4), especially when sires had few offspring. Practically, no bias was observed for estimates of any of the parameters. Using Akaike’s information criterion the true genetic model was selected as the best statistical model in at least 90% of 100 replicates when the number of offspring per sire was 100. Application of the model to lactation milk yield in dairy cattle showed that genetic variance for micro- and macro-environmental sensitivities existed. Conclusion: The algorithm and model selection criterion presented here can contribute to better understand genetic control of macro- and micro-environmental sensitivities. Designs or datasets should have at least 100 sires each with 100 offspring.

Genetic heterogeneity of residual variance : estimation of variance components using double hierarchical generalized linear models

Background: The sensitivity to microenvironmental changes varies among animals and may be under genetic control. It is essential to take this element into account when aiming at breeding robust farm animals. Here, linear mixed models with genetic effects in the residual variance part of the model can be used. Such models have previously been fitted using EM and MCMC algorithms. Results: We propose the use of double hierarchical generalized linear models (DHGLM), where the squared residuals are assumed to be gamma distributed and the residual variance is fitted using a generalized linear model. The algorithm iterates between two sets of mixed model equations, one on the level of observations and one on the level of variances. The method was validated using simulations and also by re-analyzing a data set on pig litter size that was previously analyzed using a Bayesian approach. The pig litter size data contained 10,060 records from 4,149 sows. The DHGLM was implemented using the ASReml software and the algorithm converged within three minutes on a Linux server. The estimates were similar to those previously obtained using Bayesian methodology, especially the variance components in the residual variance part of the model. Conclusions: We have shown that variance components in the residual variance part of a linear mixed model can be estimated using a DHGLM approach. The method enables analyses of animal models with large numbers of observations. An important future development of the DHGLM methodology is to include the genetic correlation between the random effects in the mean and residual variance parts of the model as a parameter of the DHGLM.

Genetic and environmental heterogeneity of residual variance of weight traits in Nellore beef cattle

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

Estimation models of variance components for farrowing interval in swine

Este trabalho teve como objetivo principal avaliar a importância da inclusão dos efeitos genético materno, comum de leitegada e de ambiente permanente no modelo de estimação de componentes de variância para a característica intervalo de parto em fêmeas suínas. Foram utilizados dados que consistiam de 1.013 observações de fêmeas Dalland (C-40), registradas em dois rebanhos. As estimativas dos componentes de variância foram realizadas pelo método da máxima verossimilhança restrita livre de derivadas. Foram testados oito modelos, que continham os efeitos fixos (grupos de contemporâneo e covariáveis) e os efeitos genético aditivo direto e residual, mas variavam quanto à inclusão dos efeitos aleatórios genético materno, ambiental comum de leitegada e ambiental permanente. O teste da razão de verossimilhança (LR) indicou a não necessidade da inclusão desses efeitos no modelo. No entanto observou-se que o efeito ambiental permanente causou mudança nas estimativas de herdabilidade, que variaram de 0,00 a 0,03. Conclui-se que os valores de herdabilidade obtidos indicam que esta característica não apresentaria ganho genético como resposta à seleção. O efeito ambiental comum de leitegada e o genético materno não apresentaram influência sobre esta característica. Já o ambiental permanente, mesmo sem ter sido significativo o seu efeito pelo LR, deve ser considerado nos modelos genéticos para essa característica, pois sua presença causou mudança nas estimativas da variância genética aditiva.

Development of a new high-performance liquid chromatographic method for the determination of ceftazidime

A rapid, accurate, and sensitive high-performance liquid chromatographic (HPLC) method was developed and validated for the determination of ceftazidime in pharmaceuticals. The method validation parameters yielded good results and included range, linearity, precision, accuracy, specificity, and recovery. The excipients in the commercial powder for injection did not interfere with the assay. Reversed-phase chromatography was used for the HPLC separation on a Waters C18 (WAT 054275; Milford, MA) column with methanol-water (70 + 30, v/v) as the mobile phase pumped isocratically at a flow rate of 1.0 mL/min. The effluent was monitored at 245 nm. The calibration graph for ceftazidime was linear from 50.0 to 300.0 mu g/mL. The values for interday and intraday precision (relative standard deviation) were < 1 %. The results obtained by the HPLC method were calculated statistically by analysis of variance. We concluded that the HPLC method is satisfactory for the determination of ceftazidime in the raw material and pharmaceuticals.

Common factors method to predict the carcass composition tissue in kid goats

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

Induced variational method from supersymmetric quantum mechanics and the screened Coulomb potential

The formalism of supersymmetric quantum mechanics supplies a trial wave function to be used in the variational method. The screened Coulomb potential is analyzed within this approach. Numerical and exact results for energy eigenvalues are compared.

Variance heterogeneity for milk yield in Brazilian and Colombian Holstein herds

Heterogeneity of variances for milk yield (MY) was determined for Brazilian and Colombian herds. The herds were grouped as high or low variability within each country, using as criterion the phenotypic standard deviation (PSD) of MY in the contemporary groups of cows, from the first to the sixth calving. Brazilian and Colombian herds with PSD greater than 1,168 kg or 1,012 kg, respectively, were classified as high variability while the herd groups with values lower than those were classified as low variability. The genetic parameters for MY within each herd group were estimated using bivariate analysis in an animal model and the restricted maximum likelihood method with a derivative free algorithm, using 72,280 first lactations of cows, daughters of 1,880 sires. Heterogeneous variances were found, and Brazilian herds with high PSD had the greatest additive and residual genetic variances and heritability coefficients for MY. MY genetic correlation coefficients between herds of high and low variability within each country were 0.96 and 0.93 while between countries they varied from 0.72 to 0.81, suggesting that there was a reclassification of animals in the two countries and also heterogeneity of variances. This phenomenon leads to the questioning of the strategy of imported semen usage and the need to do genetic evaluations to identify sires with greater genetic potential for (sub) tropical environmental conditions.

Estimation of variance components and prediction of breeding values in rubber tree breeding using the REML/BLUP procedure

The present paper deals with estimation of variance components, prediction of breeding values and selection in a population of rubber tree [Hevea brasiliensis (Willd. ex Adr. de Juss.) Müell.-Arg.] from Rio Branco, State of Acre, Brazil. The REML/BLUP (restricted maximum likelihood/best linear unbiased prediction) procedure was applied. For this purpose, 37 rubber tree families were obtained and assessed in a randomized complete block design, with three unbalanced replications. The field trial was carried out at the Experimental Station of UNESP, located in Selvíria, State of Mato Grosso do Sul, Brazil. The quantitative traits evaluated were: girth (G), bark thickness (BT), number of latex vessel rings (NR), and plant height (PH). Given the unbalanced condition of the progeny test, the REML/BLUP procedure was used for estimation. The narrow-sense individual heritability estimates were 0.43 for G, 0.18 for BT, 0.01 for NR, and 0.51 for PH. Two selection strategies were adopted: one short-term (ST - selection intensity of 8.85%) and the other long-term (LT - selection intensity of 26.56%). For G, the estimated genetic gains in relation to the population average were 26.80% and 17.94%, respectively, according to the ST and LT strategies. The effective population sizes were 22.35 and 46.03, respectively. The LT and ST strategies maintained 45.80% and 28.24%, respectively, of the original genetic diversity represented in the progeny test. So, it can be inferred that this population has potential for both breeding and ex situ genetic conservation as a supplier of genetic material for advanced rubber tree breeding programs. Copyright by the Brazilian Society of Genetics.

Influence of cavity preparation and curing method on the marginal seal of resin composite restorations: an in vitro evaluation.

OBJECTIVE: To evaluate the influence of cavity design and photocuring method on the marginal seal of resin composite restorations. METHOD AND MATERIALS: Seventy-two bovine teeth were divided into 2 groups: group 1 received box-type cavity preparations, and group 2 received plate-type preparations. Each group was divided into 3 subgroups. After etching and bonding, Z250 resin composite (3M Espe) was applied in 2 equal increments and cured with 1 of 3 techniques: (1) conventional curing for 30 seconds at 650 mW/cm2; (2) 2-step photocuring, in which the first step was performed 14 mm from the restoration for 10 seconds at 180 mW/cm2 and the second step was performed in direct contact for 20 seconds at 650 mW/cm2; or (3) progressive curing using Jetlite 4000 (J. Morita) for 8 seconds at 125 mW/cm2 and then 22 seconds at 125 mW/cm2 up to 500 mW/cm2. The specimens were thermocycled for 500 cycles and then submitted to dye penetration with a 50% silver nitrate solution. Microleakage was assessed using a stereomicroscope. Data were analyzed using analysis of variance and Tukey test (5% level of significance). RESULTS: A statistically significant difference was found between groups when a double interaction between photocuring and cavity preparation was considered (P = .029). CONCLUSIONS: No one type of cavity preparation or photocuring method prevented micro-leakage. The plate-type preparation showed the worst dye penetration when conventional and progressive photocuring methods were used. The best results were found using the 2-step photocuring with the plate-type preparation.

Genetic (co)variance components for ratio of lamb weight to ewe metabolic weight as an indicator of ewe efficiency

Data from a multibreed commercial flock located at Mid-West of Brazil, supported by Programa de Melhoramento Genético de Caprinos e Ovinos de Corte (GENECOC), were used to estimate genetic parameters of traits related to ewe productivity by Average Information Restricted Maximum Likelihood method applied to an animal model. The analyzed traits were litter weight at birth (LWB) and at weaning (LWW), ewe weight at weaning (EW) and ewe production efficiency, estimated by WEE=LWW/EW 0.75. The heritabilities were 0.26±0.05, 0.32±0.06, 0.37±0.03 and 0.10±0.02 for LWB, LWW, EW and WEE, respectively. Significant effects for direct heterosis were observed for LWW and EW. Recombination losses were important for EW and WEE. Genetic correlations of LWB with LWW, EW and WEE were 0.68, 0.37 and 0.15, respectively; of LWW with EW and WEE were 0.30 and 0.34, respectively; and between EW and WEE was -0.25. Even though it is a low heritability trait, WEE can be indicated as a selection criteria for improving the ewe productivity without increasing the mature weight of animals due to its genetic correlations with LWW and other traits. © 2011 Elsevier B.V.

A useful empirical bayesian method to analyse industrial data from saturated factorial designs

The use of saturated two-level designs is very popular, especially in industrial applications where the cost of experiments is too high. Standard classical approaches are not appropriate to analyze data from saturated designs, since we could only get the estimates of the main factor effects and we would not have degrees of freedom to estimate the variance of the error. In this paper, we propose the use of empirical Bayesian procedures to get inferences for data obtained from saturated designs. The proposed methodology is illustrated assuming a simulated data set. © 2013 Growing Science Ltd. All rights reserved.