Nothing
#' @export
plot.hfpca <- function(x, ...){
dims <- dim(x$gam_pca)
num.plot <- ceiling(dims[3]/3)
TT <- dim(x$gam_pca)[2]
time <- seq(0,1,len=TT)
cnt <- 1
for (ii in 1:num.plot){
graphics::layout(matrix(c(1,2,3), 1, 3, byrow = TRUE))
graphics::matplot(time,t(x$gam_pca[,,cnt]),type="l")
graphics::title(main=sprintf("gam: PD %d", cnt))
if (dims[3] >= cnt + 1){
graphics::matplot(time,t(x$gam_pca[,,cnt+1]),type="l")
graphics::title(main=sprintf("gam: PD %d", cnt+1))
} else {
graphics::plot.new()
}
if (dims[3] >= cnt + 2){
graphics::matplot(time,t(x$gam_pca[,,cnt+2]),type="l")
graphics::title(main=sprintf("gam: PD %d", cnt+2))
}else {
graphics::plot.new()
}
cnt <- cnt + 3
}
graphics::layout(1)
cumm_coef = 100*cumsum(x$latent)/sum(x$latent)
plot(cumm_coef,type="l",col="blue",main="Coefficient Cumulative Percentage", ylab = "Percentage")
}
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.