在使用creditmetrics软件包时提示下标出界，请高手指教
N<-4
n<-100000
r<-0.025
ead<-c(320,1300,150,34)
rc<-c(“AAA”,”AA”,”A”,”BBB”,”BB”,”B”,”CCC”,”D”)
lgd<-0.3873
rating<-c(“AAA”,”AA”,”BB”,”C”)
loadnames<-c(“Load1”,”Load2”,”Load3”,”Load4”)
alpha<-0.99
rho<-matrix(c(1,0.625,-0.7619,-0.7202,0.625,1,-0.7679,-0.7857,-0.7619,-0.7679,1,0.5893,-0.7202,-0.7857,0.5893,1),4,4,dimnames=list(loadnames,loadnames),byrow=TRUE)
M<-matrix(c(97.12,2.58,0.18,0.02,0.03,0.01,0.06,0,2.85,89.73,6.68,0.49,0.05,0.08,0.06,0.06,0.1,2.16,91.43,5.42,0.4,0.17,0.07,0.25,0.03,0.14,3.66,89.63,4.47,0.73,0.52,0.82,0.08,0.05,0.2,6.3,81.79,6.59,2.29,2.7,0,0.03,0.16,0.27,3.94,70.43,19.21,5.96,0,0,0.07,0.1,0.29,3.73,67.83,27.98,0,0,0,0,0,0,0,100)/100,8,8,dimnames=list(rc,rc),byrow=TRUE)cm.CVaR(M,lgd,ead,N,n,r,rho,alpha,rating)
错误: 下标出界
此外: 警告信息：
In qnorm(p, mean, sd, lower.tail, log.p) : 产生了NaNs