Description Usage Arguments Value Examples
After partitioning using PULS, this function can plot the functional waves and color different clusters as well as their medoids.
1 2 3 4 5 6 7 8 9 10 |
toclust.fd |
A functional data object (i.e., having class |
intervals |
A data set (or matrix) with rows are intervals and columns are the beginning and ending indexes of of the interval. |
puls.obj |
A |
xlab |
Labels for x-axis. If not provided, the labels stored in |
ylab |
Labels for y-axis. If not provided, the labels stored in |
lwd |
Linewidth of normal waves. |
alpha |
Transparency of normal waves. |
lwd.med |
Linewidth of medoid waves. |
A ggplot2 object.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 | library(fda)
# Build a simple fd object from already smoothed smoothed_arctic
data(smoothed_arctic)
NBASIS <- 300
NORDER <- 4
y <- t(as.matrix(smoothed_arctic[, -1]))
splinebasis <- create.bspline.basis(rangeval = c(1, 365),
nbasis = NBASIS,
norder = NORDER)
fdParobj <- fdPar(fdobj = splinebasis,
Lfdobj = 2,
# No need for any more smoothing
lambda = .000001)
yfd <- smooth.basis(argvals = 1:365, y = y, fdParobj = fdParobj)
Jan <- c(1, 31); Feb <- c(31, 59); Mar <- c(59, 90)
Apr <- c(90, 120); May <- c(120, 151); Jun <- c(151, 181)
Jul <- c(181, 212); Aug <- c(212, 243); Sep <- c(243, 273)
Oct <- c(273, 304); Nov <- c(304, 334); Dec <- c(334, 365)
intervals <-
rbind(Jan, Feb, Mar, Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec)
PULS4_pam <- PULS(toclust.fd = yfd$fd, intervals = intervals,
nclusters = 4, method = "pam")
ggwave(toclust.fd = yfd$fd, intervals = intervals, puls = PULS4_pam)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.