Quantiles Built-in functions
qexp(0.99,0.1)
[1] 46.0517
Quantiles with optim built-in function 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

LS0tCnRpdGxlOiAiQUNUVS00NTcgKE1hc3QgNzI0KSBSaXNrIFRoZW9yeSIKYXV0aG9yOiAiUHJvZi4gTWVsaW5hIE1haWxob3QiCmRhdGU6ICIyMDI0LTEyLTI4IgpvdXRwdXQ6CiAgaHRtbF9ub3RlYm9vazogZGVmYXVsdAogIHBkZl9kb2N1bWVudDogZGVmYXVsdAotLS0KCgpRdWFudGlsZXMgQnVpbHQtaW4gZnVuY3Rpb25zCmBgYHtyfQpxZXhwKDAuOTksMC4xKQpgYGAKClF1YW50aWxlcyB3aXRoIG9wdGltIGJ1aWx0LWluIGZ1bmN0aW9uCkRpc2NyZXRlIFNpbXBsaWZpZWQgQ2FzZQoKYGBge3J9Cm49MTAKWSA8LSBjKDEsMiwzLDQsNSw2LDcsOCw5LDEwKQpGaW52IDwtIGZ1bmN0aW9uKGFhKXsKCXNjb3JlIDwtIDAKCWZvcihjYyBpbiAxOm4pewoJCW9iaiA8LSBmdW5jdGlvbihjYyl7CgkJCXZhbCA8LSAwCgkJCWZvcihpIGluIFkpewoJCQkJdmFsW2ldIDwtIHtpZihpPGNjKXsoYWEtMSkqKGktY2MpfSBlbHNleyhhYSkqKGktY2MpfQoJCQkJfQoJCQl9CgkJcmV0dXJuKHN1bSh2YWwpKX0KCXNjb3JlW2NjXSA8LSBvYmooY2MpCgl9CnNjb3JlCnBsb3Qoc2NvcmUpCglRdWFudGlsZSA8LSBtYXRjaChtaW4oYWJzKHNjb3JlKSksc2NvcmUpCglyZXR1cm4oUXVhbnRpbGUpCn0KRmludigwLjY1KSAKYGBgCg==