stata按组别产生随机数

根据样本（var1）、时间的不同（var2），生成符合二项分布（次数n，概率p）的随机数（var3），观察个数设置为1000，种子数1234。每个样本、时间，生成新变量（var4），新变量var4=mean（var3/n-p）。举例：

var1	var2	n	p	var4
1	2018-6	123	.7
2	2017-12	76	.6
3	2017-6	34	.3
4	2016-12	22	.5
5	2016-6	10	.3

谢谢各位！