我在做一个t分布下的rolling GARCH模型，来计算VaR
因此需要自由度来查t分布表来确定分位数
再用标准差，分位数来计算VaR
问题来了，rolling GARCH产生了1000个GARCH equation，每个equation都一个自由度
所以我写了这个代码来创建一个序列dof， 收集每个GARCH equation的自由度（j从1到1000）
series dof=eq_{!j}.@df
再用
quan=@tdist（0.01,dof）
计算出1000个GARCH equation对应的1000个分位数

最后计算VaR

但是问题我卡在了series dof=eq_{!j}.@df这步上
dof里面全是995，而实际的自由度明显不是这个
是不是这个代码有问题？
求大神帮忙啊