The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.

sasfunclust

R build status Travis build status The package sasfunclust implements the sparse and smooth functional clustering (SaS-Funclust) method proposed by Centofanti et al. (2021). SaS-Funclust is a new method for clustering functional data that aims to classify a sample of curves into homogeneous groups while jointly detecting the most informative portions of domain. The method relies on a general functional Gaussian mixture model whose parameters are estimated by maximizing a log-likelihood function penalized with the functional adaptive pairwise fusion penalty (FAPFP) and a roughness penalty. The package comprises two main functions and . The former performs the SaS-Funclust for fixed number of clusters G, tuning parameter of the smoothness penalty λs, and tuning parameter of the FAPFP λl. The latter executes the K-fold cross-validation procedure described in Centofanti et al. (2021) to choose G, λs, and λl.

Installation

The development version can be installed from GitHub with:

# install.packages("devtools")
devtools::install_github("unina-sfere/sasfunclust")

Example

This is a basic example which shows you how to apply the two main functions sasfclust and sasfclust_cv on a synthetic dataset generated as described in the simulation study of Centofanti et al. (2021).

We start by loading and attaching the sasfunclust package.

library(sasfunclust)

Then, we generate the synthetic dataset as follows.

n_i=20
train<-simulate_data("Scenario I",n_i=n_i,var_e = 1,var_b = 0.5^2)

To apply sasfclust_cv, sequences of G, λs, and λl should be defined.

lambda_s_seq=10^seq(-4,-3)
lambda_l_seq=10^seq(0,3)
G_seq=2

And, then, sasfclust_cv is executed.

mod_cv<-sasfclust_cv(X=train$X,grid=train$grid,G_seq=G_seq,
lambda_l_seq = lambda_l_seq,lambda_s_seq =lambda_s_seq,maxit = 50,K_fold = 5,q=30,ncores = 7)

The results are plotted.

plot(mod_cv)

By using the model selection method described in Centofanti et al. (2021), the optimal values of G, λs, and λl, are 2, 10 − 4, and 102, respectively.

Finally, sasfclust is applied with G, λs, and λl fixed to their optimal values.

mod<-sasfclust(X=train$X,grid=train$grid,G=mod_cv$G_opt,
lambda_l = mod_cv$lambda_l_opt,lambda_s =mod_cv$lambda_s_opt,maxit = 50,q=30)

The cluster membership vector and the plot of the estimated cluster mean functions and the classified curves are obtained as follows.

print(mod$clus$classes)
plot(mod)

References

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.