#Quantiles
qexp(0.99,0.1)
## [1] 46.0517
#Quantiles with optim
#Discrete Simplified Case
n=10
Y <- c(1,2,3,4,5,6,7,8,9,10)
Finv <- function(aa){
score <- 0
for(cc in 1:n){
obj <- function(cc){
val <- 0
for(i in Y){
val[i] <- {if(i<cc){(aa-1)*(i-cc)} else{(aa)*(i-cc)}
}
}
return(sum(val))}
score[cc] <- obj(cc)
}
score
plot(score)
Quantile <- match(min(abs(score)),score)
return(Quantile)
}
Finv(0.65)
## [1] 7