200论坛币求助一个小程序

如题， 我想用stata计算如下公式：



其中bin是二项分布~ min里面的意思是取其中最小值，例如n小于等于1则取值为1 否则取值为 2/n+1
PS: n小于等于30

多谢！！！

local sum0=0
forvalues x=0(1)7 {
   local sum0 = `sum0'+binomialp(7,`x',0.5)*min(2/(`x'+1),1)
}
di 22.5-20*`sum0'

loc n=30  /*总试验次数*/
loc s=0
forv i=0/`n'{
loc s=`s'+binomialp(`n',`i',.5)*min(2/(`i'+1),1)
}
di 22.5-20*`s'

*这里楼主需要说明所提出的Bin(k;n;p)的定义。（本悬赏应给dxystata）

很好，我也学习了
谢谢

dxystata 发表于 2012-3-21 23:57
local sum0=0
forvalues x=0(1)7 {
   local sum0 = `sum0'+binomialp(7,`x',0.5)*min(2/(`x'+1),1)

另，可否解释一下local的用法，我对这个不太明白~ 多谢~

local 定义暂元，类似SAS中局部宏变量！
用local定义后 引用的话需加·‘