load_perm.py

#All required modules in the permutation pipeline:
import numpy as np
from sklearn import datasets
import scipy.stats as stats
import matplotlib.pyplot as plt
import seaborn as sns
import random
import copy

irisd = datasets.load_iris()
#Load only the first two features in each species:
#Setosa:
feat_set = irisd.data[np.where(irisd.target[:]==0)[0],:2]
#Versicolor:
feat_ver = irisd.data[np.where(irisd.target[:]==1)[0],:2]
#Virginica:
feat_vir = irisd.data[np.where(irisd.target[:]==2)[0],:2]