y <- c(
0, 0, 0, 1,
0, 1, 1, 1
)
group <- gl(2, 4)
fidsamples <- gfilogisreg(y ~ 0 + group, N = 500) # (N=500 is not serious)
fcdf <- gfiCDF(~ exp(group1) / exp(group2), fidsamples)
fcdf(1)
plot(fcdf)
Run the code above in your browser using DataLab