require(tidyr)
eem_list %>%
eem_extract(eem_names((eem_list))[1], keep = TRUE) %>%
ggeem() +
geom_eempeakloc()
# We plot only component 3 of the PARAFAC model to avoid letter overload
pf4[[1]] %>%
(\(x, comp = 3) { #please change comp to the index of the comp. to plot.
x$A <- x$A[, comp, drop = FALSE]
x$B <- x$B[, comp, drop = FALSE]
x$C <- x$C[, comp, drop = FALSE]
x
})() %>%
ggeem() +
geom_eempeakloc() +
coord_cartesian(xlim = c(250,455), ylim = c(290,578))
Run the code above in your browser using DataLab