感谢Sungmoo版主：将程序变成ado命令（问题完美解决）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请路过的老师帮个忙：

              在下经常利用如下代码处理数据，因此想将这段代码变成ado文件。
          命令格式如stata命令格式  或其它格式

        gsca  m  f  r  varlist  [if] [in]

我用的代码如下（示例数据在附件中）：

foreach vi in dans blzh rzh lch lh water {
egen `vi'std=std(`vi'), mean(0) std(1)
bys m f (r):egen `vi'e=mean(`vi'std)
bys m (f):egen `vi'mg=mean(`vi'e)
bys f (m):egen `vi'fg=mean(`vi'e)
gen `vi'sca=`vi'e-`vi'mg-`vi'fg
drop `vi' `vi'std
}
keep if r==1
drop r

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享1 收藏5 回帖

关键词：sungmoo Sung GMO Sun stata命令 程序

jzhyue 发表于 2009-12-10 15:20
请路过的老师帮个忙：

              在下经常利用如下代码处理数据，因此想将这段代码变成ado文件。
          命令格式如stata命令格式  或其它格式

        gsca  m  f  r  varlist  [if] [in]

我用的代码如下（示例数据在附件中）：

foreach vi in dans blzh rzh lch lh water {
egen `vi'std=std(`vi'), mean(0) std(1)
bys m f (r):egen `vi'e=mean(`vi'std)
bys m (f):egen `vi'mg=mean(`vi'e)
bys f (m):egen `vi'fg=mean(`vi'e)
gen `vi'sca=`vi'e-`vi'mg-`vi'fg
drop `vi' `vi'std
}
keep if r==1
drop r

试试把下面的代码保存为gsca.ado文件：

capture program drop gsca
program gsca, sortpreserve
version 10.0
args m f r
syntax varlist(min=1) [if] [in]
marksample touse
foreach vi in dans blzh rzh lch lh water {
egen `vi'std=std(`vi'), mean(0) std(1)
bys `m' `f' (r):egen `vi'e=mean(`vi'std)
bys `m' (f):egen `vi'mg=mean(`vi'e)
bys `f' (m):egen `vi'fg=mean(`vi'e)
gen `vi'sca=`vi'e-`vi'mg-`vi'fg
drop `vi' `vi'std
}
keep if `r'==1
drop `r'
end

lizhiyi 发表于 2009-12-10 16:26 foreach vi in dans blzh rzh lch lh water {

我猜作者这一步想使用变量通配符（对m f r以外的变量操作）。

楼主最好把自己想达到的结果详细说明一下。

谢谢！可能我没有把要求表达清楚。


关键每一次foreach vi in后面的varlist是不一样的，即：数据结构相似，数据和变量不一样（m  f  r 结构一样，而后面的varlist 结构相同，但变量数是可多可少可不同的）
比如我想运行
gsca   m   f  r   x1 x2 x3 x4

lizhiyi 发表于 2009-12-10 16:26
试试把下面的代码保存为gsca.ado文件：

capture program drop gsca
program gsca, sortpreserve
version 10.0
args m f r
syntax varlist(min=1) [if] [in]
marksample touse
foreach vi in dans blzh rzh lch lh water {
egen `vi'std=std(`vi'), mean(0) std(1)
bys `m' f (r):egen `vi'e=mean(`vi'std)
bys `m' (f):egen `vi'mg=mean(`vi'e)
bys `f' (m):egen `vi'fg=mean(`vi'e)
gen `vi'sca=`vi'e-`vi'mg-`vi'fg
drop `vi' `vi'std
}
keep if `r'==1
drop `r'
end

sungmoo 发表于 2009-12-10 16:33

lizhiyi 发表于 2009-12-10 16:26 foreach vi in dans blzh rzh lch lh water {

我猜作者这一步想使用变量通配符（对m f r以外的变量操作）。

同意sungmoo表达的意思

jzhyue 发表于 2009-12-10 16:46 m  f  r 结构一样，而后面的varlist 结构相同，但变量数是可多可少可不同的

m f r可以按一定规则生成吧？

也就是说，m、f、r的生成，可以是程序的一部分吧？

sungmoo 发表于 2009-12-10 17:06

jzhyue 发表于 2009-12-10 16:46 m  f  r 结构一样，而后面的varlist 结构相同，但变量数是可多可少可不同的

m f r可以按一定规则生成吧？

也就是说，m、f、r的生成，可以是程序的一部分吧？

m、f、r是试验设计的参数，m、f代表不同的因子，各有不同的水平，r代表重复数，构成m*f*r的数据结构，以后的数据为试验调查数据。
m、f、r是事先设计好的，不能由程序生成固定的。
（m、f、r不能有缺失，而调查数据可以有缺失）

cap pr drop gsca
pr gsca
syntax varlist(min=1) [if] [in]
foreach vi in `varlist' {
egen `vi'std=std(`vi') `if'
bys m f (r):egen `vi'e=mean(`vi'std) `if'
bys m (f):egen `vi'mg=mean(`vi'e) `if'
bys f (m):egen `vi'fg=mean(`vi'e) `if'
g `vi'sca=`vi'e-`vi'mg-`vi'fg `if'
drop `vi' `vi'std
}
keep if r==1
drop r
end

9# sungmoo
谢谢
还有个问题：就示例数据而言，这个ado，很好。一旦利用其它类似数据，对应于m、f、r三个变量可能名字就不一样，就需要改变变量的名字（或修改程序中的m、f、r)，也就是说现在这个ado缺乏通用性（当然在运行之前我用三个rename就可以搞定）
我想是否可以用个什么命令（比如ds或lookfor或其它）将前三个变量名字提出来，让这个ado具备通用性？这也是我发求助的初衷。

或者：事先将变量顺序排好，我们引用第一个变量、第二个变量
能否用
ds
loca m=word(r(varlist),1)
loca  f=word(r(varlist),2)
我知道regress y x1-x9
就有这个能力将两类变量分开，不知是通过什么途径？