我有一个关于并购 (M&As) 的大数据框(900k 行)。
df 有四列:date(并购完成时),目标国家(被合并/收购的国家的公司),收购方_国家(收购方是哪个国家的公司),以及big_corp(无论收购方是否是大公司,其中 TRUE 表示该公司很大)。
这是我的 df 的示例:
> df <- structure(list(date = c(2000L, 2000L, 2001L, 2001L, 2001L, 2002L,
2002L, 2002L), target_nation = c("Uganda", "Uganda", "Uganda",
"Uganda", "Uganda", "Uganda", "Uganda", "Uganda"), acquiror_nation = c("France",
"Germany", "France", "France", "Germany", "France", "France",
"Germany"), big_corp_TF = c(TRUE, FALSE, TRUE, FALSE, FALSE,
TRUE, TRUE, TRUE)), row.names = c(NA, -8L))
> df
date target_nation acquiror_nation big_corp_TF
1: 2000 Uganda France TRUE
2: 2000 Uganda Germany FALSE
3: 2001 Uganda France TRUE
4: 2001 Uganda France FALSE
5: 2001 Uganda Germany FALSE
6: 2002 Uganda France TRUE
7: 2002 Uganda France TRUE
8: 2002 Uganda Germany TRUE
根据这些数据,我想创建一个新变量,表示特定收购方国家的大公司进行的并购的份额,计算2年的平均值。(对于我的实际练习,我将计算 5 年的平均值,但让我们在这里简单一点)。因此,法国的大公司将会有一个新的变量,德国的大公司也会有一个新的变量。
到目前为止,我所做的是1)统计某一年特定目标国家的并购总数; 2)统计某一收购方国家的大公司在某一年在特定目标国家进行的并购总数。我加入了这两个dfs以方便计算我想要的平均值。这是我使用的代码和生成的新 df:
##counting total rows for target nations
df2 <- df %>%
group_by(date, target_nation) %>%
count(target_nation)
##counting total rows conducted by small or big corps for certain acquiror nations
df3 <- df %>%
group_by(date, target_nation, acquiror_nation) %>%
count(big_corp_TF)
##selecting rows that were conducted by big corps
df33 <- df3 %>%
filter(big_corp_TF == TRUE)
##merging df2 and df33
df4 <- df2 %>%
left_join(df33, by = c("date" = "date", "target_nation" = "target_nation"))
df4 <- as.data.frame(df4)
> df4
date target_nation n.x acquiror_nation big_corp_TF n.y
1 2000 Uganda 2 France TRUE 1
2 2001 Uganda 3 France TRUE 1
3 2002 Uganda 3 France TRUE 2
4 2002 Uganda 3 Germany TRUE 1
这里的 n.x 是某一年特定 target_nation 的并购总数(行); n.y 是特定收购者国家/地区的大公司在特定目标国家/地区进行的并购总数(行)。
有了这个新的数据框架 df4,我现在可以轻松计算特定收购方国家的大公司在特定年份在特定目标国家中进行的并购的份额。例如,让我们计算一下法国的这一份额:
df5 <- df4 %>%
filter(acquiror_nation == "France") %>%
mutate(France_bigcorp_share_1year = n.y / n.x)
date target_nation n.x acquiror_nation big_corp_TF n.y France_bigcorp_share_1year
1 2000 Uganda 2 France TRUE 1 0.5000000
2 2001 Uganda 3 France TRUE 1 0.3333333
3 2002 Uganda 3 France TRUE 2 0.6666667
然而,我无法弄清楚如何计算特定收购方国家的大公司进行的并购的份额,计算平均数2 years.
这就是所需变量的样子:
date target_nation n.x acquiror_nation big_corp_TF n.y France_bigcorp_share_2years
1 2000 Uganda 2 France TRUE 1 0.5000000
2 2001 Uganda 3 France TRUE 1 0.4000000
3 2002 Uganda 3 France TRUE 2 0.5000000
请注意,2000 年的份额将保持不变,因为没有前一年可以使其成为 2 年平均值; 2001 年它将变为 0.4(因为 (1+1)/(2+3) = 0.4); 2002 年它将变为 0.5(因为 (1+2)/(3+3) = 0.5)。
您知道如何编写计算两年平均份额的代码吗?我想我需要在这里使用 for 循环,但我不知道如何使用。任何建议,将不胜感激。
--
EDIT:AnilGoyal 的代码与示例数据完美配合,但我的实际数据显然更混乱,因此我想知道是否有解决我遇到的问题的方法。
我的实际数据集有时会跳过一年,或者有时不包括前几行中包含的 acquiror_nations。请查看我的实际数据的更准确样本:
> df_new <- structure(list(date = c(2000L, 2000L, 2001L, 2001L, 2001L, 2002L,
2002L, 2002L, 2003L, 2003L, 2004L, 2004L, 2004L, 2006L, 2006L
), target_nation = c("Uganda", "Uganda", "Uganda", "Uganda",
"Uganda", "Uganda", "Uganda", "Uganda", "Uganda", "Uganda", "Uganda",
"Uganda", "Uganda", "Uganda", "Uganda"), acquiror_nation = c("France",
"Germany", "France", "France", "Germany", "France", "France",
"Germany", "Germany", "Germany", "France", "France", "Germany",
"France", "France"), big_corp_TF = c(TRUE, FALSE, TRUE, FALSE, FALSE,
TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, FALSE, TRUE, TRUE, TRUE)), row.names = c(NA,
-15L))
> df_new
date target_nation acquiror_nation big_corp_TF
1: 2000 Uganda France TRUE
2: 2000 Uganda Germany FALSE
3: 2001 Uganda France TRUE
4: 2001 Uganda France FALSE
5: 2001 Uganda Germany FALSE
6: 2002 Uganda France TRUE
7: 2002 Uganda France TRUE
8: 2002 Uganda Germany TRUE
9: 2003 Uganda Germany TRUE
10: 2003 Uganda Germany FALSE
11: 2004 Uganda France TRUE
12: 2004 Uganda France FALSE
13: 2004 Uganda Germany TRUE
14: 2006 Uganda France TRUE
15: 2006 Uganda France TRUE
注意:2003 年法国没有任何行;并且没有2005年。
如果我运行 Anil 的第一个代码,结果如下:
date target_nation acquiror_nation n1 n2 share
<int> <chr> <chr> <dbl> <int> <dbl>
1 2000 Uganda France 2 1 0.5
2 2001 Uganda France 3 1 0.4
3 2002 Uganda France 3 2 0.5
4 2004 Uganda France 3 1 0.5
5 2006 Uganda France 2 2 0.6
注:法国 2003 年和 2005 年没有结果;我希望有 2003 年和 2005 年的结果(因为我们正在计算 2 年平均值,因此我们应该能够得到 2003 年和 2005 年的结果)。另外,2006 年的份额实际上是不正确的,因为它应该是 1(应该采用 2005 年的值(均为 0)而不是 2004 年的值来计算平均值)。
我希望能够收到以下内容:
date target_nation acquiror_nation n1 n2 share
<int> <chr> <chr> <dbl> <int> <dbl>
1 2000 Uganda France 2 1 0.5
2 2001 Uganda France 3 1 0.4
3 2002 Uganda France 3 2 0.5
4 2003 Uganda France 2 0 0.4
5 2004 Uganda France 3 1 0.2
6 2005 Uganda France 0 0 0.33
7 2006 Uganda France 2 2 1.0
注意:请注意 2006 年的结果也不同(因为我们现在采用 2005 年而不是 2004 年作为两年平均值)。
您认为有可能找到一种方法来输出所需的小标题吗?我知道这是原始数据的问题:它只是缺少某些数据点。然而,将它们包含到原始数据集中似乎非常不方便;最好将它们包含在中间,例如计算完n1和n2后。但最方便的方法是什么?
EDIT2:Anil 的新代码可以很好地处理上面的数据样本,但在处理更复杂的数据样本(包括多个 target_nation)时,它会遇到一个不受欢迎的问题。这是一个更短但更复杂的数据样本:
> df_new_complex <- structure(list(date = c(2000L, 2000L, 2001L, 2001L, 2001L, 2003L,
2003L, 1999L, 2001L, 2002L, 2002L), target_nation = c("Uganda",
"Uganda", "Uganda", "Uganda", "Uganda", "Uganda", "Uganda", "Mozambique",
"Mozambique", "Mozambique", "Mozambique"), acquiror_nation = c("France",
"Germany", "France", "France", "Germany", "Germany", "Germany",
"Germany", "France", "France", "Germany"), big_corp_TF = c(TRUE,
FALSE, TRUE, FALSE, FALSE, TRUE, FALSE, FALSE, TRUE, FALSE, TRUE
)), row.names = c(NA, -11L))
> df_new_complex
date target_nation acquiror_nation big_corp_TF
1: 2000 Uganda France TRUE
2: 2000 Uganda Germany FALSE
3: 2001 Uganda France TRUE
4: 2001 Uganda France FALSE
5: 2001 Uganda Germany FALSE
6: 2003 Uganda Germany TRUE
7: 2003 Uganda Germany FALSE
8: 1999 Mozambique Germany FALSE
9: 2001 Mozambique France TRUE
10: 2002 Mozambique France FALSE
11: 2002 Mozambique Germany TRUE
如您所见,此数据样本包含两个 target_nations。阿尼尔的代码,其中param <- c("France", "Germany")
,产生以下小标题:
date target_nation acquiror_nation n1 n2 share
<dbl> <chr> <chr> <dbl> <int> <dbl>
1 1999 Mozambique France 1 0 0
2 1999 Mozambique Germany 1 0 0
3 1999 Uganda France 0 0 0
4 1999 Uganda Germany 0 0 0
5 2000 Mozambique France 0 0 0
6 2000 Mozambique Germany 0 0 0
7 2000 Uganda France 2 1 0.25
8 2000 Uganda Germany 2 0 0.167
9 2001 Mozambique France 1 1 0.4
10 2001 Mozambique Germany 1 0 0.333
11 2001 Uganda France 3 1 0.333
12 2001 Uganda Germany 3 0 0.25
13 2002 Mozambique France 2 0 0.2
14 2002 Mozambique Germany 2 1 0.25
15 2002 Uganda France 0 0 0.25
16 2002 Uganda Germany 0 0 0.25
17 2003 Mozambique France 0 0 0.25
18 2003 Mozambique Germany 0 0 0.25
19 2003 Uganda France 2 0 0.167
20 2003 Uganda Germany 2 1 0.25
这里不受欢迎的是,该代码为乌干达创建了 1999 年,为莫桑比克创建了 2003 年(后者问题不大)。 1999 年,乌干达没有数据样本中所示的投资,因此拥有数值是没有意义的(它可能有 NA,或者根本不存在)。莫桑比克在2003年也没有投资,所以我不想计算莫桑比克当年的份额。
我找到了一个解决方法,我可以在代码的早期过滤特定的目标国家,就像这样:
correct1 <- df_new_complex %>%
filter(target_nation == "Mozambique") %>%
mutate(d = 1) %>% ...
#I do the same for another target_nation
correct2 <- df_new_complex %>%
filter(target_nation == "Uganda") %>%
mutate(d = 1) %>% ...
#I then use rbind
correct <- rbind(correct1, correct2)
#which produces the desired tibble (without a year 2003 for Mozambique and 1999 for Uganda).
> correct
date target_nation acquiror_nation n1 n2 share
<dbl> <chr> <chr> <dbl> <int> <dbl>
1 1999 Mozambique France 1 0 0
2 1999 Mozambique Germany 1 0 0
3 2000 Mozambique France 0 0 0
4 2000 Mozambique Germany 0 0 0
5 2001 Mozambique France 1 1 1
6 2001 Mozambique Germany 1 0 0
7 2002 Mozambique France 2 0 0.33
8 2002 Mozambique Germany 2 1 0.333
9 2000 Uganda France 2 1 0.5
10 2000 Uganda Germany 2 0 0.25
11 2001 Uganda France 3 1 0.286
12 2001 Uganda Germany 3 0 0.2
13 2002 Uganda France 0 0 0.167
14 2002 Uganda Germany 0 0 0.167
15 2003 Uganda France 2 0 0
16 2003 Uganda Germany 2 1 0.25
什么是更快的方法来做到这一点?我有一份所需的 target_nations 列表。也许可以创建一个循环,我可以通过该循环对一个 target_nation 进行计算,然后对另一个进行计算;然后将它们绑定;然后是另一个;然后rbind等。或者有更好的方法吗?