Linear regression without correspondences is the problem of performing a linear regression fit to a dataset for which the correspondences between the independent samples and the observations are ...unknown. Such a problem naturally arises in diverse domains such as computer vision, data mining, communications and biology. In its simplest form, it is tantamount to solving a linear system of equations, for which the entries of the right hand side vector have been permuted. This type of data corruption renders the linear regression task considerably harder, even in the absence of other corruptions, such as noise, outliers or missing entries. Existing methods are either applicable only to noiseless data or they are very sensitive to initialization or they work only for partially shuffled data. In this paper we address these issues via an algebraic geometric approach, which uses symmetric polynomials to extract permutation-invariant constraints that the parameters <inline-formula> <tex-math notation="LaTeX">\xi ^{*} \in \mathbb {R} ^{\text {n}} </tex-math></inline-formula> of the linear regression model must satisfy. This naturally leads to a polynomial system of n equations in n unknowns, which contains <inline-formula> <tex-math notation="LaTeX">\xi ^{*} </tex-math></inline-formula> in its root locus. Using the machinery of algebraic geometry we prove that as long as the independent samples are generic, this polynomial system is always consistent with at most n ! complex roots, regardless of any type of corruption inflicted on the observations. The algorithmic implication of this fact is that one can always solve this polynomial system and use its most suitable root as initialization to the Expectation Maximization algorithm. To the best of our knowledge, the resulting method is the first working solution for small values of n able to handle thousands of fully shuffled noisy observations in milliseconds.
•Common existing regression approaches are statistically unsound.•The regression can be conducted using residuals which uphold the mass balance.•Regression approach impacts mainly the apparent ...energetic heterogeneity, nf.•Equilibration time and applied particle size fraction impact both coefficients.•Fine particle size fractions give a reasonable estimate of the equilibrium.
An accurate understanding of equilibrium conditions is an essential starting point for dimensioning and modelling efforts in adsorption systems. However, key assumptions and methodologies in the generation and analysis of adsorption data sets might be contributing to the error in the description of the equilibrium. This study investigates the potential impacts of the equilibration time, test adsorbent particle size selection, and the regression technique applied, on the fitted coefficients of the Freundlich equilibrium model in the adsorption of arsenate on Granular Ferric Hydroxide. The choice of regression algorithm was found to impact primarily the exponent describing the energetic heterogeneity of the adsorption surface and resulting non-linearity of the isotherm, nf. Non-linear regression with a hybrid error function on the observed solute concentration rather than loading was found to be a suitable approach. Insufficient equilibration time was found to impact predominantly the strength of adsorption, given by Kf, when applying a small adsorbent size fraction, and both coefficients when applying a coarser technical size fraction. Application of the small size fractions in lieu of the larger in accelerated trials appeared to over-estimate both coefficients, but is a reasonable alternative to long contact times.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
As machine learning becomes widely used for automated decisions, attackers have strong incentives to manipulate the results and models generated by machine learning algorithms. In this paper, we ...perform the first systematic study of poisoning attacks and their countermeasures for linear regression models. In poisoning attacks, attackers deliberately influence the training data to manipulate the results of a predictive model. We propose a theoretically-grounded optimization framework specifically designed for linear regression and demonstrate its effectiveness on a range of datasets and models. We also introduce a fast statistical attack that requires limited knowledge of the training process. Finally, we design a new principled defense method that is highly resilient against all poisoning attacks. We provide formal guarantees about its convergence and an upper bound on the effect of poisoning attacks when the defense is deployed. We evaluate extensively our attacks and defenses on three realistic datasets from health care, loan assessment, and real estate domains.
In the period 1991-2015, algorithmic advances in Mixed Integer Optimization (MIO) coupled with hardware improvements have resulted in an astonishing 450 billion factor speedup in solving MIO ...problems. We present a MIO approach for solving the classical best subset selection problem of choosing k out of p features in linear regression given n observations. We develop a discrete extension of modern first-order continuous optimization methods to find high quality feasible solutions that we use as warm starts to a MIO solver that finds provably optimal solutions. The resulting algorithm (a) provides a solution with a guarantee on its suboptimality even if we terminate the algorithm early, (b) can accommodate side constraints on the coefficients of the linear regression and (c) extends to finding best subset solutions for the least absolute deviation loss function. Using a wide variety of synthetic and real datasets, we demonstrate that our approach solves problems with n in the 1000s and p in the 100s in minutes to provable optimality, and finds near optimal solutions for n in the 100s and p in the 1000s in minutes. We also establish via numerical experiments that the MIO approach performs better than Lasso and other popularly used sparse learning procedures, in terms of achieving sparse solutions with good predictive power.
Full text
Available for:
BFBNIB, INZLJ, NMLJ, NUK, PNG, SAZU, UL, UM, UPUK, ZRSKP
The linearization of adsorption equations is controversial. The estimation of fitting parameters strongly depends on the linearization method, magnitude of experimental error, and data range. ...Although many studies contrast linear versions of these equations with their non-linear counterparts, linearization is preferred due to its simplicity since a line could be represented with fewer experimental points than a curve. An in-depth analysis was carried out to compare the accuracy of linear and non-linear models. Although different transformations linearize Langmuir isotherms, only one form yields reliable fitting parameters. Linear transformations could also lead to a statistical bias, favoring a model that does not represent the experimental behavior. Similar observations are discussed regarding the pseudo-second-order kinetic model. Linearization of Freundlich isotherms, pseudo-first-order kinetic models, and fixed-bed adsorption models through logarithms implies that attention must be taken on the logarithm limits by properly selecting the data range. Linearization also promotes the incorrect interpretation of models due to oversimplification. The linearized van't Hoff equation would yield a reasonable fit with fewer experimental points than the non-linear regression, which requires more data to assure convergence. In this sense, there is convincing evidence that non-linear regression is a more robust and reliable tool for adsorption modeling.
Full text
Available for:
BFBNIB, DOBA, GIS, IJS, IZUM, KILJ, KISLJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
The COVID-19 pandemic originated from the city of Wuhan of China has highly affected the health, socio-economic and financial matters of the different countries of the world. India is one of the ...countries which is affected by the disease and thousands of people on daily basis are getting infected. In this paper, an analysis of daily statistics of people affected by the disease are taken into account to predict the next days trend in the active cases in Odisha as well as India.
A valid global data set is collected from the WHO daily statistics and correlation among the total confirmed, active, deceased, positive cases are stated in this paper. Regression model such as Linear and Multiple Linear Regression techniques are applied to the data set to visualize the trend of the affected cases.
Here a comparison of Linear Regression and Multiple Linear Regression model is performed where the score of the model R2tends to be 0.99 and 1.0 which indicates a strong prediction model to forecast the next coming days active cases. Using the Multiple Linear Regression model as on July month, the forecast value of 52,290 active cases are predicted towards the next month of 15th August in India and 9,358 active cases in Odisha if situation continues like this way.
These models acquired remarkable accuracy in COVID-19 recognition. A strong correlation factor determines the relationship among the dependent (active) with the independent variables (positive, deceased, recovered).
•Multiple linear regression model is proposed for prediction of Active cases in COVID-19 daily data.•The model predicts a value of 52,290 active cases in India and 9358 active cases in Odisha towards the 15th of August.•The ANOVA results shows a significant P value that accepts the proposed model.•Statistical results show MLR model has fair predictive potential over the LR model.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ
The depopulation of rural areas is a global problem that is particularly serious in Europe and in large areas of Spain, such as the interior of the Valencia Region. Despite the consensus in the ...diagnosis of the problem, there is still a debate about the policies and measures to solve it. For this reason, in order to help policymakers, this work carries out an analysis of the determinants that have contributed to the population growth of small rural municipalities in the Valencia Region. The analysis includes variables that have not been used in previous research, such as cultural heritage, cycling and walking routes, organic farming or renewable energy. Among the variables that contribute to the fight against depopulation, in addition to basic public services (education and health), there are variables from the primary, secondary and tertiary sectors. There are therefore activities and resources to try to diversify the economic activity to revitalise small rural municipalities.
•More than a half of Valencian small rural municipalities are losing population.•New variables tested as depopulation determinants.•Schools are the main driver of population growth in rural areas.•Primary, secondary and tertiary sector activities help to sustain rural populations.•Sectoral diversification of the economy could be the key to rural development.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Researchers examine contrasts between analysis of variance (ANOVA) effects but seldom contrasts between regression coefficients even though such coefficients are an ANOVA generalization. Regression ...weight contrasts can be analyzed by reparameterizing the linear model. Two pairwise contrast models are developed for the study of qualitative differences among predictors. One leads to tests of null hypotheses that the regression weight for a reference predictor equals each of the other weights. The second involves ordered predictors and null hypotheses that the weight for a predictor equals that for the variables just above or below in the ordering. As illustration, qualitative differences in high school math course content are related to math achievement. The models facilitate the study of qualitative differences among predictors and the allocation of resources. They also readily generalize to moderated, hierarchical, and generalized linear forms.
Full text
Available for:
NUK, OILJ, SAZU, UKNU, UL, UM, UPUK
Recent results established that EM enjoys global convergence for Gaussian Mixture Models. For Mixed Linear Regression, however, only local convergence results have been established, and those only ...for the high signal-to-noise ratio (SNR) regime. In this work, we completely characterize the global optimality of EM: we show that starting from any randomly initialized point, the EM algorithm converges to the true parameter β * at the minimax statistical rates under all SNR regimes. Toward this goal, we first show the global convergence of the EM algorithm at the population level. Then we provide a complete characterization of statistical and computational behaviors of EM under all SNR regimes with finite samples. In particular: (i) When the SNR is sufficiently large, the EM updates converge to the true parameter β * at the standard parametric convergence rate O (( d / n ) 1/2 ) after O (log( n / d )) iterations. (ii) In the regime where the SNR is above O (( d / n ) 1/4 ) and below some constant, the EM iterates converge to a O (SNR −1 ( d / n ) 1/2 ) neighborhood of the true parameter, when the number of iterations is of the order O (SNR −2 log( n / d )). (iii) In the low SNR regime where the SNR is below O (( d / n ) 1/4 ), we show that EM converges to a O (( d / n ) 1/4 ) neighborhood of the true parameters, after O (( n / d ) 1/2 ) iterations. By providing tight convergence guarantees of the EM algorithm in middle-to-low SNR regimes, we reveal that in low SNR, EM changes rate, matching the n −1/4 rate of the MLE, a behavior that previous work had been unable to show.
When the outcome is binary, psychologists often use nonlinear modeling strategies such as logit or probit. These strategies are often neither optimal nor justified when the objective is to estimate ...causal effects of experimental treatments. Researchers need to take extra steps to convert logit and probit coefficients into interpretable quantities, and when they do, these quantities often remain difficult to understand. Odds ratios, for instance, are described as obscure in many textbooks (e.g., Gelman & Hill, 2006, p. 83). I draw on econometric theory and established statistical findings to demonstrate that linear regression is generally the best strategy to estimate causal effects of treatments on binary outcomes. Linear regression coefficients are directly interpretable in terms of probabilities and, when interaction terms or fixed effects are included, linear regression is safer. I review the Neyman-Rubin causal model, which I use to prove analytically that linear regression yields unbiased estimates of treatment effects on binary outcomes. Then, I run simulations and analyze existing data on 24,191 students from 56 middle schools (Paluck, Shepherd, & Aronow, 2013) to illustrate the effectiveness of linear regression. Based on these grounds, I recommend that psychologists use linear regression to estimate treatment effects on binary outcomes.
Full text
Available for:
CEKLJ, FFLJ, NUK, ODKLJ, PEFLJ, UPUK