Sorry for that. Here it goes:
mydata <- data.frame(
stringsAsFactors = FALSE,
row.names = c("ASU_PRY","SAC_PRY","POP_PRY",
"VEN_PRY","IND_PRY","BOL2_PRY","EGY_PRY","ISR_PRY",
"POP_BRA","BOL1_URY","ASU_ARG","SAC_URY","ADH_PRY",
"PAR_PRY","EGY_URY","SAC_ARG","BOL1_ARG","EGY_BRA",
"ISR_BRA","BOL2_URY","VEN_BRA","ISR_ARG","ADH_BRA",
"IND_URY","ISR_URY","IND_ARG","IND_BRA","PAR_BRA",
"EGY_ARG","SAC2_BRA","OLI_BRA","ADH_URY","USH_BRA",
"ADH_ARG","PAR_ARG","VEN_ARG","PAR_URY","VEN_URY",
"OLI_PRY","USH_ARG","POP_ARG","OLI_ARG","SAC_BRA","USH_URY",
"USH_PRY","OLI_URY","POP_URY","ASU_URY","ASU_BRA"),
bill = c("ASU_PRY","SAC_PRY","POP_PRY",
"VEN_PRY","IND_PRY","BOL2_PRY","EGY_PRY","ISR_PRY",
"POP_BRA","BOL1_URY","ASU_ARG","SAC_URY","ADH_PRY",
"PAR_PRY","EGY_URY","SAC_ARG","BOL1_ARG","EGY_BRA",
"ISR_BRA","BOL2_URY","VEN_BRA","ISR_ARG","ADH_BRA",
"IND_URY","ISR_URY","IND_ARG","IND_BRA","PAR_BRA",
"EGY_ARG","SAC2_BRA","OLI_BRA","ADH_URY","USH_BRA",
"ADH_ARG","PAR_ARG","VEN_ARG","PAR_URY","VEN_URY",
"OLI_PRY","USH_ARG","POP_ARG","OLI_ARG","SAC_BRA","USH_URY",
"USH_PRY","OLI_URY","POP_URY","ASU_URY","ASU_BRA"),
cntry = c("Paraguay","Paraguay",
"Paraguay","Paraguay","Paraguay","Paraguay","Paraguay",
"Paraguay","Brazil","Uruguay","Argentina","Uruguay",
"Paraguay","Paraguay","Uruguay","Argentina","Argentina",
"Brazil","Brazil","Uruguay","Brazil","Argentina",
"Brazil","Uruguay","Uruguay","Argentina","Brazil",
"Brazil","Argentina","Brazil","Brazil","Uruguay","Brazil",
"Argentina","Argentina","Argentina","Uruguay",
"Uruguay","Paraguay","Argentina","Argentina","Argentina",
"Brazil","Uruguay","Paraguay","Uruguay","Uruguay",
"Uruguay","Brazil"),
yearsubm = c(1991L,2010L,1995L,2007L,
2007L,2015L,2015L,2009L,1995L,2013L,1991L,2010L,
2006L,2006L,2013L,2011L,2013L,2013L,2008L,2016L,2007L,
2009L,2006L,2008L,2008L,2007L,2005L,2006L,2017L,
2015L,2002L,2006L,1999L,2006L,2006L,2006L,2006L,
2006L,2002L,1998L,1995L,2002L,2010L,1998L,1998L,
2002L,1995L,1991L,1991L),
yearappr = c(1991L,2012L,1995L,2013L,
2007L,2016L,2016L,2010L,1995L,2014L,1991L,2011L,
2006L,2006L,2015L,2012L,2014L,2015L,2009L,2016L,2009L,
2011L,2009L,2008L,2008L,2008L,2008L,2006L,2017L,
2015L,2004L,2008L,2001L,2006L,2006L,2006L,2006L,
2006L,2002L,1999L,1995L,2002L,2010L,1999L,1998L,
2003L,1995L,1991L,1991L)
)
gdp_long <- data.frame(
stringsAsFactors = FALSE,
year = c("1990","1990","1990","1990",
"1991","1991","1991","1991","1992","1992","1992",
"1992","1993","1993","1993","1993","1994","1994",
"1994","1994","1995","1995","1995","1995","1996",
"1996","1996","1996","1997","1997","1997","1997",
"1998","1998","1998","1998","1999","1999","1999",
"1999","2000","2000","2000","2000","2001","2001",
"2001","2001","2002","2002","2002","2002","2003","2003",
"2003","2003","2004","2004","2004","2004","2005",
"2005","2005","2005","2006","2006","2006","2006",
"2007","2007","2007","2007","2008","2008","2008",
"2008","2009","2009","2009","2009","2010","2010",
"2010","2010","2011","2011","2011","2011","2012",
"2012","2012","2012","2013","2013","2013","2013",
"2014","2014","2014","2014","2015","2015","2015",
"2015","2016","2016","2016","2016","2017","2017","2017",
"2017","2018","2018","2018","2018"),
country = c("Argentina","Brazil",
"Paraguay","Uruguay","Argentina","Brazil","Paraguay",
"Uruguay","Argentina","Brazil","Paraguay","Uruguay",
"Argentina","Brazil","Paraguay","Uruguay","Argentina",
"Brazil","Paraguay","Uruguay","Argentina","Brazil",
"Paraguay","Uruguay","Argentina","Brazil","Paraguay",
"Uruguay","Argentina","Brazil","Paraguay","Uruguay",
"Argentina","Brazil","Paraguay","Uruguay","Argentina",
"Brazil","Paraguay","Uruguay","Argentina","Brazil",
"Paraguay","Uruguay","Argentina","Brazil","Paraguay",
"Uruguay","Argentina","Brazil","Paraguay","Uruguay",
"Argentina","Brazil","Paraguay","Uruguay","Argentina",
"Brazil","Paraguay","Uruguay","Argentina","Brazil",
"Paraguay","Uruguay","Argentina","Brazil","Paraguay",
"Uruguay","Argentina","Brazil","Paraguay","Uruguay",
"Argentina","Brazil","Paraguay","Uruguay","Argentina",
"Brazil","Paraguay","Uruguay","Argentina","Brazil",
"Paraguay","Uruguay","Argentina","Brazil","Paraguay",
"Uruguay","Argentina","Brazil","Paraguay","Uruguay",
"Argentina","Brazil","Paraguay","Uruguay",
"Argentina","Brazil","Paraguay","Uruguay","Argentina","Brazil",
"Paraguay","Uruguay","Argentina","Brazil",
"Paraguay","Uruguay","Argentina","Brazil","Paraguay",
"Uruguay","Argentina","Brazil","Paraguay","Uruguay"),
gdp = c(204022.709576019,
1190328.7904512,15072.4009298585,22100.9277328455,225605.750803947,
1202589.17699285,15444.4660701346,22883.0372408505,
247277.638553345,1196095.19543709,15706.4703187341,
24698.0255693715,261428.668260367,1254943.07905259,
16481.7981450403,25354.3869873348,276686.170019566,
1328357.24917717,17358.286599311,27200.5267832485,
268813.868517326,1384413.92509244,18542.6095659408,26806.7722354053,
283670.377245591,1414185.73384553,18834.4303983244,
28302.0426702529,306679.014181442,1462195.16145589,
19633.4841408321,29730.6489648803,318486.704781525,
1467138.81261922,19646.8423228185,31074.1443152037,
307704.474172928,1474004.10627881,19378.4507957219,
30471.5507314667,305276.689188273,1538682.66125963,18930.0062043638,
29883.4709399851,291817.529327736,1560068.7562445,
18772.11959474,28734.7114813469,260025.512893345,
1607704.86065684,18768.1015361058,26512.9415151065,
283004.073515091,1626046.02392186,19579.0234359506,
26726.4459698004,308558.133837269,1719705.69987883,20373.4263282599,
28063.8801838554,335870.650498235,1774772.93573377,
20808.0914770334,30157.4827229858,362898.670600149,
1845089.2393129,21808.3608199,31393.5104782487,
395587.315892942,1957083.76872605,22990.7278981641,33447.120372214,
411637.215422162,2056781.44097077,24452.736056455,
35847.3341150742,387274.363587343,2054193.76311241,
24389.6549826671,37368.513657464,426487.434874632,
2208838.10859319,27107.573694974,40284.5318602359,
452093.534425313,2296626.68016874,28259.404602775,42364.0729820816,
447453.153875727,2340748.92041541,28107.2240515083,
43862.9895982757,458215.851022343,2411084.2746306,
30473.148310721,45897.1526850981,446702.64932211,
2423235.07223542,31954.4146425446,47383.6656404033,
458902.812680434,2337312.89010631,32938.7301255886,47559.3364376997,
449356.129667934,2260744.4619857,34359.3462542358,
48362.9932304934,462021.24554028,2290651.15086124,
36061.4300076708,49616.242186666,450168.774848251,
2320824.15750549,37271.7135902413,50420.0668036221)
)
master_dataset <- fuzzy_left_join(mydata,
gdp_long,
by=c("yearsubm"="Year",
"yearappr"="Year",
"cntry"="country"),
match_fun = list(`>=`,`<=`,`==`)
)
master_dataset %>% group_by(bill,cntry) %>% summarise(mean_gdp=mean(gdp,na.rm=TRUE))