第一步中，令c等于几，很关键吗？

这里，令c=某个常数，主要的作用是，把原数据库整体当作一个单位来复制。

而生成的n，标识了每个观测值按前述复制规则复制的次数。

对n排序后，即可得到原数据库保持原样被依次复制后的结构。

多谢sungmoo版主的讲解，使我真正理解了这个函数，但还有一个问题，就是在没有对n排序前，我发现n的出现并无规律可循，例如

   +--------------+
     | year   c   n |
     |--------------|
  1. | 2001   1   2 |
  2. | 2001   1   1 |
  3. | 2001   1   3 |
  4. | 2002   1   3 |
  5. | 2002   1   1 |
     |--------------|
  6. | 2002   1   2 |
  7. | 2003   1   2 |
  8. | 2003   1   3 |
  9. | 2003   1   1 |
 10. | 2004   1   3 |
     |--------------|
 11. | 2004   1   2 |
 12. | 2004   1   1 |
 13. | 2005   1   3 |
 14. | 2005   1   1 |
 15. | 2005   1   2 |
     |--------------|
 16. | 2006   1   1 |
 17. | 2006   1   2 |
 18. | 2006   1   3 |
 19. | 2007   1   3 |
 20. | 2007   1   1 |
     |--------------|
 21. | 2007   1   2 |
 22. | 2008   1   2 |
 23. | 2008   1   1 |
 24. | 2008   1   3 |
     +--------------+

请问，n的赋值规则是什么？

这个好像是随机的，我也不太清楚具体规则。（所以，后面要补上sort n）

另一个比较直观的土办法：

//生成模拟数据
clear
set obs 10
gen id = _n
gen rnd = uniform()

//数据倍增nfold倍
local N = _N
forval i = 1/ (nfold-1) {
     expand 2 in 1/`N'
}