set.seed(1)
sim <- simulation_DATA.beta(n = 120, p = 5, s = 2)
y_low <- pmax(sim$Y - 0.05, 0)
y_high <- pmin(sim$Y + 0.05, 1)
interval_fit <- sb_beta_interval(
sim$X,
Y_low = y_low,
Y_high = y_high,
B = 5,
step.num = 0.4
)
attr(interval_fit, "interval")
Run the code above in your browser using DataLab