# library(LSTMfactors)
# set.seed(358)
#
# data(data.DAPCS)
# response <- as.matrix(data.DAPCS[, 3:22]) ## loading data
#
# PA.obj <- EFAfactors::PA(response = response, plot = FALSE)
# LSTM.obj <- LSTM(response = response, plot = FALSE)
#
#
# # CDF.obj <- EFAfactors::CDF(response = response,
# #                            nfact.max = 10,
# #                            N.pop = 10000, N.Samples = 5000, mtry = "sqrt")
#
# plot(PA.obj)
# plot(LSTM.obj)
# # plot(CDF.obj)