# \donttest{
sim <- simulate_fmi_data(N_A = 20, N_B = 20, T_points = 30)
g_fac <- factor(sim$group_vec)
idx_A <- which(g_fac == levels(g_fac)[1])
# --- FIX: Added Y= and argvals= ---
fpca_A <- refund::fpca.sc(Y = sim$Y_mat[idx_A, ], argvals = sim$argvals, npc = 3)
fpca_B <- refund::fpca.sc(Y = sim$Y_mat[-idx_A, ], argvals = sim$argvals, npc = 3)
build_scree_plot(fpca_A, fpca_B, levels(g_fac))
# }
Run the code above in your browser using DataLab