person	mill
922	3
8675	9
10634	16
10640	16
10642	16
...

以上为原始数据，反映某人在某工厂工作。
需要将其转为如下形式，用于做网络分析。

	3	9	16	...
922	1	0	0	...
8675	0	1	0	...
10634	0	0	1	...
10640	0	0	1	...
10642	0	0	1	...
...	...	...	...	...

第一行是工厂名单，第一列是人名单
数据为哑变量，1表示该人在该工厂工作，0表示不在

请同学们和各方高人帮助：

（1）如何自动在STATA中生成如下的空白矩阵？如果手动操作，显然是可以做的，但是因为同形式的数据表非常多，需要做大量转换，手动去做不现实。

	3	9	16	...
922
8675
10634
10640
10642
...

（2）应如何写语句，根据哑变量的原则，自动填写上表？

我自己弄了好几天，拼拼凑凑，用excel+stata+spss可以实现。但是非常希望能把所有语句都写在STATA中，因为数据太多，只有这样才能随时计算、更新。
真诚求教各位同学和各方高人，谢谢了！