{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Predicting on new samples\n", "\n", "The model can be loaded from disk, and used to predict on unseen samples. Here, we keep using the DustPedia and H-ATLAS datasets, but these can be replaced by other UV-NIR SED fitted datasets." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import os\n", "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "\n", "pd.options.display.max_columns = 99" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "os.chdir('..') # change to root directory" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "dict_keys(['fullbay', 'fullbayerr', 'shortbay', 'shortbayerr', 'redshift', 'observed', 'observederr', 'obserr_to_short', 'obs_to_short'])" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Load data\n", "import pickle\n", "from firenet.util import add_uncertainty_features\n", "with open('./data/d_data.pkl', 'rb') as infile:\n", " d_data = pickle.load(infile)\n", " \n", "d_data = add_uncertainty_features(d_data)\n", "d_data.keys()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Load predictor that uses all data\n", "from firenet.ml.modelstore import ModelStore\n", "\n", "pred = ModelStore().load(d_data, name='nnet_alldata')" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
GALEX_FUVGALEX_NUVSDSS_uSDSS_gSDSS_rSDSS_iSDSS_z2MASS_J2MASS_H2MASS_KsWISE_3.4WISE_4.6WISE_12WISE_22
id
NGC53581.756567e+193.544724e+193.709311e+201.694206e+213.080484e+214.145965e+215.426019e+216.893655e+218.360429e+216.609317e+213.423306e+212.013436e+215.835990e+202.759253e+20
NGC32526.295663e+191.014847e+203.032180e+201.022812e+211.527502e+211.894299e+212.336673e+212.844187e+213.403842e+212.767714e+211.838617e+211.113030e+213.448660e+214.372679e+21
ESO407-0022.785232e+195.601229e+192.264886e+208.530379e+201.336088e+211.702040e+212.127060e+212.631649e+213.178567e+212.588323e+211.667594e+219.863319e+202.688774e+212.900365e+21
NGC41944.120649e+209.560917e+202.320822e+216.091552e+219.006956e+211.119523e+221.449513e+221.840004e+222.248789e+221.933412e+221.710598e+221.566078e+221.696050e+235.295668e+23
NGC71723.195182e+204.884440e+201.914068e+218.348678e+211.554232e+222.240839e+223.119440e+224.411116e+225.804870e+224.923501e+223.200892e+221.953044e+225.200287e+227.908303e+22
\n", "
" ], "text/plain": [ " GALEX_FUV GALEX_NUV SDSS_u SDSS_g \\\n", "id \n", "NGC5358 1.756567e+19 3.544724e+19 3.709311e+20 1.694206e+21 \n", "NGC3252 6.295663e+19 1.014847e+20 3.032180e+20 1.022812e+21 \n", "ESO407-002 2.785232e+19 5.601229e+19 2.264886e+20 8.530379e+20 \n", "NGC4194 4.120649e+20 9.560917e+20 2.320822e+21 6.091552e+21 \n", "NGC7172 3.195182e+20 4.884440e+20 1.914068e+21 8.348678e+21 \n", "\n", " SDSS_r SDSS_i SDSS_z 2MASS_J \\\n", "id \n", "NGC5358 3.080484e+21 4.145965e+21 5.426019e+21 6.893655e+21 \n", "NGC3252 1.527502e+21 1.894299e+21 2.336673e+21 2.844187e+21 \n", "ESO407-002 1.336088e+21 1.702040e+21 2.127060e+21 2.631649e+21 \n", "NGC4194 9.006956e+21 1.119523e+22 1.449513e+22 1.840004e+22 \n", "NGC7172 1.554232e+22 2.240839e+22 3.119440e+22 4.411116e+22 \n", "\n", " 2MASS_H 2MASS_Ks WISE_3.4 WISE_4.6 \\\n", "id \n", "NGC5358 8.360429e+21 6.609317e+21 3.423306e+21 2.013436e+21 \n", "NGC3252 3.403842e+21 2.767714e+21 1.838617e+21 1.113030e+21 \n", "ESO407-002 3.178567e+21 2.588323e+21 1.667594e+21 9.863319e+20 \n", "NGC4194 2.248789e+22 1.933412e+22 1.710598e+22 1.566078e+22 \n", "NGC7172 5.804870e+22 4.923501e+22 3.200892e+22 1.953044e+22 \n", "\n", " WISE_12 WISE_22 \n", "id \n", "NGC5358 5.835990e+20 2.759253e+20 \n", "NGC3252 3.448660e+21 4.372679e+21 \n", "ESO407-002 2.688774e+21 2.900365e+21 \n", "NGC4194 1.696050e+23 5.295668e+23 \n", "NGC7172 5.200287e+22 7.908303e+22 " ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from firenet.ml import FeatureSelect\n", "\n", "X_reg = FeatureSelect.select_xreg(d_data)\n", "X_unc = FeatureSelect.select_xunc(d_data)\n", "X_reg.head()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
GALEX_FUVGALEX_NUVSDSS_uSDSS_gSDSS_rSDSS_iSDSS_z2MASS_J2MASS_H2MASS_KsWISE_3.4WISE_4.6WISE_12WISE_22
id
NGC5358-2.289781-1.984863-0.965152-0.305480-0.0458270.0831800.2000360.3040040.3877830.28571121.534446-0.230508-0.768331-1.093654
NGC3252-1.465450-1.258091-0.782736-0.254695-0.0805090.0129570.1041070.1894670.2674780.17763021.264491-0.2179840.2731590.376256
ESO407-002-1.777229-1.473807-0.867044-0.291122-0.0962550.0088800.1056900.1981380.2801410.19092821.222090-0.2280670.2074640.240362
NGC4194-1.618182-1.252648-0.867506-0.448420-0.278570-0.184115-0.0719260.0316710.1188010.05317622.233148-0.0383350.9962911.490773
NGC7172-2.000775-1.816456-1.223314-0.583653-0.313755-0.154860-0.0111940.1392770.2585220.18700322.505271-0.2145590.2107560.392812
\n", "
" ], "text/plain": [ " GALEX_FUV GALEX_NUV SDSS_u SDSS_g SDSS_r SDSS_i \\\n", "id \n", "NGC5358 -2.289781 -1.984863 -0.965152 -0.305480 -0.045827 0.083180 \n", "NGC3252 -1.465450 -1.258091 -0.782736 -0.254695 -0.080509 0.012957 \n", "ESO407-002 -1.777229 -1.473807 -0.867044 -0.291122 -0.096255 0.008880 \n", "NGC4194 -1.618182 -1.252648 -0.867506 -0.448420 -0.278570 -0.184115 \n", "NGC7172 -2.000775 -1.816456 -1.223314 -0.583653 -0.313755 -0.154860 \n", "\n", " SDSS_z 2MASS_J 2MASS_H 2MASS_Ks WISE_3.4 WISE_4.6 \\\n", "id \n", "NGC5358 0.200036 0.304004 0.387783 0.285711 21.534446 -0.230508 \n", "NGC3252 0.104107 0.189467 0.267478 0.177630 21.264491 -0.217984 \n", "ESO407-002 0.105690 0.198138 0.280141 0.190928 21.222090 -0.228067 \n", "NGC4194 -0.071926 0.031671 0.118801 0.053176 22.233148 -0.038335 \n", "NGC7172 -0.011194 0.139277 0.258522 0.187003 22.505271 -0.214559 \n", "\n", " WISE_12 WISE_22 \n", "id \n", "NGC5358 -0.768331 -1.093654 \n", "NGC3252 0.273159 0.376256 \n", "ESO407-002 0.207464 0.240362 \n", "NGC4194 0.996291 1.490773 \n", "NGC7172 0.210756 0.392812 " ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "X_reg = pred.reg.log_normaliser.transform(X_reg)\n", "X_unc = pred.unc.log_normaliser.transform(X_unc)\n", "X_reg.head()" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "Y_pred, Y_prederr = pred.predict(X_reg, X_unc)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PACS_70PACS_100PACS_160SPIRE_250SPIRE_350SPIRE_500
id
NGC5358-0.2704550.0309070.1612520.055492-0.166314-0.518062
NGC32521.4192601.7763221.9179391.7456211.4335451.062444
ESO407-0021.3084021.6945401.8381741.6742241.3845071.007923
NGC41942.4519062.5232832.2566911.7511991.3351910.864598
NGC71721.5738761.9180621.9547081.7095161.3921640.981167
\n", "
" ], "text/plain": [ " PACS_70 PACS_100 PACS_160 SPIRE_250 SPIRE_350 SPIRE_500\n", "id \n", "NGC5358 -0.270455 0.030907 0.161252 0.055492 -0.166314 -0.518062\n", "NGC3252 1.419260 1.776322 1.917939 1.745621 1.433545 1.062444\n", "ESO407-002 1.308402 1.694540 1.838174 1.674224 1.384507 1.007923\n", "NGC4194 2.451906 2.523283 2.256691 1.751199 1.335191 0.864598\n", "NGC7172 1.573876 1.918062 1.954708 1.709516 1.392164 0.981167" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Y_pred.head()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAQ1UlEQVR4nO3df4xlZ13H8ffH/gCDyrZ0WprdxYWwUdAEqJO6hoRgS0x/GLaJNCkauzRLNmhVCCa6YiJq/GPrH1QbTHGl6JYgtKliVyhqbSHGP1qdQik/Vu3S1HaytTv8WqhVTOHrH/MsTHfv7JyZnZl75+n7ldycc57z3LnfM3f2M88895yzqSokSX35vnEXIElafYa7JHXIcJekDhnuktQhw12SOnTmuAsAOO+882rbtm3jLkOSNpQHHnjgy1U1NWrfRIT7tm3bmJmZGXcZkrShJPnPxfY5LSNJHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR2aiCtUJfVt296PD+r36L4r17iS5w5H7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOeRGTpBUZemGSxsORuyR1aFC4J9mU5I4k/5bkUJKfSnJukruTPNyW57S+SXJTksNJHkpy0doegiTpRENH7n8M/F1V/SjwKuAQsBe4p6q2A/e0bYDLge3tsQe4eVUrliQtaclwT/JDwOuAWwCq6v+q6uvATuBA63YAuKqt7wRurXn3AZuSXLjqlUuSFjVk5P4yYA748ySfSfL+JC8ALqiqJwDa8vzWfzPw+ILnz7a2Z0myJ8lMkpm5ubnTOghJ0rMNCfczgYuAm6vqNcB/870pmFEyoq1OaqjaX1XTVTU9NTU1qFhJ0jBDwn0WmK2q+9v2HcyH/ZPHp1va8uiC/lsXPH8LcGR1ypUkDbFkuFfVfwGPJ/mR1nQp8EXgILCrte0C7mzrB4Fr21kzO4Bjx6dvJEnrY+hFTL8KfCjJ2cAjwHXM/2K4Pclu4DHg6tb3LuAK4DDwdOsrSVpHg8K9qh4EpkfsunRE3wKuP826JEmnwStUJalDhrskdchwl6QOGe6S1CFv+SvpWbyVbx8cuUtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KH/J+YJE2Mof8L1KP7rlzjSja+QSP3JI8m+VySB5PMtLZzk9yd5OG2PKe1J8lNSQ4neSjJRWt5AJKkky1nWuanq+rVVTXdtvcC91TVduCetg1wObC9PfYAN69WsZKkYU5nzn0ncKCtHwCuWtB+a827D9iU5MLTeB1J0jINDfcC/iHJA0n2tLYLquoJgLY8v7VvBh5f8NzZ1vYsSfYkmUkyMzc3t7LqJUkjDf1A9bVVdSTJ+cDdSf7tFH0zoq1OaqjaD+wHmJ6ePmm/JGnlBo3cq+pIWx4FPgpcDDx5fLqlLY+27rPA1gVP3wIcWa2CJUlLWzLck7wgyQ8eXwd+Bvg8cBDY1brtAu5s6weBa9tZMzuAY8enbyRJ62PItMwFwEeTHO//l1X1d0n+Fbg9yW7gMeDq1v8u4ArgMPA0cN2qVy1JOqUlw72qHgFeNaL9K8ClI9oLuH5VqpMkrYi3H5CkDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4PDPckZST6T5GNt+6VJ7k/ycJLbkpzd2p/Xtg+3/dvWpnRJ0mKWM3J/O3BowfYNwI1VtR34GrC7te8GvlZVLwdubP0kSetoULgn2QJcCby/bQe4BLijdTkAXNXWd7Zt2v5LW39J0joZOnL/I+A3gO+07RcBX6+qZ9r2LLC5rW8GHgdo+4+1/s+SZE+SmSQzc3NzKyxfkjTKkuGe5GeBo1X1wMLmEV1rwL7vNVTtr6rpqpqempoaVKwkaZgzB/R5LfDGJFcAzwd+iPmR/KYkZ7bR+RbgSOs/C2wFZpOcCbwQ+OqqVy5pWbbt/fi4S9A6WnLkXlW/VVVbqmobcA1wb1X9AvBJ4E2t2y7gzrZ+sG3T9t9bVSeN3CVJa+d0znP/TeCdSQ4zP6d+S2u/BXhRa38nsPf0SpQkLdeQaZnvqqpPAZ9q648AF4/o87/A1atQmyRphbxCVZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHVrWf5AtSZNg296PD+r36L4r17iSyeXIXZI6ZLhLUoeWDPckz0/yL0k+m+QLSX6vtb80yf1JHk5yW5KzW/vz2vbhtn/b2h6CJOlEQ0bu3wIuqapXAa8GLkuyA7gBuLGqtgNfA3a3/ruBr1XVy4EbWz9J0jpa8gPVqirgqbZ5VnsUcAnw8639APC7wM3AzrYOcAfw3iRpX0fSKhv64aKeWwbNuSc5I8mDwFHgbuBLwNer6pnWZRbY3NY3A48DtP3HgBeN+Jp7kswkmZmbmzu9o5AkPcugcK+qb1fVq4EtwMXAK0Z1a8ucYt/Cr7m/qqaranpqampovZKkAZZ1tkxVfR34FLAD2JTk+LTOFuBIW58FtgK0/S8EvroaxUqShhlytsxUkk1t/fuBNwCHgE8Cb2rddgF3tvWDbZu2/17n2yVpfQ25QvVC4ECSM5j/ZXB7VX0syReBjyT5A+AzwC2t/y3AB5McZn7Efs0a1C1JOoUhZ8s8BLxmRPsjzM+/n9j+v8DVq1KdJGlFvEJVkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4tGe5Jtib5ZJJDSb6Q5O2t/dwkdyd5uC3Pae1JclOSw0keSnLRWh+EJOnZhozcnwF+vapeAewArk/ySmAvcE9VbQfuadsAlwPb22MPcPOqVy1JOqUlw72qnqiqT7f1bwKHgM3ATuBA63YAuKqt7wRurXn3AZuSXLjqlUuSFrWsOfck24DXAPcDF1TVEzD/CwA4v3XbDDy+4GmzrU2StE4Gh3uSHwD+CnhHVX3jVF1HtNWIr7cnyUySmbm5uaFlSJIGGBTuSc5iPtg/VFV/3ZqfPD7d0pZHW/sssHXB07cAR078mlW1v6qmq2p6ampqpfVLkkYYcrZMgFuAQ1X1ngW7DgK72vou4M4F7de2s2Z2AMeOT99IktbHmQP6vBb4ReBzSR5sbe8C9gG3J9kNPAZc3fbdBVwBHAaeBq5b1YolSUtaMtyr6p8ZPY8OcOmI/gVcf5p1SZJOg1eoSlKHDHdJ6pDhLkkdGvKBqqQx2Lb34+MuQRuYI3dJ6pDhLkkdMtwlqUOGuyR1yA9UJXVr6IfSj+67co0rWX+O3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkd8iImaZ15t0etB0fuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUNLhnuSDyQ5muTzC9rOTXJ3kofb8pzWniQ3JTmc5KEkF61l8ZKk0YaM3P8CuOyEtr3APVW1HbinbQNcDmxvjz3AzatTpiRpOZYM96r6J+CrJzTvBA609QPAVQvab6159wGbkly4WsVKkoZZ6Zz7BVX1BEBbnt/aNwOPL+g329pOkmRPkpkkM3NzcyssQ5I0ymp/oJoRbTWqY1Xtr6rpqpqemppa5TIk6bltpeH+5PHplrY82tpnga0L+m0Bjqy8PEnSSqw03A8Cu9r6LuDOBe3XtrNmdgDHjk/fSJLWz5J3hUzyYeD1wHlJZoF3A/uA25PsBh4Drm7d7wKuAA4DTwPXrUHN0kTybo+aJEuGe1W9eZFdl47oW8D1p1uUJOn0eIWqJHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6tCSp0JKUu+GXqPw6L4r17iS1ePIXZI6ZLhLUoeclpFOwVsKaKNy5C5JHXLkruckR+TqnSN3SeqQ4S5JHTLcJalDhrskdcgPVLUh+AGotDyGu5alx8u0pR45LSNJHXLkLqc8pIGW829l3H+9Gu7r4Lk4leEvDGm8DPfTYIBJmlSG+wTxl4Wk1bIm4Z7kMuCPgTOA91fVvrV4nbViyEra6Fb9bJkkZwB/AlwOvBJ4c5JXrvbrSJIWtxYj94uBw1X1CECSjwA7gS+uwWs5ypY0kcZ9IsVahPtm4PEF27PAT57YKckeYE/bfCrJvy/jNc4DvrziCsfP+sdvox+D9Y/fqhxDbjitp//wYjvWItwzoq1OaqjaD+xf0QskM1U1vZLnTgLrH7+NfgzWP36TfgxrcYXqLLB1wfYW4MgavI4kaRFrEe7/CmxP8tIkZwPXAAfX4HUkSYtY9WmZqnomya8Af8/8qZAfqKovrPLLrGg6Z4JY//ht9GOw/vGb6GNI1UnT4ZKkDc67QkpShwx3SerQhgj3JOcmuTvJw215ziL9vp3kwfYY+4e4SS5L8u9JDifZO2L/85Lc1vbfn2Tb+le5uAH1vyXJ3ILv+VvHUediknwgydEkn19kf5Lc1I7voSQXrXeNpzKg/tcnObbg+/87613jqSTZmuSTSQ4l+UKSt4/oM7HvwcD6J/c9qKqJfwB/COxt63uBGxbp99S4a11QyxnAl4CXAWcDnwVeeUKfXwbe19avAW4bd93LrP8twHvHXespjuF1wEXA5xfZfwXwCeavzdgB3D/umpdZ/+uBj427zlPUfyFwUVv/QeA/RvwMTex7MLD+iX0PNsTInfnbFxxo6weAq8ZYy1DfvQ1DVf0fcPw2DAstPK47gEuTjLoIbByG1D/RquqfgK+eostO4Naadx+wKcmF61Pd0gbUP9Gq6omq+nRb/yZwiPkr2Bea2PdgYP0Ta6OE+wVV9QTMf8OB8xfp9/wkM0nuSzLuXwCjbsNw4g/Gd/tU1TPAMeBF61Ld0obUD/Bz7c/pO5JsHbF/kg09xkn2U0k+m+QTSX5s3MUspk05vga4/4RdG+I9OEX9MKHvwcTczz3JPwIvHrHrt5fxZV5SVUeSvAy4N8nnqupLq1Phsg25DcOgWzWMyZDa/hb4cFV9K8nbmP8r5JI1r2z1TPL3f4hPAz9cVU8luQL4G2D7mGs6SZIfAP4KeEdVfePE3SOeMlHvwRL1T+x7MDEj96p6Q1X9+IjHncCTx/9Ua8uji3yNI235CPAp5n/TjsuQ2zB8t0+SM4EXMjl/hi9Zf1V9paq+1Tb/DPiJdapttWzoW2VU1Teq6qm2fhdwVpLzxlzWsyQ5i/lg/FBV/fWILhP9HixV/yS/BxMT7ks4COxq67uAO0/skOScJM9r6+cBr2WNbjM80JDbMCw8rjcB91b7lGYCLFn/CXOjb2R+TnIjOQhc287Y2AEcOz79txEkefHxz2iSXMz8v+evjLeq72m13QIcqqr3LNJtYt+DIfVP8nswMdMyS9gH3J5kN/AYcDVAkmngbVX1VuAVwJ8m+Q7z3+B9VTW2cK9FbsOQ5PeBmao6yPwPzgeTHGZ+xH7NuOo90cD6fy3JG4FnmK//LWMreIQkH2b+bIbzkswC7wbOAqiq9wF3MX+2xmHgaeC68VQ62oD63wT8UpJngP8BrpmgwQHMD7B+Efhckgdb27uAl8CGeA+G1D+x74G3H5CkDm2UaRlJ0jIY7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalD/w+Ui5ZKTlYKcgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.hist(Y_pred['PACS_100'], bins=30);" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PACS_70PACS_100PACS_160SPIRE_250SPIRE_350SPIRE_500
id
NGC53581.836500e+213.675806e+214.962474e+213.889896e+212.334165e+211.038445e+21
NGC32524.827820e+221.098534e+231.522056e+231.023558e+234.989255e+222.122928e+22
ESO407-0023.392286e+228.253350e+221.148853e+237.876157e+224.042008e+221.698296e+22
NGC41944.842326e+245.707301e+243.089153e+249.646004e+233.701173e+231.252408e+23
NGC71721.199904e+242.650531e+242.883889e+241.639784e+247.896506e+233.065052e+23
\n", "
" ], "text/plain": [ " PACS_70 PACS_100 PACS_160 SPIRE_250 \\\n", "id \n", "NGC5358 1.836500e+21 3.675806e+21 4.962474e+21 3.889896e+21 \n", "NGC3252 4.827820e+22 1.098534e+23 1.522056e+23 1.023558e+23 \n", "ESO407-002 3.392286e+22 8.253350e+22 1.148853e+23 7.876157e+22 \n", "NGC4194 4.842326e+24 5.707301e+24 3.089153e+24 9.646004e+23 \n", "NGC7172 1.199904e+24 2.650531e+24 2.883889e+24 1.639784e+24 \n", "\n", " SPIRE_350 SPIRE_500 \n", "id \n", "NGC5358 2.334165e+21 1.038445e+21 \n", "NGC3252 4.989255e+22 2.122928e+22 \n", "ESO407-002 4.042008e+22 1.698296e+22 \n", "NGC4194 3.701173e+23 1.252408e+23 \n", "NGC7172 7.896506e+23 3.065052e+23 " ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "F_input, F_pred = pred.reg.log_normaliser.inverse_transform(X_reg, Y_pred)\n", "F_pred.head()" ] } ], "metadata": { "kernelspec": { "display_name": "firenet", "language": "python", "name": "firenet" }, "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.7.4" } }, "nbformat": 4, "nbformat_minor": 4 }