data("abide_9002_23")
res <- lapply(list(1:100, 101:200, 201:261), function(part) {
set.seed(123) # When using partial_forder, all parts must use the same seed.
fset <- frank.flm(nsim=99, formula.full = Y ~ Group + Sex + Age,
formula.reduced = Y ~ Group + Sex,
curve_sets = list(Y = abide_9002_23$curve_set[part,]),
factors = abide_9002_23$factors, savefuns = "return")
partial_forder(fset, measure="erl")
})
combine_forder(res)
Run the code above in your browser using DataLab