在Stata在金融场景里使用的两个例子里，我们讨论了如何下载美联储加息前后两天的美国国债利率期限结构的数据。自从2017年以来，美国总共加过7次息， 2017年3月16， 6月15，12月14；2018年3月20，6月12，9月25和12月18日。

当时有位读者留言，说Stata因为没有列表格式，所以循环起来有点麻烦。最近想到了另外一个办法来解决这个问题。而这个办法的核心思想就是从定义的暂元集里取出子元素，这和列表(list)的思维有异曲同工之妙。

最初看到这个例子是在Stata官网，他们需要做一个平行(parallel)的循环，即从暂元集A取出的元素和暂元集B取出的元素要并列地展示在一起。下面这个例子：

local agrp "cat dog cow pig"local bgrp "meow woof moo oinkoink "local n : word count `agrp'local m: word count `bgrp'forvalues i = 1/`n' {     local a : word `i' of `agrp'     local b : word `i' of `bgrp'     di "`a' says `b'"}

我们得到的结果是：

cat says meowdog says woofcow says moopig says oinkoink

如果我们掌握了这种从暂元集里取元素的循环方式，对于上面寻找加息前后两天数据的例子，我们可以使用如下的程序：

clearfreduse DGS30 DGS20 DGS10 DGS7 DGS5 DGS3 DGS2 DGS1 DGS6MO DGS3MO DGS1MOsave fred_data,replace

clearinput year month day2017 3 162017 6 152017 12 142018 3 202018 6 122018 9 252018 12 18 end

gen daten=mdy(month,day,year)format daten %tdkeep datenlevelsof daten,local(fdate)save daten,replace

use fred_data,clearkeep if year(daten)>=2017

levelsof daten,local(fulldate)

local n : word count `fdate'local m: word count `fulldate'tempname memholdpostfile `memhold' day using date.dta,replace

forvalue i=1/`n'{forvalue j=1/`m'{local a : word `i' of `fdate'local b : word `j' of `fulldate'

if `a'==`b'{local b1 : word `=`j'-2' of `fulldate'post  `memhold' (`b1') local b2 : word `=`j'-1' of `fulldate'post  `memhold' (`b2') local b3 : word `=`j'+1' of `fulldate'post  `memhold' (`b3') local b4 : word `=`j'+2' of `fulldate'post  `memhold' (`b4') }}}postclose `memhold'

use date,clearformat day %tdsort day

感兴趣的读者可以动手操作一下，以充分理解这种从暂元里取子元素的方式。

--

我们受大家欢迎的文章还包括：

Stata画图难看？matplotlib来帮忙！

Python和Stata的对话

Stata和阿尔法狗的距离，只隔着一个TensorFlow

如何在Stata里下载百度指数？

新学期特辑 | MIT华人教授的金融学第一课

AER论文复制 | 货币政策、金融危机与信贷繁荣的破灭

如何找到顶尖期刊论文可供复制的数据和程序

使用Jupyter Notebook提高实证研究的可复制性

如果想从Stata转到Python，也许你可以这样做

Stata应用 | 中金的研究报告一定就可信吗？

周末特辑 | 更适合金融背景学生的Stata学习案例

Econometrica论文复制 | 历史位置、集体行为和经济发展：越南的经验