{
"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",
" GALEX_FUV | \n",
" GALEX_NUV | \n",
" SDSS_u | \n",
" SDSS_g | \n",
" SDSS_r | \n",
" SDSS_i | \n",
" SDSS_z | \n",
" 2MASS_J | \n",
" 2MASS_H | \n",
" 2MASS_Ks | \n",
" WISE_3.4 | \n",
" WISE_4.6 | \n",
" WISE_12 | \n",
" WISE_22 | \n",
"
\n",
" \n",
" id | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" NGC5358 | \n",
" 1.756567e+19 | \n",
" 3.544724e+19 | \n",
" 3.709311e+20 | \n",
" 1.694206e+21 | \n",
" 3.080484e+21 | \n",
" 4.145965e+21 | \n",
" 5.426019e+21 | \n",
" 6.893655e+21 | \n",
" 8.360429e+21 | \n",
" 6.609317e+21 | \n",
" 3.423306e+21 | \n",
" 2.013436e+21 | \n",
" 5.835990e+20 | \n",
" 2.759253e+20 | \n",
"
\n",
" \n",
" NGC3252 | \n",
" 6.295663e+19 | \n",
" 1.014847e+20 | \n",
" 3.032180e+20 | \n",
" 1.022812e+21 | \n",
" 1.527502e+21 | \n",
" 1.894299e+21 | \n",
" 2.336673e+21 | \n",
" 2.844187e+21 | \n",
" 3.403842e+21 | \n",
" 2.767714e+21 | \n",
" 1.838617e+21 | \n",
" 1.113030e+21 | \n",
" 3.448660e+21 | \n",
" 4.372679e+21 | \n",
"
\n",
" \n",
" ESO407-002 | \n",
" 2.785232e+19 | \n",
" 5.601229e+19 | \n",
" 2.264886e+20 | \n",
" 8.530379e+20 | \n",
" 1.336088e+21 | \n",
" 1.702040e+21 | \n",
" 2.127060e+21 | \n",
" 2.631649e+21 | \n",
" 3.178567e+21 | \n",
" 2.588323e+21 | \n",
" 1.667594e+21 | \n",
" 9.863319e+20 | \n",
" 2.688774e+21 | \n",
" 2.900365e+21 | \n",
"
\n",
" \n",
" NGC4194 | \n",
" 4.120649e+20 | \n",
" 9.560917e+20 | \n",
" 2.320822e+21 | \n",
" 6.091552e+21 | \n",
" 9.006956e+21 | \n",
" 1.119523e+22 | \n",
" 1.449513e+22 | \n",
" 1.840004e+22 | \n",
" 2.248789e+22 | \n",
" 1.933412e+22 | \n",
" 1.710598e+22 | \n",
" 1.566078e+22 | \n",
" 1.696050e+23 | \n",
" 5.295668e+23 | \n",
"
\n",
" \n",
" NGC7172 | \n",
" 3.195182e+20 | \n",
" 4.884440e+20 | \n",
" 1.914068e+21 | \n",
" 8.348678e+21 | \n",
" 1.554232e+22 | \n",
" 2.240839e+22 | \n",
" 3.119440e+22 | \n",
" 4.411116e+22 | \n",
" 5.804870e+22 | \n",
" 4.923501e+22 | \n",
" 3.200892e+22 | \n",
" 1.953044e+22 | \n",
" 5.200287e+22 | \n",
" 7.908303e+22 | \n",
"
\n",
" \n",
"
\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",
" GALEX_FUV | \n",
" GALEX_NUV | \n",
" SDSS_u | \n",
" SDSS_g | \n",
" SDSS_r | \n",
" SDSS_i | \n",
" SDSS_z | \n",
" 2MASS_J | \n",
" 2MASS_H | \n",
" 2MASS_Ks | \n",
" WISE_3.4 | \n",
" WISE_4.6 | \n",
" WISE_12 | \n",
" WISE_22 | \n",
"
\n",
" \n",
" id | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" NGC5358 | \n",
" -2.289781 | \n",
" -1.984863 | \n",
" -0.965152 | \n",
" -0.305480 | \n",
" -0.045827 | \n",
" 0.083180 | \n",
" 0.200036 | \n",
" 0.304004 | \n",
" 0.387783 | \n",
" 0.285711 | \n",
" 21.534446 | \n",
" -0.230508 | \n",
" -0.768331 | \n",
" -1.093654 | \n",
"
\n",
" \n",
" NGC3252 | \n",
" -1.465450 | \n",
" -1.258091 | \n",
" -0.782736 | \n",
" -0.254695 | \n",
" -0.080509 | \n",
" 0.012957 | \n",
" 0.104107 | \n",
" 0.189467 | \n",
" 0.267478 | \n",
" 0.177630 | \n",
" 21.264491 | \n",
" -0.217984 | \n",
" 0.273159 | \n",
" 0.376256 | \n",
"
\n",
" \n",
" ESO407-002 | \n",
" -1.777229 | \n",
" -1.473807 | \n",
" -0.867044 | \n",
" -0.291122 | \n",
" -0.096255 | \n",
" 0.008880 | \n",
" 0.105690 | \n",
" 0.198138 | \n",
" 0.280141 | \n",
" 0.190928 | \n",
" 21.222090 | \n",
" -0.228067 | \n",
" 0.207464 | \n",
" 0.240362 | \n",
"
\n",
" \n",
" NGC4194 | \n",
" -1.618182 | \n",
" -1.252648 | \n",
" -0.867506 | \n",
" -0.448420 | \n",
" -0.278570 | \n",
" -0.184115 | \n",
" -0.071926 | \n",
" 0.031671 | \n",
" 0.118801 | \n",
" 0.053176 | \n",
" 22.233148 | \n",
" -0.038335 | \n",
" 0.996291 | \n",
" 1.490773 | \n",
"
\n",
" \n",
" NGC7172 | \n",
" -2.000775 | \n",
" -1.816456 | \n",
" -1.223314 | \n",
" -0.583653 | \n",
" -0.313755 | \n",
" -0.154860 | \n",
" -0.011194 | \n",
" 0.139277 | \n",
" 0.258522 | \n",
" 0.187003 | \n",
" 22.505271 | \n",
" -0.214559 | \n",
" 0.210756 | \n",
" 0.392812 | \n",
"
\n",
" \n",
"
\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",
" PACS_70 | \n",
" PACS_100 | \n",
" PACS_160 | \n",
" SPIRE_250 | \n",
" SPIRE_350 | \n",
" SPIRE_500 | \n",
"
\n",
" \n",
" id | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" NGC5358 | \n",
" -0.270455 | \n",
" 0.030907 | \n",
" 0.161252 | \n",
" 0.055492 | \n",
" -0.166314 | \n",
" -0.518062 | \n",
"
\n",
" \n",
" NGC3252 | \n",
" 1.419260 | \n",
" 1.776322 | \n",
" 1.917939 | \n",
" 1.745621 | \n",
" 1.433545 | \n",
" 1.062444 | \n",
"
\n",
" \n",
" ESO407-002 | \n",
" 1.308402 | \n",
" 1.694540 | \n",
" 1.838174 | \n",
" 1.674224 | \n",
" 1.384507 | \n",
" 1.007923 | \n",
"
\n",
" \n",
" NGC4194 | \n",
" 2.451906 | \n",
" 2.523283 | \n",
" 2.256691 | \n",
" 1.751199 | \n",
" 1.335191 | \n",
" 0.864598 | \n",
"
\n",
" \n",
" NGC7172 | \n",
" 1.573876 | \n",
" 1.918062 | \n",
" 1.954708 | \n",
" 1.709516 | \n",
" 1.392164 | \n",
" 0.981167 | \n",
"
\n",
" \n",
"
\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",
" PACS_70 | \n",
" PACS_100 | \n",
" PACS_160 | \n",
" SPIRE_250 | \n",
" SPIRE_350 | \n",
" SPIRE_500 | \n",
"
\n",
" \n",
" id | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" NGC5358 | \n",
" 1.836500e+21 | \n",
" 3.675806e+21 | \n",
" 4.962474e+21 | \n",
" 3.889896e+21 | \n",
" 2.334165e+21 | \n",
" 1.038445e+21 | \n",
"
\n",
" \n",
" NGC3252 | \n",
" 4.827820e+22 | \n",
" 1.098534e+23 | \n",
" 1.522056e+23 | \n",
" 1.023558e+23 | \n",
" 4.989255e+22 | \n",
" 2.122928e+22 | \n",
"
\n",
" \n",
" ESO407-002 | \n",
" 3.392286e+22 | \n",
" 8.253350e+22 | \n",
" 1.148853e+23 | \n",
" 7.876157e+22 | \n",
" 4.042008e+22 | \n",
" 1.698296e+22 | \n",
"
\n",
" \n",
" NGC4194 | \n",
" 4.842326e+24 | \n",
" 5.707301e+24 | \n",
" 3.089153e+24 | \n",
" 9.646004e+23 | \n",
" 3.701173e+23 | \n",
" 1.252408e+23 | \n",
"
\n",
" \n",
" NGC7172 | \n",
" 1.199904e+24 | \n",
" 2.650531e+24 | \n",
" 2.883889e+24 | \n",
" 1.639784e+24 | \n",
" 7.896506e+23 | \n",
" 3.065052e+23 | \n",
"
\n",
" \n",
"
\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
}