R语言计算贸易公式

各位大神，初学R语言，求助编写程序，计算RCA指标，不胜感激！数据结构见图片。[attachimg]2232654[/attachimg
第一栏的year代表年份，从1995年到2013年；ISO代表国家，一共126国家；sitc4代表产品，共774种；export_value代表出口额。即数据中包含了1995年到2013年126个国家各国的774产品的出口数据。利用这个数据计算RCA指标，公式为：
RCA= (Eij/Eit)/(Wej/Wet)
Eij是某年度i国j产品的出口额，
Eit是某年度i国出口总额，
Wej是某年度全球j产品的出口总额（即126国家出口总额），
Wet是某年度全球产品出口总额（即126国家出口总额）
分子的含义是某年度i国j产品占i国的出口份额，分母的含义是某年度全球j产品占全球出口份额。

1. library(foreign)
   
2. mydata <- read.dta("data2.dta")
   
3. 
   
4. dat1 <- as.array(xtabs(export_value ~ iso + sitc4 + year, data = mydata))
   
5. dat2 <- apply(dat1, c(1, 3), sum)
   
6. dat3 <- apply(dat1, c(2, 3), sum)
   
7. dat4 <- apply(dat1, 3, sum)
   
8. 
   
9. n.iso <- length(unique(mydata$iso))
   
10. n.sitc4 <- length(unique(mydata$sitc4))
    
11. n.year <- length(unique(mydata$year))
    
12. RCA.array <- array(dim = dim(dat1), dimnames = dimnames(dat1))
    
13. for (i in 1:n.iso) {
    
14.         for (j in 1:n.sitc4) {
    
15.                 for (t in 1:n.year) {
    
16.                         RCA.array[i, j, t] <-
    
17.                                 (dat1[i, j, t] / dat2[i, t]) / (dat3[j, t] / dat4[t])
    
18.                 }
    
19.         }
    
20. }
    
21. 
    
22. library(reshape2)
    
23. RCA.df <- melt(RCA.array)

复制代码

有数据没？是每年每个国家每个产品算一个RCA吗？

stzhao 发表于 2017-4-16 13:14
有数据没？是每年每个国家每个产品算一个RCA吗？

有的，请留下您的邮箱，我发给您。

majormei 发表于 2017-4-16 14:19
有的，请留下您的邮箱，我发给您。

19*126*774 = 1852956，你发的图片的数据行数超过这个数，发给我看看吧 427511@qq.com

stzhao 发表于 2017-4-16 16:08
19*126*774 = 1852956，你发的图片的数据行数超过这个数，发给我看看吧 427511@qq.com

谢谢热心帮助！赞！