toy-example

TOY EXAMPLE

This script shows a toy example of usage of SparseFunClust (without alignment).

Generate the data

library(SparseFunClust)
set.seed(24032023)
n <- 50
x <- seq(0,1,len=500)
out <- generate.data.FV17(n, x)
data <- out$data
trueClust <- out$true.partition
matplot(x, t(data), type='l', col=trueClust,
        xlab = 'x', ylab = 'data', main = 'Simulated data')

Run Sparse Functional Clustering (no alignment)

K <- 2            # run with 2 groups only
method <- 'kmea'  # version with K-means clustering
tuning.m <- FALSE # don't perform tuning of the sparsity parameter (faster)
result <- SparseFunClust(data, x, K = K, do.alignment = FALSE,
                         clust.method = method, tuning.m = tuning.m)

Plot / explore results

table(trueClust,result$labels)

##          
## trueClust  1  2
##         1 50  0
##         2 10 40

cer(trueClust,result$labels)

## [1] 0.1818182

matplot(x,t(data),type='l',lty=1,col=result$labels+1,ylab='',
        main='clustering results')
lines(x,colMeans(data[which(result$labels==1),]),lwd=2)
lines(x,colMeans(data[which(result$labels==2),]),lwd=2)

plot(x,result$w,type='l',lty=1,lwd=2,ylab='',
     main='estimated weighting function')
abline(v=0.5)