detailed results

fabio_SPA.R

@@ -182,28 +182,36 @@ rm(B); gc()
 #----------------------------------------
 country <- "BRA"
 product <- "Cattle"
-product <- "Soyabeans"
-cutoff <- 0.0001
+#----------------------------------------
+country <- "BRA"
+roduct <- "Soyabeans"
 #----------------------------------------
 country <- "IDN"
 product <- "Wood fuel"
 product <- "Industrial roundwood, coniferous"
 product <- "Industrial roundwood, non-coniferous"
-cutoff <- 0.001
 #----------------------------------------
 country <- "IDN"
 product <- "Oil, palm fruit"
-cutoff <- 0.0001
 #----------------------------------------
 
-
+cutoff <- 0.00001
 p <- which(index$ISO==country & index$item==product)
 
 results <- spa5(p,x,D,cutoff)
 
+results$fd <- results$rest <- results$value
+results$fd[apply(results[,1:6], 1, max, na.rm = T) < nrow(D)] <- 0
+results$rest[apply(results[,1:6], 1, min, na.rm = T) > 0] <- 0
 data.table::fwrite(results, paste0("./output/results_spa_",year,"_",country,"_",product,".csv"))
 
 
+sum(results$rest)
+sum(results$fd)
+
+
+# Prepare index for sankeys ----------------------------------------
+
 index <- rbind(index,
                data.frame(country = countries$Country[countries$ISO == colnames(Y_fabio)[1:192]],
                           ISO = colnames(Y_fabio)[1:192],

output/results_spa_2013_BRA_Cattle.csv

@@ -5,73 +5,241 @@ L0,L1,L2,L3,L4,L5,value
 2177,2201,,,,,1851.029392567
 2177,2202,,,,,3499.62558027717
 2177,2204,,,,,45.6189793946491
+2177,5835,,,,,11.0500680812195
+2177,10645,,,,,14.9804489657131
 2177,11815,,,,,97.4383446373217
+2177,11822,,,,,9.6827596770616
+2177,21825,,,,,19.0966298829048
 2177,23125,,,,,413.457842072229
+2177,23130,,,,,21.9334201413904
+2177,23131,,,,,25.1597486226068
 2177,23132,,,,,46.4573179867071
-2177,0,,,,,120.640510658745
+2177,0,,,,,18.7374352878484
+2177,2195,33850,,,,20.2566775951514
+2177,2195,34762,,,,7.62742440431271
 2177,2195,34764,,,,92.9496301150254
 2177,2195,34765,,,,121.035514641531
+2177,2195,34769,,,,10.6953636761621
+2177,2195,34771,,,,6.18571038932216
+2177,2195,34772,,,,5.24437475136561
 2177,2195,34777,,,,28733.8235873238
+2177,2195,34785,,,,18.8235082990542
 2177,2195,34790,,,,479.4817069689
 2177,2195,34805,,,,622.551789883603
+2177,2195,34806,,,,9.85903948899163
+2177,2195,34810,,,,10.7476499087502
+2177,2195,34811,,,,6.24577924838515
+2177,2195,34812,,,,7.92573646957704
+2177,2195,34814,,,,17.7628688842138
+2177,2195,34816,,,,33.319069693683
+2177,2195,34818,,,,5.8914216372973
 2177,2195,34827,,,,381.368665897359
+2177,2195,34828,,,,12.427102564267
+2177,2195,34832,,,,35.8627197413181
 2177,2195,34833,,,,224.104373278093
 2177,2195,34836,,,,67.3470889317983
 2177,2195,34837,,,,165.473742041314
+2177,2195,34839,,,,15.1508258435432
+2177,2195,34840,,,,14.107212398741
 2177,2195,34842,,,,44.4082256589384
 2177,2195,34847,,,,271.287738126841
+2177,2195,34848,,,,18.6423003238532
 2177,2195,34851,,,,64.0553911641312
 2177,2195,34854,,,,71.4350767059517
 2177,2195,34865,,,,120.939368030355
 2177,2195,34872,,,,51.693490487605
+2177,2195,34875,,,,5.63525530129437
+2177,2195,34884,,,,13.4735393451078
 2177,2195,34888,,,,70.9526542171588
+2177,2195,34897,,,,8.76091104540292
 2177,2195,34899,,,,1475.39028136829
+2177,2195,34905,,,,22.2421395008677
+2177,2195,34909,,,,9.89674980770705
 2177,2195,34910,,,,70.7247941331419
+2177,2195,34913,,,,42.263385036879
+2177,2195,34917,,,,21.3166618265154
+2177,2195,34918,,,,10.4206408251127
 2177,2195,34921,,,,70.6827465006803
+2177,2195,34925,,,,7.82901292128413
+2177,2195,34926,,,,25.012880089912
+2177,2195,34927,,,,8.55367020836541
 2177,2195,34929,,,,78.1911501208396
 2177,2195,34932,,,,216.910641885374
+2177,2195,34933,,,,10.8439499061687
 2177,2195,34934,,,,223.381673205747
 2177,2195,34938,,,,564.88621892509
+2177,2195,34939,,,,29.075006715327
+2177,2195,34946,,,,14.7704878198166
 2177,2195,34952,,,,304.024966137228
-2177,2195,0,,,,629.104550387085
+2177,2195,0,,,,142.235474719335
 2177,2200,31617,,,,160.079300472507
+2177,2200,31650,,,,33.8392546940995
 2177,2200,31652,,,,118.874792305189
+2177,2200,31678,,,,5.21919533962734
+2177,2200,33850,,,,8.81340178978898
 2177,2200,34777,,,,511.485498115203
+2177,2200,34827,,,,28.6553613652362
 2177,2200,34952,,,,761.407923377177
-2177,2200,0,,,,109.378534814702
+2177,2200,0,,,,32.8513216259495
 2177,2201,31617,,,,516.841597990013
+2177,2201,31646,,,,4.60424339350127
 2177,2201,31650,,,,109.255440392765
 2177,2201,31652,,,,383.806259987355
+2177,2201,31655,,,,4.75230166346549
+2177,2201,31678,,,,16.8510060425852
 2177,2201,34777,,,,782.749695635928
-2177,2201,0,,,,58.3763985609376
+2177,2201,0,,,,32.1688474613857
 2177,2202,31617,,,,3393.97643119373
-2177,2202,0,,,,105.649149083446
-2177,2204,0,,,,45.6189793946491
+2177,2202,33817,,,,16.5166371207697
+2177,2202,0,,,,89.1325119626764
+2177,2204,34777,,,,5.64868298211608
+2177,2204,34886,,,,4.87597571542902
+2177,2204,34936,,,,16.108229945368
+2177,2204,0,,,,18.986090751736
+2177,5835,34805,,,,11.0433508668393
+2177,5835,0,,,,0.00671721438020526
+2177,10645,34834,,,,4.64186131736256
+2177,10645,0,,,,10.3385876483505
 2177,11815,34851,,,,83.4610726521721
 2177,11815,0,,,,13.9772719851496
+2177,11822,0,,,,9.6827596770616
+2177,21825,34928,,,,19.00433244953
+2177,21825,0,,,,0.0922974333748265
 2177,23125,34938,,,,413.385047043481
 2177,23125,0,,,,0.0727950287474136
-2177,23132,0,,,,46.4573179867071
+2177,23130,34938,,,,21.9334201413904
+2177,23130,0,,,,0
+2177,23131,23107,,,,4.97227025103442
+2177,23131,34938,,,,9.00035113996334
+2177,23131,0,,,,11.187127231609
+2177,23132,34017,,,,15.6127695954839
+2177,23132,0,,,,30.8445483912232
+2177,2195,33850,0,,,20.2566775951514
+2177,2200,31617,31617,,,9.34583902812346
+2177,2200,31617,33850,,,5.356002926059
 2177,2200,31617,34777,,,113.518638106146
-2177,2200,31617,0,,,46.5606623663606
+2177,2200,31617,34954,,,5.85735881631526
+2177,2200,31617,0,,,26.0014615958629
+2177,2200,31650,34777,,,17.1130952454422
+2177,2200,31650,0,,,16.7261594486572
+2177,2200,31652,31627,,,10.069991981969
+2177,2200,31652,31632,,,6.02849973975814
+2177,2200,31652,31650,,,4.41414097829032
 2177,2200,31652,34777,,,84.708213731093
-2177,2200,31652,0,,,34.1665785740961
+2177,2200,31652,0,,,13.6539458740787
+2177,2200,31678,34777,,,4.36282986076863
+2177,2200,31678,0,,,0.856365478858716
+2177,2200,33850,0,,,8.81340178978898
+2177,2201,31617,31017,,,12.5801480250604
+2177,2201,31617,31617,,,30.1745345187984
+2177,2201,31617,33850,,,17.2927111935934
 2177,2201,31617,34777,,,366.513060384763
-2177,2201,31617,0,,,150.32853760525
+2177,2201,31617,34791,,,9.86234346715417
+2177,2201,31617,34837,,,8.5567896836287
+2177,2201,31617,34934,,,10.7706085073698
+2177,2201,31617,34954,,,18.9114187886222
+2177,2201,31617,0,,,42.1799834210225
+2177,2201,31646,0,,,4.60424339350128
+2177,2201,31650,31565,,,8.10976507937402
+2177,2201,31650,31650,,,7.54479387241504
+2177,2201,31650,31735,,,5.8146729401307
 2177,2201,31650,34777,,,55.2523622173669
-2177,2201,31650,0,,,54.0030781753975
+2177,2201,31650,0,,,32.5338462834777
+2177,2201,31652,31627,,,32.51257803067
+2177,2201,31652,31632,,,19.4639746037252
+2177,2201,31652,31633,,,4.5458677869137
+2177,2201,31652,31635,,,4.30304903424403
+2177,2201,31652,31649,,,11.2385994284325
+2177,2201,31652,31650,,,14.251759410734
 2177,2201,31652,34777,,,273.494002150371
-2177,2201,31652,0,,,110.312257836984
+2177,2201,31652,0,,,23.9964295422642
+2177,2201,31655,0,,,4.75230166346549
+2177,2201,31678,34777,,,14.0860932696634
+2177,2201,31678,0,,,2.7649127729218
+2177,2202,31617,28017,,,21.6706315567501
 2177,2202,31617,31017,,,82.6108542037444
+2177,2202,31617,31083,,,5.21891426484218
+2177,2202,31617,31085,,,10.4807705758539
+2177,2202,31617,31133,,,5.32722060134877
 2177,2202,31617,31617,,,198.149025498955
+2177,2202,31617,33017,,,5.45264818070275
 2177,2202,31617,33850,,,113.55714100944
+2177,2202,31617,34017,,,8.92390373258146
 2177,2202,31617,34777,,,2406.80450936654
 2177,2202,31617,34791,,,64.7636750099699
+2177,2202,31617,34810,,,10.6808412308513
+2177,2202,31617,34816,,,21.8383330763754
+2177,2202,31617,34828,,,4.97157476070078
 2177,2202,31617,34837,,,56.1904123543064
+2177,2202,31617,34847,,,10.7253515485365
+2177,2202,31617,34865,,,6.5756461229211
+2177,2202,31617,34899,,,5.45131670451146
+2177,2202,31617,34921,,,21.163134968435
+2177,2202,31617,34932,,,5.2021143980951
 2177,2202,31617,34934,,,70.7280365314831
 2177,2202,31617,34954,,,124.186810617861
-2177,2202,31617,0,,,276.985966601432
-2177,2202,31617,31017,0,,82.6108542037443
+2177,2202,31617,34956,,,8.96333168244435
+2177,2202,31617,34957,,,8.90525158169501
+2177,2202,31617,0,,,115.434981614787
+2177,2202,33817,34953,,,10.3409047359317
+2177,2202,33817,0,,,6.17573238483801
+2177,23131,23107,0,,,4.97227025103442
+2177,23132,34017,34954,,,7.8919723241407
+2177,23132,34017,0,,,7.72079727134318
+2177,2200,31617,31617,34777,,6.62750846174553
+2177,2200,31617,31617,0,,2.71833056637793
+2177,2200,31617,33850,0,,5.356002926059
+2177,2200,31652,31627,34777,,7.94860155224191
+2177,2200,31652,31627,0,,2.12139042972704
+2177,2200,31652,31632,0,,6.02849973975814
+2177,2200,31652,31650,0,,4.41414097829031
+2177,2201,31617,31017,0,,12.5801480250604
+2177,2201,31617,31617,34777,,21.3979699683232
+2177,2201,31617,31617,0,,8.77656455047518
+2177,2201,31617,33850,0,,17.2927111935934
+2177,2201,31650,31565,0,,8.10976507937406
+2177,2201,31650,31650,0,,7.54479387241502
+2177,2201,31650,31735,34777,,5.81190128376094
+2177,2201,31650,31735,0,,0.00277165636976609
+2177,2201,31652,31627,34777,,25.663330086529
+2177,2201,31652,31627,0,,6.84924794414099
+2177,2201,31652,31632,0,,19.4639746037252
+2177,2201,31652,31633,0,,4.54586778691369
+2177,2201,31652,31635,0,,4.30304903424403
+2177,2201,31652,31649,0,,11.2385994284325
+2177,2201,31652,31650,34777,,7.20736075353179
+2177,2201,31652,31650,0,,7.0443986572022
+2177,2202,31617,28017,34837,,6.59954729898556
+2177,2202,31617,28017,0,,15.0710842577645
+2177,2202,31617,31017,31017,,20.690779168441
+2177,2202,31617,31017,33850,,5.47422650089129
+2177,2202,31617,31017,34791,,13.9580417822427
+2177,2202,31617,31017,34934,,9.223501226804
+2177,2202,31617,31017,0,,33.2643055253653
+2177,2202,31617,31083,0,,5.21891426484218
+2177,2202,31617,31085,0,,10.4807705758539
+2177,2202,31617,31133,34791,,5.31497979206039
+2177,2202,31617,31133,0,,0.0122408092883779
+2177,2202,31617,31617,31017,,4.8230329785616
+2177,2202,31617,31617,31617,,11.5684469536448
+2177,2202,31617,31617,33850,,6.62975635972631
 2177,2202,31617,31617,34777,,140.5154036175
-2177,2202,31617,31617,0,,57.6336218814546
-2177,2202,31617,33850,0,,113.55714100944
+2177,2202,31617,31617,34954,,7.25034366107763
+2177,2202,31617,31617,0,,27.3620419284443
+2177,2202,31617,33017,0,,5.45264818070275
+2177,2202,31617,33850,31050,,6.56184148979299
+2177,2202,31617,33850,31135,,5.13718168279578
+2177,2202,31617,33850,34934,,5.46867766397125
+2177,2202,31617,33850,0,,96.3894401728796
+2177,2202,31617,34017,34954,,4.51087174828996
+2177,2202,31617,34017,0,,4.41303198429149
+2177,2202,31617,31017,31017,31017,5.18222873644839
+2177,2202,31617,31017,31017,0,15.5085504319926
+2177,2202,31617,31017,33850,0,5.4742265008913
+2177,2202,31617,31617,31017,0,4.8230329785616
+2177,2202,31617,31617,31617,34777,8.20364868727358
+2177,2202,31617,31617,31617,0,3.36479826637123
+2177,2202,31617,31617,33850,0,6.62975635972631
+2177,2202,31617,33850,31050,0,6.56184148979299
+2177,2202,31617,33850,31135,34791,4.99569802656549
+2177,2202,31617,33850,31135,0,0.141483656230292