example_diplomaticexchange.py

__author__ = 'lisette.espin'

#######################################################################
# References
# - http://www.albany.edu/faculty/kretheme/PAD637/ClassNotes/Spring%202013/Lab8.pdf
#######################################################################

#######################################################################
# Dependencies
#######################################################################
import numpy as np
from libs.mrqap import MRQAP

#######################################################################
# Data
# Source: http://vlado.fmf.uni-lj.si/pub/networks/data/ucinet/ucidata.htm
#######################################################################
X1 = np.loadtxt('data/crudematerials.dat')
X2 = np.loadtxt('data/foods.dat')
X3 = np.loadtxt('data/manufacturedgoods.dat')
X4 = np.loadtxt('data/minerals.dat')
Y = np.loadtxt('data/diplomatic.dat')
X = {'CRUDEMATERIALS':X1, 'FOODS':X2, 'MANUFACTUREDGOODS':X3, 'MINERALS':X4}
Y = {'DIPLOMATIC':Y}
np.random.seed(473)

#######################################################################
# QAP
#######################################################################
mrqap = MRQAP(Y=Y, X=X, npermutations=2000, diagonal=False, directed=True)
mrqap.mrqap()
mrqap.summary()
mrqap.plot('betas')
mrqap.plot('tvalues')