Nothing
plot.fpcat <-
function(x, nscore=c(1, 2), main="PCA of probability density functions", sub.title=NULL, ...)
{
if (length(nscore) > 2) {
warning(paste("Since dad-4, nscore must be a length 2 numeric vector. The scores number", nscore[1], "and", nscore[2], "are plotted."))
nscore <- nscore[1:2]
}
times <- x$times
numtimes <- as.numeric(times)
inertia=x$inertia$inertia
coor=x$scores[-1]
ngroup <- nrow(coor)
if (max(nscore)>ncol(coor))
stop("The components of nscore must be smaller than the number of score columns in the x$scores data frame")
group=x$scores[, 1]
# ind=combn(nscore, 2)
if (length(nscore) > 1) {
i1 = nscore[1]; i2 = nscore[2]
coorij <- data.frame(coor[-nrow(coor), c(i1, i2)],
coor[-1, c(i1, i2)], stringsAsFactors = TRUE)
names(coorij) <- c("x1", "y1", "x2", "y2")
graph <- ggplot(coorij)
graph <- graph + aes_q(x = as.name("x1"), y = as.name("y1"),
xend = as.name("x2"), yend = as.name("y2"))
graph <- graph + geom_segment(arrow = arrow(unit(0.1, "inches")))
graph <- graph + labs(title = main, subtitle = sub.title,
x = paste0(names(coor)[i1], " (", inertia[i1], "%)"),
y = paste0(names(coor)[i2], " (", inertia[i2], "%)"))
print(graph)
}
if (length(nscore) == 1) {
if ("Date" %in% class(times))
times <- as.POSIXct(times)
graph <- ggplot(data.frame(time = times, coor = coor[, nscore]), stringsAsFactors = TRUE)
graph <- graph + aes_q(x = as.name("time"), y = as.name("coor"))
graph <- graph + geom_line()
graph <- graph + labs(title = main, subtitle = sub.title,
x = "time", y = paste0(names(coor)[nscore], " (", inertia[nscore], "%)"))
print(graph)
}
return(invisible(NULL))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.