求助关于数据处理的

id year close_ind _freq
33 2002 8.6142308 234
33 2003 6.6412448 241
33 2004 5.021358 243
33 2005 3.3047863 234
33 2006 3.2140332 181
33 2007 6.8583966 237
33 2008 4.5016327 245
33 2009 5.308983 236
33 2010 6.77075 240
33 2011 5.5315385 26

63 1997 22.609063 32
63 1998 30.843699 246
63 1999 32.157197 239
63 2000 34.388025 238
63 2001 31.044914 232
63 2002 20.079785 233
63 2003 18.002647 238
63 2004 24.637866 239
63 2005 26.674387 212
63 2006 30.023718 234
63 2007 50.723602 236
63 2008 47.497366 243
63 2009 34.937303 241
63 2010 32.367437 238
63 2011 26.131539 26
有以上的金融数据，如果我要bysort id,通过id分组来处理，比如：id为33的，产生一个变量var，2004年var的值等于[2004年的close_ind -[（2003年的close_ind*2003年的_freq +2002年的close_ind*2002年的_freq ）/（2002年的_freq +2003年的_freq ）]]/（2003年的close_ind*2003年的_freq +2002年的close_ind*2002年的_freq ），以此类推，2005年的var值=2005年值与前两年（03，04年）平均值的差再除以前两年（03，04年）平均值。
该如何实现。
请各位热心的朋友不吝赐教，在下不胜感激。

bys id (year): g var=       ///
(close_ind-(close_ind[_n-1]*_freq[_n-1]+close_ind[_n-2]*_freq[_n-2])/(_freq[_n-1]+_freq[_n-2]))       ///
/(close_ind[_n-1]*_freq[_n-1]+close_ind[_n-2]*_freq[_n-2])

一直在等着答案，谢谢sungmoo