在LSＶ模型中，
                                     H(i,t) = | P(i,t) - P(t) | - AF(i,t)
        其中，p(i,t)=B(i,t)/(B(i,t)+S(i,t))，即表示在时期t内买入股票i的基金占买卖股票i的所有基金的比例，p(t)为时期t内预期买入股票i的基金占买卖股票i的所有基金的比例，调整因子AF(i,t)=E[|p(i,t)-p(t)|]，它是在没有羊群行为存在的零假设条件下的|B(i,t)/(B(i,t)+S(i,t))-p(t)|的期望值。因为B(i,t)是服从参数为p(t)的二项分布，所以调整因子AF如下:

                       
问题：
这个涉及循环，每一个样本都要循环，还要嵌套，搞了半天也不会，只好求助大家了，用stata怎么编程呢？

给出的数据变量分别是季度tq、股票代码stk、买入B、卖出S、总交易T、均值P（t）
----------------------- copy starting from the next line -----------------------

------------------ copy up to and including the previous line ------------------


我的命令这样的是错的，
[/CODE]
bys tq:n=_n
forval i=1/n{
forval k=0/T{
local sum=0
local P5=`k'/T
local c=comb(T,k)
local a=(`p5'-m1)*`c'*m1^k*(1-m1)^(T-k)
local sum=`sum'+`a'
}
}

[/CODE]



谢谢回答！