Augmented Lagrangian method with nonmonotone penalty parameters for constrained optimization

At each outer iteration of standard Augmented Lagrangian methods one tries to solve a box-constrained optimization problem with some prescribed tolerance. In the continuous world, using exact arithmetic, this subproblem is always solvable. Therefore, the possibility of finishing the subproblem resolution without satisfying the theoretical stopping conditions is not contemplated in usual convergence theories. However, in practice, one might not be able to solve the subproblem up to the required precision. This may be due to different reasons. One of them is that the presence of an excessively large penalty parameter could impair the performance of the box-constraint optimization solver. In this paper a practical strategy for decreasing the penalty parameter in situations like the one mentioned above is proposed. More generally, the different decisions that may be taken when, in practice, one is not able to solve the Augmented Lagrangian subproblem will be discussed. As a result, an improved Augmented Lagrangian method is presented, which takes into account numerical difficulties in a satisfactory way, preserving suitable convergence theory. Numerical experiments are presented involving all the CUTEr collection test problems.

STATISTICAL TEST FOR GENOTYPE AND ENVIRONMENT CONTRIBUTION IN THE GENOTYPES x ENVIRONMENTS INTERACTION MATRIX

The objective of the present work was to propose a method for testing the contribution of each level of the factors in a genotypes x environments (GxE) interaction using multi-environment trials analyses by means of an F test. The study evaluated a data set, with twenty genotypes and thirty-four environments, in a block design with four replications. The sum of squares within rows (genotypes) and columns (environments) of the GxE matrix was simulated, generating 10000 experiments to verify the empirical distribution. Results indicate a noncentral chi-square distribution for rows and columns of the GxE interaction matrix, which was also verified by the Kolmogorov-Smirnov test and Q-Q plot. Application of the F test identified the genotypes and environments that contributed the most to the GxE interaction. In this way, geneticists can select good genotypes in their studies.

A modified Primal-Dual Logarithmic-Barrier Method for solving the Optimal Power Flow problem with discrete and continuous control variables

The aim of solving the Optimal Power Flow problem is to determine the optimal state of an electric power transmission system, that is, the voltage magnitude and phase angles and the tap ratios of the transformers that optimize the performance of a given system, while satisfying its physical and operating constraints. The Optimal Power Flow problem is modeled as a large-scale mixed-discrete nonlinear programming problem. This paper proposes a method for handling the discrete variables of the Optimal Power Flow problem. A penalty function is presented. Due to the inclusion of the penalty function into the objective function, a sequence of nonlinear programming problems with only continuous variables is obtained and the solutions of these problems converge to a solution of the mixed problem. The obtained nonlinear programming problems are solved by a Primal-Dual Logarithmic-Barrier Method. Numerical tests using the IEEE 14, 30, 118 and 300-Bus test systems indicate that the method is efficient. (C) 2012 Elsevier B.V. All rights reserved.

A method for harvesting unfermented pollen from stingless bees (Hymenoptera, Apidae, Meliponini)

Pollen traps used for harvesting pollen from Apis mellifera do not work for stingless bees, as most species have small entrances and rapidly deposit large quantities of propolis at any barrier in front of the nest. Some stingless beekeepers harvest pollen by removing it directly from pollen pots, but this pollen is normally fermented and unpalatable. The aim of this study was to test a new method for harvesting pollen from stingless bee colonies before it begins to ferment. Colonies of Scaptotrigona depilis were removed and replaced by empty hives, which were occupied by the returning foragers and used for storing pollen and nectar. After one week, the pollen and honey were harvested directly from the storing pots and weighed. On average, the colonies produced 8.7 g of honey and 54.2 g of unfermented pollen (n = 10). This method is a viable option for harvesting unfermented pollen from stingless bees, especially with species that harvest large amounts of pollen. The unfermented pollen of S. depilis was well received in taste tests, receiving higher scores than fermented pollen, and similar scores to A. mellifera pollen, so could have great commercial possibilities. It is also a good method for studying the foraging of stingless bees because the amount of harvested food can be easily and precisely quantified.

On reducing test length for FSMs with extra states

A long-standing problem when testing from a deterministic finite state machine is to guarantee full fault coverage even if the faults introduce extra states in the implementations. It is well known that such tests should include the sequences in a traversal set which contains all input sequences of length defined by the number of extra states. This paper suggests the SPY method, which helps reduce the length of tests by distributing sequences of the traversal set and reducing test branching. It is also demonstrated that an additional assumption about the implementation under test relaxes the requirement of the complete traversal set. The results of the experimental comparison of the proposed method with an existing method indicate that the resulting reduction can reach 40%. Experimental results suggest that the additional assumption about the implementation can help in further reducing the test suite length. Copyright (C) 2011 John Wiley & Sons, Ltd.

GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

MODEL PARTITIONS AND COMPACT TEST CASE SUITES

We present a generalized test case generation method, called the G method. Although inspired by the W method, the G method, in contrast, allows for test case suite generation even in the absence of characterization sets for the specification models. Instead, the G method relies on knowledge about the index of certain equivalences induced at the implementation models. We show that the W method can be derived from the G method as a particular case. Moreover, we discuss some naturally occurring infinite classes of FSM models over which the G method generates test suites that are exponentially more compact than those produced by the W method.

Short communication: Principal components and factor analytic models for test-day milk yield in Brazilian Holstein cattle

A total of 46,089 individual monthly test-day (TD) milk yields (10 test-days), from 7,331 complete first lactations of Holstein cattle were analyzed. A standard multivariate analysis (MV), reduced rank analyses fitting the first 2, 3, and 4 genetic principal components (PC2, PC3, PC4), and analyses that fitted a factor analytic structure considering 2, 3, and 4 factors (FAS2, FAS3, FAS4), were carried out. The models included the random animal genetic effect and fixed effects of the contemporary groups (herd-year-month of test-day), age of cow (linear and quadratic effects), and days in milk (linear effect). The residual covariance matrix was assumed to have full rank. Moreover, 2 random regression models were applied. Variance components were estimated by restricted maximum likelihood method. The heritability estimates ranged from 0.11 to 0.24. The genetic correlation estimates between TD obtained with the PC2 model were higher than those obtained with the MV model, especially on adjacent test-days at the end of lactation close to unity. The results indicate that for the data considered in this study, only 2 principal components are required to summarize the bulk of genetic variation among the 10 traits.

The modified Hodge test is a useful tool for ruling out klebsiella pneumoniae carbapenemase

OBJECTIVE: Enterobacteriaceae bacteria harboring Klebsiella pneumoniae carbapenemase are a serious worldwide threat. The molecular identification of these pathogens is not routine in Brazilian hospitals, and a rapid phenotypic screening test is desirable. This study aims to evaluate the modified Hodge test as a phenotypic screening test for Klebsiella pneumoniae carbapenemase. METHOD: From April 2009 to July 2011, all Enterobacteriaceae bacteria that were not susceptible to ertapenem according to Vitek2 analysis were analyzed with the modified Hodge test. All positive isolates and a random subset of negative isolates were also assayed for the presence of blaKPC. Isolates that were positive in modified Hodge tests were sub-classified as true-positives (E. coli touched the ertapenem disk) or inconclusive (distortion of the inhibition zone of E. coli, but growth did not reach the ertapenem disk). Negative results were defined as samples with no distortion of the inhibition zone around the ertapenem disk. RESULTS: Among the 1521 isolates of Enterobacteriaceae bacteria that were not susceptible to ertapenem, 30% were positive for blaKPC, and 35% were positive according to the modified Hodge test (81% specificity). Under the proposed sub-classification, true positives showed a 98% agreement with the blaKPC results. The negative predictive value of the modified Hodge test for detection was 100%. KPC producers showed high antimicrobial resistance rates, but 90% and 77% of these isolates were susceptible to aminoglycoside and tigecycline, respectively. CONCLUSION: Standardizing the modified Hodge test interpretation may improve the specificity of KPC detection. In this study, negative test results ruled out 100% of the isolates harboring Klebsiella pneumoniae carbapenemase-2. The test may therefore be regarded as a good epidemiological tool.

Evaluation of Inductively Coupled Plasma Mass Spectrometry for Determining Ca, Cu, Fe, Mg, Mn, Se and Zn in Bovine Semen Samples using a Simple Sample Dilution Method

A simple and fast method for the determination of Ca, Cu, Fe, Mg, Mn, Se and Zn in bovine semen by quadrupole inductively coupled plasma spectrometry (q-ICP-MS) is described. Prior to analysis, samples (200 mu L) were diluted 1:50 in a solution containing 0.01% v/v Triton (R) X-100 and 0.5% v/v nitric acid and directly analyzed by ICP-MS. The limits of detection of the method are 0.3, 0.03, 0.2, 0.04, 0.04, 0.03 and 0.03 mu g L-1 for Ca-44, Cu-63, Fe-57, Mg-24, Zn-64, Se-82 and Mn-55, respectively. For purposes of comparison and method validation, four ordinary bovine semen samples were directly analyzed by ICP-MS and by flame atomic absorption spectrometry (FAAS) or graphite furnace atomic absorption spectrometry (GF AAS), with no statistical difference between the techniques at the 95% level when applying the t-test. Then, the proposed method was applied in the determinations of Ca, Cu, Fe, Mg, Mn, Se and Zn in collected samples of bovine semen from different breeds, which are used in reproduction programs and artificial insemination.

A new method for calculating the convergent point of a moving group

Context. Convergent point (CP) search methods are important tools for studying the kinematic properties of open clusters and young associations whose members share the same spatial motion. Aims. We present a new CP search strategy based on proper motion data. We test the new algorithm on synthetic data and compare it with previous versions of the CP search method. As an illustration and validation of the new method we also present an application to the Hyades open cluster and a comparison with independent results. Methods. The new algorithm rests on the idea of representing the stellar proper motions by great circles over the celestial sphere and visualizing their intersections as the CP of the moving group. The new strategy combines a maximum-likelihood analysis for simultaneously determining the CP and selecting the most likely group members and a minimization procedure that returns a refined CP position and its uncertainties. The method allows one to correct for internal motions within the group and takes into account that the stars in the group lie at different distances. Results. Based on Monte Carlo simulations, we find that the new CP search method in many cases returns a more precise solution than its previous versions. The new method is able to find and eliminate more field stars in the sample and is not biased towards distant stars. The CP solution for the Hyades open cluster is in excellent agreement with previous determinations.

Assessment of Pharmaceutical Equivalence: Difference Test or Equivalence Test?

Pharmaceutical equivalence is an important step towards the confirmation of similarity and Interchangeability among pharmaceutical products, particularly regarding those that win not be tested for bioequivalence. The aim of this paper is to compare traditional difference testing to two one-side equivalence tests in the assessment of pharmaceutical equivalence, by means of equivalence studies between similar, generic and reference products of acyclovir cream, atropine sulfate injection, meropenem for injection, and metronidazole injection. All tests were performed in accordance with the Brazilian Pharmacopeia or the United States Pharmacopeia. All four possible combinations of results arise in these comparisons of difference testing and equivalence testing. Most of the former did not show significant difference, whereas the latter presented similarity. We concluded that equivalence testing is more appropriate than difference testing, what can make it a useful tool to assess pharmaceutical equivalence in products that will not be tested for bioequivalence.

Towards a precise test for malaria diagnosis in the Brazilian Amazon: comparison among field microscopy, a rapid diagnostic test, nested PCR, and a computational expert system based on artificial neural networks

Abstract Background Accurate malaria diagnosis is mandatory for the treatment and management of severe cases. Moreover, individuals with asymptomatic malaria are not usually screened by health care facilities, which further complicates disease control efforts. The present study compared the performances of a malaria rapid diagnosis test (RDT), the thick blood smear method and nested PCR for the diagnosis of symptomatic malaria in the Brazilian Amazon. In addition, an innovative computational approach was tested for the diagnosis of asymptomatic malaria. Methods The study was divided in two parts. For the first part, passive case detection was performed in 311 individuals with malaria-related symptoms from a recently urbanized community in the Brazilian Amazon. A cross-sectional investigation compared the diagnostic performance of the RDT Optimal-IT, nested PCR and light microscopy. The second part of the study involved active case detection of asymptomatic malaria in 380 individuals from riverine communities in Rondônia, Brazil. The performances of microscopy, nested PCR and an expert computational system based on artificial neural networks (MalDANN) using epidemiological data were compared. Results Nested PCR was shown to be the gold standard for diagnosis of both symptomatic and asymptomatic malaria because it detected the major number of cases and presented the maximum specificity. Surprisingly, the RDT was superior to microscopy in the diagnosis of cases with low parasitaemia. Nevertheless, RDT could not discriminate the Plasmodium species in 12 cases of mixed infections (Plasmodium vivax + Plasmodium falciparum). Moreover, the microscopy presented low performance in the detection of asymptomatic cases (61.25% of correct diagnoses). The MalDANN system using epidemiological data was worse that the light microscopy (56% of correct diagnoses). However, when information regarding plasma levels of interleukin-10 and interferon-gamma were inputted, the MalDANN performance sensibly increased (80% correct diagnoses). Conclusions An RDT for malaria diagnosis may find a promising use in the Brazilian Amazon integrating a rational diagnostic approach. Despite the low performance of the MalDANN test using solely epidemiological data, an approach based on neural networks may be feasible in cases where simpler methods for discriminating individuals below and above threshold cytokine levels are available.

Cost effectiveness of OptiMal® rapid diagnostic test for malaria in remote areas of the Amazon Region, Brazil

Abstract Background In areas with limited structure in place for microscopy diagnosis, rapid diagnostic tests (RDT) have been demonstrated to be effective. Method The cost-effectiveness of the Optimal® and thick smear microscopy was estimated and compared. Data were collected on remote areas of 12 municipalities in the Brazilian Amazon. Data sources included the National Malaria Control Programme of the Ministry of Health, the National Healthcare System reimbursement table, hospitalization records, primary data collected from the municipalities, and scientific literature. The perspective was that of the Brazilian public health system, the analytical horizon was from the start of fever until the diagnostic results provided to patient and the temporal reference was that of year 2006. The results were expressed in costs per adequately diagnosed cases in 2006 U.S. dollars. Sensitivity analysis was performed considering key model parameters. Results In the case base scenario, considering 92% and 95% sensitivity for thick smear microscopy to Plasmodium falciparum and Plasmodium vivax, respectively, and 100% specificity for both species, thick smear microscopy is more costly and more effective, with an incremental cost estimated at US$549.9 per adequately diagnosed case. In sensitivity analysis, when sensitivity and specificity of microscopy for P. vivax were 0.90 and 0.98, respectively, and when its sensitivity for P. falciparum was 0.83, the RDT was more cost-effective than microscopy. Conclusion Microscopy is more cost-effective than OptiMal® in these remote areas if high accuracy of microscopy is maintained in the field. Decision regarding use of rapid tests for diagnosis of malaria in these areas depends on current microscopy accuracy in the field.

Comparison between the six-minute walk test and the six-minute step test in post stroke patients

Abstract Background The Stroke remains one of the major chronic diseases worldwide, and is considered a major cause of disability, which results not only in persistent neurological deficits, but also in the high physical deconditioning, nevertheless there are not many forms of assessing functional capacity in this population. We aimed to investigate the feasibility of the Six Minute Walk Teste and the Six-Minute Step Test (6MST) in post-stroke patients and compare the behavior of physiological variables during the 6MST and the Six-Minute Walk Test (6MWT), by correlating the functional performance obtained in both tests. Method The 6MWT was carried out according to the American Thoracic Society (ATS) and the 6MST was performed in six minutes in order to compare it to the 6MWT in a 20 cm step. Was included post-stroke individuals able to walk without aid. All of them did the 6MWT and the 6MST. Results 12 patients participated in the study. There was no statistical difference in the parameters analyzed when tests were compared. There was poor correlation between the functional performance in both tests. Conclusion The 6MWT and the 6MST is feasible for post-stroke patients and physiological responses are equal during the performance of both tests. However, there was no correlation with respect to functional performance, which was assessed by the distance walked in the 6MWT and by the number of steps climbed in the 6MST.