These are the codes used in the article Pizzol M, Scotti M, Identifying marginal supplying countries of wood products via trade network analysis, The International Journal of Life Cycle Assessment (2016) DOI: 10.1007/s11367-016-1222-6
The publisher version and supplementary material (SM) are available at the IJLCA website, a free preprint version is available on my Researchgate profile. All material is available on email request ([email protected]).
- Forestry_Trade_Flows_E_All_Data_reduced.csv is a reduced sample of the trade data, i.e. trade of Sawnwood (C). If you want to reproduce all results, you need to download the full data from FAOSTAT here (select everything and press the .csv icon) and check that they are in the same format as in the example provided here
- Forestry_E_All_Data_reduced.csv is a reduced sample of the production data, i.e. production of Sawnwood (C). If you want to reproduce all results, you need to download the full data from FAOSTAT here (select everything and press the .csv icon) and check that they are in the same format as in the example provided here
- Demo.R is a R script demo where you can play around with igraph package on the wood trade data
- Historical_increments.R is a R script for calculating historical increments from time series of production data
- Countries_selection.txt a list of countries
- Countries_selection_ROW.txt a list of world regions
The next three scripts has to be run in order. The first two take quite some time to run...so if you are only interested in their output I suggest you retrieve the article's SM.
- 1_Wood_trade_network_loop_100.R is a R script for identifying communities and looping them. You will need the full FAOSTAT trade data see above to make this work.
- 2_Wood_trade_network_clusters_to_contingency.R is the R script to build the contingency tables.
- 3_Wood trade network_figures.R is a script used to generate graphics. You can load it using the SM thus skipping 1 and 2.