有一套调查数据，样本为个人。变量中已有年龄、预期每年的养老金收入，也估计出各年龄死亡概率曲线，如果公式如下，应该怎么编程序计算每个人生命周期的预期养老金收入现值综合？
公式：

S（下标a，n）=1-exp(-6.5749+0 .1095487*n-1.146528*ln(n+0.4)
其中，p代表每年养老金收入，已有每个人的收入变量；S（下标a，n）代表a岁的人活到n岁的概率；g为养老金增长率，以固定值0.1计算，d为贴现率，以固定值0.01计算。
假定每个样本编号为ID。

希望有过类似处理经验的朋友帮忙！谢谢！