Next Article in Journal
Vibration Prediction of the Robotic Arm Based on Elastic Joint Dynamics Modeling
Next Article in Special Issue
Comparison of the Usability of Apple M1 Processors for Various Machine Learning Tasks
Previous Article in Journal
Cycling through 360° Virtual Reality Tourism for Senior Citizens: Empirical Analysis of an Assistive Technology
Previous Article in Special Issue
Sensor Screening Methodology for Virtually Sensing Transmission Input Loads of a Wind Turbine Using Machine Learning Techniques and Drivetrain Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast Weighted Fuzzy C-Medoids Clustering for Time Series Data Based on P-Splines

1
College of Computer and Information Engineering, Henan Normal University, Xinxiang 453007, China
2
Engineering Lab of Intelligence Business & Internet of Things, Xinxiang 453007, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(16), 6163; https://doi.org/10.3390/s22166163
Submission received: 26 June 2022 / Revised: 10 August 2022 / Accepted: 15 August 2022 / Published: 17 August 2022
(This article belongs to the Special Issue Intelligent Sensors and Machine Learning)

Abstract

:
The rapid growth of digital information has produced massive amounts of time series data on rich features and most time series data are noisy and contain some outlier samples, which leads to a decline in the clustering effect. To efficiently discover the hidden statistical information about the data, a fast weighted fuzzy C-medoids clustering algorithm based on P-splines (PS-WFCMdd) is proposed for time series datasets in this study. Specifically, the P-spline method is used to fit the functional data related to the original time series data, and the obtained smooth-fitting data is used as the input of the clustering algorithm to enhance the ability to process the data set during the clustering process. Then, we define a new weighted method to further avoid the influence of outlier sample points in the weighted fuzzy C-medoids clustering process, to improve the robustness of our algorithm. We propose using the third version of mueen’s algorithm for similarity search (MASS 3) to measure the similarity between time series quickly and accurately, to further improve the clustering efficiency. Our new algorithm is compared with several other time series clustering algorithms, and the performance of the algorithm is evaluated experimentally on different types of time series examples. The experimental results show that our new method can speed up data processing and the comprehensive performance of each clustering evaluation index are relatively good.

1. Introduction

With the rapid development of computer information technology, the research on time series clustering is gradually rising in various fields such as finance, biology, medicine, meteorology, electricity, industry, and agriculture [1,2,3]. To better summarize and discover the effective information on these time series data, time series clustering technology has received extensive attention from researchers [4]. Through clustering, we can mine the key hidden information in time series data and find some regularities between time series sample points. It provides indispensable preparation for forecasting and anomaly detection of time series data. Consequently, the study of time series clustering is of great significance.
The data is often imperfect. Due to some human errors, machine failures, or unavoidable natural influences, noise points and outlier samples may appear in the data, this phenomenon also exists in time series data.To better avoid the impact of abnormal fluctuations in time series data onto clustering, some researchers try to smooth the data before clustering. Abraha et al. (2003) [5] proposed using B-spline fitting function data and using a K-means algorithm to divide the estimated model coefficients, which proved the strong consistency and effectiveness of the clustering method. Iorio et al. (2016) [6] used a P-spline smoother to model the raw data and then used the k-means algorithm to cluster the fitted data objects according to the best spline coefficients. D’Urso et al. (2021) [7,8] applied the B-spline fitting coefficient to the robust fuzzy clustering method and proposed four different robust fuzzy clustering methods of time series data, then applied them to the clustering of COVID-19 time series. The above studies all smooth the data to reduce the negative effects of noise. In this paper, we fit the time series in the original data using P-splines [9]. P-spline shrinks the difference penalty based on B-spline. This method will be described in detail in Section 2.1.
With the extensive research on time series clustering, many targeted methods have been proposed, such as model-based clustering, feature-based clustering, deep learning-based clustering, and fuzzy theory-based clustering. Among them, the Fuzzy C-Means (FCM) algorithm is the classical representative of fuzzy clustering [10], which adds the fuzzy concept to the k-means algorithm, and is widely used in the clustering of time series [11]. The membership value between the sample point and the cluster center is calculated to indicate the closeness of each sample point in the cluster center [12]. The FCM algorithm provides a more flexible clustering criterion so that time series samples with greater similarity can also be classified into the class to which they belong [13,14]. Compared with the K-means algorithm, FCM provides a better clustering effect [15,16]. In addition, different fuzzy systems are also applied to time series data. Pattanayak et al. (2020) [17] used FCM to determine intervals of unequal lengths, and the support vector machine also considered the membership value when modeling fuzzy logical relationships. Kumar et al. (2021) [18] proposed the nested Particle Swarm Optimization algorithm and exhaustive search Particle Swarm Optimization algorithm for fuzzy time series prediction. Xian et al. (2021) [19] proposed a fuzzy c-means clustering method based on n-Pythagorean fuzzy sets, and effectively improved the prediction accuracy by using the clustering results. Al-qaness et al. (2022) [20,21] proposed optimizing the adaptive neuro-fuzzy inference system using the enhanced marine predator algorithm, and then again proposed combining the improved Aquila optimizer with adversarial-based learning techniques to optimize the adaptive neuro-fuzzy inference system. These fuzzy systems have excellent performance on time series data, which shows that the research on time series based on fuzzy theory has far-reaching significance.
Due to the high-dimension and complexity of time series data, there will inevitably be some noise sample points in the dataset, which affects the updating of fuzzy clustering centers. To solve this problem, a new clustering algorithm is derived based on FCM, called Fuzzy C-medoids (FCMdd) [22,23]. The center of each cluster in the FCMdd algorithm is a real object in the selected data set, which greatly reduces the influence of noise data onto the new medoids, and can better deal with the problem that the FCM algorithm is sensitive to noise data [24,25,26]. Therefore, FCMdd can obtain higher quality clustering results, which makes many researchers apply FCMdd to time series clustering [27]. Izakian H et al. (2015) [28] proposed a hybrid technique using the advantages of FCM and FCMdd and used dynamic time warping distance to cluster time series. Liu et al. (2018) [29] proposed two incremental FCMdd algorithms based on dynamic time warping distance, which can easily obtain higher-quality clustering results in the face of time-shifted time series data. D’Urso P et al. (2021) [30] clustered multivariate financial time series through an appropriate cutting process, in which the Medoids surrounding partition strategy combined with dynamic time warping distances was used in the clustering process. The above studies have improved the performance of the FCMdd algorithm in different aspects, but they did not consider the impact of outlier samples on the FCMdd algorithm. In this paper, We define a new method to update the weights of each sample to improve the FCMdd algorithm, which reduces the impact of outlier samples on the clustering process and improves the robustness of our algorithm.
A robust time series clustering method not only handles the adverse effects of noise in the dataset, but also classifies sample points quickly. One of the main influencing factors of time series clustering is the complexity of the similarity metric for time series [31]. So far, many effective similarity methods have been studied for time series data [32], among which two main types are used. One is lock-step measures, its main representative method is Euclidean Distance (ED), the other is elastic measures, its main representative methods are Dynamic Time Warping (DTW) [33] and Time Warp Edit Distance (TWED) [34]. ED is simple and common, but it is susceptible to the time offset and has a small range of applications; DTW is a well-known time series similarity measurement criterion [35], and it is also the most widely used and improved algorithm so far, but its computational cost is O N 2 when the sample object length is N. Among the edit distances, there is the TWED algorithm that satisfies the triangular inequality relationship, which can deal with time series of different sampling rates, and is also an excellent time series similarity measurement method, but the time complexity is the same as DTW. Among these similarity measures, DTW and TWED are popular, but their high time complexity leads to low clustering efficiency. This paper chooses the MASS 3 algorithm [36] in the matrix configuration file to measure the similarity between time series, which has certain robustness and can improve the accuracy of time series clustering. Moreover, the matrix configuration file has some outstanding advantages, such as simplicity, no parameters, time and space complexity are relatively low [37].
In response to the problems discussed above, we propose a fast weighted fuzzy C-medoids clustering for time series data based on P-splines, named the PS-WFCMdd. PS-WFCMdd can well deal with the negative problems caused by noise and outliers in time series data, and has fast and efficient clustering performance. The main contributions of this paper include the following three aspects:
  • Abnormal fluctuations and noises in time series data usually make it challenging to define a suitable clustering model and increase the burden of related measurements. To solve these problems, we use the P-spline to fit clustering data. Smoothing the fitted data reduces abnormal fluctuations in the data compared to the original data. At the same time, the influence of noise on fuzzy clustering is avoided. Taking the fitting data as the input of the clustering algorithm, it is easier to find similar sample points in the clustering process;
  • In the traditional fuzzy clustering method, no matter how far the sample point is from the center of the data cluster, the membership degree of the sample point belonging to this class is not 0. This makes FCM vulnerable to boundary samples. In the proposed PS-WFCMdd, we design a new weighted fuzzy C-medoids, which can reduce the influence of outlier sample points on cluster center update by defining a new weighting method;
  • Among the time series similarity metrics, ED is simple and fast, but it is vulnerable to time offset. DTW algorithm has high time complexity and low efficiency in the face of a large-scale time-series data model. As a time-series measurement standard, MASS 3 has the characteristics of fast computing speed and small space overhead due to its unique segmented computing method. Therefore, we propose to apply MASS 3 to weighted fuzzy C-medoids to measure the similarity between time series data and medoids more accurately and quickly.
The rest of the paper is organized as follows: in Section 2, we introduce some basic knowledge used in our algorithm. Section 3 presents our method. The experimental results are presented in Section 4, which also includes a discussion. Finally, the conclusions are presented in Section 5.

2. Background and Related Work

2.1. Time Series: A P-Spline Representation

B-splines are formed by connecting polynomial segments in a specific way. As an excellent curve representation method, it can maintain the continuity and local shape controllability of the fitted data after fitting the time series [38].
Suppose the set containing N time series data is T = T 1 , T 2 , T i , T N , where T i = T i 1 , T i 2 , , T i j , let B h T i j ; p denotes the value at T i j of the h-th B-spline of degree p for a given equidistant grid of knots. Based on the above conditions, the vector set L T i of the j fitting coefficients corresponding to the i-th time series can be obtained:
L T i = h b h B h T i ; p + θ , i = 1 , , N ,
where θ is the error vector and b h is the B-spline coefficient estimated through simple linear least-squares.
B-splines perform well for nonparametric modeling problems but have limited control over smoothness and data fitting. In this regard, Eilers and Marx [9] use a relatively large number of nodes and different penalties for the coefficients of adjacent B-splines, showing the connection of the common spline penalty on the integral of the squared second derivative. The difference penalty is integrated into the B-spline basis, which created P-splines. If a dense set of splines is used, the curve is obtained by minimizing L T i B b 2 . An element of B is B h T i j ; p , compared with the original sequence data, the fitting curve results obtained by b are quite different. To avoid overfitting, b can be estimated in a penalized regression setting:
b ^ = argmin b L B b 2 + λ D d b 2 ,
where D d is a d-th order difference penalty matrix such that Δ d . The parameter λ continuously controls the smoothness of the final fit. Overfitting occurs when the smoothing parameter approaches zero, while for high values of λ , the result of the fitted function will be close to the polynomial model. The optimal spline coefficients follow as:
b ^ = B B + λ D D 1 B y
Based on the above, the linear least-squares fitting model based on P-spline can be obtained by using differential penalty for the coefficients of the B-spline curve, which is as follows:
L T i = h b ^ h B h T i ; p + θ , i = 1 , , N .
In a broad sense, the P-spline curve is limited by polynomial curve fitting and extended on the basis of a linear regression model. P-splines have many useful properties: P-splines show no boundary effects, just like many types of kernel smoothers; P-splines can fit polynomial data exactly; P-splines can conserve moments of the data.
We show a comparison of the two time series results with different fluctuations after P-spline fitting in Figure 1. Here, we can see that whether it is noisy time series A or time series B with slight irregularity, the fitting curve obtained by fitting the time series with P-spline is smoother than the original curve, but it reasonably retains the key characteristics of the time series. Therefore, this paper proposes to use P-spline to fit the time series first, which can eliminate the noise peak contained in the time series, but retain the contour and trend of the sample. Then, the fitting data is used as the input of clustering to improve the final clustering effect.

2.2. Weighted FCMdd

The FCMdd algorithm can objectively assign each data object to a certain cluster with a certain value of membership, thereby generating fuzzy clusters. Compared with the clustering method of the k-medoids algorithm that directly divides the sample categories, the FCMdd algorithm can reasonably perform uncertainty analysis on the categories of time series. In fact, the objective function of FCMdd is the same as that of FCM, except that V in the clustering process is selected from the real sample space, rather than the points in the C continuous spaces generated by the algorithm as the new clustering center.
Suppose the set containing N time series data is T = T 1 , T 2 , T i , , T N , V = v 1 , v 2 , , v C , v i T represents the cluster centers obtained during the clustering process. Then, the objective function of FCMdd is defined as:
J m U , V = k = 1 C i = 1 N u i k m d i k T i , v k s . t . 0 u i j 1 , k = 1 C u i j = 1
where C represents the number of clusters, respectively; u i k represents the membership degree of the sample to the k-th clustering; m is a fuzzification parameter used to control the degree of clustering fuzziness. The larger the value of m, the wider the fuzzy decision boundary, its default value is 2; U is the numbership matrix; V is the set of cluster centers. The definition of d i k T i , v k represents the similarity between the object T i and the cluster center point v k . Both the cluster assignment results in k-medoids and the update results of medoids are discrete, while the cluster assignment results in FCMdd are continuous variables. For a randomly initialized medoid group or an updated medoid group, FCMdd can update u i k for the existing C medoid members. The update method is as follows:
u i k = l = 1 C d i k T i , v k d l k T i , v l 1 m 1 1 .
Any clustering algorithm needs to select the best cluster center. For FCMdd, selecting the best instance object as the medoid is one of the cores of the algorithm. The commonest approach is to select the object with the smallest distance from all objects in the dataset as the new medoid, but this approach will make the time complexity of the algorithm very high. We consider that only q sample points that maximize the membership in each cluster should be considered as new medoids. In this way, the update formula of cluster center v k can be defined as:
v c = min x ξ j = 1 N u c j m d c j T j , T
where ξ is the set of q medoid candidates.
The weighted Fuzzy C-Medoids clustering algorithm (WFCMdd) obtains w i by weighting each data object T i and adding them to the objective function. The objective function of the minimization of WFCMdd can be defined as follows:
J m U , V = k = 1 C i = 1 N u i k m w i d i k T i , v k s . t . 0 u i j 1 , k = 1 C u i j = 1 .
Under constraint i = 1 c u i k = 1 , the value of u i k are the same as Equation (6). The new medoid v k of cluster C are defined as follows:
v k = min x ξ j = 1 N u j k m w j d j k T j , T .

2.3. MASS 3

MASS is an algorithm for creating distance distributions for queries on time series. Using it to compute data requires only a small space overhead, with a small constant factor, allowing large datasets to be processed in main memory. In the continuous development of MASS, MASS 2 and MASS 3 came into being. Among them, MASS 3 is the best, it uses fast Fourier transform (FFT) to achieve semi-convolution, which is used to calculate all sliding dot products between query and sliding window. In order to avoid the problem of numerical error in convolution calculation when the data sample length n is very large, MASS 3 processes the data samples in segments instead of all n sample points at one time. The segments must overlap by exactly m-1 samples to maintain the continuity of the distance profile. Where the size of the block is a power of 2, and the size of the last block may not be a power of 2. MASS 3 does not depend on the query length (m), nor does it have the problem of dimension disaster, and its calculation cost is O N log ( N ) .
The specific calculation steps are as follows:
  • Scan the time series with a sliding window.
  • Z-Normalize the window to get Q.
  • The standard deviations of moving windows of a fixed size can be calculated in one linear scan. We assume T is the time series. Calculate cumulative sums of T and T 2 and store. Subtract the two cumulative sums to obtain the sum over any window, and finally calculate the standard deviation for all windows. We can obtain the expression for the standard deviation as:
    σ i = C i + m 2 C i 2 m C i + m C i m 2
    where: C = T , C 2 = T 2 .
  • Use a half convolution to compute all sliding dot products between the query and the sliding window. If x and y are vectors of polynomial coefficients, the expression for this half convolution is:
    Half conf x , y = i f f t f f t x f f t y .
  • The final vectorized working formula is:
    M x ^ , y ^ = 2 m Half conf T , Q σ T .

3. Weighted Fuzzy C-Medoids Clustering Based on P-Splines

In the process of FCMdd algorithm clustering, after classifying an outlier sample into a certain cluster, in the next iterative update, because of the constraint condition A, the outlier sample will obtain a higher value of membership. In turn, the selected medoids are biased, and finally affect the clustering results. Although the traditional WFCMdd assigns different weights to each data sample point in the clustering process to obtain a better clustering effect. However, this method does not solve the negative impact of outlier samples on updating medoids. To solve this problem, we propose to judge the importance between different data and new medoids according to the distance. This increases the importance of data in more dense areas in the data space and reduces the importance of peripheral data (outlier data). To solve this problem, after the initial medoids V is obtained, the update weight W i of the sample point x i and each cluster center v k are obtained as:
W i = k = 1 C exp x k v i s t d x i , v k
where s t d is the standard deviation between the calculated input data. By updating the weights in this way, so that the data points are closer to the v k point, then W i will be larger, and W i will be smaller on the contrary.
Based on the above research contents, we propose a fast weighted fuzzy C-medoids algorithm. The specific process of PS-WFCMdd is shown in Figure 2. The input to the algorithm is the original time series dataset T = T 1 , T 2 , T i , T N . Then, use P-spline to fit each sample point in T to obtain the fitted data set L = L T 1 , L T 2 , L T i , L T N . The P-spline will smooth the sample while preserving the main characteristics of the sample. The obtained fitting data set L is used as the processing object of clustering, and initial medoids are randomly selected from it. MASS 3 is used to efficiently calculate the similarity between each fitting data object and medoids, and convert it into the corresponding fuzzy membership matrix. The C N fuzzy membership matrix is generated after each iteration, and the intensity of the square color in the matrix represents the membership of each sample point relative to the cluster center. In the future medoids update, to solve the problem that FCMdd is sensitive to outlier samples, we propose a new weighting method to update the medoids. To redefine the weight relationship between different data objects and new medoids according to the distance, increase the weight of data objects close to medoids, and reduce the weight of outlier data objects farther from medoids. Thus, the clustering efficiency and robustness of PS-WFCMdd are improved.
In simple terms, after transforming the time series into projection coefficients on the basis of the P-spline function, it is clustered with new weights based on WFCMdd, and MASS 3 is used to measure the similarity between time series data and medoids in the clustering process. Then, the objective function of PS-WFCMdd is defined as:
J m U , V = k = 1 C i = 1 N u i k m W i M i k T i , v k s . t . 0 u i j 1 , k = 1 C u i j = 1
The iterative expression to obtain the final membership matrix and medoids is:
u i k = l = 1 C M i k L T i , v k M l k L T i , v l 1 m + 1 1 ,
v k = min x ξ j = 1 N u j k m W j M j k L T j , T .
According to Equations (4), (13), (15) and (16), the steps of PS-WFCMdd can be summarized as Algorithm 1.
Algorithm 1 Weighted fuzzy C-medoids clustering algorithm based on P-splines(PS-WFCMdd).
Input: 
Time series data set T, determine the values of parameters C, λ and max.iter;
Output: 
Clustered sample class label Y.
  1:
Fit each sample in the cluster dataset T using Equation (4), the new fitted data is obtained as L = L T 1 , L T 2 , L T i , L T N ;
  2:
Set i t e r = 0;
  3:
Randomly set C initial medoids: L T i 1 , L T i 2 , L T i C .
  4:
Repeat
  5:
    Store the current centroids L T i 1 , L T i 2 , L T i C old = L T i 1 , L T i 2 , L T i C ;
  6:
    Compute the fuzzy membership matrix U of T i in clusters (C) using Equation (15);
  7:
    Update the the new centroids L T i 1 , L T i 2 , L T i C using Equation (16);
  8:
    Update the weight W k using Equation (13);
  9:
     i t e r = i t e r + 1;
10:
Until L T i 1 , L T i 2 , L T i C old = L T i 1 , L T i 2 , L T i C or iter=max.iter;
11:
Calculate the classification label (Y) to which the sample belongs according to the final fuzzy membership U.

4. Experimental Results

4.1. Experiment Preparation and Evaluation Index

To demonstrate the effectiveness of our proposed method on time series datasets of different dimensions and data lengths, in the experiments, we selected 18 real-world datasets from the open time series database UCR for evaluation. Among them, 10 common time series data, whose detailed description is shown in Table 1, and the detailed description of 8 large-scale time series data in Table 2. The length of the time series in all datasets is between 96 and 2844, the size of the dataset is between 40 and 8236, and the number of classes in the dataset is between 2 and 8.
To make the experiment more complete, the data objects we choose are diversified, and a brief introduction is given below. Car recorded the contours of four different types of cars extracted from the traffic video. Lightning2 is a transient electromagnetic event recorded by a suite of optical and radio frequency instruments. Wafer is a database of wafers recorded during the processing of silicon wafers for semiconductor manufacturing. StarLightCurves are recorded celestial brightness data over time. Car, Lightning2, Wafer, and StarLightCurves are all time-series databases constructed from collections of process control measurements recorded from various sensors. Beef, Meat, Strawberry, and Rock are the classification records of items made by various spectrometers. Beefly, Facefour, Fish, and Herring record image databases for contour problems. Both the ECG200 and ECG5000 track the electrical activity recorded during a heartbeat. ChlorineConcentration and Mallat are simulated datasets constructed by software applications. Gunpoint recorded the movement trajectory data of two actors with their hands. Hourstwenty is the dataset recorded using smart home devices.
A statement is required here, PS-WFCMdd, PS-k-means, and PS-WFCMdd-sD set the corresponding parameters of P-splines for each experimental dataset are consistent. Each series has been modeled by P-splines taking cubic bases and third-order penalties. For the P-spline modeling, a fixed smoothing parameter λ is chosen for each dataset, and λ ranges from 0.3 to 0.6. The distance metric criterion for both FCM and PS-k-means is ED. Note that the number of centers for each dataset is also fixed, here we specify the number of classes in Table 1 and Table 2 as the C value.
In this experiment, our proposed PS-WFCMdd algorithm will be compared with six other time series clustering algorithms: traditional FCM algorithm; K-means based on P-splines [6] (PS-k-means); based on soft-DTW [39] (it is an improvement on the DTW algorithm as a differentiable loss function which can calculate all alignment costs Soft minimum) K-means algorithm (K-means-sD algorithm) [40]; soft-DTW-based K-medoids algorithm (K-medoids-sD); K-shape algorithm based on cross-correlation distance metric [41]; to better demonstrate the superiority of MASS 3, replace MASS 3 in PS-WFCMdd with softDTW, denoted as PS-WFCMdd-sD.
In the validation step, we use three widely used criteria to measure the performance of the clustering results. Due to the insufficient penalty of the Rand coefficient, the scores of the clustering algorithms are generally relatively high, and the performance of the algorithm cannot be significantly compared. Therefore, we use the adjusted Rand coefficient (ARI) [42], which is a measure of the Rand coefficient. In an improved version, the purpose is to remove the influence of random labels on the Rand coefficient evaluation results. The Fowlkes–Mallows index (FMI) [43] can be regarded as the geometric mean between precision and recall. This index shows similarity even if the number of clusters of the noisy dataset is different from that of the source data, making FMI one of the key metrics for evaluating the effectiveness of clustering. The amount of mutual information can be understood as the amount by which the uncertainty of the other decreases when guiding one. Normalized Mutual Information (NMI) is one of the most popular information-theoretic metrics in community detection methods. NMI is a normalization of the mutual information score used to scale the results between 0 (no mutual information) and 1 (perfect correlation). These standards are briefly described below.
The ARI score range is [−1, 1], negative scores indicate poor results, and a score equal to 1 indicates that the clustering results are completely consistent with the true labels. Where TP represents the number of true positive samples (the clustering results can group similar samples into the same cluster), TN represents the number of true negative samples (clustering results can group dissimilar samples into different clusters), FP represents the number of false positive samples (that is, the clustering results classify dissimilar samples into the same cluster), FN represents the number of false negative samples (clustering results group similar samples into different clusters). The definition of the ARI is:
ARI = 2 · ( T P · T N F N · F P ) ( T P + F N ) ( F N + T N ) + ( T P + F P ) ( F P + T N ) .
The parameter in the expression of FMI has the same meaning as the parameter in Equation (17) and can be defined as:
FMI = T P ( ( T P + F P ) · ( T P + F N ) ) .
The minimum score of the Fowlkes–Mallows index is 0, indicating that the clustering results are different from the actual class labels; the highest score is 1, indicating that all data have been correctly classified.
NMI value range is between [0, 1], and the optimal value is 1. Given class labels η , and a result cluster labels Y obtained by a clustering algorithm. The NMI is defined as:
NMI η , Y = 2 MI η ; Y H η + H Y
where H · means entropy; MI is given by
MI η , Y = i = 1 η j = 1 Y η i Y j N log N η i Y j η i Y j
where η i is the number of the samples in cluster η i .

4.2. Experimental Comparison and Analysis

4.2.1. Results of PS-WFCMdd Method on Ten Common Time Series Datasets

To confirm that the PS-WFCMdd method is effective on common low-dimensional datasets, in the first part of this subsection, we compare the performance of PS-WFCMdd and six other clustering methods on the 10 low-dimensional datasets in Table 1. Table 3, Table 4 and Table 5, they show the ARI index scores, FMI index scores, and NMI index scores of the seven clustering methods on each dataset, where the bold content represents the best value in the comparison results.
Table 3 shows the clustering results under ARI. We can see that on the Beef, GunPoint, and Meat datasets, the PS-WFCMdd algorithm has an obvious score advantage. On the BeetleFly, Car, ECG200, FaceFour, and Lightning2 datasets, the PS-WFCMdd algorithm also has the best ARI scores, but the performance is not significant. On the Fish and Herring datasets, the ARI score of the K-means-sD algorithm is the best, followed by the SRFCM algorithm, PS-WFCMdd performs the next best. There is no significant difference in the performance of PS-WFCMdd-sD and PS-WFCMdd under ARI on the four datasets FaceFour, Fish, Lightning2, and Meat datasets.
Table 4 shows the clustering results under FMI. We can see that on the Beef and Meat datasets, the PS-WFCMdd algorithm has an obvious score advantage. On the BeetleFly, Fish, and GunPoint datasets, the FMI score of the PS-WFCMdd algorithm is also the best. On the Car, FaceFour, and Herring datasets, the FMI score of the PS-WFCMdd-sD algorithm is the best, and the SRFCM algorithm has a certain advantage over the remaining algorithms. On the ECG200 dataset, the FMI score of the K-means-sD algorithm is the best. On the Lightning2 dataset, the K-shape algorithm has the best score.
Table 5 shows the clustering results under NMI. We can see that on the Beef, BeetleFly, GunPoint, Car, ECG200, FaceFour, and Meat datasets, the PS-WFCMdd algorithm has obvious score advantages. On the Fish dataset, the NMI score of the K-means-sD algorithm is also the best value. On the Car, FaceFour, and Herring datasets, the NMI score of the K-means-sD algorithm is the best, and the PS-WFCMdd algorithm has a certain advantage over the remaining algorithms.
By analyzing the results of the above three tables, it can be concluded that our algorithm can obtain relatively excellent performance under the three evaluation indicators when faced with ordinary low-dimensional time series data sets. The improved algorithm based on the K-means algorithm also achieved good results in some indicators, indicating that the performance of the algorithm should not be underestimated. Although the PS-WFCMdd-sD algorithm achieves good performance under the FMI evaluation index, the overall performance is still inferior to the PS-WFCMdd algorithm. This shows that the MASS 3 we used has good performance on low-dimensional data in terms of clustering effect and the time complexity of the PS-WFCMdd-sD algorithm is high, we will further analyze it in Section 4.4.

4.2.2. Results of PS-WFCMdd Method on Nine Large-Scale Time Series Datasets

To confirm that the PS-WFCMdd method is also effective on large-scale datasets in Table 2, in this section, we compare the clustering effects of PS-WFCMdd and six other clustering methods on the eight datasets in Figure 3, Figure 4 and Figure 5 show the ARI scores, FMI scores, NMI scores of the seven clustering methods on each dataset.
By looking at Figure 3. We can see that on the ChlorineConcentration, Rock, and Strawberry datasets, the PS-WFCMdd algorithm has obvious score advantages. On the ECG5000, HouseTwenty, and Wafer datasets, the ARI score of the PS-WFCMdd algorithm is also the best. On the StarLightCurves dataset, the ARI score of the PS-WFCMdd-sD algorithm is the best. On the Mallat dataset, the ARI score of the K-shape algorithm is the best. The ARI scores of the PS-WFCMdd-sD and PS-WFCMdd algorithm on the ChlorineConcentration and HouseTwenty data have little difference.
According to the results in Figure 4. On the Rock, StarLightCurves, and Strawberry datasets, the PS-WFCMdd algorithm has obvious advantages. On the ECG5000 and Mallat datasets, the PS-WFCMdd algorithm is barely optimal. The FCM algorithm achieves excellent performance on the ChlorineConcentration dataset, and the PS-K-means algorithm achieves the highest score on the HouseTwenty dataset. On the Wafer dataset, the FMI score of the PS-WFCMdd-sD algorithm is outstanding.
It can be seen from Figure 5 that on the ChlorineConcentration, Rock, StarLightCurves, Strawberry, and Wafer datasets, the PS-WFCMdd algorithm has obvious score advantages; on the HouseTwenty dataset, the PS-WFCMdd algorithm achieves only a small advantage; on the ECG5000 dataset, the NMI score of the FCM algorithm is the optimal value; on the Mallat dataset, the K-shape algorithm has the best FMI score.
From the above analysis, it can be concluded that our algorithm can achieve excellent performance under three evaluation indicators when faced with large-scale time series datasets. On some datasets, the FCM algorithm and the K-shape algorithm have also achieved good results, indicating that these two algorithms have certain advantages in processing large-scale time series data. The PS-WFCMdd-sD algorithm has achieved good performance under the three evaluation indicators, but the overall performance is still better than the PS-WFCMdd algorithm.

4.3. Experimental Analysis from the Perspective of Data Characteristics

To more intuitively explain the robustness of the PS-WFCMdd algorithm from the data features, the trend of all experimental data and the cluster center results obtained by the PS-WFCMdd algorithm are shown in Figure 6. The K-means-sD algorithm, which is generally used in time series clustering, is selected for the corresponding visual comparison, as shown in Figure 7. In both figures, the colored lines represent cluster centers, and the light-colored lines represent the original sample set. In addition, the comprehensive performance of K-means-sD in the experiments in this paper is also relatively good. Let us first observe the clustering centers in the two visualization charts. Through comparison, we can see that since k-means-sd is the clustering center gradually generated in the clustering process, the trend of the clustering centers of some data is quite different from the original data, and the difference is more obvious in BeetleFly, HouseTwenty, and Wafer. However, PS-WFCMdd is the cluster center found in the corresponding original data, and the trend of the obtained cluster center is the same as that of the original data.
We then evaluate the data from the trend in the time series dynamic analysis. From the display of the source data set alone, it can be seen that BeetleFly, Car, Fish, Herring, and ChlorineConcentration all exhibit periodic changes, among which BeetleFly has a larger time series amplitude. In terms of the experimental results on Fish and Herring, PS-WFCMdd performs generally in the three evaluation indicators, while K-means-sD has better experimental results. In the experimental results of Car and ChlorineConcentration, the advantage of PS-WFCMdd is small. PS-WFCMdd performs well in the experimental results on BeetleFly. This shows that when faced with data with obvious periodicity and stable data trend, because there are few noise and abnormal samples in these data, PS-WFCMdd has no advantage in experimental results. The remaining 13 datasets contain a large number of irregular dynamic time series, which inevitably contain a lot of noise and edge data. From the experimental results of these data, the performance of PS-WFCMdd is relatively good. However, on Housetwenty and Wafer data with mutation frequencies, this type of data often presents more difficult modeling problems. PS-WFCMdd performs best among all comparison algorithms, but is not outstanding and does not make a practical solution for mutation frequency time series.
Based on the above experimental analysis, P-splines are used to fit time series data and use the fitted data as the input of clustering, which can effectively improve the final clustering effect. Under the influence of the new weight, the PS-WFCMdd algorithm can still obtain a higher clustering effectiveness index score under the condition of noise and outlier samples in the sample data. However, our algorithm also has some limitations. For example, compared with most algorithms, PS-WFCMdd needs to adjust more smoothing parameters, so the algorithm is not very concise. When faced with time series with obvious periodicity and stable data trends and irregularly changing data with mutation, our algorithm performs in general and needs further optimization and improvement.

4.4. Time Complexity Analysis

The time complexity of each comparison algorithm and PS-WFCMdd algorithm in this experiment is shown in Table 6, where N represents the number of sample data, C represents the number of clusters, and K represents the number of iterations. It can be seen that the time complexity of PS-WFCMdd is the same as that of K-shape, and it is second only to FCM and PS-K-means algorithm. The time complexity of the remaining clustering algorithms is higher. We can see that because of the difference in time complexity between soft-DTW and MASS 3, the running time consumption of PS-WFCMdd-sD and PS-WFCMdd is one order of magnitude different. Especially in the face of large time series, the advantages of the PS-WFCMdd algorithm are more obvious.

5. Conclusions

This paper proposes a clustering method called PS-WFCMdd suitable for time series characterized by being fast, parsimonious, and noise-resistant. The idea behind our proposal is simple and effective. It can solve the effect of noisy points and outlier samples on fuzzy clustering. We recommend fitting each series with a P-spline smoother before clustering, and using the fitted data as input to the clustering analysis. We also define a new weighting method to avoid the impact of outlier sample points on medoids update. The new weight allocation makes the selected new medoids more appropriate. In the experiment, compared with soft-DTW, MASS 3 has a better advantage in our new weighted fuzzy C-medoids method, MASS 3 can measure the similarity between time series data and medoids more accurately and quickly. In addition, the performance of our weighted fuzzy C-medoids method based on P-spline has been analyzed through three clustering evaluation indexes in Section 4. Our experimental results show that our algorithm can ensure high-quality clustering performance within a reasonable computing time, whether on ordinary time series data sets or large time series data.
In future work, we will try to study a better cluster initialization method to reduce the number of iterations. The distance measurement criterion is further improved to reduce the time complexity of the algorithm, to better improve the comprehensive performance of the algorithm. There are still many open-ended problems that need to be solved. The characteristics of each time series dataset are different. To obtain higher clustering accuracy, it is necessary to analyze the characteristics of the data more carefully and study corresponding characteristics. The corresponding clustering method is obtained, and a more appropriate similarity measure criterion is selected. These considerations can be gradually realized in future research, which must be very interesting. Finally, the code URL for PS-WFCMdd is: https://github.com/houqinchen/PS-WFCMdd (accessed on 1 February 2022).

Author Contributions

All the authors devoted a lot of effort to this paper. J.X. and Q.H. put forward the concept and carried out experiments and statistics, Y.S. and X.M. actively participated in the writing of the paper, when revising the manuscript, K.Q. gave many suggestions. The specific contributions of each author are as follows: Conceptualization, J.X. and Q.H.; methodology, Q.H.; software, Q.H.; validation, Q.H. and K.Q.; formal analysis, J.X. and K.Q.; investigation, J.X. and Q.H.; resources, J.X. and Q.H.; data curation, Y.S.; writing—original draft preparation, J.X.; writing—review and editing, J.X. and Q.H.; visualization, Y.S. and X.M.; supervision, J.X. and K.Q.; project administration, J.X.; funding acquisition, J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 61976082, 62076089, and 62002103.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The dataset that support the findings of this study are available in https://www.cs.ucr.edu/~eamonn/time_series_data_2018/ (accessed on 1 February 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liao, T.W. Clustering of time series data—A survey. Pattern Recognit. 2005, 38, 1857–1874. Available online: https://www.sciencedirect.com/science/article/abs/pii/S0031320305001305 (accessed on 1 January 2022). [CrossRef]
  2. Fu, T.A. Review on time series data mining. Eng. Appl. Artif. Intell. 2011, 24, 164–181. Available online: https://www.sciencedirect.com/science/article/pii/S0952197610001727 (accessed on 1 January 2022). [CrossRef]
  3. Aghabozorgi, S.; Shirkhorshidi, A.S.; Wah, T.Y. Time-series clustering—A decade review. Inf. Syst. 2015, 53, 16–38. Available online: https://www.sciencedirect.com/science/article/abs/pii/S0306437915000733 (accessed on 1 January 2022). [CrossRef]
  4. Caiado, J.; Maharaj, E.A.; D’Urso, P. Time series clustering. In Handbook of Cluster Analysis; Chapman and Hall/CRC: Boca Raton, FL, USA, 2015; pp. 241–263. Available online: https://www.taylorfrancis.com/chapters/edit/10.1201/b19706-18/time-series-clustering-jorge-caiado-elizabeth-ann-maharaj-pierpaolo-urso (accessed on 1 June 2022).
  5. Abraham, C.; Cornillon, P.A.; Matzner-Løber, E.; Molinari, N. Unsupervised curve clustering using B-splines. Scand. J. Stat. 2003, 30, 581–595. Available online: https://onlinelibrary.wiley.com/doi/abs/10.1111/1467-9469.00350 (accessed on 1 June 2022). [CrossRef]
  6. Iorio, C.; Frasso, G.; D’Ambrosio, A.; Siciliano, R. Parsimonious time series clustering using p-splines. Expert Syst. Appl. 2016, 52, 26–38. Available online: https://www.sciencedirect.com/science/article/pii/S0957417416000063 (accessed on 1 June 2022). [CrossRef]
  7. D’Urso, P.; Giovanni, L.D.; Vitale, V. Spatial robust fuzzy clustering of COVID 19 time series based on B-splines. Spat. Stat. 2021, 38, 100518. Available online: https://www.sciencedirect.com/science/article/pii/S2211675321000282 (accessed on 1 March 2022). [CrossRef]
  8. D’Urso, P.; García-Escudero, L.A.; De Giovanni, L.; Vincenzina, V.; Agustín, M. Robust fuzzy clustering of time series based on B-splines. Int. J. Approx. Reason. 2021, 136, 223–246. Available online: https://www.sciencedirect.com/science/article/abs/pii/S0888613X21000918 (accessed on 1 March 2022). [CrossRef]
  9. Eilers, P.; Rijnmond, D.M.; Marx, B.D. Flexible smoothing with B-splines and penalties. Stat. Sci. 1996, 11, 89–121. Available online: https://projecteuclid.org/journals/statistical-science/volume-11/issue-2/Flexible-smoothing-with-B-splines-and-penalties/10.1214/ss/1038425655.full (accessed on 1 June 2022). [CrossRef]
  10. Bezdek, J.C.; Ehrlich, R.; Full, W. FCM: The fuzzy c-means clustering algorithm. Comput. Geosci. 1984, 10, 191–203. Available online: https://www.sciencedirect.com/science/article/pii/0098300484900207 (accessed on 1 June 2022). [CrossRef]
  11. Aghabozorgi, S.; Wah, T.Y. Effective clustering of time-series data using FCM. Int. J. Mach. Learn. Comput. 2014, 4, 170–176. Available online: http://www.ijmlc.org/index.php?m=content&c=index&a=show&catid=44&id=445 (accessed on 1 June 2022). [CrossRef]
  12. Liu, B.; Zhang, T.; Li, Y.; Liu, Z.; Zhang, Z. Kernel Probabilistic K-Means Clustering. Sensors 2021, 21, 1892. Available online: https://www.mdpi.com/1424-8220/21/5/1892 (accessed on 1 June 2022). [CrossRef]
  13. Bei, H.; Mao, Y.; Wang, W.; Zhang, X. Fuzzy clustering method based on improved weighted distance. Math. Probl. Eng. 2021, 2021, 6687202. Available online: https://www.hindawi.com/journals/mpe/2021/6687202/ (accessed on 1 June 2022). [CrossRef]
  14. Chen, H.; Xie, Z.; Huang, Y.; Gai, D. Intuitionistic fuzzy C-means algorithm based on membership information transfer-ring and similarity measurement. Sensors 2021, 21, 696. Available online: https://www.mdpi.com/1424-8220/21/3/696 (accessed on 1 June 2022). [CrossRef]
  15. Mahmoudi, M.R.; Baleanu, D.; Mansor, Z.; Bui, A.T.; Kim, P. Fuzzy clustering method to compare the spread rate of Covid-19 in the high risks countries. Chaos Solitons Fractals 2020, 140, 110230. Available online: https://www.sciencedirect.com/science/article/pii/S0960077920306263 (accessed on 1 June 2022). [CrossRef]
  16. Kannan, S.R.; Ramathilagam, S.; Chung, P.C. Effective fuzzy c-means clustering algorithms for data clustering problems. Expert Syst. Appl. 2012, 39, 6292–6300. Available online: https://www.sciencedirect.com/science/article/abs/pii/S0957417411016095 (accessed on 1 June 2022). [CrossRef]
  17. Pattanayak, R.M.; Panigrahi, S.; Behera, H.S. High-order fuzzy time series forecasting by using membership values along with data and support vector machine. Arab. J. Sci. Eng. 2020, 45, 10311–10325. Available online: https://linkspringer.53yu.com/article/10.1007/s13369-020-04721-1 (accessed on 1 February 2022). [CrossRef]
  18. Kumar, N.; Susan, S. Particle swarm optimization of partitions and fuzzy order for fuzzy time series forecasting of COVID-19. Appl. Soft Comput. 2021, 110, 107611. Available online: https://www.sciencedirect.com/science/article/pii/S1568494621005329 (accessed on 1 February 2022). [CrossRef]
  19. Xian, S.; Cheng, Y. Pythagorean fuzzy time series model based on Pythagorean fuzzy c-means and improved Markov weighted in the prediction of the new COVID-19 cases. Soft Comput. 2021, 25, 13881–13896. Available online: https://linkspringer.53yu.com/article/10.1007/s00500-021-06259-2 (accessed on 1 February 2022). [CrossRef]
  20. Al-qaness, M.A.A.; Ewees, A.A.; Fan, H.; AlRassas, A.M.; Elaziz, M.A. Modified aquila optimizer for forecasting oil production. Geo-Spat. Inf. Sci. 2022, 1–17. Available online: https://wwwtandfonline.53yu.com/doi/abs/10.1080/10095020.2022.2068385 (accessed on 1 February 2022). [CrossRef]
  21. Al-qaness, M.A.A.; Ewees, A.A.; Fan, H.; Abualigah, L.; Elaziz, M.A. Boosted ANFIS model using augmented marine predator algorithm with mutation operators for wind power forecasting. Appl. Energy. 2022, 314, 118851. Available online: https://www.sciencedirect.com/science/article/abs/pii/S0306261922002872 (accessed on 1 June 2022). [CrossRef]
  22. Haga, N.; Honda, K.; Ichihashi, H.; Notsu, A. Linear fuzzy clustering of relational data based on extended fuzzy c-medoids. In Proceedings of the IEEE International Conference on Fuzzy Systems (IEEE World Congress on Computational Intelligence), Hong Kong, China, 1–6 June 2008; pp. 366–371. Available online: https://ieeexplore.ieee.org/abstract/document/4630392 (accessed on 1 June 2022).
  23. Nguyen, A.C.; Ngo, T.L.; Mai, D.S. A hybrid fuzzy c-medoids clustering using the where optimization algorithm. J. Sci. Tech. Sect. Inf. Commun. Technol. 2021, 10, 1. Available online: https://jst.lqdtu.edu.vn/index.php/ict/article/view/287 (accessed on 1 June 2022).
  24. Krishnapuram, R.; Joshi, A.; Yi, L. A fuzzy relative of the k-medoids algorithm with application to web document and snippet clustering. IEEE Int. Fuzzy Syst. Conf. Proc. 1999, 3, 1281–1286. Available online: https://ieeexplore.ieee.org/abstract/document/790086 (accessed on 1 June 2022).
  25. Sardar, T.H.; Ansari, Z. Distributed big data clustering using mapreduce-based fuzzy C-medoids. J. Inst. Eng. 2022, 103, 73–82. Available online: https://link.springer.com/article/10.1007/s40031-021-00647-w (accessed on 1 June 2022). [CrossRef]
  26. Labroche, N. New incremental fuzzy c medoids clustering algorithms. In Proceedings of the Annual Meeting of the North American Fuzzy Information Processing Society, Toronto, ON, Canada, 12–14 July 2010; pp. 1–6. Available online: https://ieeexplore.ieee.org/abstract/document/5548263 (accessed on 1 June 2022).
  27. Coppi, R.; D’Urso, P.; Giordani, P. Fuzzy c-medoids clustering models for time-varying data. In Modern Information Processing; Elsevier Science: Amsterdam, The Netherlands, 2006; pp. 195–206. Available online: https://www.sciencedirect.com/science/article/pii/B9780444520753500170 (accessed on 1 June 2022).
  28. Izakian, H.; Pedrycz, W.; Jamal, I. Fuzzy clustering of time series data using dynamic time warping distance. Eng. Appl. Artif. Intell. 2015, 39, 235–244. Available online: https://www.sciencedirect.com/science/article/abs/pii/S0952197614003078 (accessed on 1 June 2022). [CrossRef]
  29. Liu, Y.; Chen, J.; Wu, S.; Liu, Z.; Chao, H. Incremental fuzzy C medoids clustering of time series data using dynamic time warping distance. PLoS ONE 2018, 13, e0197499. Available online: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0197499 (accessed on 1 June 2022). [CrossRef]
  30. D’Urso, P.; Giovanni, L.; Massari, R. Trimmed fuzzy clustering of financial time series based on dynamic time warping. Ann. Oper. Res. 2021, 299, 235–244. Available online: https://link.springer.com/article/10.1007/s10479-019-03284-1 (accessed on 1 June 2022). [CrossRef]
  31. He, G.; Pan, Y.; Xia, X.; He, J.; Peng, R.; Xiong, N.N. A Fast Semi-Supervised Clustering Framework for Large-Scale Time Series Data. IEEE Trans. 2019, 51, 4201–4216. Available online: https://link.springer.com/article/ (accessed on 1 June 2022). [CrossRef]
  32. Serra, J.; Arcos, J.L. An empirical evaluation of similarity measures for time series classification. Knowl.-Based Syst. 2014, 67, 305–314. Available online: https://www.sciencedirect.com/science/article/abs/pii/S0950705114001658 (accessed on 1 June 2022). [CrossRef]
  33. Petitjean, F.; Ketterlin, A.; Gançarski, P. A global averaging method for dynamic time warping with applications to clustering. Pattern Recognit. 2011, 44, 678–693. Available online: https://www.sciencedirect.com/science/article/abs/pii/S003132031000453X (accessed on 1 June 2022). [CrossRef]
  34. Marteau, P.F. Time warp edit distance with stiffness adjustment for time series matching. IEEE Trans. Pattern Anal. Mach. Intell. 2008, 31, 306–318. Available online: https://ieeexplore.ieee.org/abstract/document/4479483 (accessed on 1 June 2022). [CrossRef]
  35. Suris, F.N.A.; Bakar, M.A.A.; Ariff, N.M.; Nadzir, M.S.M.; Ibrahim, K. Malaysia PM10 Air Quality Time Series Clustering Based on Dynamic Time Warping. Atmosphere 2022, 13, 503. Available online: https://www.mdpi.com/2073-4433/13/4/503 (accessed on 1 June 2022). [CrossRef]
  36. Yeh, C.M.; Zhu, Y.; Ulanova, L.; Begum, N.; Ding, Y.; Dau, H.A.; Silva, D.F.; Mueen, A.; Keogh, E. Matrix profile I: All pairs similarity joins for time series: A unifying view that includes motifs, discords and shapelets. In Proceedings of the IEEE 16th International Conference on Data Mining (ICDM), Barcelona, Spain, 12–15 December 2016; pp. 1317–1322. Available online: https://ieeexplore.ieee.org/abstract/document/7837992 (accessed on 1 June 2022).
  37. Abdullah, M.; Sheng, Z.; Yan, Z.; Michael, Y.; Kaveh, K.; Krishnamurthy, V.; Chetan, K.G.; Eamonn, K. The Fastest Similarity Search Algorithm for Time Series Subsequences under Euclidean Distance. Available online: http://www.cs.unm.edu/~mueen/FastestSimilaritySearch.html (accessed on 15 March 2022).
  38. Eilers, P.; Marx, B.D. Splines, knots, and penalties. Wiley Interdiscip. Rev. Comput. Stat. 2010, 2, 637–653. Available online: https://wires.onlinelibrary.wiley.com/doi/abs/10.1002/wics.125 (accessed on 1 June 2022). [CrossRef]
  39. Cuturi, M.; Blondel, M. Soft-dtw: A differentiable loss function for time-series. Int. Conf. Mach. Learn. PMLR 2017, 70, 894–903. Available online: http://proceedings.mlr.press/v70/cuturi17a.html?source=post_page628e4799533c (accessed on 1 June 2022).
  40. Kobylin, O.; Lyashenko, V. Time series clustering based on the k-means algorithm. J. Multiapp 2020, 1, 1–7. Available online: http://newinera.com/index.php/JournalLaMultiapp/article/view/191 (accessed on 1 May 2022). [CrossRef]
  41. Paparrizos, J.; Gravano, L. k-Shape: Efficient and Accurate Clustering of Time Series. ACM Sigmod Rec. 2015, 45, 69–76. Available online: https://dl.acm.org/doi/abs/10.1145/2723372.2737793 (accessed on 1 May 2022). [CrossRef]
  42. Steinley, D. Properties of the Hubert-Arable Adjusted Rand Index. Psychol. Methods 2004, 9, 386–396. Available online: https://psycnet.apa.org/record/2004-17801-007 (accessed on 1 May 2022). [CrossRef]
  43. Fowlkes, E.B.; Mallows, C.L. A method for comparing two hierarchical clusterings. J. Am. Stat. Assoc. 1983, 78, 553–569. Available online: https://www.tandfonline.com/doi/abs/10.1080/01621459.1983.10478008 (accessed on 1 May 2022). [CrossRef]
Figure 1. Comparison of P−spline fitting processing results.
Figure 1. Comparison of P−spline fitting processing results.
Sensors 22 06163 g001
Figure 2. Flow chart of PS−WFCMdd.
Figure 2. Flow chart of PS−WFCMdd.
Sensors 22 06163 g002
Figure 3. ARI of clustering results with different clustering algorithms.
Figure 3. ARI of clustering results with different clustering algorithms.
Sensors 22 06163 g003
Figure 4. FMI of clustering results with different clustering algorithms.
Figure 4. FMI of clustering results with different clustering algorithms.
Sensors 22 06163 g004
Figure 5. NMI of clustering results with different clustering algorithms.
Figure 5. NMI of clustering results with different clustering algorithms.
Sensors 22 06163 g005
Figure 6. Display the trend of the original data and the cluster centers obtained by the PS−WFCMdd.
Figure 6. Display the trend of the original data and the cluster centers obtained by the PS−WFCMdd.
Sensors 22 06163 g006
Figure 7. Display the trend of the original data and the cluster centers obtained by the K−means−sD.
Figure 7. Display the trend of the original data and the cluster centers obtained by the K−means−sD.
Sensors 22 06163 g007
Table 1. Descriptions of ten common time series data.
Table 1. Descriptions of ten common time series data.
DatasetNumber of ClassesTime Series LengthSize of Set
Beef547060
BeetleFly251240
Car4577120
ECG200296200
FaceFour435088
Fish7463175
GunPoint2150200
Herring2512128
Lightning22637121
Meat3448120
Table 2. Descriptions of eight large-scale time series data.
Table 2. Descriptions of eight large-scale time series data.
DatasetNumber of ClassesTime Series LengthSize of Set
ChlorineConcentration31663840
ECG500051404500
HouseTwenty22000159
Mallat810242345
StarLightCurves310248236
Strawberry5235983
Rock4284470
Wafer21526174
Table 3. ARI index scores of different clustering algorithms.
Table 3. ARI index scores of different clustering algorithms.
DatasetFCMPS-K-MeansK-ShapeK-Means-sDPS-WFCMdd-sDK-Medoids-sDPS-WFCMdd
Beef0.032290.068460.015370.068460.05139−0.026050.12876
BeetleFly0.11333−0.044420.11333−0.044420.035370.161970.18711
Car0.087530.143380.091260.133450.00699−0.004750.16309
ECG2000.007570.242080.249190.283070.066600.204030.34122
FaceFour0.318350.263860.269570.122860.360960.370890.37165
Fish0.219650.170130.200090.288030.190320.174310.25790
GunPoint−0.00512−0.00512−0.00512−0.005120.015940.006000.15795
Herring−0.01435−0.01488−0.008520.078760.039360.008930.04904
Lightning20.011940.030000.000000.011940.037640.037370.03801
Meat0.562960.495990.562960.562960.757340.587450.83168
Table 4. FMI index scores of different clustering algorithms.
Table 4. FMI index scores of different clustering algorithms.
DatasetFCMPS-K-MeansK-ShapeK-Means-sDPS-WFCMdd-sDK-Medoids-sDPS-WFCMdd
Beef0.253180.264230.241210.264230.408050.204120.42276
BeetleFly0.533330.453050.533330.453050.735360.718240.74662
Car0.317880.370090.342380.418420.465970.254240.43560
ECG2000.540560.685260.675560.709170.489460.627840.66386
FaceFour0.520850.481060.477050.379070.575720.561730.56499
Fish0.368780.306370.331060.419140.377030.303040.38720
GunPoint0.495240.495240.495240.495240.420980.495290.51423
Herring0.496950.500800.515990.581250.603800.512510.54077
Lightning20.527410.602820.740290.527410.526370.544140.56615
Meat0.766570.663580.766570.766570.756380.720420.88439
Table 5. NMI index scores of different clustering algorithms.
Table 5. NMI index scores of different clustering algorithms.
DatasetFCMPS-K-MeansK-ShapeK-Means-sDPS-WFCMdd-sDK-Medoids-sDPS-WFCMdd
Beef0.253650.319470.247970.319470.255520.210940.39062
BeetleFly0.118710.007330.118710.007330.331750.353210.38343
Car0.166040.267080.222200.266040.114720.066760.31387
ECG2000.005760.140240.150850.171600.149980.173600.23618
FaceFour0.374400.443300.349800.247020.136120.502330.54600
Fish0.380360.343830.371700.402860.237390.270900.35590
GunPoint0.001110.001110.001110.001110.087020.047700.20012
Herring0.000750.002710.001280.181050.133130.130730.18397
Lightning20.015230.010260.000000.015230.083810.076390.08742
Meat0.733680.590920.733680.733680.715120.643520.83955
Table 6. Time complexity comparison.
Table 6. Time complexity comparison.
AlgorithmTime Complexity
FCM O ( N × C × K )
PS-K-means O ( N × C × K )
K-shape O N × log ( N ) × C × K
K-means-sD O N 2 × C × K
PS-WFCMdd-sD O N 2 × C × K
K-Medoids-sD O N 2 × C × K
PS-WFCMdd O N × log ( N ) × C × K
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, J.; Hou, Q.; Qu, K.; Sun, Y.; Meng, X. A Fast Weighted Fuzzy C-Medoids Clustering for Time Series Data Based on P-Splines. Sensors 2022, 22, 6163. https://doi.org/10.3390/s22166163

AMA Style

Xu J, Hou Q, Qu K, Sun Y, Meng X. A Fast Weighted Fuzzy C-Medoids Clustering for Time Series Data Based on P-Splines. Sensors. 2022; 22(16):6163. https://doi.org/10.3390/s22166163

Chicago/Turabian Style

Xu, Jiucheng, Qinchen Hou, Kanglin Qu, Yuanhao Sun, and Xiangru Meng. 2022. "A Fast Weighted Fuzzy C-Medoids Clustering for Time Series Data Based on P-Splines" Sensors 22, no. 16: 6163. https://doi.org/10.3390/s22166163

APA Style

Xu, J., Hou, Q., Qu, K., Sun, Y., & Meng, X. (2022). A Fast Weighted Fuzzy C-Medoids Clustering for Time Series Data Based on P-Splines. Sensors, 22(16), 6163. https://doi.org/10.3390/s22166163

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop