该公式中研究对象的指标ab、bc、cd为已知变量，ref1、ref2、ref3是随着age(0到90岁)的变化而变化。
假设第一研究对象age1=20。则Y等于上述公式从0至age1（20）的累加
第二个观测对象age1=30, 则Y等于上述公式从0至30的累加。
不知道这样描述问题是不是更合适，这个问题已经困扰我很久了，但是编程0基础，始终没有得到解决，我运行了以下命令只能得到每个观测对象的0岁的y，而不能累加：
gen Y=.
local N=`age1'+1
forvalues i =1(1)`N' {
           replace Y=((ab-0.4)*ref1[`N]'+(bc-0.19)*ref2[`N']+(cd-2.0.045)*ref3[`N'])/1000*366-0.79
               }


显然这个程序是不对的，但是真不知道如何实现这个公式，望对于stata编程比较精通的精英们能解此惑，或给点建议，多谢了！