The ELMVIS method 
is an interesting Machine Learning method that optimize a cost function by changing assignment between two sets of samples, or by changing the order of samples in one set which is the same. The cost function is learned by an Extreme Learning Machine (ELM)[13, 12, 10]
, a fast method for training feed-forward neural networks with convenient mathematical properties[11, 14]. Such optimization problem is found in various applications like open-loop Traveling Salesman problem  or clustering  (mapping between samples and clusters), but not in Neural Networks. ELMVIS is unique in a sense that it combines the optimal assignment task with neural network optimization problem; the latter is optimized at each step of ELMVIS.
A recent advance in ELMVIS+ method  set its runtime speed comparable or faster than other state-of-the-art methods in visualization application. However there are unresolved problems like a greedy optimization leading to a local optimum. Also ELMVIS+ can be applied to a much wider range of problems than a simple visualization or a visualization accounting for the class information , which have not been tested or reported yet. This paper addresses the aforementioned drawbacks, and presents the most recent research advances in the family of ELMVIS methods.
The proposed incremental ELMVIS allows for iterative growth of dataset size and model complexity. Incremental ELMVIS learns an approximate global data structure with a few data samples and a simple ELM model, because at a very small scale global and local optimums are similar or the same. Then more data samples are added to the model, choosing the ones that better fit an existing ELM. After adding a batch of new samples, the current model is refined by ELMVIS+ method. This refinement keeps the global optimum due to the greedy optimization and only small changes. More neurons are added to ELM as the dataset size grows, to better separate the data.
Iterative ELMVIS is useful for semi-supervised learning, starting from the data with known outputs and adding more data with unknown outputs, simultaneously updating the model. It can even apply to completely unsupervised datasets, where it finds an input-output dependency, learns it with ELM model, and then simultaneously expands the supervised part of a dataset and updates an ELM model that encodes the input-output dependency.
The experiments have shown the ability of iterative ELMVIS to improve global optimum, successfully perform semi-supervised and unsupervised learning with complex tasks. Current version of the method is found limited to good separation between only two classes in data (which it learns first), ignoring samples of the additional classes until the first two ones are exhausted, and poorly fitting the additional classes into the learned two-class model. Solution to this problem will be considered in further works on the topic.
An iterative extension of ELMVIS+ methodology is rather straight-forward, as explained below. ELMVIS methods start with a list of (visualization) samples and an unordered set of data samples; and it finds an optimal order of data samples in the set by a greedy search with changing positions of many samples at once (ELMVIS) or only two samples at once (ELMVIS+).
Iterative ELMVIS splits data samples into fixed, candidate and available ones. Fixed samples have their order fixed and cannot be moved by an iterative ELMVIS. Candidate samples are a small number of samples which are chosen from candidate+available ones to maximize the cost function. This cost function takes into account fixed and current candidate samples, but ignores the available samples. Once current candidate samples are chosen optimally, they are added to the fixed ones, and the method is repeated with a few more candidate samples - until the available data samples are exhausted.
2.1 Extreme Learning Machine
Extreme Learning Machine is a way of training feedforward neural networks  with a single hidden layer that features randomly assigned input weights , explicit non-iterative solution for output weights and an extreme computation speed and scalability . This model is used as a non-linear cost function in all ELMVIS methods. This short summary introduces the notations to the reader.
The goal of ELM model is to approximate the projection function using a representative training dataset. As ELM is a deterministic model, the function is assumed to be deterministic, and a noise is added to cover the deviation of true outputs from the predictions by a deterministic function
The Extreme Learning Machine  (ELM) is a neural network with input, hidden and output neurons. The hidden layer weights and biases are initialized randomly and are fixed. The hidden layer neurons apply a transformation function
to their outputs that is usually a non-linear transformation function with bounded output like sigmoid or hyperbolic tangent.
The output of the hidden layer is denoted by with an expression
where the function is applied element-wise, and can also be gathered in a matrix for convenience.
2.2 ELMVIS+ method
ELMVIS+ method  approximates a relation between visualization (or input) space and data space by an ELM model. The task has representative samples and , but their order is unknown. The method assumes fixed order of samples joined in matrix , and finds a suitable order of samples joined in matrix by exchanging pairs of rows in . Contrary to a common use of ELM, data samples are the outputs of ELM and visualization coordinates are the inputs (thus ELM predicts original data ). Visualization coordinates
are chosen arbitrary and fixed – they can be distributed randomly with normal or uniform distribution, or initialized on a regular grid.
The optimization criterion is a cosine similarity betweenand predicted by ELM. A low error means that its possible to reconstruct data from the given visualization points, thus the visualization points keep information about the data. The reconstruction is approximated by the ELM model in ELMVIS.
There is an explicit formula for a change of error (negative cosine similarity) for swapping two rows in and re-training ELM with this new dataset. The readers can refer to the original paper  for the full formula. It is based on the expression for the change on error in case a row in is changed by amount:
2.3 Incremental ELMVIS
Incremental ELMVIS splits all data samples in three groups: fixed, candidate and available samples. The separation is maintained with two indexes: is the number of fixed samples, and is the number of fixed+candidate ones.
Incremental ELMVIS works similar to ELMVIS+. First, the initial numbers of fixed and candidate samples are given by and . Then swap indexes and are selected randomly to replace one candidate sample with an available one. The change of error is computed by the formula (3) for the change in the candidate row of the data matrix . Compared to ELMVIS+, the change in the available sample is ignored. Then for the negative , an update step is performed for sample as in equation (4) and for both and as in equation (5).
If there is no improvement during a large number of swaps, the current candidate samples are added to the fixed ones (), and more samples are added as candidates (). Candidate samples are already initialized with the data samples at indexes . The method then repeats for another iteration. Iterations stop when no available samples are left.
In the original ELMVIS+, matrix took the most space and limited the maximum amount of processed samples (its memory size is ). In incremental ELMVIS, only a part of the whole matrix is needed. That relaxed memory requirements of the method, and allows to use a very large pool of available samples. The memory constraints of incremental ELMVIS apply only to the number of optimized data samples.
3 Experimental Results
Incremental ELMVIS method is developed for two main applications. The first one is achieving a better global optimum in ELMVIS+. The original EMLVIS+ is a good visualization method, however it has an unwanted feature: with a large number of neurons in ELM it fragments clusters in the visualized data. This happens with large amount of data and a complex ELM model. An iterative ELMVIS that starts with small amount of data and a simple model keeps all similar data together; then more data samples are gradually added while the total picture changes little due to local minimum in ELMVIS+ optimization.
The second application is finding unknown relations in datasets. This is an unsupervised learning field relevant to the current Big Data trends, when a large amount of interesting data is available - but there pre-processing like manual labeling or classification. It is possible to extract relations inside data automatically by iteratively growing an ELMVIS+ model between two sets of data samples (they don’t have to be related to visualization). Results for both applications are presented below.
3.1 Better Optimum with ELMVIS
ELMVIS+ method is fast and works with large datasets, but it has a greedy optimization approach that leads to local optimality of the solution. Such local optimum is close to a global one for small datasets and simple ELM models, but with a large dataset and many neurons in ELM model the visualization data is split into multiple small clusters with local optimality, non-representative of a global picture.
A better global optimality is achievable with an incremental ELMVIS. This experiment uses MNIST digits  with their original features (grayscale pixels), with 500 digits for each of the classes 0-4. It starts by seeding several cluster as shown on Figure 1 by bold samples, and a simple ELM model. Then gradually added data fits into the existing model (Figure 1, left). An ELM learns an easy separation between two clusters, and an incremental ELMVIS prefers to add samples of these clusters until they are available (Figure 1, right).
The incremental ELMVIS learns a good model that separated between two different kinds of data. Then there is no data samples of these two types left, it begins adding more types, starting at the boundary (Figure 2, left). These additional classes are mapped to a single area, although they go over the two previously learned cluster as there is no space left on the visualization (Figure 2, right). The ELMVIS still ignores the last available class (digits 4) because it is not represented on the visualization space.
3.2 Better Optimum with Semi-supervised ELMVIS
In the previous experiment, there were no sharp borders between clusters because ELMVIS used all the visualization space to show only two clusters, and then has to map additional data clusters other them. Sharper borders can be obtained by running the original ELMVIS+ on the fixed set of data points after each iteration of the incremental ELMVIS. That will make space for more clusters by compacting the existing ones; while still preserving the global structure as ELMVIS+ optimization goes to the local optimum only.
In addition, a semi-supervised approach is tested where the clusters are initialized with a larger number of samples. This experiment uses 5 classes of digits with 1000 samples per class, initialized with 20 samples per class as shown on Figure 3 (left). A large initialization set forces ELM to learn all the classes, and add samples from all of them instead of only two (Figure 3, right).
When the method need space to map more samples from a particular class, an additional ELMVIS+ step moves existing clusters to give that space while keeping sharp boundaries between the classes. The effect is shown on Figure 4 where all clusters are moved to give more space for digits 2. The global structure is well preserved, with large clusters keeping their place.
The added ELMVIS+ step refines the candidate samples placement within the fixed ones, that is necessary towards the end of visualization when there may be no samples of the desired class left, or no spaces left within the desired class area. An improvement of ELMVIS+ step fitting the candidate samples is presented on Figure 5.
3.3 Data Structure Discovery with Unsupervised ELMVIS
The inputs to ELMVIS are not limited to visualization coordinates; they can be arbitrary data. Thus ELMVIS is a feasible method for finding structure in the data in an unsupervised manner. This experiments takes a number of MNIST digits in random (undefined) order, and maps them to the same number of classes (in zero/one encoding), or to the same number of different MNIST digits of the same classes.
First dataset has zero-one classes as inputs, and MNIST digits as outpus. However it cannot be used to train a supervised model because it is unknown which input corresponds to which output – a common situation in the analysis of large automatically acquired data corpora that has not been manually labeled. The goal of an incremental ELMVIS is to reconstruct the correct input-output pairing. The experiment uses 100 samples per class with 2 or 5 classes of digits. Note that the best pairing across all permutations of classes is reported, as in the unsupervised setup with equal amount of samples there is no way to tell ELMVIS which class should go to which digit.
The resulting confusion matrices are shown on Figure 6. An incremental ELMVIS successfully paired classes with pictures of MNIST digits. The method mapped classes arbitrary (i.e. class is mapped to digit 4), but this is to be expected from a purely unsupervised method.
Another experiment is performed in a similar setup, but instead of binary class representations the incremental ELMVIS tried to map MNIST digits to other MNIST digits (of the same classes). The mapping is successful with two classes. With more than two classes the same feature always appears: two random classes are separated well while other classes are randomly mixed with them. This outcome is in line with the results observed in section 3.1 where two classes are clearly separated at the beginning, followed by other classes mapped over them.
An iterative extension to the original ELMVIS+ method is proposed in this paper. It iteratively selects a small number of best fitting samples from all the available ones, and adds them to the model. It allows for a much larger set of potential samples than ELMVIS+ by limiting memory requirements to already fitted samples rather than to all available ones, keeping the high speed of the ELMVIS+ at the same time.
The method improves global structure of ELMVIS+ visualization by starting with a small dataset, and gradually adding more data or increasing the complexity of the model. It preserves the global structure, sharp boundaries between classes, and has a possible semi-supervised extension where the samples are mapped to the specified places on the visualization space.
The proposed method is capable of unsupervised data structure detection. It excels in reconstructing a randomly shuffled dataset with unknown pairing between inputs and outputs. It is also capable of finding a mapping between two complex data spaces as shown on MNIST digits example.
The methodology needs further investigation and improvement to counter the observed drawbacks, specifically the tendency of learning an easy model first leading to problems in incorporating more complex parts to the global picture.
-  Anton Akusok, Kaj-Mikael Björk, Yoan Miche, and Amaury Lendasse. High-Performance Extreme Learning Machines: A Complete Toolbox for Big Data Applications. IEEE Access, 3:1011–1025, July 2015.
-  Anton Akusok, Yoan Miche, Kaj-Mikael Björk, Rui Nian, Paula Lauren, and Amaury Lendasse. ELMVIS+: Improved Nonlinear Visualization Technique Using Cosine Distance and Extreme Learning Machines. In Jiuwen Cao, Kezhi Mao, Jonathan Wu, and Amaury Lendasse, editors, Proceedings of ELM-2015 Volume 2: Theory, Algorithms and Applications (II), pages 357–369. Springer International Publishing, Cham, 2016.
-  Anton Akusok, Yoan Miche, Kaj-Mikael Björk, Rui Nian, Paula Lauren, and Amaury Lendasse. ELMVIS+: Fast Nonlinear Visualization Technique based on Cosine Distance and Extreme Learning Machines. Neurocomputing, forthcoming 2016.
-  Ethem Alpaydin. Introduction to machine learning. MIT press, 2014.
-  Eric Cambria and others. Extreme Learning Machines [Trends & Controversies]. IEEE Intelligent Systems, 28(6):30–59, Nov.-Dec. 2013.
-  Andrey Gritsenko, Anton Akusok, Yoan Miche, Kaj-Mikael Björk, Baek, Stephen, and Amaury Lendasse. Combined Nonlinear Visualization and Classification: ELMVIS++C. Vancouver, Canada, July 2016. Accepted for publication.
-  Gregory Gutin and Abraham P. Punnen, editors. The traveling salesman problem and its variations. Combinatorial optimization. Kluwer Academic, Dordrecht, London, 2002.
-  Simon Haykin. Neural Networks: A Comprehensive Foundation (2nd Edition). Prentice Hall, 2nd edition, July 1998.
-  Guang-Bin Huang. An Insight into Extreme Learning Machines: Random Neurons, Random Features and Kernels. Cognitive Computation, 6(3):376–390, 2014.
-  Guang-Bin Huang. What are Extreme Learning Machines? Filling the Gap Between Frank Rosenblatt’s Dream and John von Neumann’s Puzzle. Cognitive Computation, 7(3):263–278, 2015.
-  Guang-Bin Huang, Lei Chen, and Chee-Kheong Siew. Universal approximation using incremental constructive feedforward networks with random hidden nodes. IEEE Transactions on Neural Networks, 17(4):879–892, July 2006.
-  Guang-Bin Huang, Hongming Zhou, Xiaojian Ding, and Rui Zhang. Extreme learning machine for regression and multiclass classification. Systems, Man, and Cybernetics, Part B: Cybernetics, IEEE Transactions on, 42(2):513–529, April 2012.
-  Guang-Bin Huang, Qin-Yu Zhu, and Chee-Kheong Siew. Extreme learning machine: Theory and applications. Neural Networks Selected Papers from the 7th Brazilian Symposium on Neural Networks (SBRN ’04)7th Brazilian Symposium on Neural Networks, 70(1–3):489–501, December 2006.
-  Yoan Miche, Antti Sorjamaa, Patrick Bas, Olli Simula, Christian Jutten, and Amaury Lendasse. OP-ELM: Optimally-Pruned Extreme Learning Machine. IEEE Transactions on Neural Networks, 21(1):158–162, January 2010.
C. Radhakrishna Rao and Sujit Kumar Mitra.
Generalized inverse of a matrix and its applications.
Proceedings of the Sixth Berkeley Symposium on Mathematical Statistics and Probability, Volume 1: Theory of Statistics, pages 601–620, Berkeley, CA, 1972. University of California Press.
-  Y. Lecun, L. Bottou, Y. Bengio, and P. Haffner. Gradient-based learning applied to document recognition. Proceedings of the IEEE, 86(11):2278–2324, November 1998.