```{r} #Quantiles

qexp(0.99,0.1)

#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)

```