"""
from ep_bioinfo_load105 import *
from Bio import Cluster
from matplotlib import pyplot as plt
try:
mPD
PD
ProtNames
except:
mPD = load_mPD()
PD,ProtNames = load_PD()
R = Cluster.somcluster(mPD,transpose=1,nxgrid=40,nygrid=40,niter=1)
from minisom import MiniSom
### Initialization and training ###
som = MiniSom(40,40,15,sigma=1.0,learning_rate=0.5)
#som.random_weights_init(mPD)
som.weights
som.random_weights_init(transpose(mPD))
print("Training...")
som.train_random(transpose(mPD),100) # training with 100 iterations
print("\n...ready!")
timg = np.zeros(shape=(40,40))
for c in R[0]:
timg[c[0],c[1]]=timg[c[0],c[1]]+1
