x <- list(E1 <- list(A1 <- c(0.14, 0.22, 0.28),
A2 <- c(130000, 150000, 200000),
X <- c(350000, 400000, 525000)),
E2 <- list(A1 <- c(0.2, 0.3, 0.4),
A2 <- c(165000, 205000, 250000),
X <- c(550000, 600000, 650000)),
E3 <- list(A1 <- c(0.2, 0.4, 0.52),
A2 <- c(200000, 400000, 500000),
X <- c(625000, 700000, 800000)))
probs <- c(0.1, 0.5, 0.9)
true.seed <- c(0.27, 210000)
fit <- expert(x, "cooke", probs, true.seed, 0.03)
Fn <- cdf(fit)
Fn
knots(Fn) # the group boundaries
Fn(knots(Fn)) # true values of the cdf
plot(Fn) # graphic
Run the code above in your browser using DataLab