We revisit the subject of Be star candidates towards the Magellanic Clouds, previously studied by the authors using SPM4 proper motions. We obtain GAIA DR2 parallaxes and proper motions for 2357 and 994 LMC and SMC Be candidates, respectively. Parallaxes and proper motions vs. color V-I easily reveal the presence of the redder galactic contaminant foreground, as concluded in our previous work, but this time we do find a few red Be stars candidates consistent with being true Magellanic objects. A membership assessment to each Magellanic Cloud is done for each Be candidate, based on the distribution of their parallaxes and proper motions. From a compilation of published catalogues of spectroscopically confirmed Be stars, we found that 40 (LMC) and 64 (SMC) of these Be candidates, are in fact Be stars. Near infrared IRSF JHKs magnitudes were obtained for about 70% the Be stars candidates with GAIA DR2 astrometric data. Mid-infrared SAGE IRAC 3.6, 4.5, 5.8 and 8.0 $\mu$m magnitudes were obtained for about 85% as well. 6 LMC and 7 SMC confirmed Be stars show optical, near- or mid-infrared colours redder than what has been typically measured for Classical Be stars. Several of the Be candidates follow those redder-than-expected colours distributions suggesting the existence of more red Magellanic Be stars.

Statistical pattern recognition methods have provided competitive solutions for variable star classification at a relatively low computational cost. In order to perform supervised classification, a set of features is proposed and used to train an automatic classification system. Quantities related to the magnitude density of the light curves and their Fourier coefficients have been chosen as features in previous studies. However, some of these features are not robust to the presence of outliers and the calculation of Fourier coefficients is computationally expensive for large data sets. We propose and evaluate the performance of a new robust set of features using supervised classifiers in order to look for new Be star candidates in the OGLE-IV Gaia south ecliptic pole field. We calculated the proposed set of features on six types of variable stars and on a set of Be star candidates reported in the literature. We evaluated the performance of these features using classification trees and random forests along with K-nearest neighbours, support vector machines, and gradient boosted trees methods. We tuned the classifiers with a 10-fold cross-validation and grid search. We validated the performance of the best classifier on a set of OGLE-IV light curves and applied this to find new Be star candidates. The random forest classifier outperformed the others. By using the random forest classifier and colour criteria we found 50 Be star candidates in the direction of the Gaia south ecliptic pole field, four of which have infrared colours consistent with Herbig Ae/Be stars. Supervised methods are very useful in order to obtain preliminary samples of variable stars extracted from large databases. As usual, the stars classified as Be stars candidates must be checked for the colours and spectroscopic characteristics expected for them.

We present a proper motion investigation of a sample of Be star candidates towards the Magellanic Clouds, which has resulted in the identification of separate populations, in the Galactic foreground and in the Magellanic background. Be stars are broadly speaking B-type stars that have shown emission lines in their spectra. In this work, we studied a sample of 2446 and 1019 Be star candidates towards the LMC and SMC respectively, taken from the literature and proposed as possible Be stars due to their variability behaviour in the OGLE-II I band. JHKs magnitudes from the IRSF catalog and proper motions from the SPM4 catalog, were obtained for 1188 and 619 LMC and SMC Be stars candidates, respectively. Color-color and vector-point diagrams were used to identify different populations among the Be star candidates. In the LMC sample, two populations with distinctive infrared colours and kinematics were found, the bluer sample is consistent with being in the LMC and the redder one with belonging to the Milky Way disk. This settles the nature of the redder sample which had been described in previous publications as a possible unknown subclass of stars among the Be candidates in the LMC. In the SMC sample, a similar but less evident result was obtained, since this apparent unknown subclass was not seen in this galaxy. We confirm that in the selection of Be stars by their variability, although generally successful, there is a higher risk of contamination by Milky Way objects towards redder B$-$V and V$-$I colors.

The investigation of the non-linearity of the Leavitt law is a topic that began more than seven decades ago, when some of the studies in this field found that the Leavitt law has a break at about ten days. The goal of this work is to investigate a possible statistical cause of this non-linearity. By applying linear regressions to OGLE-II and OGLE-IV data, we find that, in order to obtain the Leavitt law by using linear regression, robust techniques to deal with influential points and/or outliers are needed instead of the ordinary least-squares regression traditionally used. In particular, by using $M$- and $MM$-regressions we establish firmly and without doubts the linearity of the Leavitt law in the Large Magellanic Cloud, without rejecting or excluding Cepheid data from the analysis. This implies that light curves of Cepheids suggesting blending, bumps, eclipses or period changes, do not affect the Leavitt law for this galaxy. For the SMC, including this kind of Cepheids, it is not possible to find an adequate model, probably due to the geometry of the galaxy. In that case, a possible influence of these stars could exist.

The universality and linearity of the Leavitt law are hypotheses commonly adopted in studies of galaxy distances using Cepheid variables as standard candles. In order to test these hypotheses, we obtain slopes of the Leavitt law using linear regressions of fundamental-mode Cepheids observed by the Optical Gravitational Lensing Experiment project in the Magellanic Clouds. We find that the slopes in $VI$-bands and in the Wesenheit index behave exponentially, indicating non-linearity. We also find that the slopes obtained using long-period Cepheids can be considered as universal in the $VI$-bands, but not in the Wesenheit index.