如何得到样本个体得分及排名？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请问：如何得到样本个体得分及排名？对应的stata命令是什么？是在对面板数据做主成分分析中遇到此问题，希望得到解答。谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：stata命令 Stata 主成分分析 面板数据 tata 如何 样本

这样的问法，很难让人了解问题所在！具体一点（有程序、例子最好）！

黃河泉 发表于 2016-10-23 11:15
这样的问法，很难让人了解问题所在！具体一点（有程序、例子最好）！

您好，目前在做关于七个省市的区域创新能力测算，使用其1998-2013年间的13项指标，已经进行了主成分分析，得到了三项公因子。下一步想要测算每个省份每年（1998-2013）的区域创新能力得分，如何实现呢？或者是只处理2013年的数据，得出每个省份在各个公因子上的得分及总得分，用以比较。希望得到您的解答。谢谢

赟瞾 发表于 2016-10-23 11:39
您好，目前在做关于七个省市的区域创新能力测算，使用其1998-2013年间的13项指标，已经进行了主成分分析， ...

这样清楚多了，我有一个类似的程式（你对应改一下）：

1. // Principal Component
   
2. // define a wrap program
   
3. capture program drop pca_new
   
4. program pca_new, eclass
   
5. pca `0'
   
6. tempname L
   
7. matrix `L' = e(L)
   
8. matrix b = `L'[1..., 1]'
   
9. matlist b
   
10. ereturn post b
    
11. end
    
12. 
    
13. // run -pca_new-
    
14. statsby _b, by(ID) saving("data\3-mycomp.dta", replace): pca_new ncskew duvol
    
15. 
    
16. // -merge- the result to the data set in memory
    
17. merge m:1 ID using "data\3-mycomp.dta", nogenerate
    
18. 
    
19. // standardize the variables
    
20. foreach x of varlist ncskew duvol {
    
21.   bysort ID: egen mean_`x'= mean(`x')
    
22.   bysort ID: egen sd_`x' = sd(`x')
    
23.   gen std_`x' = (`x'-mean_`x')/sd_`x'
    
24.   drop mean_`x' sd_`x'
    
25. }
    
26. 
    
27. // apply the formula
    
28. gen pc = std_ncskew*_b_ncskew + std_duvol*_b_duvol
    
29. 
    
30. drop std_* _b_*

复制代码

请做个参考！

赟瞾 发表于 2016-10-23 11:39
您好，目前在做关于七个省市的区域创新能力测算，使用其1998-2013年间的13项指标，已经进行了主成分分析， ...

或许你也可看看 help pca postestimation 的 predict 功能（比较简单，但我还没时间去仔细研究！）。

黃河泉 发表于 2016-10-23 12:07
这样清楚多了，我有一个类似的程式（你对应改一下）：请做个参考！

非常感谢您的解答，我会好好学习下的