Given two samples of size n_1 and n_2 from a data set where each sample consists of K functional observations (channels), each recorded on T grid points,  the function energy method implements a hypothesis test of equality of channel-wise mean at each channel using the bootstrapped distribution of maximum energy to control family wise error. The function energy_method_complex accomodates complex valued functional observations.

Provides a function for distribution free control chart based on the change point model, for multivariate statistical process control. The main constituent of the chart is the energy test that focuses on the discrepancy between empirical characteristic functions of two random vectors. This new control chart highlights in three aspects. Firstly, it is distribution free, requiring no knowledge of the random processes. Secondly, this control chart can monitor mean and variance simultaneously. Thirdly it is devised for multivariate time series which is more practical in real data application. Fourthly, it is designed for online detection (Phase II), which is central for real time surveillance of stream data. For more information please refer to O. Okhrin and Y.F. Xu (2017) <https://github.com/YafeiXu/working_paper/raw/master/CPM102.pdf>.

Data published by the United States Federal Energy Regulatory Commission including
    electric company financial data, natural gas company financial data, 
    hydropower plant data, liquified natural gas plant data, oil company financial data
    natural gas company financial data, and natural gas storage field data.

Compute energy landscapes using a digital elevation model
 and body mass of animals.

Fully robust versions of the elastic net estimator are introduced for linear and binary and multinomial regression, in particular high dimensional data. The algorithm searches for outlier free subsets on which the classical elastic net estimators can be applied. A reweighting step is added to improve the statistical efficiency of the proposed estimators. Selecting appropriate tuning parameters for elastic net penalties are done via cross-validation.