里面的p我是直接按[x,j]调用矩阵里面的数据的,但是写循环就会出问题。。。我写了一个。。。但是怎么都出问题。。。想问问为什么?怎样改正呢?
I<-function(x,y,z,r1,r2,A1,A2,A3){
C2<-matrix()
for(m in 1:(55-j-x)){
for(j in (21-z):105){
C2[m,j]=(c[z,j]-c[z,j+1])*(b[y,m+j]-b[y,m+j+1])*a[x,j]*A2/((1+r2)^(j+m+1))
}
}
sum(C2)}
I(30,30,5,0.025,0.035,120000,80000,200000)
a、b、c都是生命表矩阵


雷达卡







京公网安备 11010802022788号







