目前，尽管试图进行绘制，但试图编写函数以计算保证金收益和准备金/保单价值，但仍会出现相同的错误，并且对我的解决方法颇为迷惑。任何帮助将不胜感激。
#Annuity-due function
AnU = function(x, n, i)
{
An = 0
for (k in 0:(n-1))
{
An = An + (1/(1+i))^k*lx[x+k]/lx[x]
}
return(An)
}
#Term/whole life assurance function
TWa <- function(x, n, i){
if (missing(n))
n = 120 - x - 1
Ta = 0
for(k in 0:(n-1))
{
Ta = Ta + (1/(1+i))^(k+1)*(lx[x+k]-lx[x+k+1])/lx[x]
}
return(Ta)
}
#Setting Premium
P <- 100000*TWa(x = 25, n = 40, i = 0.04)/AnU(x = 25, n = 40, i = 0.04)
#Policy value function
age <- c(1:120)
term <- c(1:120)
time <- c(1:120)
policy_TW <- function(time, age, term, rate)
return(100000*TWa(x = age + time, n = term - time, i = rate) -
P * AnU(x = age + time, n = term - time, i = rate))
for (time in 0:40)
cat("At time", time, " policy value is ", policy_TW(time, 25, 40, 0.04), "\n")
#Plot
plot(0:40, policy_TW(0:40, 25, 40, 0.04), ylim = c(0,6000), type = "b")
> plot(0:40, policy_TW(0:40, 25, 40, 0.04), ylim = c(0,6000), type = "b")
Warning messages:
1: In 0:(n - 1) :
numerical expression has 41 elements: only the first used
2: In 0:(n - 1) :
numerical expression has 41 elements: only the first used
编辑：为了重现性，我忘记添加lx代码。这是从链接的xls文件计算出来的，该文件随后被命名为data https://www.actuaries.org.uk/documents/am92-permanent-assurances-males
# Calculate the values for lx
qx = c(rep(NA,16),data$qx)
lx = vector()
lx[1:17] = c(rep(NA,16),10000)
for (i in 18:120)
{
lx[i]=lx[i-1]*(1-qx[i-1])
}
该语法0:(n-1)创建一个序列，即一个从零到的整数整数n-1。但是，您n本身就是一个向量，因为它被定义为term - time并且time是所有0到40的整数的向量（因此它包含41个元素）。因此，仅使用的第一个元素n来创建序列。例如，如果n是c(4,5,10)，则仅使用4，而忽略5和10，然后0:(n-1)返回c(0,1,2,3)。
您将需要重新定义AnU()和TWa()功能。我不知道正是你需要什么，但我想你只需更换线路for (k in 0:(n-1))用for (k in 0:(max(n)-1))。我想测试一下，但由于lx您的代码中未定义对象，而这两个函数的公式中都没有定义，因此我无法重现您的错误。
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句