Efficient clustering in collaborative filtering recommender system: Hybrid method based on genetic algorithm and gravitational emulation local search algorithm

The ultimate goal of the Recommender System (RS) is to offer a proposal that is very close to the user's real opinion. Data clustering can be effective in increasing the accuracy of production proposals by the RS. In this paper, single-objective hybrid evolutionary approach is proposed for clustering items in the offline collaborative filtering RS. This method, after generating a population of randomized solutions, at each iteration, improves the population of solutions first by Genetic Algorithm (GA) and then by using the Gravitational Emulation Local Search (GELS) algorithm. Simulation results on standard datasets indicate that although the proposed hybrid meta-heuristic algorithm requires a relatively high run time, it can lead to more appropriate clustering of existing data and thus improvement of the Mean Absolute Error (MAE), Root Mean Square Error (RMSE) and Coverage criteria.     

Determination of the structural ensemble of the molten globule state of a protein by computer simulations

We have used computer simulations to investigate the structural nature of the molten globule (MG) state of canine milk lysozyme. To sample the conformational space efficiently, we performed replica-exchange umbrella sampling simulations with the radius of gyration as a reaction coordinate. We applied the Weighted Histogram Analysis Method to the trajectory of the simulations to obtain the potential of mean force, from which we identified representative structures corresponding to local minima in the free energy surface. The representative structures obtained in this way are in accord with the characteristics of the MG state reported previously by experimental studies. We conjecture that the MG state comprises a series of partially structured states undergoing relatively fast conformational interchange.     

Modelling time‐varying low‐temperature‐induced mortality rates for pupae of Tuta absoluta (Gelechiidae,Lepidoptera)

In a series of laboratory experiments, acclimated pupae of Tuta absoluta were exposed to various constant low temperatures in order to estimate their maximum survival times (Kaplan–Meier, Lt_99.99). A Weibull function was fitted to the data points, describing maximum survival time as a function of temperature. In another experiment at ?6°C, the progress of mortality increasing with exposure time was identified. These values were fitted by a sigmoidal function converging asymptotically to 100% mortality for very long exposure times. Analysing mortality data from the maximum survival experiment by a generalized linear model showed a significant common slope parameter (p < .001) that reveals parallelism of the survival curves at each temperature if a log time axis is used. These curves appear stretched (time scaled) if plotted with a nonlogarithmic time axis. By combining these mathematical relations, it was possible to calculate a species‐specific ‘mortality surface’ which exhibits mortalities, depending on temperature and duration of exposure. In order to accumulate hourly mortalities for courses of varying temperatures, an algorithm was developed which yields mortality values from that surface taking into account the attained mortality level. In validation experiments, recorded mortalities were compared against modelled mortalities. Prediction of mortality was partially supported by the model, but pupae experiencing intensely fluctuating temperatures showed decreased mortality, probably caused by rapid cold hardening during exposure. Despite this observation, mortality data converged to distinct levels very close to 100% depending on the intensity of temperature fluctuations that were characteristic for different types of experiments. The highest mortality limit occurred at intensely fluctuating temperatures in laboratory experiments. This constituted a benchmark that was not reached under various field conditions. Thus, it was possible to identify temperature limits for the extinction of field populations of Tuta absoluta pupae.     

Multiple imputation for multivariate data with missing and below-threshold measurements: time-series concentrations of pollutants in the Arctic

Many chemical and environmental data sets are complicated by the existence of fully missing values or censored values known to lie below detection thresholds. For example, week-long samples of airborne particulate matter were obtained at Alert, NWT, Canada, between 1980 and 1991, where some of the concentrations of 24 particulate constituents were coarsened in the sense of being either fully missing or below detection limits. To facilitate scientific analysis, it is appealing to create complete data by filling in missing values so that standard complete-data methods can be applied. We briefly review commonly used strategies for handling missing values and focus on the multiple-imputation approach, which generally leads to valid inferences when faced with missing data. Three statistical models are developed for multiply imputing the missing values of airborne particulate matter. We expect that these models are useful for creating multiple imputations in a variety of incomplete multivariate time series data sets.     

Maximum likelihood estimation of two-level latent variable models with mixed continuous and polytomous data

Two-level data with hierarchical structure and mixed continuous and polytomous data are very common in biomedical research. In this article, we propose a maximum likelihood approach for analyzing a latent variable model with these data. The maximum likelihood estimates are obtained by a Monte Carlo EM algorithm that involves the Gibbs sampler for approximating the E-step and the M-step and the bridge sampling for monitoring the convergence. The approach is illustrated by a two-level data set concerning the development and preliminary findings from an AIDS preventative intervention for Filipina commercial sex workers where the relationship between some latent quantities is investigated.     

A study of deleterious gene structure in plants using Markov chain Monte Carlo

The characteristics of deleterious genes have been of great interest in both theory and practice in genetics. Because of the complex genetic mechanism of these deleterious genes, most current studies try to estimate the overall magnitude of mortality effects on a population, which is characterized classically by the number of lethal equivalents. This number is a combination of several parameters, each of which has a distinct biological effect on genetic mortality. In conservation and breeding programs, it is important to be able to distinguish among different combinations of these parameters that lead to the same number of lethal equivalents, such as a large number of mildly deleterious genes or a few lethal genes, The ability to distinguish such parameter combinations requires more than one generation of mating. We propose a model for survival data from a two-generation mating experiment on the plant species Brassica rapa, and we enable inference with Markov chain Monte Carlo. This computational strategy is effective because a vast amount of missing genotype information must be accounted for. In addition to the lethal equivalents, the two-generation data provide separate information on the average intensity of mortality and the average number of deleterious genes carried by an individual. In our Markov chain Monte Carlo algorithm, we use a vector proposal distribution to overcome inefficiency of a single-site Gibbs sampler. Information about environmental effects is obtained from an outcrossing experiment conducted in parallel with the two-generation mating experiments.     

Statistical inference for serial dilution assay data

Serial dilution assays are widely employed for estimating substance concentrations and minimum inhibitory concentrations. The Poisson-Bernoulli model for such assays is appropriate for count data but not for continuous measurements that are encountered in applications involving substance concentrations. This paper presents practical inference methods based on a log-normal model and illustrates these methods using a case application involving bacterial toxins.     

Approximate standard errors in semiparametric models

We consider semiparametric models with p regressor terms and q smooth terms. We obtain an explicit expression for the estimate of the regression coefficients given by the back-fitting algorithm. The calculation of the standard errors of these estimates based on this expression is a considerable computational exercise. We present an alternative, approximate method of calculation that is less demanding. With smoothing splines, the method is exact, while with loess, it gives good estimates of standard errors. We assess the adequacy of our approximation and of another approximation with the help of two examples.     

Beyond the rotamer library: genetic algorithm combined with the disturbing mutation process for upbuilding protein side-chains

The disturbing genetic algorithm, incorporating the disturbing mutation process into the genetic algorithm flow, has been developed to extend the searching space of side-chain conformations and to improve the quality of the rotamer library. Moreover, the growing generation amount idea, simulating the real situation of the natural evolution, is introduced to improve the searching speed. In the calculations using the pseudo energy scoring function of the root mean squared deviation, the disturbing genetic algorithm method has been shown to be highly efficient. With the real energy function based on AMBER force field, the program has been applied to rebuilding side-chain conformations of 25 high-quality crystallographic structures of single-protein and protein-protein complexes. The averaged root mean standard deviation of atom coordinates in side-chains and veracities of the torsion angles of chi(1) and chi(1) + chi(2) are 1.165 A, 88.2 and 72.9% for the buried residues, respectively, and 1.493 A, 79.2 and 64.7% for all residues, showing that the method has equal precision to the program SCWRL, whereas it performs better in the prediction of buried residues and protein-protein interfaces. This method has been successfully used in redesigning the interface of the Basnase-Barstar complex, indicating that it will have extensive application in protein design, protein sequence and structure relationship studies, and research on protein-protein interaction.