# NOT RUN {
set.seed(666L)
dat <- data.frame(
group = gl(2, 5),
y = c(2*rlogis(5L), 10 + 2*rlogis(5L))
)
gfi <- gfilinreg(
y ~ 0 + group, distr = "logistic", data = dat, L = 25L, nthreads = 2L
)
fcdf <- gfiCDF(~ group1 - group2, gfi)
fcdf(0)
plot(fcdf)
# }
Run the code above in your browser using DataLab