{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "708c9745",
   "metadata": {},
   "source": [
    "### Imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "53b10294",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "from pathlib import Path\n",
    "import sys\n",
    "import os\n",
    "\n",
    "base_dir = os.path.abspath(os.path.join(os.getcwd(), \"..\"))\n",
    "sys.path.append(base_dir)\n",
    "print(base_dir)\n",
    "\n",
    "from Fahrsimulator_MSY2526_AI.model_training.tools import evaluation_tools, scaler, mad_outlier_removal\n",
    "from sklearn.preprocessing import StandardScaler, MinMaxScaler\n",
    "from sklearn.svm import OneClassSVM\n",
    "from sklearn.model_selection import GridSearchCV, KFold, ParameterGrid, train_test_split\n",
    "import matplotlib.pyplot as plt\n",
    "import tensorflow as tf\n",
    "import pickle\n",
    "from sklearn.metrics import (roc_auc_score, accuracy_score, precision_score, \n",
    "                                recall_score, f1_score, confusion_matrix, classification_report)  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "68101229",
   "metadata": {},
   "source": [
    "### load Dataset"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "24a765e8",
   "metadata": {},
   "outputs": [],
   "source": [
    "dataset_path = Path(r\"/home/jovyan/data-paulusjafahrsimulator-gpu/first_AU_dataset/output_windowed.parquet\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "471001b0",
   "metadata": {},
   "outputs": [],
   "source": [
    "df = pd.read_parquet(path=dataset_path)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0fdecdaa",
   "metadata": {},
   "source": [
    "### Load Performance data and Subject Split"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "692d1b47",
   "metadata": {},
   "outputs": [],
   "source": [
    "performance_path =  Path(r\"/home/jovyan/data-paulusjafahrsimulator-gpu/subject_performance/3new_au_performance.csv\")\n",
    "performance_df = pd.read_csv(performance_path)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ea617e3f",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Subject IDs aus dem Haupt-Dataset nehmen\n",
    "subjects_from_df = df[\"subjectID\"].unique()\n",
    "\n",
    "# Performance-Subset nur für vorhandene Subjects\n",
    "perf_filtered = performance_df[\n",
    "    performance_df[\"subjectID\"].isin(subjects_from_df)\n",
    "][[\"subjectID\", \"overall_score\"]]\n",
    "\n",
    "# Merge: nur Subjects, die sowohl im df als auch im Performance-CSV vorkommen\n",
    "merged = (\n",
    "    pd.DataFrame({\"subjectID\": subjects_from_df})\n",
    "    .merge(perf_filtered, on=\"subjectID\", how=\"inner\")\n",
    ")\n",
    "\n",
    "# Sicherstellen, dass keine Scores fehlen\n",
    "if merged[\"overall_score\"].isna().any():\n",
    "    raise ValueError(\"Es fehlen Score-Werte für manche Subjects.\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ae43df8d",
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_sorted = merged.sort_values(\"overall_score\", ascending=False).reset_index(drop=True)\n",
    "\n",
    "scores = merged_sorted[\"overall_score\"].values\n",
    "n_total = len(merged_sorted)\n",
    "n_small = n_total // 3\n",
    "n_large = n_total - n_small\n",
    "\n",
    "# Schritt 1: zufällige Start-Aufteilung\n",
    "idx = np.arange(n_total)\n",
    "np.random.shuffle(idx)\n",
    "\n",
    "small_idx = idx[:n_small]\n",
    "large_idx = idx[n_small:]\n",
    "\n",
    "def score_diff(small_idx, large_idx):\n",
    "    return abs(scores[small_idx].mean() - scores[large_idx].mean())\n",
    "\n",
    "diff = score_diff(small_idx, large_idx)\n",
    "threshold = 0.01\n",
    "max_iter = 100\n",
    "count = 0\n",
    "\n",
    "# Schritt 2: random swaps bis Differenz klein genug\n",
    "while diff > threshold and count < max_iter:\n",
    "    # Zwei zufällige Elemente auswählen\n",
    "    si = np.random.choice(small_idx)\n",
    "    li = np.random.choice(large_idx)\n",
    "    \n",
    "    # Tausch durchführen\n",
    "    new_small_idx = small_idx.copy()\n",
    "    new_large_idx = large_idx.copy()\n",
    "    \n",
    "    new_small_idx[new_small_idx == si] = li\n",
    "    new_large_idx[new_large_idx == li] = si\n",
    "\n",
    "    # neue Differenz berechnen\n",
    "    new_diff = score_diff(new_small_idx, new_large_idx)\n",
    "\n",
    "    # Swap akzeptieren, wenn es besser wird\n",
    "    if new_diff < diff:\n",
    "        small_idx = new_small_idx\n",
    "        large_idx = new_large_idx\n",
    "        diff = new_diff\n",
    "\n",
    "    count += 1\n",
    "\n",
    "# Finalgruppen\n",
    "group_small = merged_sorted.loc[small_idx].reset_index(drop=True)\n",
    "group_large = merged_sorted.loc[large_idx].reset_index(drop=True)\n",
    "\n",
    "print(\"Finale Score-Differenz:\", diff)\n",
    "print(\"Größe Gruppe 1:\", len(group_small))\n",
    "print(\"Größe Gruppe 2:\", len(group_large))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9d1b414e",
   "metadata": {},
   "outputs": [],
   "source": [
    "group_large['overall_score'].mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "fa71f9a5",
   "metadata": {},
   "outputs": [],
   "source": [
    "group_small['overall_score'].mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "79ecb4a2",
   "metadata": {},
   "outputs": [],
   "source": [
    "training_subjects = group_large['subjectID'].values\n",
    "test_subjects = group_small['subjectID'].values\n",
    "print(training_subjects)\n",
    "print(test_subjects)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4353f87c",
   "metadata": {},
   "source": [
    "### Data cleaning with mad"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "76610052",
   "metadata": {},
   "outputs": [],
   "source": [
    "# SET\n",
    "threshold_mad = 5\n",
    "column_praefix ='AU'\n",
    "\n",
    "au_columns = [col for col in df.columns if col.startswith(column_praefix)]\n",
    "cleaned_df = mad_outlier_removal.mad_outlier_removal(df,columns=au_columns, threshold=threshold_mad)\n",
    "print(cleaned_df.shape)\n",
    "print(df.shape)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9a6c1732",
   "metadata": {},
   "source": [
    "#### TO DO\n",
    "    * pipeline aus Autoencoder und SVM\n",
    "    * group k fold\n",
    "    * AE überpüfen, loss dokumentieren"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "877309d9",
   "metadata": {},
   "outputs": [],
   "source": [
    "### Variational Autoencoder with Classifier Head\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import tensorflow as tf\n",
    "from tensorflow import keras\n",
    "from tensorflow.keras import layers, Model\n",
    "from sklearn.model_selection import GroupKFold\n",
    "from sklearn.preprocessing import StandardScaler\n",
    "from sklearn.metrics import (\n",
    "    accuracy_score, precision_score, recall_score, f1_score, \n",
    "    roc_auc_score, confusion_matrix, classification_report\n",
    ")\n",
    "import matplotlib.pyplot as plt\n",
    "from collections import defaultdict\n",
    "\n",
    "# ============================================================================\n",
    "# 1. CREATE LABELS\n",
    "# ============================================================================\n",
    "\n",
    "# Low workload: baseline + n-back level 1,4\n",
    "low_all = cleaned_df[\n",
    "    ((cleaned_df[\"PHASE\"] == \"baseline\") |\n",
    "     ((cleaned_df[\"STUDY\"] == \"n-back\") & (cleaned_df[\"PHASE\"] != \"baseline\") & (cleaned_df[\"LEVEL\"].isin([1,4]))))\n",
    "].copy()\n",
    "low_all['label'] = 0\n",
    "print(f\"Low workload samples: {low_all.shape[0]}\")\n",
    "\n",
    "# High workload n-back: level 2,3,5,6\n",
    "high_nback = cleaned_df[\n",
    "    (cleaned_df[\"STUDY\"]==\"n-back\") &\n",
    "    (cleaned_df[\"LEVEL\"].isin([2, 3, 5, 6])) &\n",
    "    (cleaned_df[\"PHASE\"].isin([\"train\", \"test\"]))\n",
    "].copy()\n",
    "high_nback['label'] = 1\n",
    "print(f\"High n-back samples: {high_nback.shape[0]}\")\n",
    "\n",
    "# High workload k-drive\n",
    "high_kdrive = cleaned_df[\n",
    "    (cleaned_df[\"STUDY\"] == \"k-drive\") & (cleaned_df[\"PHASE\"] != \"baseline\")\n",
    "].copy()\n",
    "high_kdrive['label'] = 1\n",
    "print(f\"High k-drive samples: {high_kdrive.shape[0]}\")\n",
    "\n",
    "# Combine all high workload\n",
    "high_all = pd.concat([high_nback, high_kdrive])\n",
    "print(f\"Total high workload samples: {high_all.shape[0]}\")\n",
    "\n",
    "# Complete labeled dataset\n",
    "labeled_df = pd.concat([low_all, high_all]).reset_index(drop=True)\n",
    "print(f\"\\nTotal labeled samples: {labeled_df.shape[0]}\")\n",
    "print(f\"Class distribution:\\n{labeled_df['label'].value_counts()}\")\n",
    "\n",
    "# ============================================================================\n",
    "# 2. TRAIN/TEST SPLIT BY SUBJECTS\n",
    "# ============================================================================\n",
    "\n",
    "train_df = labeled_df[labeled_df['subjectID'].isin(training_subjects)].copy()\n",
    "test_df = labeled_df[labeled_df['subjectID'].isin(test_subjects)].copy()\n",
    "\n",
    "print(f\"\\nTraining subjects: {training_subjects}\")\n",
    "print(f\"Test subjects: {test_subjects}\")\n",
    "print(f\"Train samples: {train_df.shape[0]}, Test samples: {test_df.shape[0]}\")\n",
    "\n",
    "# Extract features and labels\n",
    "au_columns = [col for col in labeled_df.columns if col.startswith('AU')]\n",
    "print(f\"\\nUsing {len(au_columns)} AU features: {au_columns}\")\n",
    "\n",
    "X_train = train_df[au_columns].values\n",
    "y_train = train_df['label'].values\n",
    "groups_train = train_df['subjectID'].values\n",
    "\n",
    "X_test = test_df[au_columns].values\n",
    "y_test = test_df['label'].values\n",
    "\n",
    "# Normalize features\n",
    "scaler = StandardScaler()\n",
    "X_train_scaled = scaler.fit_transform(X_train)\n",
    "X_test_scaled = scaler.transform(X_test)\n",
    "\n",
    "print(f\"\\nTrain class distribution: {np.bincount(y_train)}\")\n",
    "print(f\"Test class distribution: {np.bincount(y_test)}\")\n",
    "\n",
    "# ============================================================================\n",
    "# 3. VAE WITH CLASSIFIER HEAD MODEL\n",
    "# ============================================================================\n",
    "\n",
    "class Sampling(layers.Layer):\n",
    "    \"\"\"Reparameterization trick for VAE\"\"\"\n",
    "    def call(self, inputs):\n",
    "        z_mean, z_log_var = inputs\n",
    "        batch = tf.shape(z_mean)[0]\n",
    "        dim = tf.shape(z_mean)[1]\n",
    "        epsilon = tf.random.normal(shape=(batch, dim))\n",
    "        return z_mean + tf.exp(0.5 * z_log_var) * epsilon\n",
    "\n",
    "def build_vae_classifier(input_dim, latent_dim, encoder_dims=[32, 16], \n",
    "                         decoder_dims=[16, 32], classifier_dims=[16]):\n",
    "    \"\"\"\n",
    "    Build VAE with classifier head\n",
    "    \n",
    "    Args:\n",
    "        input_dim: Number of input features (20 AUs)\n",
    "        latent_dim: Dimension of latent space (2-5)\n",
    "        encoder_dims: Hidden layer sizes for encoder\n",
    "        decoder_dims: Hidden layer sizes for decoder\n",
    "        classifier_dims: Hidden layer sizes for classifier\n",
    "    \"\"\"\n",
    "    \n",
    "    # ---- ENCODER ----\n",
    "    encoder_inputs = keras.Input(shape=(input_dim,), name='encoder_input')\n",
    "    x = encoder_inputs\n",
    "    \n",
    "    for i, dim in enumerate(encoder_dims):\n",
    "        x = layers.Dense(dim, activation='relu', name=f'encoder_dense_{i}')(x)\n",
    "        x = layers.BatchNormalization(name=f'encoder_bn_{i}')(x)\n",
    "        x = layers.Dropout(0.2, name=f'encoder_dropout_{i}')(x)\n",
    "    \n",
    "    z_mean = layers.Dense(latent_dim, name='z_mean')(x)\n",
    "    z_log_var = layers.Dense(latent_dim, name='z_log_var')(x)\n",
    "    z = Sampling()([z_mean, z_log_var])\n",
    "    \n",
    "    encoder = Model(encoder_inputs, [z_mean, z_log_var, z], name='encoder')\n",
    "    \n",
    "    # ---- DECODER ----\n",
    "    latent_inputs = keras.Input(shape=(latent_dim,), name='latent_input')\n",
    "    x = latent_inputs\n",
    "    \n",
    "    for i, dim in enumerate(decoder_dims):\n",
    "        x = layers.Dense(dim, activation='relu', name=f'decoder_dense_{i}')(x)\n",
    "        x = layers.BatchNormalization(name=f'decoder_bn_{i}')(x)\n",
    "    \n",
    "    decoder_outputs = layers.Dense(input_dim, activation='linear', name='decoder_output')(x)\n",
    "    decoder = Model(latent_inputs, decoder_outputs, name='decoder')\n",
    "    \n",
    "    # ---- CLASSIFIER HEAD ----\n",
    "    x = latent_inputs\n",
    "    for i, dim in enumerate(classifier_dims):\n",
    "        x = layers.Dense(dim, activation='relu', name=f'classifier_dense_{i}')(x)\n",
    "        x = layers.Dropout(0.3, name=f'classifier_dropout_{i}')(x)\n",
    "    \n",
    "    classifier_output = layers.Dense(1, activation='sigmoid', name='classifier_output')(x)\n",
    "    classifier = Model(latent_inputs, classifier_output, name='classifier')\n",
    "    \n",
    "    # ---- FULL MODEL ----\n",
    "    inputs = keras.Input(shape=(input_dim,), name='vae_input')\n",
    "    z_mean, z_log_var, z = encoder(inputs)\n",
    "    reconstructed = decoder(z)\n",
    "    classification = classifier(z)\n",
    "    \n",
    "    model = Model(inputs, [reconstructed, classification], name='vae_classifier')\n",
    "    \n",
    "    return model, encoder, decoder, classifier\n",
    "\n",
    "# ============================================================================\n",
    "# 4. CUSTOM TRAINING LOOP WITH COMBINED LOSS\n",
    "# ============================================================================\n",
    "\n",
    "class VAEClassifier(keras.Model):\n",
    "    def __init__(self, encoder, decoder, classifier, **kwargs):\n",
    "        super().__init__(**kwargs)\n",
    "        self.encoder = encoder\n",
    "        self.decoder = decoder\n",
    "        self.classifier = classifier\n",
    "        self.total_loss_tracker = keras.metrics.Mean(name=\"total_loss\")\n",
    "        self.reconstruction_loss_tracker = keras.metrics.Mean(name=\"reconstruction_loss\")\n",
    "        self.kl_loss_tracker = keras.metrics.Mean(name=\"kl_loss\")\n",
    "        self.classification_loss_tracker = keras.metrics.Mean(name=\"classification_loss\")\n",
    "        self.accuracy_tracker = keras.metrics.BinaryAccuracy(name=\"accuracy\")\n",
    "        \n",
    "    @property\n",
    "    def metrics(self):\n",
    "        return [\n",
    "            self.total_loss_tracker,\n",
    "            self.reconstruction_loss_tracker,\n",
    "            self.kl_loss_tracker,\n",
    "            self.classification_loss_tracker,\n",
    "            self.accuracy_tracker,\n",
    "        ]\n",
    "    \n",
    "    def train_step(self, data):\n",
    "        x, y = data\n",
    "        \n",
    "        with tf.GradientTape() as tape:\n",
    "            # Forward pass\n",
    "            z_mean, z_log_var, z = self.encoder(x, training=True)\n",
    "            reconstruction = self.decoder(z, training=True)\n",
    "            classification = self.classifier(z, training=True)\n",
    "            \n",
    "            # Reconstruction loss (MSE)\n",
    "            reconstruction_loss = tf.reduce_mean(\n",
    "            keras.losses.mse(x, reconstruction))\n",
    "            \n",
    "            # KL divergence loss\n",
    "            kl_loss = -0.5 * tf.reduce_mean(\n",
    "                tf.reduce_sum(\n",
    "                    1 + z_log_var - tf.square(z_mean) - tf.exp(z_log_var),\n",
    "                    axis=1\n",
    "                )\n",
    "            )\n",
    "            \n",
    "            # Classification loss (binary crossentropy)\n",
    "            # Classification loss (binary crossentropy)\n",
    "            classification_loss = tf.reduce_mean(\n",
    "                keras.losses.binary_crossentropy(tf.expand_dims(y, -1), classification)\n",
    "            )\n",
    "            \n",
    "            # Combined loss with weights\n",
    "            total_loss = reconstruction_loss + kl_loss + classification_loss\n",
    "        \n",
    "        # Backpropagation\n",
    "        grads = tape.gradient(total_loss, self.trainable_weights)\n",
    "        self.optimizer.apply_gradients(zip(grads, self.trainable_weights))\n",
    "        \n",
    "        # Update metrics\n",
    "        self.total_loss_tracker.update_state(total_loss)\n",
    "        self.reconstruction_loss_tracker.update_state(reconstruction_loss)\n",
    "        self.kl_loss_tracker.update_state(kl_loss)\n",
    "        self.classification_loss_tracker.update_state(classification_loss)\n",
    "        self.accuracy_tracker.update_state(y, classification)\n",
    "        \n",
    "        return {\n",
    "            \"total_loss\": self.total_loss_tracker.result(),\n",
    "            \"reconstruction_loss\": self.reconstruction_loss_tracker.result(),\n",
    "            \"kl_loss\": self.kl_loss_tracker.result(),\n",
    "            \"classification_loss\": self.classification_loss_tracker.result(),\n",
    "            \"accuracy\": self.accuracy_tracker.result(),\n",
    "        }\n",
    "    \n",
    "    def test_step(self, data):\n",
    "        x, y = data\n",
    "        \n",
    "        z_mean, z_log_var, z = self.encoder(x, training=False)\n",
    "        reconstruction = self.decoder(z, training=False)\n",
    "        classification = self.classifier(z, training=False)\n",
    "        \n",
    "        # Reconstruction loss (MSE)\n",
    "        reconstruction_loss = tf.reduce_mean(\n",
    "            keras.losses.mse(x, reconstruction))\n",
    "        kl_loss = -0.5 * tf.reduce_mean(\n",
    "            tf.reduce_sum(1 + z_log_var - tf.square(z_mean) - tf.exp(z_log_var), axis=1)\n",
    "        )\n",
    "        # Classification loss (binary crossentropy)\n",
    "        classification_loss = tf.reduce_mean(\n",
    "            keras.losses.binary_crossentropy(tf.expand_dims(y, -1), classification)\n",
    "        )\n",
    "        total_loss = reconstruction_loss + kl_loss + classification_loss\n",
    "        \n",
    "        self.total_loss_tracker.update_state(total_loss)\n",
    "        self.reconstruction_loss_tracker.update_state(reconstruction_loss)\n",
    "        self.kl_loss_tracker.update_state(kl_loss)\n",
    "        self.classification_loss_tracker.update_state(classification_loss)\n",
    "        self.accuracy_tracker.update_state(y, classification)\n",
    "        \n",
    "        return {\n",
    "            \"total_loss\": self.total_loss_tracker.result(),\n",
    "            \"reconstruction_loss\": self.reconstruction_loss_tracker.result(),\n",
    "            \"kl_loss\": self.kl_loss_tracker.result(),\n",
    "            \"classification_loss\": self.classification_loss_tracker.result(),\n",
    "            \"accuracy\": self.accuracy_tracker.result(),\n",
    "        }\n",
    "\n",
    "# ============================================================================\n",
    "# 5. GROUP K-FOLD CROSS-VALIDATION WITH GRID SEARCH\n",
    "# ============================================================================\n",
    "\n",
    "# Hyperparameter grid\n",
    "param_grid = {\n",
    "    'latent_dim': [2, 5],\n",
    "    'encoder_dims': [[32, 16], [64, 32]],\n",
    "    'learning_rate': [0.001, 0.005],\n",
    "    'batch_size': [32, 64],\n",
    "}\n",
    "\n",
    "# Generate all combinations\n",
    "from itertools import product\n",
    "keys = param_grid.keys()\n",
    "values = param_grid.values()\n",
    "param_combinations = [dict(zip(keys, v)) for v in product(*values)]\n",
    "\n",
    "print(f\"\\nTotal hyperparameter combinations: {len(param_combinations)}\")\n",
    "\n",
    "# Group K-Fold setup\n",
    "n_splits = 5\n",
    "gkf = GroupKFold(n_splits=n_splits)\n",
    "\n",
    "# Store results\n",
    "cv_results = []\n",
    "\n",
    "# Grid search with cross-validation\n",
    "for idx, params in enumerate(param_combinations):\n",
    "    print(f\"\\n{'='*80}\")\n",
    "    print(f\"Testing combination {idx+1}/{len(param_combinations)}: {params}\")\n",
    "    print(f\"{'='*80}\")\n",
    "    \n",
    "    fold_results = []\n",
    "    \n",
    "    for fold, (train_idx, val_idx) in enumerate(gkf.split(X_train_scaled, y_train, groups_train)):\n",
    "        print(f\"\\nFold {fold+1}/{n_splits}\")\n",
    "        \n",
    "        X_fold_train, X_fold_val = X_train_scaled[train_idx], X_train_scaled[val_idx]\n",
    "        y_fold_train, y_fold_val = y_train[train_idx], y_train[val_idx]\n",
    "        \n",
    "        # Build model\n",
    "        model, encoder, decoder, classifier = build_vae_classifier(\n",
    "            input_dim=len(au_columns),\n",
    "            latent_dim=params['latent_dim'],\n",
    "            encoder_dims=params['encoder_dims'],\n",
    "            decoder_dims=list(reversed(params['encoder_dims'])),\n",
    "            classifier_dims=[16]\n",
    "        )\n",
    "        \n",
    "        vae_classifier = VAEClassifier(encoder, decoder, classifier)\n",
    "        vae_classifier.compile(optimizer=keras.optimizers.Adam(params['learning_rate']))\n",
    "        \n",
    "        # Early stopping\n",
    "        early_stop = keras.callbacks.EarlyStopping(\n",
    "            monitor='val_total_loss',\n",
    "            patience=10,\n",
    "            restore_best_weights=True,\n",
    "            mode='min'\n",
    "        )\n",
    "        \n",
    "        # Train\n",
    "        history = vae_classifier.fit(\n",
    "            X_fold_train, y_fold_train,\n",
    "            validation_data=(X_fold_val, y_fold_val),\n",
    "            epochs=60,\n",
    "            batch_size=params['batch_size'],\n",
    "            callbacks=[early_stop],\n",
    "            verbose=0\n",
    "        )\n",
    "        \n",
    "        # Evaluate on validation fold\n",
    "        z_mean_val, _, _ = encoder.predict(X_fold_val, verbose=0)\n",
    "        y_pred_proba = classifier.predict(z_mean_val, verbose=0).flatten()\n",
    "        y_pred = (y_pred_proba > 0.5).astype(int)\n",
    "        \n",
    "        fold_metrics = {\n",
    "            'accuracy': accuracy_score(y_fold_val, y_pred),\n",
    "            'precision': precision_score(y_fold_val, y_pred, zero_division=0),\n",
    "            'recall': recall_score(y_fold_val, y_pred, zero_division=0),\n",
    "            'f1': f1_score(y_fold_val, y_pred, zero_division=0),\n",
    "            'roc_auc': roc_auc_score(y_fold_val, y_pred_proba),\n",
    "            'final_recon_loss': history.history['val_reconstruction_loss'][-1],\n",
    "            'final_kl_loss': history.history['val_kl_loss'][-1],\n",
    "            'final_class_loss': history.history['val_classification_loss'][-1],\n",
    "        }\n",
    "        \n",
    "        fold_results.append(fold_metrics)\n",
    "        print(f\"  Accuracy: {fold_metrics['accuracy']:.4f}, F1: {fold_metrics['f1']:.4f}, AUC: {fold_metrics['roc_auc']:.4f}\")\n",
    "        \n",
    "        # Clear session to free memory\n",
    "        keras.backend.clear_session()\n",
    "    \n",
    "    # Average across folds\n",
    "    avg_results = {\n",
    "        'params': params,\n",
    "        'mean_accuracy': np.mean([r['accuracy'] for r in fold_results]),\n",
    "        'std_accuracy': np.std([r['accuracy'] for r in fold_results]),\n",
    "        'mean_f1': np.mean([r['f1'] for r in fold_results]),\n",
    "        'std_f1': np.std([r['f1'] for r in fold_results]),\n",
    "        'mean_roc_auc': np.mean([r['roc_auc'] for r in fold_results]),\n",
    "        'std_roc_auc': np.std([r['roc_auc'] for r in fold_results]),\n",
    "        'mean_recon_loss': np.mean([r['final_recon_loss'] for r in fold_results]),\n",
    "        'mean_kl_loss': np.mean([r['final_kl_loss'] for r in fold_results]),\n",
    "        'mean_class_loss': np.mean([r['final_class_loss'] for r in fold_results]),\n",
    "        'fold_results': fold_results\n",
    "    }\n",
    "    \n",
    "    cv_results.append(avg_results)\n",
    "    \n",
    "    print(f\"\\nMean CV Accuracy: {avg_results['mean_accuracy']:.4f} ± {avg_results['std_accuracy']:.4f}\")\n",
    "    print(f\"Mean CV F1: {avg_results['mean_f1']:.4f} ± {avg_results['std_f1']:.4f}\")\n",
    "    print(f\"Mean CV AUC: {avg_results['mean_roc_auc']:.4f} ± {avg_results['std_roc_auc']:.4f}\")\n",
    "\n",
    "# ============================================================================\n",
    "# 6. SELECT BEST MODEL AND EVALUATE ON TEST SET\n",
    "# ============================================================================\n",
    "\n",
    "# Find best hyperparameters based on mean F1 score\n",
    "best_idx = np.argmax([r['mean_f1'] for r in cv_results])\n",
    "best_params = cv_results[best_idx]['params']\n",
    "\n",
    "print(f\"\\n{'='*80}\")\n",
    "print(\"BEST HYPERPARAMETERS (based on CV F1 score):\")\n",
    "print(f\"{'='*80}\")\n",
    "for key, value in best_params.items():\n",
    "    print(f\"{key}: {value}\")\n",
    "print(f\"\\nCV Performance:\")\n",
    "print(f\"  Accuracy: {cv_results[best_idx]['mean_accuracy']:.4f} ± {cv_results[best_idx]['std_accuracy']:.4f}\")\n",
    "print(f\"  F1 Score: {cv_results[best_idx]['mean_f1']:.4f} ± {cv_results[best_idx]['std_f1']:.4f}\")\n",
    "print(f\"  ROC-AUC: {cv_results[best_idx]['mean_roc_auc']:.4f} ± {cv_results[best_idx]['std_roc_auc']:.4f}\")\n",
    "\n",
    "# Train final model on all training data\n",
    "print(f\"\\n{'='*80}\")\n",
    "print(\"TRAINING FINAL MODEL ON ALL TRAINING DATA\")\n",
    "print(f\"{'='*80}\")\n",
    "\n",
    "final_model, final_encoder, final_decoder, final_classifier = build_vae_classifier(\n",
    "    input_dim=len(au_columns),\n",
    "    latent_dim=best_params['latent_dim'],\n",
    "    encoder_dims=best_params['encoder_dims'],\n",
    "    decoder_dims=list(reversed(best_params['encoder_dims'])),\n",
    "    classifier_dims=[16]\n",
    ")\n",
    "\n",
    "final_vae_classifier = VAEClassifier(final_encoder, final_decoder, final_classifier)\n",
    "final_vae_classifier.compile(optimizer=keras.optimizers.Adam(best_params['learning_rate']))\n",
    "\n",
    "final_history = final_vae_classifier.fit(\n",
    "    X_train_scaled, y_train,\n",
    "    validation_split=0.2,\n",
    "    epochs=100,\n",
    "    batch_size=best_params['batch_size'],\n",
    "    callbacks=[keras.callbacks.EarlyStopping(monitor='val_total_loss', patience=15, restore_best_weights=True, mode='min')],\n",
    "    verbose=1\n",
    ")\n",
    "\n",
    "# Evaluate on held-out test set\n",
    "print(f\"\\n{'='*80}\")\n",
    "print(\"EVALUATION ON HELD-OUT TEST SET\")\n",
    "print(f\"{'='*80}\")\n",
    "\n",
    "z_mean_test, _, _ = final_encoder.predict(X_test_scaled, verbose=0)\n",
    "y_test_pred_proba = final_classifier.predict(z_mean_test, verbose=0).flatten()\n",
    "y_test_pred = (y_test_pred_proba > 0.5).astype(int)\n",
    "\n",
    "test_metrics = {\n",
    "    'accuracy': accuracy_score(y_test, y_test_pred),\n",
    "    'precision': precision_score(y_test, y_test_pred),\n",
    "    'recall': recall_score(y_test, y_test_pred),\n",
    "    'f1': f1_score(y_test, y_test_pred),\n",
    "    'roc_auc': roc_auc_score(y_test, y_test_pred_proba),\n",
    "}\n",
    "\n",
    "print(\"\\nTest Set Performance:\")\n",
    "for metric, value in test_metrics.items():\n",
    "    print(f\"  {metric.capitalize()}: {value:.4f}\")\n",
    "\n",
    "print(\"\\nConfusion Matrix:\")\n",
    "print(confusion_matrix(y_test, y_test_pred))\n",
    "\n",
    "print(\"\\nClassification Report:\")\n",
    "print(classification_report(y_test, y_test_pred, target_names=['Low Workload', 'High Workload']))\n",
    "\n",
    "# ============================================================================\n",
    "# 7. VISUALIZATION\n",
    "# ============================================================================\n",
    "\n",
    "# Plot training history\n",
    "fig, axes = plt.subplots(2, 2, figsize=(15, 10))\n",
    "\n",
    "axes[0, 0].plot(final_history.history['reconstruction_loss'], label='Train')\n",
    "axes[0, 0].plot(final_history.history['val_reconstruction_loss'], label='Val')\n",
    "axes[0, 0].set_title('Reconstruction Loss')\n",
    "axes[0, 0].set_xlabel('Epoch')\n",
    "axes[0, 0].set_ylabel('Loss')\n",
    "axes[0, 0].legend()\n",
    "axes[0, 0].grid(True)\n",
    "\n",
    "axes[0, 1].plot(final_history.history['kl_loss'], label='Train')\n",
    "axes[0, 1].plot(final_history.history['val_kl_loss'], label='Val')\n",
    "axes[0, 1].set_title('KL Divergence Loss')\n",
    "axes[0, 1].set_xlabel('Epoch')\n",
    "axes[0, 1].set_ylabel('Loss')\n",
    "axes[0, 1].legend()\n",
    "axes[0, 1].grid(True)\n",
    "\n",
    "axes[1, 0].plot(final_history.history['classification_loss'], label='Train')\n",
    "axes[1, 0].plot(final_history.history['val_classification_loss'], label='Val')\n",
    "axes[1, 0].set_title('Classification Loss')\n",
    "axes[1, 0].set_xlabel('Epoch')\n",
    "axes[1, 0].set_ylabel('Loss')\n",
    "axes[1, 0].legend()\n",
    "axes[1, 0].grid(True)\n",
    "\n",
    "axes[1, 1].plot(final_history.history['accuracy'], label='Train')\n",
    "axes[1, 1].plot(final_history.history['val_accuracy'], label='Val')\n",
    "axes[1, 1].set_title('Classification Accuracy')\n",
    "axes[1, 1].set_xlabel('Epoch')\n",
    "axes[1, 1].set_ylabel('Accuracy')\n",
    "axes[1, 1].legend()\n",
    "axes[1, 1].grid(True)\n",
    "\n",
    "plt.tight_layout()\n",
    "plt.show()\n",
    "\n",
    "# Visualize latent space (if 2D or 3D)\n",
    "if best_params['latent_dim'] == 2:\n",
    "    z_mean_train, _, _ = final_encoder.predict(X_train_scaled, verbose=0)\n",
    "    \n",
    "    plt.figure(figsize=(10, 8))\n",
    "    scatter = plt.scatter(z_mean_train[:, 0], z_mean_train[:, 1], \n",
    "                         c=y_train, cmap='RdYlBu', alpha=0.6, edgecolors='k')\n",
    "    plt.colorbar(scatter, label='Workload (0=Low, 1=High)')\n",
    "    plt.xlabel('Latent Dimension 1')\n",
    "    plt.ylabel('Latent Dimension 2')\n",
    "    plt.title('2D Latent Space Representation (Training Data)')\n",
    "    plt.grid(True, alpha=0.3)\n",
    "    plt.show()\n",
    "    \n",
    "    # Test set latent space\n",
    "    plt.figure(figsize=(10, 8))\n",
    "    scatter = plt.scatter(z_mean_test[:, 0], z_mean_test[:, 1], \n",
    "                         c=y_test, cmap='RdYlBu', alpha=0.6, edgecolors='k')\n",
    "    plt.colorbar(scatter, label='Workload (0=Low, 1=High)')\n",
    "    plt.xlabel('Latent Dimension 1')\n",
    "    plt.ylabel('Latent Dimension 2')\n",
    "    plt.title('2D Latent Space Representation (Test Data)')\n",
    "    plt.grid(True, alpha=0.3)\n",
    "    plt.show()\n",
    "\n",
    "print(\"\\n\" + \"=\"*80)\n",
    "print(\"TRAINING COMPLETE!\")\n",
    "print(\"=\"*80)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "79bcfc58",
   "metadata": {},
   "outputs": [],
   "source": [
    "### Save Trained VAE Classifier Model\n",
    "from pathlib import Path\n",
    "from datetime import datetime\n",
    "\n",
    "# Define save path\n",
    "model_dir = Path(\"/home/jovyan/data-paulusjafahrsimulator-gpu/trained_models\")\n",
    "model_dir.mkdir(parents=True, exist_ok=True)\n",
    "\n",
    "timestamp = datetime.now().strftime(\"%Y%m%d_%H%M%S\")\n",
    "model_path = model_dir / f\"vae_classifier_{timestamp}.keras\"\n",
    "\n",
    "# Save the complete model\n",
    "final_vae_classifier.save(model_path)\n",
    "\n",
    "print(f\"Model saved to: {model_path}\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d700e517",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "30d8d100",
   "metadata": {},
   "outputs": [],
   "source": [
    "### Plot Confusion Matrix for Final Model\n",
    "from sklearn.metrics import ConfusionMatrixDisplay\n",
    "x = Path(\"/home/jovyan/data-paulusjafahrsimulator-gpu/trained_models/vae_classifier_20251210_230121.keras\")\n",
    "# Load the saved model\n",
    "print(f\"Loading model from: {x}\")\n",
    "# loaded_vae_classifier = tf.keras.models.load_model(x)\n",
    "loaded_vae_classifier = final_vae_classifier\n",
    "print(\"✓ Model loaded successfully!\")\n",
    "\n",
    "# Extract encoder and classifier from loaded model\n",
    "loaded_encoder = loaded_vae_classifier.encoder\n",
    "loaded_classifier = loaded_vae_classifier.classifier\n",
    "\n",
    "# Get predictions on test set\n",
    "z_mean_test, _, _ = loaded_encoder.predict(X_test_scaled, verbose=0)\n",
    "y_test_pred_proba = loaded_classifier.predict(z_mean_test, verbose=0).flatten()\n",
    "y_test_pred = (y_test_pred_proba > 0.5).astype(int)\n",
    "\n",
    "# Create and plot confusion matrix\n",
    "cm = confusion_matrix(y_test, y_test_pred)\n",
    "disp = ConfusionMatrixDisplay(confusion_matrix=cm, \n",
    "                               display_labels=['Low Workload', 'High Workload'])\n",
    "\n",
    "fig, ax = plt.subplots(figsize=(8, 6))\n",
    "disp.plot(ax=ax, cmap='Blues', values_format='d')\n",
    "plt.title('Confusion Matrix - Test Set (Loaded Model)')\n",
    "plt.tight_layout()\n",
    "plt.show()\n",
    "\n",
    "# Print metrics\n",
    "print(f\"\\nTest Set Performance (Loaded Model):\")\n",
    "print(f\"  Accuracy: {accuracy_score(y_test, y_test_pred):.4f}\")\n",
    "print(f\"  Precision: {precision_score(y_test, y_test_pred):.4f}\")\n",
    "print(f\"  Recall: {recall_score(y_test, y_test_pred):.4f}\")\n",
    "print(f\"  F1 Score: {f1_score(y_test, y_test_pred):.4f}\")\n",
    "print(f\"  ROC-AUC: {roc_auc_score(y_test, y_test_pred_proba):.4f}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e826a998",
   "metadata": {},
   "source": [
    "TO DO\n",
    "    * autoencoder langsam anfangen mit 19 schichten\n",
    "    * dann AE und SVM mit hybridem training wie bei claude?!\n",
    "    * dataset aus eyetracking verwenden?"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
