from PyQt5.QtGui import *
from matplotlib.backends.backend_qt5 import FigureCanvasQT
from PyQt5 import QtGui, QtCore, QtWidgets, uic
import joblib
import sys, os, time
import seaborn as sns
import numpy as np
import pandas as pd
from NiBAx.core.plotcanvas import PlotCanvas
from NiBAx.core.baseplugin import BasePlugin
from NiBAx.core.gui.SearchableQComboBox import SearchableQComboBox
[documentos]class SPARE(QtCore.QObject):
"""This class is an adapter to the worker class."""
#constructor
def __init__(self,datamodel, model):
super(SPARE,self).__init__()
self.datamodel = datamodel
self.sparemodel = model
self.SPAREs = None
# Initialize thread
self.thread = None
[documentos] def InterruptComputation(self):
self.thread.requestInterruption()
[documentos] def DoSPAREsComputation(self,asynchronous = False):
if(asynchronous):
# Using this example: https://realpython.com/python-pyqt-qthread/
self.thread = QtCore.QThread()
self.worker = BrainAgeWorker(self.datamodel.data, self.sparemodel)
self.worker.moveToThread(self.thread)
self.thread.started.connect(self.worker.run)
self.worker.done.connect(self.thread.quit)
self.worker.done.connect(self.worker.deleteLater)
self.thread.finished.connect(self.thread.deleteLater)
self.worker.progress.connect(self.updateProgress)
self.worker.done.connect(lambda y_hat: self.OnComputationDone(y_hat))
self.thread.start()
else:
worker = BrainAgeWorker(self.datamodel.data,self.sparemodel)
worker.progress.connect(self.updateProgress)
worker.done.connect(lambda y_hat: self.OnComputationDone(y_hat))
worker.run()
[documentos] def OnComputationDone(self, y_hat):
self.SPAREs = y_hat
self.done.emit()
[documentos] def updateProgress(self, txt, vl):
msg = txt + " | progress: " + str(vl) + str("\\") + str(len(self.sparemodel['BrainAge']['scaler'])-1)
print(msg,flush=True)
self.sendprogress.emit(txt, vl)
[documentos]class BrainAgeWorker(QtCore.QObject):
"""Class for doing spare computation."""
#constructor
def __init__(self, data, model):
super(BrainAgeWorker, self).__init__()
self.data = data
self.model = model
[documentos] def run(self):
y_hat = pd.DataFrame.from_dict({'SPARE_BA': np.full((self.data.shape[0],),np.nan),
'SPARE_AD': np.full((self.data.shape[0],),np.nan)})
# SPARE-BA
idx = ~self.data[self.model['BrainAge']['predictors'][0]].isnull()
y_hat_test = np.zeros((np.sum(idx),))
n_ensembles = np.zeros((np.sum(idx),))
for i,_ in enumerate(self.model['BrainAge']['scaler']):
# Predict validation (fold) and test
self.progress.emit('Computing SPARE-BA | Task 1 of 2', i)
if QtCore.QThread.currentThread().isInterruptionRequested():
self.progress.emit('Cancelled.', 0)
self.done.emit(pd.DataFrame())
return
test = np.logical_not(self.data[idx]['participant_id'].isin(np.concatenate(self.model['BrainAge']['train']))) | self.data[idx]['participant_id'].isin(self.model['BrainAge']['validation'][i])
X = self.data[idx].loc[test, self.model['BrainAge']['predictors']].values
X = self.model['BrainAge']['scaler'][i].transform(X)
y_hat_test[test] += (self.model['BrainAge']['svm'][i].predict(X) - self.model['BrainAge']['bias_ints'][i]) / self.model['BrainAge']['bias_slopes'][i]
n_ensembles[test] += 1.
y_hat_test /= n_ensembles
y_hat.loc[idx.values, 'SPARE_BA'] = y_hat_test
idx = ~self.data[self.model['AD']['predictors'][0]].isnull()
y_hat_test = np.zeros((np.sum(idx),))
n_ensembles = np.zeros((np.sum(idx),))
for i,_ in enumerate(self.model['AD']['scaler']):
# Predict validation (fold) and test
self.progress.emit('Computing SPARE-AD | Task 2 of 2', i)
if QtCore.QThread.currentThread().isInterruptionRequested():
self.progress.emit('Cancelled.', 0)
# Emit the result
self.done.emit(pd.DataFrame())
return
test = np.logical_not(self.data[idx]['participant_id'].isin(np.concatenate(self.model['AD']['train']))) | self.data[idx]['participant_id'].isin(self.model['AD']['validation'][i])
X = self.data[idx].loc[test, self.model['AD']['predictors']].values
X = self.model['AD']['scaler'][i].transform(X)
y_hat_test[test] += self.model['AD']['svm'][i].decision_function(X)
n_ensembles[test] += 1.
y_hat_test /= n_ensembles
y_hat.loc[idx.values, 'SPARE_AD'] = y_hat_test
self.progress.emit('All done.', i)
# Emit the result
self.done.emit(y_hat)