celldetective 1.2.1__py3-none-any.whl → 1.2.2.post1__py3-none-any.whl

celldetective/gui/plot_signals_ui.py

@@ -1,40 +1,24 @@
-from PyQt5.QtWidgets import QMainWindow, QApplication, QMessageBox, QScrollArea, QButtonGroup, QComboBox, QFrame, \
-	QCheckBox, QFileDialog, QGridLayout, QTextEdit, QLineEdit, QVBoxLayout, QWidget, QLabel, QHBoxLayout, QPushButton, \
-	QRadioButton, QSpacerItem, QSizePolicy
-from PyQt5.QtCore import Qt, QSize, QRect
+from PyQt5.QtWidgets import QMessageBox, QComboBox, \
+	QCheckBox, QLineEdit, QVBoxLayout, QWidget, QLabel, QHBoxLayout, QPushButton
+from PyQt5.QtCore import Qt
 from PyQt5.QtGui import QIcon, QDoubleValidator
-from sklearn.preprocessing import MinMaxScaler
-
-from celldetective.gui.gui_utils import center_window, FeatureChoice, ListWidget, QHSeperationLine, FigureCanvas, GeometryChoice, OperationChoice
-from superqt import QLabeledSlider, QColormapComboBox
-from superqt.fonticon import icon
-from fonticon_mdi6 import MDI6
-from celldetective.utils import extract_experiment_channels, get_software_location, _extract_labels_from_config
-from celldetective.io import interpret_tracking_configuration, load_frames, auto_load_number_of_frames, load_experiment_tables
-from celldetective.measure import compute_haralick_features, contour_of_instance_segmentation
-from celldetective.signals import columnwise_mean, mean_signal
+
+from celldetective.gui.gui_utils import center_window
+from celldetective.gui.generic_signal_plot import GenericSignalPlotWidget
+from superqt import QLabeledSlider, QColormapComboBox, QSearchableComboBox
+from celldetective.utils import get_software_location, _extract_labels_from_config
+from celldetective.io import load_experiment_tables
+from celldetective.signals import mean_signal
 import numpy as np
-from tifffile import imread
-import json
-from shutil import copyfile
 import os
 import matplotlib.pyplot as plt
 plt.rcParams['svg.fonttype'] = 'none'
-from mpl_toolkits.axes_grid1 import make_axes_locatable
 from glob import glob
 from natsort import natsorted
-from tifffile import imread
-from pathlib import Path, PurePath
-import gc
 import pandas as pd
-from tqdm import tqdm
-from lifelines import KaplanMeierFitter
-from matplotlib.cm import viridis, tab10
 import math
 from celldetective.gui import Styles
 from matplotlib import colormaps
-import matplotlib.cm as mcm
-
 
 
 class ConfigSignalPlot(QWidget, Styles):
@@ -63,24 +47,17 @@ class ConfigSignalPlot(QWidget, Styles):
 		self.ax2=None
 		self.auto_close = False
 
-
-		print('Parent wells: ', self.wells)
-
-
 		self.well_option = self.parent_window.parent_window.well_list.currentIndex()
 		self.position_option = self.parent_window.parent_window.position_list.currentIndex()
 		self.interpret_pos_location()
-		#self.load_available_tables()
-		#self.config_path = self.exp_dir + self.config_name
 
 		self.screen_height = self.parent_window.parent_window.parent_window.screen_height
 		center_window(self)
-		#self.setMinimumHeight(int(0.8*self.screen_height))
-		#self.setMaximumHeight(int(0.8*self.screen_height))
 		self.populate_widget()
-		#self.load_previous_measurement_instructions()
+
 		if self.auto_close:
 			self.close()
+		self.setAttribute(Qt.WA_DeleteOnClose)
 
 	def interpret_pos_location(self):
 		
@@ -119,7 +96,7 @@ class ConfigSignalPlot(QWidget, Styles):
 		main_layout.addWidget(panel_title, alignment=Qt.AlignCenter)
 
 		labels = [QLabel('population: '), QLabel('class: '), QLabel('time of\ninterest: '), QLabel('cmap: ')]
-		self.cb_options = [['targets','effectors'],['class'], ['t0'], list(plt.colormaps())]
+		self.cb_options = [['targets','effectors'],['class'], ['t0'], []]
 		self.cbs = [QComboBox() for i in range(len(labels))]
 		self.cbs[-1] = QColormapComboBox()
 		
@@ -135,7 +112,8 @@ class ConfigSignalPlot(QWidget, Styles):
 				self.cbs[i].addItems(self.cb_options[i])
 			choice_layout.addLayout(hbox)
 		
-		for cm in list(colormaps):
+		all_cms = list(colormaps)
+		for cm in all_cms:
 			try:
 				self.cbs[-1].addColormap(cm)
 			except:
@@ -244,42 +222,26 @@ class ConfigSignalPlot(QWidget, Styles):
 		self.feature_choice_widget.setWindowTitle("Select numeric feature")
 		layout = QVBoxLayout()
 		self.feature_choice_widget.setLayout(layout)
-		self.feature_cb = QComboBox()
+		self.feature_cb = QSearchableComboBox()
 		self.feature_cb.addItems(feats)
 		hbox = QHBoxLayout()
 		hbox.addWidget(QLabel('feature: '), 33)
 		hbox.addWidget(self.feature_cb, 66)
 		#hbox.addWidget((QLabel('Plot two features')))
 		layout.addLayout(hbox)
-		self.checkBox_feature = QCheckBox(self.feature_choice_widget)
-		self.checkBox_feature.setGeometry(QRect(170, 120, 81, 20))
-		self.checkBox_feature.setText('Plot two features')
-		#hbox.addWidget(self.checkBox_feature)
-		layout.addWidget(self.checkBox_feature, alignment=Qt.AlignCenter)
-		self.checkBox_feature.stateChanged.connect(self.enable_second_feature)
-
-		#layout = QVBoxLayout()
-		#self.feature_two_choice_widget.setLayout(layout)
-		self.feature_two_cb = QComboBox()
-		self.feature_two_cb.addItems(feats)
-		self.feature_two_cb.setEnabled(False)
-		hbox_two = QHBoxLayout()
-		hbox_two.addWidget(QLabel('feature two: '), 33)
-		hbox_two.addWidget(self.feature_two_cb, 66)
-		layout.addLayout(hbox_two)
-
 
 		self.set_feature_btn = QPushButton('set')
+		self.set_feature_btn.setStyleSheet(self.button_style_sheet)
 		self.set_feature_btn.clicked.connect(self.compute_signals)
 		layout.addWidget(self.set_feature_btn)
 		self.feature_choice_widget.show()
 		center_window(self.feature_choice_widget)
 
-	def enable_second_feature(self):
-		if self.checkBox_feature.isChecked():
-			self.feature_two_cb.setEnabled(True)
-		else:
-			self.feature_two_cb.setEnabled(False)
+	# def enable_second_feature(self):
+	# 	if self.checkBox_feature.isChecked():
+	# 		self.feature_two_cb.setEnabled(True)
+	# 	else:
+	# 		self.feature_two_cb.setEnabled(False)
 	
 	def compute_signals(self):
 
@@ -293,245 +255,38 @@ class ConfigSignalPlot(QWidget, Styles):
 				print(e, ' The query is misunderstood and will not be applied...')
 
 			self.feature_selected = self.feature_cb.currentText()
-			if self.checkBox_feature.isChecked():
-				self.second_feature_selected=self.feature_two_cb.currentText()
-				print(self.second_feature_selected)
 			self.feature_choice_widget.close()
 			self.compute_signal_functions()
-			# prepare survival
-
-			# plot survival
-			self.survivalWidget = QWidget()
-			self.scroll = QScrollArea()
-			self.survivalWidget.setMinimumHeight(int(0.8*self.screen_height))
-			self.survivalWidget.setWindowTitle('signals')
-			self.plotvbox = QVBoxLayout(self.survivalWidget)
-			self.plotvbox.setContentsMargins(30,30,30,30)
-			self.survival_title = QLabel('Signal function')
-			self.survival_title.setStyleSheet("""
-				font-weight: bold;
-				padding: 0px;
-				""")
-			self.plotvbox.addWidget(self.survival_title, alignment=Qt.AlignCenter)
-
-			plot_buttons_hbox = QHBoxLayout()
-			plot_buttons_hbox.addWidget(QLabel(''),80, alignment=Qt.AlignLeft)
-
-			self.legend_btn = QPushButton('')
-			self.legend_btn.setIcon(icon(MDI6.text_box,color="black"))
-			self.legend_btn.setStyleSheet(self.button_select_all)
-			self.legend_btn.setToolTip('Show or hide the legend')
-			self.legend_visible = True
-			self.legend_btn.clicked.connect(self.show_hide_legend)
-			plot_buttons_hbox.addWidget(self.legend_btn, 5,alignment=Qt.AlignRight)
-
-
-			self.log_btn = QPushButton('')
-			self.log_btn.setIcon(icon(MDI6.math_log,color="black"))
-			self.log_btn.setStyleSheet(self.button_select_all)
-			self.log_btn.clicked.connect(self.switch_to_log)
-			self.log_btn.setToolTip('Enable or disable log scale')
-			plot_buttons_hbox.addWidget(self.log_btn, 5, alignment=Qt.AlignRight)
-
-
-			self.ci_btn = QPushButton('')
-			self.ci_btn.setIcon(icon(MDI6.arrow_expand_horizontal,color="blue"))
-			self.ci_btn.setStyleSheet(self.button_select_all)
-			self.ci_btn.clicked.connect(self.switch_ci)
-			self.ci_btn.setToolTip('Show or hide confidence intervals.')
-			plot_buttons_hbox.addWidget(self.ci_btn, 5, alignment=Qt.AlignRight)
-
-			self.cell_lines_btn = QPushButton('')
-			self.cell_lines_btn.setIcon(icon(MDI6.view_headline,color="black"))
-			self.cell_lines_btn.setStyleSheet(self.button_select_all)
-			self.cell_lines_btn.clicked.connect(self.switch_cell_lines)
-			self.cell_lines_btn.setToolTip('Show or hide individual cell signals.')
-			plot_buttons_hbox.addWidget(self.cell_lines_btn, 5, alignment=Qt.AlignRight)
-
-
-			self.fig, self.ax = plt.subplots(1,1,figsize=(4,3))
-			#self.setMinimumHeight(100)
-			self.survival_window = FigureCanvas(self.fig, title="Survival")
-			self.survival_window.setContentsMargins(0,0,0,0)
-			if self.df is not None:
-				self.initialize_axis()
-			plt.tight_layout()
-
-			cmap_lbl = self.cbs[-1].currentText()
-			self.cmap = getattr(mcm, cmap_lbl)
-			self.ax.set_prop_cycle('color',[self.cmap(i) for i in np.linspace(0, 1, len(self.well_indices))])
-
-
-			self.fig.set_facecolor('none')  # or 'None'
-			self.fig.canvas.setStyleSheet("background-color: transparent;")
-			self.survival_window.canvas.draw()
-
-			#self.survival_window.layout.addWidget(QLabel('WHAAAAATTT???'))
-
-			self.plot_options = [QRadioButton() for i in range(3)]
-			self.radio_labels = ['well', 'pos', 'both']
-			radio_hbox = QHBoxLayout()
-			radio_hbox.setContentsMargins(30,30,30,30)
-			self.plot_btn_group = QButtonGroup()
-			for i in range(3):
-				self.plot_options[i].setText(self.radio_labels[i])
-				#self.plot_options[i].toggled.connect(self.plot_survivals)
-				self.plot_btn_group.addButton(self.plot_options[i])
-				radio_hbox.addWidget(self.plot_options[i], 33, alignment=Qt.AlignCenter)
-			self.plot_btn_group.buttonClicked[int].connect(self.plot_survivals)
-
-			if self.position_indices is not None:
-				if len(self.well_indices)>1 and len(self.position_indices)==1:
-					self.plot_btn_group.buttons()[0].click()
-					for i in [1,2]:
-						self.plot_options[i].setEnabled(False)
-				elif len(self.well_indices)>1:
-					self.plot_btn_group.buttons()[0].click()
-				elif len(self.well_indices)==1 and len(self.position_indices)==1:
-					self.plot_btn_group.buttons()[1].click()
-					for i in [0,2]:
-						self.plot_options[i].setEnabled(False)
-			else:
-				if len(self.well_indices)>1:
-					self.plot_btn_group.buttons()[0].click()
-				elif len(self.well_indices)==1:
-					self.plot_btn_group.buttons()[2].click()
-
-
-			# elif len(self.well_indices)>1:		
-			# 	self.plot_btn_group.buttons()[0].click()
-			# else:
-			# 	self.plot_btn_group.buttons()[1].click()
-
-			# if self.position_indices is not None:
-			# 	for i in [0,2]:
-			# 		self.plot_options[i].setEnabled(False)
-
-
-			#self.plot_options[0].setChecked(True)
-			self.plotvbox.addLayout(radio_hbox)
-			self.plotvbox.addLayout(plot_buttons_hbox)
-			self.plotvbox.addWidget(self.survival_window)
-			self.class_selection_lbl = QLabel('Select class')
-			self.class_selection_lbl.setStyleSheet("""
-				font-weight: bold;
-				padding: 0px;
-				""")
-			self.plotvbox.addWidget(self.class_selection_lbl, alignment=Qt.AlignCenter)
-			class_selection_hbox = QHBoxLayout()
-			class_selection_hbox.setContentsMargins(30,30,30,30)
-			self.all_btn = QRadioButton('*')
-			self.all_btn.setChecked(True)
-			self.event_btn = QRadioButton('event')
-			self.no_event_btn = QRadioButton('no event')
-			self.class_btn_group = QButtonGroup()
-			for btn in [self.all_btn, self.event_btn, self.no_event_btn]:
-				self.class_btn_group.addButton(btn)
-
-			self.class_btn_group.buttonClicked[int].connect(self.set_class_to_plot)
-
-			class_selection_hbox.addWidget(self.all_btn, 33, alignment=Qt.AlignLeft)
-			class_selection_hbox.addWidget(self.event_btn, 33, alignment=Qt.AlignCenter)
-			class_selection_hbox.addWidget(self.no_event_btn, 33, alignment=Qt.AlignRight)
-			self.plotvbox.addLayout(class_selection_hbox)
-
-
-			self.select_pos_label = QLabel('Select positions')
-			self.select_pos_label.setStyleSheet("""
-				font-weight: bold;
-				padding: 0px;
-				""")
-			self.plotvbox.addWidget(self.select_pos_label, alignment=Qt.AlignCenter)
-
-			self.select_option = [QRadioButton() for i in range(2)]
-			self.select_label = ['name', 'spatial']
-			select_hbox = QHBoxLayout()
-			select_hbox.setContentsMargins(30,30,30,30)
-			self.select_btn_group = QButtonGroup()
-			for i in range(2):
-				self.select_option[i].setText(self.select_label[i])
-				#self.select_option[i].toggled.connect(self.switch_selection_mode)
-				self.select_btn_group.addButton(self.select_option[i])
-				select_hbox.addWidget(self.select_option[i],33, alignment=Qt.AlignCenter)
-			self.select_btn_group.buttonClicked[int].connect(self.switch_selection_mode)
-			self.plotvbox.addLayout(select_hbox)
-
-			self.look_for_metadata()
-			if self.metadata_found:
-				self.fig_scatter, self.ax_scatter = plt.subplots(1,1,figsize=(4,3))
-				self.position_scatter = FigureCanvas(self.fig_scatter)
-				self.load_coordinates()
-				self.plot_spatial_location()
-				#self.plot_positions()
-				self.ax_scatter.spines['top'].set_visible(False)
-				self.ax_scatter.spines['right'].set_visible(False)
-				self.ax_scatter.set_aspect('equal')
-				self.ax_scatter.set_xticks([])
-				self.ax_scatter.set_yticks([])
-				plt.tight_layout()
-
-				self.fig_scatter.set_facecolor('none')  # or 'None'
-				self.fig_scatter.canvas.setStyleSheet("background-color: transparent;")
-				self.plotvbox.addWidget(self.position_scatter)
-
-
-
-			self.generate_pos_selection_widget()
-
-			# if self.df is not None and len(self.ks_estimators_per_position)>0:
-			# 	self.plot_survivals()
-			self.select_btn_group.buttons()[0].click()
-			self.scroll.setVerticalScrollBarPolicy(Qt.ScrollBarAlwaysOn)
-			self.scroll.setHorizontalScrollBarPolicy(Qt.ScrollBarAlwaysOff)
-			#self.scroll.setWidgetResizable(True)
-			self.scroll.setWidget(self.survivalWidget)
-
-			self.scroll.setMinimumHeight(int(0.8*self.screen_height))
-			self.survivalWidget.setSizePolicy(QSizePolicy.Expanding, QSizePolicy.Expanding)
-			self.scroll.show()
-
+			self.interpret_pos_location()
+			self.plot_window = GenericSignalPlotWidget(parent_window=self, df=self.df, df_pos_info = self.df_pos_info, df_well_info = self.df_well_info, feature_selected=self.feature_selected, title='plot signals')
+			self.plot_window.show()
 
 	def process_signal(self):
 
-		print('you clicked!!')
 		self.FrameToMin = float(self.time_calibration_le.text().replace(',','.'))
-		print(self.FrameToMin, 'set')
+		print(f'Time calibration set to 1 frame =  {self.FrameToMin} min...')
 
 		# read instructions from combobox options
 		self.load_available_tables()
+		class_col = self.cbs[1].currentText()	
+
 		if self.df is not None:
-			self.ask_for_features()
-		else:
-			return None
 
-	def generate_pos_selection_widget(self):
-
-		self.well_names = self.df['well_name'].unique()
-		self.pos_names = self.df_pos_info['pos_name'].unique() #pd.DataFrame(self.ks_estimators_per_position)['position_name'].unique()
-		print(f'POSITION NAMES: ',self.pos_names)
-		self.usable_well_labels = []
-		for name in self.well_names:
-			for lbl in self.well_labels:
-				if name+':' in lbl:
-					self.usable_well_labels.append(lbl)
-
-		self.line_choice_widget = QWidget()
-		self.line_check_vbox = QVBoxLayout()
-		self.line_choice_widget.setLayout(self.line_check_vbox)
-		if len(self.well_indices)>1:
-			self.well_display_options = [QCheckBox(self.usable_well_labels[i]) for i in range(len(self.usable_well_labels))]
-			for i in range(len(self.well_names)):
-				self.line_check_vbox.addWidget(self.well_display_options[i], alignment=Qt.AlignLeft)
-				self.well_display_options[i].setChecked(True)
-				self.well_display_options[i].toggled.connect(self.select_survival_lines)
+			if class_col not in list(self.df.columns):
+				msgBox = QMessageBox()
+				msgBox.setIcon(QMessageBox.Warning)
+				msgBox.setText("The class of interest could not be found in the data. Abort.")
+				msgBox.setWindowTitle("Warning")
+				msgBox.setStandardButtons(QMessageBox.Ok)
+				returnValue = msgBox.exec()
+				if returnValue == QMessageBox.Ok:
+					return None
+			else:
+				self.ask_for_features()
 		else:
-			self.pos_display_options = [QCheckBox(self.pos_names[i]) for i in range(len(self.pos_names))]
-			for i in range(len(self.pos_names)):
-				self.line_check_vbox.addWidget(self.pos_display_options[i], alignment=Qt.AlignLeft)
-				self.pos_display_options[i].setChecked(True)
-				self.pos_display_options[i].toggled.connect(self.select_survival_lines)
+			return None
 
-		self.plotvbox.addWidget(self.line_choice_widget, alignment=Qt.AlignCenter)
+		#self.plotvbox.addWidget(self.line_choice_widget, alignment=Qt.AlignCenter)
 
 	def load_available_tables(self):
 
@@ -579,7 +334,8 @@ class ConfigSignalPlot(QWidget, Styles):
 		max_time = int(self.df.FRAME.max()) + 1
 		class_col = self.cbs[1].currentText()
 		time_col = self.cbs[2].currentText()
-
+		if self.abs_time_checkbox.isChecked():
+			time_col = self.frame_slider.value()
 
 		for block,movie_group in self.df.groupby(['well','position']):
 
@@ -587,28 +343,13 @@ class ConfigSignalPlot(QWidget, Styles):
 			well_signal_event, well_std_event, timeline_event, matrix_event = mean_signal(movie_group, self.feature_selected, class_col, time_col=time_col, class_value=[0], return_matrix=True, forced_max_duration=max_time)		
 			well_signal_no_event, well_std_no_event, timeline_no_event, matrix_no_event = mean_signal(movie_group, self.feature_selected, class_col, time_col=time_col, class_value=[1], return_matrix=True, forced_max_duration=max_time)
 			self.mean_plots_timeline = timeline_all
-			# self.df_pos_info.loc[self.df_pos_info['pos_path']==block[1],'signal'] = [{'mean_all': well_signal_mean, 'std_all': well_std_mean,'matrix_all': matrix_all,'mean_event': well_signal_event, 'std_event': well_std_event,
-			# 																		'matrix_event': matrix_event,'mean_no_event': well_signal_no_event, 'std_no_event': well_std_no_event, 'matrix_no_event': matrix_no_event, 'timeline':  self.mean_plots_timeline}]
-
-			if self.checkBox_feature.isChecked():
-				second_well_signal_mean, second_well_std_mean, second_timeline_all, second_matrix_all = mean_signal(movie_group, self.second_feature_selected, class_col, time_col=time_col, class_value=[0,1], return_matrix=True, forced_max_duration=max_time)
-				second_well_signal_event, second_well_std_event, second_timeline_event, second_matrix_event = mean_signal(movie_group, self.second_feature_selected, class_col, time_col=time_col, class_value=[0], return_matrix=True, forced_max_duration=max_time)
-				second_well_signal_no_event, second_well_std_no_event, second_timeline_no_event, second_matrix_no_event = mean_signal(movie_group, self.second_feature_selected, class_col, time_col=time_col, class_value=[1], return_matrix=True, forced_max_duration=max_time)
-				self.second_mean_plots_timeline = second_timeline_all
-				self.df_pos_info.loc[self.df_pos_info['pos_path'] == block[1], 'signal'] = [
-					{'mean_all': well_signal_mean, 'std_all': well_std_mean,'matrix_all': matrix_all,'mean_event': well_signal_event, 'std_event': well_std_event,
-					'matrix_event': matrix_event,'mean_no_event': well_signal_no_event, 'std_no_event': well_std_no_event, 'matrix_no_event': matrix_no_event, 'timeline':  self.mean_plots_timeline,'second_mean_all': second_well_signal_mean, 'second_std_all': second_well_std_mean, 'second_matrix_all': second_matrix_all,
-					'second_mean_event': second_well_signal_event, 'second_std_event': second_well_std_event,
-					'second_matrix_event': second_matrix_event, 'second_mean_no_event': second_well_signal_no_event,
-					'second_std_no_event': second_well_std_no_event, 'second_matrix_no_event': second_matrix_no_event,
-					'second_timeline': self.second_mean_plots_timeline}]
-			else:
-				self.df_pos_info.loc[self.df_pos_info['pos_path'] == block[1], 'signal'] = [
-					{'mean_all': well_signal_mean, 'std_all': well_std_mean, 'matrix_all': matrix_all,
-					 'mean_event': well_signal_event, 'std_event': well_std_event,
-					 'matrix_event': matrix_event, 'mean_no_event': well_signal_no_event,
-					 'std_no_event': well_std_no_event, 'matrix_no_event': matrix_no_event,
-					 'timeline': self.mean_plots_timeline}]
+
+			self.df_pos_info.loc[self.df_pos_info['pos_path'] == block[1], 'signal'] = [
+				{'mean_all': well_signal_mean, 'std_all': well_std_mean, 'matrix_all': matrix_all,
+				 'mean_event': well_signal_event, 'std_event': well_std_event,
+				 'matrix_event': matrix_event, 'mean_no_event': well_signal_no_event,
+				 'std_no_event': well_std_no_event, 'matrix_no_event': matrix_no_event,
+				 'timeline': self.mean_plots_timeline}]
 
 		# Per well
 		for well,well_group in self.df.groupby('well'):
@@ -616,24 +357,12 @@ class ConfigSignalPlot(QWidget, Styles):
 			well_signal_mean, well_std_mean, timeline_all, matrix_all = mean_signal(well_group, self.feature_selected, class_col, time_col=time_col, class_value=[0,1], return_matrix=True, forced_max_duration=max_time)
 			well_signal_event, well_std_event, timeline_event, matrix_event = mean_signal(well_group, self.feature_selected, class_col, time_col=time_col, class_value=[0], return_matrix=True, forced_max_duration=max_time)			
 			well_signal_no_event, well_std_no_event, timeline_no_event, matrix_no_event = mean_signal(well_group, self.feature_selected, class_col, time_col=time_col, class_value=[1], return_matrix=True, forced_max_duration=max_time)
-			if self.checkBox_feature.isChecked():
-				second_well_signal_mean, second_well_std_mean, second_timeline_all, second_matrix_all = mean_signal(movie_group, self.second_feature_selected, class_col, time_col=time_col, class_value=[0,1], return_matrix=True, forced_max_duration=max_time)
-				second_well_signal_event, second_well_std_event, second_timeline_event, second_matrix_event = mean_signal(movie_group, self.second_feature_selected, class_col, time_col=time_col, class_value=[0], return_matrix=True, forced_max_duration=max_time)
-				second_well_signal_no_event, second_well_std_no_event, second_timeline_no_event, second_matrix_no_event = mean_signal(movie_group, self.second_feature_selected, class_col, time_col=time_col, class_value=[1], return_matrix=True, forced_max_duration=max_time)
-				self.df_well_info.loc[self.df_well_info['well_path'] == well, 'signal'] = [
-					{'mean_all': well_signal_mean, 'std_all': well_std_mean,'matrix_all': matrix_all,'mean_event': well_signal_event, 'std_event': well_std_event,
-																					'matrix_event': matrix_event,'mean_no_event': well_signal_no_event, 'std_no_event': well_std_no_event, 'matrix_no_event': matrix_no_event, 'timeline':  self.mean_plots_timeline,'second_mean_all': second_well_signal_mean, 'second_std_all': second_well_std_mean, 'second_matrix_all': second_matrix_all,
-					'second_mean_event': second_well_signal_event, 'second_std_event': second_well_std_event,
-					 'second_matrix_event': second_matrix_event, 'second_mean_no_event': second_well_signal_no_event,
-					 'second_std_no_event': second_well_std_no_event, 'second_matrix_no_event': second_matrix_no_event,
-					 'second_timeline': self.second_mean_plots_timeline}]
-			else:
-				self.df_well_info.loc[self.df_well_info['well_path']==well,'signal'] = [{'mean_all': well_signal_mean, 'std_all': well_std_mean,'matrix_all': matrix_all,'mean_event': well_signal_event, 'std_event': well_std_event,
+			
+			self.df_well_info.loc[self.df_well_info['well_path']==well,'signal'] = [{'mean_all': well_signal_mean, 'std_all': well_std_mean,'matrix_all': matrix_all,'mean_event': well_signal_event, 'std_event': well_std_event,
 																					'matrix_event': matrix_event,'mean_no_event': well_signal_no_event, 'std_no_event': well_std_no_event, 'matrix_no_event': matrix_no_event, 'timeline':  self.mean_plots_timeline}]
 
 		self.df_pos_info.loc[:,'select'] = True
 		self.df_well_info.loc[:,'select'] = True
-		self.df_pos_info.to_csv('pos_info.csv')
 
 
 	def generate_synchronized_matrix(self, well_group, feature_selected, cclass, max_time):
@@ -698,337 +427,6 @@ class ConfigSignalPlot(QWidget, Styles):
 
 		return mean_line, std_line
 
-	def initialize_axis(self):
-
-		self.ax.clear()
-		self.ax.plot([],[])
-		self.ax.spines['top'].set_visible(False)
-		self.ax.spines['right'].set_visible(False)
-		self.ax.set_ylim(self.df[self.feature_selected].min(),self.df[self.feature_selected].max())
-		self.ax.set_xlim(-(self.df['FRAME'].max()+2)*self.FrameToMin,(self.df['FRAME'].max()+2)*self.FrameToMin)
-		self.ax.set_xlabel('time [min]')
-		self.ax.set_ylabel(self.feature_selected)
-		if self.checkBox_feature.isChecked():
-			if self.ax2:
-				self.ax2.clear()
-				self.ax2.get_yaxis().set_visible(False)
-			self.ax2 = self.ax.twinx()
-			#self.ax2.plot([], color='red', label='Feature 2')
-			self.ax2.set_ylim(self.df[self.second_feature_selected].min(), self.df[self.second_feature_selected].max())
-			self.ax2.set_ylabel(self.second_feature_selected, color='red')
-			self.ax2.get_yaxis().set_visible(True)
-			#self.ax2.tick_params('y', colors='red')
-
-	def plot_survivals(self, id):
-
-		for i in range(3):
-			if self.plot_options[i].isChecked():
-				self.plot_mode = self.radio_labels[i]
-
-		if self.target_class==[0,1]:
-			mean_signal = 'mean_all'
-			std_signal = 'std_all'
-			matrix = 'matrix_all'
-			if self.checkBox_feature.isChecked():
-				second_mean_signal='second_mean_all'
-				second_std_signal='second_std_all'
-				second_matrix='second_matrix_all'
-		elif self.target_class==[0]:
-			mean_signal = 'mean_event'
-			std_signal = 'std_event'
-			matrix = 'matrix_event'
-			if self.checkBox_feature.isChecked():
-				second_mean_signal='second_mean_event'
-				second_std_signal='second_std_event'
-				second_matrix='second_matrix_event'
-		else:
-			mean_signal = 'mean_no_event'
-			std_signal = 'std_no_event'
-			matrix = 'matrix_no_event'
-			if self.checkBox_feature.isChecked():
-				second_mean_signal='second_mean_no_event'
-				second_std_signal='second_std_no_event'
-				second_matrix='second_matrix_no_event'
-
-
-		colors = np.array([self.cmap(i / len(self.df_pos_info)) for i in range(len(self.df_pos_info))])
-		if self.checkBox_feature.isChecked():
-			second_colors = self.cmap(np.linspace(0.5, 1.5, len(self.df_pos_info)))
-		well_color = [self.cmap(i / len(self.df_well_info)) for i in range(len(self.df_well_info))]
-
-		if self.plot_mode=='pos':
-			self.initialize_axis()
-			lines = self.df_pos_info.loc[self.df_pos_info['select'],'signal'].values
-			pos_labels = self.df_pos_info.loc[self.df_pos_info['select'],'pos_name'].values
-			pos_indices = self.df_pos_info.loc[self.df_pos_info['select'],'pos_index'].values
-			well_index = self.df_pos_info.loc[self.df_pos_info['select'],'well_index'].values
-			for i in range(len(lines)):
-				if len(self.well_indices)<=1:
-					self.plot_line(lines[i], colors[pos_indices[i]], pos_labels[i], mean_signal, std_signal=std_signal, ci_option=self.show_ci, cell_lines_option=self.show_cell_lines, matrix=matrix)
-					if self.checkBox_feature.isChecked():
-						self.plot_line(lines[i],second_colors[pos_indices[i]],pos_labels[i],second_mean_signal,std_signal=second_std_signal,ci_option=self.show_ci, cell_lines_option=self.show_cell_lines, matrix=second_matrix)
-
-
-				else:
-					self.plot_line(lines[i], well_color[well_index[i]], pos_labels[i], mean_signal, std_signal=std_signal, ci_option=self.show_ci, cell_lines_option=self.show_cell_lines, matrix=matrix)
-					if self.checkBox_feature.isChecked():
-						self.plot_line(lines[i],second_colors[pos_indices[i]],pos_labels[i],second_mean_signal,std_signal=second_std_signal,ci_option=self.show_ci, cell_lines_option=self.show_cell_lines, matrix=second_matrix)
-
-		elif self.plot_mode=='well':
-			self.initialize_axis()
-			lines = self.df_well_info.loc[self.df_well_info['select'],'signal'].values	
-			well_index = self.df_well_info.loc[self.df_well_info['select'],'well_index'].values
-			well_labels = self.df_well_info.loc[self.df_well_info['select'],'well_name'].values
-			for i in range(len(lines)):	
-				if len(self.well_indices)<=1:
-					self.plot_line(lines[i], 'k', well_labels[i], mean_signal, std_signal=std_signal, ci_option=self.show_ci, cell_lines_option=self.show_cell_lines,matrix=matrix)
-					if self.checkBox_feature.isChecked():
-						self.plot_line(lines[i], 'k', well_labels[i], second_mean_signal, std_signal=second_std_signal,
-									   ci_option=self.show_ci, cell_lines_option=self.show_cell_lines, matrix=second_matrix)
-				else:
-					self.plot_line(lines[i], well_color[well_index[i]], well_labels[i], mean_signal, std_signal=std_signal, ci_option=self.show_ci, cell_lines_option=self.show_cell_lines,matrix=matrix)
-					if self.checkBox_feature.isChecked():
-						self.plot_line(lines[i], well_color[well_index[i]], well_labels[i], second_mean_signal, std_signal=second_std_signal, ci_option=self.show_ci, cell_lines_option=self.show_cell_lines,matrix=second_matrix)
-
-		elif self.plot_mode=='both':
-			self.initialize_axis()
-			lines_pos = self.df_pos_info.loc[self.df_pos_info['select'],'signal'].values
-			lines_well = self.df_well_info.loc[self.df_well_info['select'],'signal'].values	
-
-			pos_indices = self.df_pos_info.loc[self.df_pos_info['select'],'pos_index'].values
-			well_index_pos = self.df_pos_info.loc[self.df_pos_info['select'],'well_index'].values
-			well_index = self.df_well_info.loc[self.df_well_info['select'],'well_index'].values
-			well_labels = self.df_well_info.loc[self.df_well_info['select'],'well_name'].values
-			pos_labels = self.df_pos_info.loc[self.df_pos_info['select'],'pos_name'].values
-
-			for i in range(len(lines_pos)):
-				if len(self.well_indices)<=1:
-					self.plot_line(lines_pos[i], colors[pos_indices[i]], pos_labels[i], mean_signal, std_signal=std_signal, ci_option=self.show_ci, cell_lines_option=self.show_cell_lines,matrix=matrix)
-					if self.checkBox_feature.isChecked():
-						self.plot_line(lines_pos[i], second_colors[pos_indices[i]], pos_labels[i], second_mean_signal,
-									   std_signal=second_std_signal, ci_option=self.show_ci,
-									   cell_lines_option=self.show_cell_lines, matrix=second_matrix)
-
-				else:
-					self.plot_line(lines_pos[i], well_color[well_index_pos[i]], None, mean_signal, std_signal=std_signal, ci_option=False)
-					if self.checkBox_feature.isChecked():
-						self.plot_line(lines_pos[i], well_color[well_index_pos[i]], None, second_mean_signal,
-									   std_signal=second_std_signal, ci_option=False)
-
-			for i in range(len(lines_well)):
-				if len(self.well_indices)<=1:
-					self.plot_line(lines_well[i], 'k', 'pool', mean_signal, std_signal=std_signal, ci_option=False)
-				else:
-					self.plot_line(lines_well[i], well_color[well_index[i]], well_labels[i], mean_signal, std_signal=std_signal, ci_option=False)
-		self.survival_window.setMinimumHeight(int(0.5*self.screen_height))
-		self.survival_window.setMinimumWidth(int(0.8 * self.survivalWidget.width()))
-		self.survival_window.setSizePolicy(QSizePolicy.Expanding, QSizePolicy.Expanding)
-
-		self.survival_window.canvas.draw()
-
-	def plot_line(self, line, color, label, mean_signal, ci_option=True, cell_lines_option=False, alpha_ci=0.5, alpha_cell_lines=0.5, std_signal=None, matrix=None):
-		
-
-		try:
-			if 'second' in str(mean_signal):
-				self.ax2.plot(line['timeline'] * self.FrameToMin, line[mean_signal], color=color, label=label)
-			else:
-				self.ax.plot(line['timeline']*self.FrameToMin, line[mean_signal], color=color, label=label)
-			if ci_option and std_signal is not None:
-				if 'second' in str(mean_signal):
-					self.ax2.fill_between(line['timeline'] * self.FrameToMin,
-										 [a - b for a, b in zip(line[mean_signal], line[std_signal])],
-										 [a + b for a, b in zip(line[mean_signal], line[std_signal])],
-										 color=color,
-										 alpha=alpha_ci,
-										 )
-				else:
-					self.ax.fill_between(line['timeline']*self.FrameToMin,
-									[a-b for a,b in zip(line[mean_signal], line[std_signal])],
-									[a+b for a,b in zip(line[mean_signal], line[std_signal])],
-									color=color,
-									alpha=alpha_ci,
-									)
-			if cell_lines_option and matrix is not None:
-				print(mean_signal)
-				mat = line[matrix]
-				if 'second' in str(mean_signal):
-					for i in range(mat.shape[0]):
-						self.ax2.plot(line['timeline'] * self.FrameToMin, mat[i, :], color=color, alpha=alpha_cell_lines)
-				else:
-					for i in range(mat.shape[0]):
-						self.ax.plot(line['timeline']*self.FrameToMin, mat[i,:], color=color, alpha=alpha_cell_lines)
-
-		except Exception as e:
-			print(f'Exception {e}')
-
-	def switch_to_log(self):
-
-		"""
-		Switch threshold histogram to log scale. Auto adjust.
-		"""
-
-		if self.ax.get_yscale()=='linear':
-			self.ax.set_yscale('log')
-			#self.ax.set_ylim(0.01,1.05)
-		else:
-			self.ax.set_yscale('linear')
-			#self.ax.set_ylim(0.01,1.05)
-
-		#self.ax.autoscale()
-		self.survival_window.canvas.draw_idle()
-
-	def show_hide_legend(self):
-		if self.legend_visible:
-			self.ax.legend().set_visible(False)
-			self.legend_visible = False
-			self.legend_btn.setIcon(icon(MDI6.text_box_outline,color="black"))
-		else:
-			self.ax.legend().set_visible(True)
-			self.legend_visible = True
-			self.legend_btn.setIcon(icon(MDI6.text_box,color="black"))
-
-		self.survival_window.canvas.draw_idle()		
-
-	def look_for_metadata(self):
-
-		self.metadata_found = False
-		self.metafiles = glob(self.exp_dir+os.sep.join([f'W*','*','movie','*metadata.txt'])) \
-					+ glob(self.exp_dir+os.sep.join([f'W*','*','*metadata.txt'])) \
-					+ glob(self.exp_dir+os.sep.join([f'W*','*metadata.txt'])) \
-					+ glob(self.exp_dir+'*metadata.txt')
-		print(f'Found {len(self.metafiles)} metadata files...')
-		if len(self.metafiles)>0:
-			self.metadata_found = True
-
-	def switch_selection_mode(self, id):
-		print(f'button {id} was clicked')
-		for i in range(2):
-			if self.select_option[i].isChecked():
-				self.selection_mode = self.select_label[i]
-		if self.selection_mode=='name':
-			if len(self.metafiles)>0:
-				self.position_scatter.hide()
-			self.line_choice_widget.show()
-		else:
-			if len(self.metafiles)>0:
-				self.position_scatter.show()
-			self.line_choice_widget.hide()
-
-
-	def load_coordinates(self):
-
-		"""
-		Read metadata and try to extract position coordinates
-		"""
-
-		self.no_meta = False
-		try:
-			with open(self.metafiles[0], 'r') as f:
-				data = json.load(f)
-				positions = data['Summary']['InitialPositionList']
-		except Exception as e:
-			print(f'Trouble loading metadata: error {e}...')
-			return None
-
-		for k in range(len(positions)):
-			pos_label = positions[k]['Label']
-			try:
-				coords = positions[k]['DeviceCoordinatesUm']['XYStage']
-			except:
-				try:
-					coords = positions[k]['DeviceCoordinatesUm']['PIXYStage']
-				except:
-					self.no_meta = True
-
-			if not self.no_meta:
-				files = self.df_pos_info['stack_path'].values
-				pos_loc = [pos_label in f for f in files]
-				self.df_pos_info.loc[pos_loc, 'x'] = coords[0]
-				self.df_pos_info.loc[pos_loc, 'y'] = coords[1]
-				self.df_pos_info.loc[pos_loc, 'metadata_tag'] = pos_label
-
-
-	def update_annot(self, ind):
-		
-		pos = self.sc.get_offsets()[ind["ind"][0]]
-		self.annot.xy = pos
-		text = self.scat_labels[ind["ind"][0]]
-		self.annot.set_text(text)
-		self.annot.get_bbox_patch().set_facecolor('k')
-		self.annot.get_bbox_patch().set_alpha(0.4)
-
-	def hover(self, event):
-		vis = self.annot.get_visible()
-		if event.inaxes == self.ax_scatter:
-			cont, ind = self.sc.contains(event)
-			if cont:
-				self.update_annot(ind)
-				self.annot.set_visible(True)
-				self.fig_scatter.canvas.draw_idle()
-			else:
-				if vis:
-					self.annot.set_visible(False)
-					self.fig_scatter.canvas.draw_idle()
-
-	def unselect_position(self, event):
-		
-		ind = event.ind # index of selected position
-		well_idx = self.df_pos_info.iloc[ind]['well_index'].values[0]
-		selectedPos = self.df_pos_info.iloc[ind]['pos_path'].values[0]
-		currentSelState = self.df_pos_info.iloc[ind]['select'].values[0]
-		if self.plot_options[0].isChecked() or self.plot_options[2].isChecked():
-			self.df_pos_info.loc[self.df_pos_info['well_index']==well_idx,'select'] = not currentSelState
-			self.df_well_info.loc[self.df_well_info['well_index']==well_idx, 'select'] = not currentSelState
-			if len(self.well_indices)>1:
-				self.well_display_options[well_idx].setChecked(not currentSelState)
-			else:
-				for p in self.pos_display_options:
-					p.setChecked(not currentSelState)
-		else:
-			self.df_pos_info.loc[self.df_pos_info['pos_path']==selectedPos,'select'] = not currentSelState
-			if len(self.well_indices)<=1:
-				self.pos_display_options[ind[0]].setChecked(not currentSelState)
-
-		self.sc.set_color(self.select_color(self.df_pos_info["select"].values))
-		self.position_scatter.canvas.draw_idle()
-		self.plot_survivals(0)
-
-	def select_survival_lines(self):
-		
-		if len(self.well_indices)>1:
-			for i in range(len(self.well_display_options)):
-				self.df_well_info.loc[self.df_well_info['well_index']==i,'select'] = self.well_display_options[i].isChecked()
-				self.df_pos_info.loc[self.df_pos_info['well_index']==i,'select'] = self.well_display_options[i].isChecked()
-		else:
-			for i in range(len(self.pos_display_options)):
-				self.df_pos_info.loc[self.df_pos_info['pos_index']==i,'select'] = self.pos_display_options[i].isChecked()
-
-		if len(self.metafiles)>0:
-			self.sc.set_color(self.select_color(self.df_pos_info["select"].values))
-			self.position_scatter.canvas.draw_idle()
-		self.plot_survivals(0)				
-
-
-	def select_color(self, selection):
-		colors = [tab10(0) if s else tab10(0.1) for s in selection]
-		return colors
-
-	def plot_spatial_location(self):
-
-		try:
-			self.sc = self.ax_scatter.scatter(self.df_pos_info["x"].values, self.df_pos_info["y"].values, picker=True, pickradius=1, color=self.select_color(self.df_pos_info["select"].values))
-			self.scat_labels = self.df_pos_info['metadata_tag'].values
-			self.ax_scatter.invert_xaxis()
-			self.annot = self.ax_scatter.annotate("", xy=(0,0), xytext=(10,10),textcoords="offset points",
-								bbox=dict(boxstyle="round", fc="w"),
-								arrowprops=dict(arrowstyle="->"))
-			self.annot.set_visible(False)
-			self.fig_scatter.canvas.mpl_connect("motion_notify_event", self.hover)
-			self.fig_scatter.canvas.mpl_connect("pick_event", self.unselect_position)
-		except Exception as e:
-			pass
-
 
 	def switch_ref_time_mode(self):
 		if self.abs_time_checkbox.isChecked():
@@ -1036,36 +434,4 @@ class ConfigSignalPlot(QWidget, Styles):
 			self.cbs[-2].setEnabled(False)
 		else:
 			self.frame_slider.setEnabled(False)
-			self.cbs[-2].setEnabled(True)
-
-	def switch_ci(self):
-		
-		if self.show_ci:
-			self.ci_btn.setIcon(icon(MDI6.arrow_expand_horizontal,color="black"))
-		else:
-			self.ci_btn.setIcon(icon(MDI6.arrow_expand_horizontal,color="blue"))
-		self.show_ci = not self.show_ci
-		self.plot_survivals(0)
-
-	def switch_cell_lines(self):
-		
-		if self.show_cell_lines:
-			self.cell_lines_btn.setIcon(icon(MDI6.view_headline,color="black"))
-		else:
-			self.cell_lines_btn.setIcon(icon(MDI6.view_headline,color="blue"))
-		self.show_cell_lines = not self.show_cell_lines
-		self.plot_survivals(0)
-
-
-	def set_class_to_plot(self):
-
-		if self.all_btn.isChecked():
-			self.target_class=[0,1]
-		elif self.event_btn.isChecked():
-			self.target_class = [0]
-		else:
-			self.target_class = [1]
-
-		self.plot_survivals(0)
-
-
+			self.cbs[-2].setEnabled(True)