原QI1函数不能计算形如QI1(c(0.5,0.6,0.7))的形式，而在QI2函数会用到的QI1(x)，x为向量的计算形式，所以程序改为，
QI0<-function(a) {
qnorm(a)
}
#1阶上分位数积分
QI1<-function(a){
n = length(a)
result = 1:n
for (i in 1:n){
temp = integrate(QI0,lower=0, upper=a)
result = temp$value
}
result
}
QI2<-function(a){
result = integrate(QI1,lower=0, upper=a)
rusult = result$value
}
QI2(0.5)
扫码加好友，拉您进群
全部版块
我的主页
收藏
