library(factoextra)
## Warning: package 'factoextra' was built under R version 4.0.5
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 4.0.5
## Warning: replacing previous import 'lifecycle::last_warnings' by
## 'rlang::last_warnings' when loading 'pillar'
## Welcome! Want to learn more? See two factoextra-related books at https://goo.gl/ve3WBa
library(cluster)
## Warning: package 'cluster' was built under R version 4.0.5
#Import dataset
library(readr)
## Warning: replacing previous import 'lifecycle::last_warnings' by
## 'rlang::last_warnings' when loading 'hms'
data <- read_csv("C:/Users/Desktop/Downloads/raw.csv")
## 
## -- Column specification --------------------------------------------------------
## cols(
##   No = col_double(),
##   year = col_double(),
##   month = col_double(),
##   day = col_double(),
##   hour = col_double(),
##   pm2.5 = col_double(),
##   DEWP = col_double(),
##   TEMP = col_double(),
##   PRES = col_double(),
##   cbwd = col_character(),
##   Iws = col_double(),
##   Is = col_double(),
##   Ir = col_double()
## )
data
#Remove Any rows with missing values
df <- na.omit(data)
#Remove any Columns that is not needed
df1 = subset(df, select = -c(cbwd,Is,Ir) )
df1
#Scale Data
df2 <- scale(df1)
df3 <- round(df2,6)
#Read the data with only first 100 rows
df4 <- head(df3,100)
df4
##               No      year     month       day      hour     pm2.5      DEWP
##   [1,] -1.758104 -1.443338 -1.596248 -1.557732 -1.661020  0.330110 -1.229777
##   [2,] -1.758025 -1.443338 -1.596248 -1.557732 -1.516613  0.536519 -1.160494
##   [3,] -1.757946 -1.443338 -1.596248 -1.557732 -1.372205  0.656019 -0.883364
##   [4,] -1.757867 -1.443338 -1.596248 -1.557732 -1.227798  0.895018 -0.606234
##   [5,] -1.757788 -1.443338 -1.596248 -1.557732 -1.083390  0.427883 -0.606234
##   [6,] -1.757709 -1.443338 -1.596248 -1.557732 -0.938983  0.112838 -0.606234
##   [7,] -1.757630 -1.443338 -1.596248 -1.557732 -0.794575  0.069384 -0.606234
##   [8,] -1.757551 -1.443338 -1.596248 -1.557732 -0.650168  0.275792 -0.606234
##   [9,] -1.757472 -1.443338 -1.596248 -1.557732 -0.505760  0.232338 -0.675516
##  [10,] -1.757393 -1.443338 -1.596248 -1.557732 -0.361352  0.362701 -0.606234
##  [11,] -1.757314 -1.443338 -1.596248 -1.557732 -0.216945  0.449610 -0.606234
##  [12,] -1.757235 -1.443338 -1.596248 -1.557732 -0.072537  0.579974 -0.675516
##  [13,] -1.757156 -1.443338 -1.596248 -1.557732  0.071870  0.536519 -0.675516
##  [14,] -1.757077 -1.443338 -1.596248 -1.557732  0.216278  0.710337 -0.675516
##  [15,] -1.756998 -1.443338 -1.596248 -1.557732  0.360685  0.645155 -0.744799
##  [16,] -1.756919 -1.443338 -1.596248 -1.557732  0.505093  0.601701 -0.744799
##  [17,] -1.756840 -1.443338 -1.596248 -1.557732  0.649500  0.656019 -0.744799
##  [18,] -1.756761 -1.443338 -1.596248 -1.557732  0.793908  0.710337 -0.675516
##  [19,] -1.756682 -1.443338 -1.596248 -1.557732  0.938315  0.775519 -0.675516
##  [20,] -1.756603 -1.443338 -1.596248 -1.557732  1.082723  0.547383 -0.675516
##  [21,] -1.756524 -1.443338 -1.596248 -1.557732  1.227130  0.601701 -0.606234
##  [22,] -1.756445 -1.443338 -1.596248 -1.557732  1.371538  0.710337 -0.606234
##  [23,] -1.756366 -1.443338 -1.596248 -1.557732  1.515945  0.623428 -0.675516
##  [24,] -1.756287 -1.443338 -1.596248 -1.557732  1.660353  0.297520 -0.675516
##  [25,] -1.756208 -1.443338 -1.596248 -1.443908 -1.661020 -0.093571 -0.606234
##  [26,] -1.756129 -1.443338 -1.596248 -1.443908 -1.516613 -0.386888 -0.675516
##  [27,] -1.756050 -1.443338 -1.596248 -1.443908 -1.372205 -0.365161 -0.675516
##  [28,] -1.755971 -1.443338 -1.596248 -1.443908 -1.227798 -0.473797 -0.675516
##  [29,] -1.755892 -1.443338 -1.596248 -1.443908 -1.083390 -0.365161 -0.675516
##  [30,] -1.755813 -1.443338 -1.596248 -1.443908 -0.938983 -0.169616 -0.744799
##  [31,] -1.755734 -1.443338 -1.596248 -1.443908 -0.794575 -0.082707 -0.814081
##  [32,] -1.755655 -1.443338 -1.596248 -1.443908 -0.650168 -0.137025 -0.814081
##  [33,] -1.755576 -1.443338 -1.596248 -1.443908 -0.505760 -0.180480 -0.814081
##  [34,] -1.755497 -1.443338 -1.596248 -1.443908 -0.361352 -0.137025 -0.883364
##  [35,] -1.755418 -1.443338 -1.596248 -1.443908 -0.216945 -0.223934 -0.883364
##  [36,] -1.755339 -1.443338 -1.596248 -1.443908 -0.072537 -0.006662 -0.883364
##  [37,] -1.755260 -1.443338 -1.596248 -1.443908  0.071870  0.091111 -0.883364
##  [38,] -1.755181 -1.443338 -1.596248 -1.443908  0.216278 -0.093571 -0.883364
##  [39,] -1.755102 -1.443338 -1.596248 -1.443908  0.360685 -0.028389 -0.883364
##  [40,] -1.755023 -1.443338 -1.596248 -1.443908  0.505093 -0.039253 -0.883364
##  [41,] -1.754944 -1.443338 -1.596248 -1.443908  0.649500 -0.137025 -0.883364
##  [42,] -1.754865 -1.443338 -1.596248 -1.443908  0.793908 -0.310843 -0.883364
##  [43,] -1.754786 -1.443338 -1.596248 -1.443908  0.938315 -0.408615 -0.883364
##  [44,] -1.754707 -1.443338 -1.596248 -1.443908  1.082723 -0.495524 -0.883364
##  [45,] -1.754628 -1.443338 -1.596248 -1.443908  1.227130 -0.299979 -0.814081
##  [46,] -1.754549 -1.443338 -1.596248 -1.443908  1.371538 -0.289116 -0.883364
##  [47,] -1.754470 -1.443338 -1.596248 -1.443908  1.515945 -0.245661 -0.883364
##  [48,] -1.754391 -1.443338 -1.596248 -1.443908  1.660353 -0.278252 -0.952646
##  [49,] -1.754312 -1.443338 -1.596248 -1.330085 -1.661020 -0.213070 -1.091211
##  [50,] -1.754233 -1.443338 -1.596248 -1.330085 -1.516613 -0.441206 -1.229777
##  [51,] -1.754154 -1.443338 -1.596248 -1.330085 -1.372205 -0.799706 -1.299059
##  [52,] -1.754075 -1.443338 -1.596248 -1.330085 -1.227798 -0.788842 -1.368342
##  [53,] -1.753996 -1.443338 -1.596248 -1.330085 -1.083390 -0.767115 -1.437624
##  [54,] -1.753917 -1.443338 -1.596248 -1.330085 -0.938983 -0.788842 -1.506907
##  [55,] -1.753838 -1.443338 -1.596248 -1.330085 -0.794575 -0.854024 -1.576189
##  [56,] -1.753759 -1.443338 -1.596248 -1.330085 -0.650168 -0.756251 -1.576189
##  [57,] -1.753680 -1.443338 -1.596248 -1.330085 -0.505760 -0.788842 -1.645472
##  [58,] -1.753601 -1.443338 -1.596248 -1.330085 -0.361352 -0.777978 -1.645472
##  [59,] -1.753522 -1.443338 -1.596248 -1.330085 -0.216945 -0.777978 -1.645472
##  [60,] -1.753443 -1.443338 -1.596248 -1.330085 -0.072537 -0.799706 -1.714754
##  [61,] -1.753364 -1.443338 -1.596248 -1.330085  0.071870 -0.756251 -1.576189
##  [62,] -1.753285 -1.443338 -1.596248 -1.330085  0.216278 -0.723660 -1.506907
##  [63,] -1.753206 -1.443338 -1.596248 -1.330085  0.360685 -0.767115 -1.576189
##  [64,] -1.753127 -1.443338 -1.596248 -1.330085  0.505093 -0.756251 -1.576189
##  [65,] -1.753048 -1.443338 -1.596248 -1.330085  0.649500 -0.745388 -1.576189
##  [66,] -1.752969 -1.443338 -1.596248 -1.330085  0.793908 -0.745388 -1.506907
##  [67,] -1.752890 -1.443338 -1.596248 -1.330085  0.938315 -0.767115 -1.714754
##  [68,] -1.752811 -1.443338 -1.596248 -1.330085  1.082723 -0.788842 -1.576189
##  [69,] -1.752732 -1.443338 -1.596248 -1.330085  1.227130 -0.734524 -1.784037
##  [70,] -1.752653 -1.443338 -1.596248 -1.330085  1.371538 -0.712797 -1.784037
##  [71,] -1.752574 -1.443338 -1.596248 -1.330085  1.515945 -0.756251 -1.784037
##  [72,] -1.752495 -1.443338 -1.596248 -1.330085  1.660353 -0.734524 -1.922602
##  [73,] -1.752416 -1.443338 -1.596248 -1.216262 -1.661020 -0.745388 -1.922602
##  [74,] -1.752337 -1.443338 -1.596248 -1.216262 -1.516613 -0.701933 -1.922602
##  [75,] -1.752258 -1.443338 -1.596248 -1.216262 -1.372205 -0.777978 -1.922602
##  [76,] -1.752179 -1.443338 -1.596248 -1.216262 -1.227798 -0.799706 -1.991884
##  [77,] -1.752100 -1.443338 -1.596248 -1.216262 -1.083390 -0.767115 -1.991884
##  [78,] -1.752021 -1.443338 -1.596248 -1.216262 -0.938983 -0.767115 -1.991884
##  [79,] -1.751942 -1.443338 -1.596248 -1.216262 -0.794575 -0.777978 -1.922602
##  [80,] -1.751863 -1.443338 -1.596248 -1.216262 -0.650168 -0.777978 -1.991884
##  [81,] -1.751784 -1.443338 -1.596248 -1.216262 -0.505760 -0.777978 -1.922602
##  [82,] -1.751705 -1.443338 -1.596248 -1.216262 -0.361352 -0.756251 -1.922602
##  [83,] -1.751626 -1.443338 -1.596248 -1.216262 -0.216945 -0.680206 -1.853319
##  [84,] -1.751547 -1.443338 -1.596248 -1.216262 -0.072537 -0.745388 -1.853319
##  [85,] -1.751468 -1.443338 -1.596248 -1.216262  0.071870 -0.777978 -1.853319
##  [86,] -1.751389 -1.443338 -1.596248 -1.216262  0.216278 -0.647615 -1.784037
##  [87,] -1.751310 -1.443338 -1.596248 -1.216262  0.360685 -0.625888 -1.645472
##  [88,] -1.751231 -1.443338 -1.596248 -1.216262  0.505093 -0.712797 -1.714754
##  [89,] -1.751152 -1.443338 -1.596248 -1.216262  0.649500 -0.528115 -1.784037
##  [90,] -1.751073 -1.443338 -1.596248 -1.216262  0.793908 -0.462934 -1.714754
##  [91,] -1.750994 -1.443338 -1.596248 -1.216262  0.938315 -0.430343 -1.714754
##  [92,] -1.750915 -1.443338 -1.596248 -1.216262  1.082723 -0.419479 -1.645472
##  [93,] -1.750836 -1.443338 -1.596248 -1.216262  1.227130 -0.158752 -1.645472
##  [94,] -1.750757 -1.443338 -1.596248 -1.216262  1.371538  0.080247 -1.784037
##  [95,] -1.750678 -1.443338 -1.596248 -1.216262  1.515945 -0.354297 -1.645472
##  [96,] -1.750599 -1.443338 -1.596248 -1.216262  1.660353 -0.528115 -1.645472
##  [97,] -1.750520 -1.443338 -1.596248 -1.102438 -1.661020 -0.462934 -1.853319
##  [98,] -1.750441 -1.443338 -1.596248 -1.102438 -1.516613 -0.234798 -1.853319
##  [99,] -1.750362 -1.443338 -1.596248 -1.102438 -1.372205 -0.528115 -1.922602
## [100,] -1.750283 -1.443338 -1.596248 -1.102438 -1.227798 -0.593297 -1.922602
##             TEMP     PRES       Iws
##   [1,] -1.347127 0.345325 -0.444939
##   [2,] -1.347127 0.345325 -0.427002
##   [3,] -1.429261 0.442406 -0.409064
##   [4,] -1.429261 0.539486 -0.372988
##   [5,] -1.429261 0.539486 -0.355051
##   [6,] -1.511395 0.539486 -0.337114
##   [7,] -1.511395 0.636567 -0.301038
##   [8,] -1.429261 0.733647 -0.264962
##   [9,] -1.511395 0.733647 -0.228886
##  [10,] -1.429261 0.830728 -0.192810
##  [11,] -1.429261 0.927808 -0.129727
##  [12,] -1.429261 0.927808 -0.066645
##  [13,] -1.429261 0.927808 -0.003562
##  [14,] -1.429261 0.830728  0.077458
##  [15,] -1.429261 0.830728  0.158477
##  [16,] -1.429261 0.830728  0.239497
##  [17,] -1.429261 0.927808  0.275573
##  [18,] -1.429261 1.024889  0.311649
##  [19,] -1.429261 1.024889  0.374732
##  [20,] -1.429261 1.121969  0.410808
##  [21,] -1.429261 1.121969  0.446884
##  [22,] -1.429261 1.024889  0.509966
##  [23,] -1.511395 1.121969  0.573049
##  [24,] -1.511395 1.024889  0.636132
##  [25,] -1.511395 1.024889  0.699214
##  [26,] -1.511395 0.927808  0.762297
##  [27,] -1.593529 0.927808  0.843317
##  [28,] -1.593529 0.830728  0.906399
##  [29,] -1.593529 0.733647  0.987419
##  [30,] -1.675663 0.733647  1.068439
##  [31,] -1.675663 0.733647  1.149459
##  [32,] -1.757797 0.733647  1.230478
##  [33,] -1.757797 0.733647  1.329637
##  [34,] -1.757797 0.636567  1.410657
##  [35,] -1.757797 0.636567  1.491677
##  [36,] -1.757797 0.539486  1.590835
##  [37,] -1.757797 0.442406  1.653918
##  [38,] -1.757797 0.345325  1.771013
##  [39,] -1.757797 0.345325  1.870172
##  [40,] -1.757797 0.345325  1.969331
##  [41,] -1.757797 0.345325  2.032413
##  [42,] -1.757797 0.345325  2.095496
##  [43,] -1.757797 0.442406 -0.463078
##  [44,] -1.757797 0.539486 -0.445140
##  [45,] -1.757797 0.539486 -0.399995
##  [46,] -1.839931 0.636567 -0.336912
##  [47,] -1.757797 0.636567 -0.255893
##  [48,] -1.922066 0.636567 -0.192810
##  [49,] -2.004200 0.636567 -0.156734
##  [50,] -1.757797 0.636567 -0.039638
##  [51,] -1.839931 0.733647  0.122603
##  [52,] -1.922066 0.733647  0.320719
##  [53,] -1.922066 0.830728  0.500897
##  [54,] -2.004200 0.927808  0.644999
##  [55,] -2.004200 1.024889  0.825178
##  [56,] -2.086334 1.024889  1.005356
##  [57,] -2.086334 1.121969  1.149459
##  [58,] -2.086334 1.219050  1.347574
##  [59,] -2.004200 1.316130  1.545690
##  [60,] -2.004200 1.413211  1.707931
##  [61,] -1.922066 1.316130  1.888109
##  [62,] -1.839931 1.316130  2.086225
##  [63,] -1.839931 1.316130  2.266403
##  [64,] -1.757797 1.316130  2.446582
##  [65,] -1.757797 1.413211  2.590684
##  [66,] -1.922066 1.510291  2.734786
##  [67,] -1.922066 1.510291  2.851882
##  [68,] -2.004200 1.607372  2.968978
##  [69,] -2.004200 1.704452  3.149156
##  [70,] -2.086334 1.704452  3.293259
##  [71,] -2.086334 1.801532  3.455500
##  [72,] -2.250602 1.801532  3.518583
##  [73,] -2.414870 1.801532  3.581665
##  [74,] -2.497004 1.801532  3.662685
##  [75,] -2.579138 1.801532  3.725768
##  [76,] -2.497004 1.801532  3.824926
##  [77,] -2.579138 1.801532  3.924085
##  [78,] -2.332736 1.704452 -0.381856
##  [79,] -2.332736 1.801532 -0.318774
##  [80,] -2.332736 1.704452 -0.201678
##  [81,] -2.332736 1.801532 -0.102519
##  [82,] -2.250602 1.801532  0.014577
##  [83,] -2.168468 1.801532  0.113735
##  [84,] -2.086334 1.801532  0.212894
##  [85,] -2.004200 1.704452  0.312052
##  [86,] -1.922066 1.510291  0.348128
##  [87,] -1.922066 1.510291 -0.463078
##  [88,] -1.922066 1.413211 -0.444939
##  [89,] -1.922066 1.413211 -0.408863
##  [90,] -1.922066 1.413211 -0.372787
##  [91,] -1.922066 1.510291 -0.336711
##  [92,] -2.086334 1.607372 -0.273628
##  [93,] -2.004200 1.607372 -0.210546
##  [94,] -2.497004 1.607372 -0.147463
##  [95,] -2.086334 1.704452 -0.066443
##  [96,] -2.332736 1.607372 -0.003361
##  [97,] -2.414870 1.607372  0.059722
##  [98,] -2.168468 1.607372 -0.399995
##  [99,] -2.168468 1.704452 -0.318975
## [100,] -2.168468 1.607372 -0.201879
#Calculate gap statistic based on number of cluster
gap_stat <- clusGap(df4, FUN=kmeans, nstart = 25, K.max = 10, B= 50)
#Plot number of clusters vs. gap statistic
fviz_gap_stat(gap_stat)

#perform K-Means clustering with optimal K

#Make this example reproducible
set.seed(2341)

#Perform K-means clustering with k=2 clusters
km <- kmeans(df4, centers=2,nstart=25)
km
## K-means clustering with 2 clusters of sizes 24, 76
## 
## Cluster means:
##          No      year     month       day       hour      pm2.5      DEWP
## 1 -1.753272 -1.443338 -1.596248 -1.330085  0.2403456 -0.6236246 -1.558869
## 2 -1.754485 -1.443338 -1.596248 -1.389992 -0.1523416 -0.1558934 -1.162317
##        TEMP     PRES       Iws
## 1 -2.031578 1.324220 2.7067217
## 2 -1.790219 1.037662 0.2057682
## 
## Clustering vector:
##   [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
##  [38] 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [75] 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
## 
## Within cluster sum of squares by cluster:
## [1]  49.57878 169.37811
##  (between_SS / total_SS =  36.6 %)
## 
## Available components:
## 
## [1] "cluster"      "centers"      "totss"        "withinss"     "tot.withinss"
## [6] "betweenss"    "size"         "iter"         "ifault"
#Find mean of each cluster
df5 <- aggregate(df4,by=list(cluster=km$cluster),mean)
round(df5,4)
#add cluster assignment to original data
final_data <- cbind(df4,cluster = km$cluster)
final_data
##               No      year     month       day      hour     pm2.5      DEWP
##   [1,] -1.758104 -1.443338 -1.596248 -1.557732 -1.661020  0.330110 -1.229777
##   [2,] -1.758025 -1.443338 -1.596248 -1.557732 -1.516613  0.536519 -1.160494
##   [3,] -1.757946 -1.443338 -1.596248 -1.557732 -1.372205  0.656019 -0.883364
##   [4,] -1.757867 -1.443338 -1.596248 -1.557732 -1.227798  0.895018 -0.606234
##   [5,] -1.757788 -1.443338 -1.596248 -1.557732 -1.083390  0.427883 -0.606234
##   [6,] -1.757709 -1.443338 -1.596248 -1.557732 -0.938983  0.112838 -0.606234
##   [7,] -1.757630 -1.443338 -1.596248 -1.557732 -0.794575  0.069384 -0.606234
##   [8,] -1.757551 -1.443338 -1.596248 -1.557732 -0.650168  0.275792 -0.606234
##   [9,] -1.757472 -1.443338 -1.596248 -1.557732 -0.505760  0.232338 -0.675516
##  [10,] -1.757393 -1.443338 -1.596248 -1.557732 -0.361352  0.362701 -0.606234
##  [11,] -1.757314 -1.443338 -1.596248 -1.557732 -0.216945  0.449610 -0.606234
##  [12,] -1.757235 -1.443338 -1.596248 -1.557732 -0.072537  0.579974 -0.675516
##  [13,] -1.757156 -1.443338 -1.596248 -1.557732  0.071870  0.536519 -0.675516
##  [14,] -1.757077 -1.443338 -1.596248 -1.557732  0.216278  0.710337 -0.675516
##  [15,] -1.756998 -1.443338 -1.596248 -1.557732  0.360685  0.645155 -0.744799
##  [16,] -1.756919 -1.443338 -1.596248 -1.557732  0.505093  0.601701 -0.744799
##  [17,] -1.756840 -1.443338 -1.596248 -1.557732  0.649500  0.656019 -0.744799
##  [18,] -1.756761 -1.443338 -1.596248 -1.557732  0.793908  0.710337 -0.675516
##  [19,] -1.756682 -1.443338 -1.596248 -1.557732  0.938315  0.775519 -0.675516
##  [20,] -1.756603 -1.443338 -1.596248 -1.557732  1.082723  0.547383 -0.675516
##  [21,] -1.756524 -1.443338 -1.596248 -1.557732  1.227130  0.601701 -0.606234
##  [22,] -1.756445 -1.443338 -1.596248 -1.557732  1.371538  0.710337 -0.606234
##  [23,] -1.756366 -1.443338 -1.596248 -1.557732  1.515945  0.623428 -0.675516
##  [24,] -1.756287 -1.443338 -1.596248 -1.557732  1.660353  0.297520 -0.675516
##  [25,] -1.756208 -1.443338 -1.596248 -1.443908 -1.661020 -0.093571 -0.606234
##  [26,] -1.756129 -1.443338 -1.596248 -1.443908 -1.516613 -0.386888 -0.675516
##  [27,] -1.756050 -1.443338 -1.596248 -1.443908 -1.372205 -0.365161 -0.675516
##  [28,] -1.755971 -1.443338 -1.596248 -1.443908 -1.227798 -0.473797 -0.675516
##  [29,] -1.755892 -1.443338 -1.596248 -1.443908 -1.083390 -0.365161 -0.675516
##  [30,] -1.755813 -1.443338 -1.596248 -1.443908 -0.938983 -0.169616 -0.744799
##  [31,] -1.755734 -1.443338 -1.596248 -1.443908 -0.794575 -0.082707 -0.814081
##  [32,] -1.755655 -1.443338 -1.596248 -1.443908 -0.650168 -0.137025 -0.814081
##  [33,] -1.755576 -1.443338 -1.596248 -1.443908 -0.505760 -0.180480 -0.814081
##  [34,] -1.755497 -1.443338 -1.596248 -1.443908 -0.361352 -0.137025 -0.883364
##  [35,] -1.755418 -1.443338 -1.596248 -1.443908 -0.216945 -0.223934 -0.883364
##  [36,] -1.755339 -1.443338 -1.596248 -1.443908 -0.072537 -0.006662 -0.883364
##  [37,] -1.755260 -1.443338 -1.596248 -1.443908  0.071870  0.091111 -0.883364
##  [38,] -1.755181 -1.443338 -1.596248 -1.443908  0.216278 -0.093571 -0.883364
##  [39,] -1.755102 -1.443338 -1.596248 -1.443908  0.360685 -0.028389 -0.883364
##  [40,] -1.755023 -1.443338 -1.596248 -1.443908  0.505093 -0.039253 -0.883364
##  [41,] -1.754944 -1.443338 -1.596248 -1.443908  0.649500 -0.137025 -0.883364
##  [42,] -1.754865 -1.443338 -1.596248 -1.443908  0.793908 -0.310843 -0.883364
##  [43,] -1.754786 -1.443338 -1.596248 -1.443908  0.938315 -0.408615 -0.883364
##  [44,] -1.754707 -1.443338 -1.596248 -1.443908  1.082723 -0.495524 -0.883364
##  [45,] -1.754628 -1.443338 -1.596248 -1.443908  1.227130 -0.299979 -0.814081
##  [46,] -1.754549 -1.443338 -1.596248 -1.443908  1.371538 -0.289116 -0.883364
##  [47,] -1.754470 -1.443338 -1.596248 -1.443908  1.515945 -0.245661 -0.883364
##  [48,] -1.754391 -1.443338 -1.596248 -1.443908  1.660353 -0.278252 -0.952646
##  [49,] -1.754312 -1.443338 -1.596248 -1.330085 -1.661020 -0.213070 -1.091211
##  [50,] -1.754233 -1.443338 -1.596248 -1.330085 -1.516613 -0.441206 -1.229777
##  [51,] -1.754154 -1.443338 -1.596248 -1.330085 -1.372205 -0.799706 -1.299059
##  [52,] -1.754075 -1.443338 -1.596248 -1.330085 -1.227798 -0.788842 -1.368342
##  [53,] -1.753996 -1.443338 -1.596248 -1.330085 -1.083390 -0.767115 -1.437624
##  [54,] -1.753917 -1.443338 -1.596248 -1.330085 -0.938983 -0.788842 -1.506907
##  [55,] -1.753838 -1.443338 -1.596248 -1.330085 -0.794575 -0.854024 -1.576189
##  [56,] -1.753759 -1.443338 -1.596248 -1.330085 -0.650168 -0.756251 -1.576189
##  [57,] -1.753680 -1.443338 -1.596248 -1.330085 -0.505760 -0.788842 -1.645472
##  [58,] -1.753601 -1.443338 -1.596248 -1.330085 -0.361352 -0.777978 -1.645472
##  [59,] -1.753522 -1.443338 -1.596248 -1.330085 -0.216945 -0.777978 -1.645472
##  [60,] -1.753443 -1.443338 -1.596248 -1.330085 -0.072537 -0.799706 -1.714754
##  [61,] -1.753364 -1.443338 -1.596248 -1.330085  0.071870 -0.756251 -1.576189
##  [62,] -1.753285 -1.443338 -1.596248 -1.330085  0.216278 -0.723660 -1.506907
##  [63,] -1.753206 -1.443338 -1.596248 -1.330085  0.360685 -0.767115 -1.576189
##  [64,] -1.753127 -1.443338 -1.596248 -1.330085  0.505093 -0.756251 -1.576189
##  [65,] -1.753048 -1.443338 -1.596248 -1.330085  0.649500 -0.745388 -1.576189
##  [66,] -1.752969 -1.443338 -1.596248 -1.330085  0.793908 -0.745388 -1.506907
##  [67,] -1.752890 -1.443338 -1.596248 -1.330085  0.938315 -0.767115 -1.714754
##  [68,] -1.752811 -1.443338 -1.596248 -1.330085  1.082723 -0.788842 -1.576189
##  [69,] -1.752732 -1.443338 -1.596248 -1.330085  1.227130 -0.734524 -1.784037
##  [70,] -1.752653 -1.443338 -1.596248 -1.330085  1.371538 -0.712797 -1.784037
##  [71,] -1.752574 -1.443338 -1.596248 -1.330085  1.515945 -0.756251 -1.784037
##  [72,] -1.752495 -1.443338 -1.596248 -1.330085  1.660353 -0.734524 -1.922602
##  [73,] -1.752416 -1.443338 -1.596248 -1.216262 -1.661020 -0.745388 -1.922602
##  [74,] -1.752337 -1.443338 -1.596248 -1.216262 -1.516613 -0.701933 -1.922602
##  [75,] -1.752258 -1.443338 -1.596248 -1.216262 -1.372205 -0.777978 -1.922602
##  [76,] -1.752179 -1.443338 -1.596248 -1.216262 -1.227798 -0.799706 -1.991884
##  [77,] -1.752100 -1.443338 -1.596248 -1.216262 -1.083390 -0.767115 -1.991884
##  [78,] -1.752021 -1.443338 -1.596248 -1.216262 -0.938983 -0.767115 -1.991884
##  [79,] -1.751942 -1.443338 -1.596248 -1.216262 -0.794575 -0.777978 -1.922602
##  [80,] -1.751863 -1.443338 -1.596248 -1.216262 -0.650168 -0.777978 -1.991884
##  [81,] -1.751784 -1.443338 -1.596248 -1.216262 -0.505760 -0.777978 -1.922602
##  [82,] -1.751705 -1.443338 -1.596248 -1.216262 -0.361352 -0.756251 -1.922602
##  [83,] -1.751626 -1.443338 -1.596248 -1.216262 -0.216945 -0.680206 -1.853319
##  [84,] -1.751547 -1.443338 -1.596248 -1.216262 -0.072537 -0.745388 -1.853319
##  [85,] -1.751468 -1.443338 -1.596248 -1.216262  0.071870 -0.777978 -1.853319
##  [86,] -1.751389 -1.443338 -1.596248 -1.216262  0.216278 -0.647615 -1.784037
##  [87,] -1.751310 -1.443338 -1.596248 -1.216262  0.360685 -0.625888 -1.645472
##  [88,] -1.751231 -1.443338 -1.596248 -1.216262  0.505093 -0.712797 -1.714754
##  [89,] -1.751152 -1.443338 -1.596248 -1.216262  0.649500 -0.528115 -1.784037
##  [90,] -1.751073 -1.443338 -1.596248 -1.216262  0.793908 -0.462934 -1.714754
##  [91,] -1.750994 -1.443338 -1.596248 -1.216262  0.938315 -0.430343 -1.714754
##  [92,] -1.750915 -1.443338 -1.596248 -1.216262  1.082723 -0.419479 -1.645472
##  [93,] -1.750836 -1.443338 -1.596248 -1.216262  1.227130 -0.158752 -1.645472
##  [94,] -1.750757 -1.443338 -1.596248 -1.216262  1.371538  0.080247 -1.784037
##  [95,] -1.750678 -1.443338 -1.596248 -1.216262  1.515945 -0.354297 -1.645472
##  [96,] -1.750599 -1.443338 -1.596248 -1.216262  1.660353 -0.528115 -1.645472
##  [97,] -1.750520 -1.443338 -1.596248 -1.102438 -1.661020 -0.462934 -1.853319
##  [98,] -1.750441 -1.443338 -1.596248 -1.102438 -1.516613 -0.234798 -1.853319
##  [99,] -1.750362 -1.443338 -1.596248 -1.102438 -1.372205 -0.528115 -1.922602
## [100,] -1.750283 -1.443338 -1.596248 -1.102438 -1.227798 -0.593297 -1.922602
##             TEMP     PRES       Iws cluster
##   [1,] -1.347127 0.345325 -0.444939       2
##   [2,] -1.347127 0.345325 -0.427002       2
##   [3,] -1.429261 0.442406 -0.409064       2
##   [4,] -1.429261 0.539486 -0.372988       2
##   [5,] -1.429261 0.539486 -0.355051       2
##   [6,] -1.511395 0.539486 -0.337114       2
##   [7,] -1.511395 0.636567 -0.301038       2
##   [8,] -1.429261 0.733647 -0.264962       2
##   [9,] -1.511395 0.733647 -0.228886       2
##  [10,] -1.429261 0.830728 -0.192810       2
##  [11,] -1.429261 0.927808 -0.129727       2
##  [12,] -1.429261 0.927808 -0.066645       2
##  [13,] -1.429261 0.927808 -0.003562       2
##  [14,] -1.429261 0.830728  0.077458       2
##  [15,] -1.429261 0.830728  0.158477       2
##  [16,] -1.429261 0.830728  0.239497       2
##  [17,] -1.429261 0.927808  0.275573       2
##  [18,] -1.429261 1.024889  0.311649       2
##  [19,] -1.429261 1.024889  0.374732       2
##  [20,] -1.429261 1.121969  0.410808       2
##  [21,] -1.429261 1.121969  0.446884       2
##  [22,] -1.429261 1.024889  0.509966       2
##  [23,] -1.511395 1.121969  0.573049       2
##  [24,] -1.511395 1.024889  0.636132       2
##  [25,] -1.511395 1.024889  0.699214       2
##  [26,] -1.511395 0.927808  0.762297       2
##  [27,] -1.593529 0.927808  0.843317       2
##  [28,] -1.593529 0.830728  0.906399       2
##  [29,] -1.593529 0.733647  0.987419       2
##  [30,] -1.675663 0.733647  1.068439       2
##  [31,] -1.675663 0.733647  1.149459       2
##  [32,] -1.757797 0.733647  1.230478       2
##  [33,] -1.757797 0.733647  1.329637       2
##  [34,] -1.757797 0.636567  1.410657       2
##  [35,] -1.757797 0.636567  1.491677       2
##  [36,] -1.757797 0.539486  1.590835       2
##  [37,] -1.757797 0.442406  1.653918       2
##  [38,] -1.757797 0.345325  1.771013       1
##  [39,] -1.757797 0.345325  1.870172       1
##  [40,] -1.757797 0.345325  1.969331       1
##  [41,] -1.757797 0.345325  2.032413       1
##  [42,] -1.757797 0.345325  2.095496       1
##  [43,] -1.757797 0.442406 -0.463078       2
##  [44,] -1.757797 0.539486 -0.445140       2
##  [45,] -1.757797 0.539486 -0.399995       2
##  [46,] -1.839931 0.636567 -0.336912       2
##  [47,] -1.757797 0.636567 -0.255893       2
##  [48,] -1.922066 0.636567 -0.192810       2
##  [49,] -2.004200 0.636567 -0.156734       2
##  [50,] -1.757797 0.636567 -0.039638       2
##  [51,] -1.839931 0.733647  0.122603       2
##  [52,] -1.922066 0.733647  0.320719       2
##  [53,] -1.922066 0.830728  0.500897       2
##  [54,] -2.004200 0.927808  0.644999       2
##  [55,] -2.004200 1.024889  0.825178       2
##  [56,] -2.086334 1.024889  1.005356       2
##  [57,] -2.086334 1.121969  1.149459       2
##  [58,] -2.086334 1.219050  1.347574       2
##  [59,] -2.004200 1.316130  1.545690       1
##  [60,] -2.004200 1.413211  1.707931       1
##  [61,] -1.922066 1.316130  1.888109       1
##  [62,] -1.839931 1.316130  2.086225       1
##  [63,] -1.839931 1.316130  2.266403       1
##  [64,] -1.757797 1.316130  2.446582       1
##  [65,] -1.757797 1.413211  2.590684       1
##  [66,] -1.922066 1.510291  2.734786       1
##  [67,] -1.922066 1.510291  2.851882       1
##  [68,] -2.004200 1.607372  2.968978       1
##  [69,] -2.004200 1.704452  3.149156       1
##  [70,] -2.086334 1.704452  3.293259       1
##  [71,] -2.086334 1.801532  3.455500       1
##  [72,] -2.250602 1.801532  3.518583       1
##  [73,] -2.414870 1.801532  3.581665       1
##  [74,] -2.497004 1.801532  3.662685       1
##  [75,] -2.579138 1.801532  3.725768       1
##  [76,] -2.497004 1.801532  3.824926       1
##  [77,] -2.579138 1.801532  3.924085       1
##  [78,] -2.332736 1.704452 -0.381856       2
##  [79,] -2.332736 1.801532 -0.318774       2
##  [80,] -2.332736 1.704452 -0.201678       2
##  [81,] -2.332736 1.801532 -0.102519       2
##  [82,] -2.250602 1.801532  0.014577       2
##  [83,] -2.168468 1.801532  0.113735       2
##  [84,] -2.086334 1.801532  0.212894       2
##  [85,] -2.004200 1.704452  0.312052       2
##  [86,] -1.922066 1.510291  0.348128       2
##  [87,] -1.922066 1.510291 -0.463078       2
##  [88,] -1.922066 1.413211 -0.444939       2
##  [89,] -1.922066 1.413211 -0.408863       2
##  [90,] -1.922066 1.413211 -0.372787       2
##  [91,] -1.922066 1.510291 -0.336711       2
##  [92,] -2.086334 1.607372 -0.273628       2
##  [93,] -2.004200 1.607372 -0.210546       2
##  [94,] -2.497004 1.607372 -0.147463       2
##  [95,] -2.086334 1.704452 -0.066443       2
##  [96,] -2.332736 1.607372 -0.003361       2
##  [97,] -2.414870 1.607372  0.059722       2
##  [98,] -2.168468 1.607372 -0.399995       2
##  [99,] -2.168468 1.704452 -0.318975       2
## [100,] -2.168468 1.607372 -0.201879       2