final_output.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "cardiovascular-seeker",
   "metadata": {},
   "source": [
    "## External imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "moving-northeast",
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import pandas as pd\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "banned-musician",
   "metadata": {},
   "source": [
    "## Internal libraries imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "american-reverse",
   "metadata": {},
   "outputs": [],
   "source": [
    "import wntr_WSN\n",
    "import networkx_graph as ng\n",
    "import linear_regression as lr"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "headed-witness",
   "metadata": {},
   "source": [
    "## Input data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "german-injury",
   "metadata": {},
   "outputs": [],
   "source": [
    "method = 'method'                              # for file naming\n",
    "leak_node = 'J64'                             # node that's leaking\n",
    "comment = ''                                   # last word(s) in output filename\n",
    "leak_start = 48                                # in hours\n",
    "leak_end = 73                                  # in hours\n",
    "leak_area = 0.0005\n",
    "sim_duration = 130                             # total duration of the simulation"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "amateur-southwest",
   "metadata": {},
   "source": [
    "## Model each node with linear regression"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "derived-message",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Modeling the network with linreg (linear_regression.py - model_network_with_linregression())\n",
      "Creating a NetworkX graph from .inp file (networkx_graph.py - create_graph())\n",
      "Running the simulation (wntr_WSN.py - get_sim_results())\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\users\\pc\\anaconda3\\envs\\wodociagi\\lib\\site-packages\\scipy\\optimize\\minpack.py:828: OptimizeWarning: Covariance of the parameters could not be estimated\n",
      "  warnings.warn('Covariance of the parameters could not be estimated',\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\tModeling junction J1: (1 out of 309)\n",
      "\tModeling junction J2: (2 out of 309)\n",
      "\tModeling junction J3: (3 out of 309)\n",
      "\tModeling junction J4: (4 out of 309)\n",
      "\tModeling junction J5: (5 out of 309)\n",
      "\tModeling junction J6: (6 out of 309)\n",
      "\tModeling junction J7: (7 out of 309)\n",
      "\tModeling junction J8: (8 out of 309)\n",
      "\tModeling junction J9: (9 out of 309)\n",
      "\tModeling junction J10: (10 out of 309)\n",
      "\tModeling junction J11: (11 out of 309)\n",
      "\tModeling junction J12: (12 out of 309)\n",
      "\tModeling junction J13: (13 out of 309)\n",
      "\tModeling junction J14: (14 out of 309)\n",
      "\tModeling junction J15: (15 out of 309)\n",
      "\tModeling junction J16: (16 out of 309)\n",
      "\tModeling junction J17: (17 out of 309)\n",
      "\tModeling junction J18: (18 out of 309)\n",
      "\tModeling junction J19: (19 out of 309)\n",
      "\tModeling junction J20: (20 out of 309)\n",
      "\tModeling junction J21: (21 out of 309)\n",
      "\tModeling junction J22: (22 out of 309)\n",
      "\tModeling junction J24: (23 out of 309)\n",
      "\tModeling junction J25: (24 out of 309)\n",
      "\tModeling junction J26: (25 out of 309)\n",
      "\tModeling junction J27: (26 out of 309)\n",
      "\tModeling junction J28: (27 out of 309)\n",
      "\tModeling junction J29: (28 out of 309)\n",
      "\tModeling junction J30: (29 out of 309)\n",
      "\tModeling junction J31: (30 out of 309)\n",
      "\tModeling junction J32: (31 out of 309)\n",
      "\tModeling junction J33: (32 out of 309)\n",
      "\tModeling junction J34: (33 out of 309)\n",
      "\tModeling junction J35: (34 out of 309)\n",
      "\tModeling junction J36: (35 out of 309)\n",
      "\tModeling junction J37: (36 out of 309)\n",
      "\tModeling junction J38: (37 out of 309)\n",
      "\tModeling junction J39: (38 out of 309)\n",
      "\tModeling junction J40: (39 out of 309)\n",
      "\tModeling junction J41: (40 out of 309)\n",
      "\tModeling junction J42: (41 out of 309)\n",
      "\tModeling junction J43: (42 out of 309)\n",
      "\tModeling junction J44: (43 out of 309)\n",
      "\tModeling junction J45: (44 out of 309)\n",
      "\tModeling junction J46: (45 out of 309)\n",
      "\tModeling junction J47: (46 out of 309)\n",
      "\tModeling junction J48: (47 out of 309)\n",
      "\tModeling junction J49: (48 out of 309)\n",
      "\tModeling junction J50: (49 out of 309)\n",
      "\tModeling junction J52: (50 out of 309)\n",
      "\tModeling junction J53: (51 out of 309)\n",
      "\tModeling junction J54: (52 out of 309)\n",
      "\tModeling junction J56: (53 out of 309)\n",
      "\tModeling junction J57: (54 out of 309)\n",
      "\tModeling junction J58: (55 out of 309)\n",
      "\tModeling junction J59: (56 out of 309)\n",
      "\tModeling junction J60: (57 out of 309)\n",
      "\tModeling junction J61: (58 out of 309)\n",
      "\tModeling junction J62: (59 out of 309)\n",
      "\tModeling junction J63: (60 out of 309)\n",
      "\tModeling junction J64: (61 out of 309)\n",
      "\tModeling junction J65: (62 out of 309)\n",
      "\tModeling junction J66: (63 out of 309)\n",
      "\tModeling junction J67: (64 out of 309)\n",
      "\tModeling junction J68: (65 out of 309)\n",
      "\tModeling junction J69: (66 out of 309)\n",
      "\tModeling junction J70: (67 out of 309)\n",
      "\tModeling junction J71: (68 out of 309)\n",
      "\tModeling junction J72: (69 out of 309)\n",
      "\tModeling junction J73: (70 out of 309)\n",
      "\tModeling junction J75: (71 out of 309)\n",
      "\tModeling junction J76: (72 out of 309)\n",
      "\tModeling junction J77: (73 out of 309)\n",
      "\tModeling junction J78: (74 out of 309)\n",
      "\tModeling junction J79: (75 out of 309)\n",
      "\tModeling junction J80: (76 out of 309)\n",
      "\tModeling junction J81: (77 out of 309)\n",
      "\tModeling junction J82: (78 out of 309)\n",
      "\tModeling junction J83: (79 out of 309)\n",
      "\tModeling junction J84: (80 out of 309)\n",
      "\tModeling junction J85: (81 out of 309)\n",
      "\tModeling junction J86: (82 out of 309)\n",
      "\tModeling junction J87: (83 out of 309)\n",
      "\tModeling junction J88: (84 out of 309)\n",
      "\tModeling junction J89: (85 out of 309)\n",
      "\tModeling junction J90: (86 out of 309)\n",
      "\tModeling junction J91: (87 out of 309)\n",
      "\tModeling junction J92: (88 out of 309)\n",
      "\tModeling junction J93: (89 out of 309)\n",
      "\tModeling junction J94: (90 out of 309)\n",
      "\tModeling junction J95: (91 out of 309)\n",
      "\tModeling junction J96: (92 out of 309)\n",
      "\tModeling junction J97: (93 out of 309)\n",
      "\tModeling junction J98: (94 out of 309)\n",
      "\tModeling junction J99: (95 out of 309)\n",
      "\tModeling junction J100: (96 out of 309)\n",
      "\tModeling junction J102: (97 out of 309)\n",
      "\tModeling junction J103: (98 out of 309)\n",
      "\tModeling junction J104: (99 out of 309)\n",
      "\tModeling junction J105: (100 out of 309)\n",
      "\tModeling junction J106: (101 out of 309)\n",
      "\tModeling junction J107: (102 out of 309)\n",
      "\tModeling junction J108: (103 out of 309)\n",
      "\tModeling junction J109: (104 out of 309)\n",
      "\tModeling junction J110: (105 out of 309)\n",
      "\tModeling junction J111: (106 out of 309)\n",
      "\tModeling junction J112: (107 out of 309)\n",
      "\tModeling junction J113: (108 out of 309)\n",
      "\tModeling junction J114: (109 out of 309)\n",
      "\tModeling junction J115: (110 out of 309)\n",
      "\tModeling junction J116: (111 out of 309)\n",
      "\tModeling junction J117: (112 out of 309)\n",
      "\tModeling junction J118: (113 out of 309)\n",
      "\tModeling junction J119: (114 out of 309)\n",
      "\tModeling junction J120: (115 out of 309)\n",
      "\tModeling junction J121: (116 out of 309)\n",
      "\tModeling junction J122: (117 out of 309)\n",
      "\tModeling junction J123: (118 out of 309)\n",
      "\tModeling junction J124: (119 out of 309)\n",
      "\tModeling junction J125: (120 out of 309)\n",
      "\tModeling junction J126: (121 out of 309)\n",
      "\tModeling junction J127: (122 out of 309)\n",
      "\tModeling junction J128: (123 out of 309)\n",
      "\tModeling junction J129: (124 out of 309)\n",
      "\tModeling junction J130: (125 out of 309)\n",
      "\tModeling junction J131: (126 out of 309)\n",
      "\tModeling junction J132: (127 out of 309)\n",
      "\tModeling junction J133: (128 out of 309)\n",
      "\tModeling junction J134: (129 out of 309)\n",
      "\tModeling junction J135: (130 out of 309)\n",
      "\tModeling junction J136: (131 out of 309)\n",
      "\tModeling junction J137: (132 out of 309)\n",
      "\tModeling junction J138: (133 out of 309)\n",
      "\tModeling junction J139: (134 out of 309)\n",
      "\tModeling junction J140: (135 out of 309)\n",
      "\tModeling junction J141: (136 out of 309)\n",
      "\tModeling junction J142: (137 out of 309)\n",
      "\tModeling junction J143: (138 out of 309)\n",
      "\tModeling junction J144: (139 out of 309)\n",
      "\tModeling junction J145: (140 out of 309)\n",
      "\tModeling junction J146: (141 out of 309)\n",
      "\tModeling junction J147: (142 out of 309)\n",
      "\tModeling junction J148: (143 out of 309)\n",
      "\tModeling junction J149: (144 out of 309)\n",
      "\tModeling junction J150: (145 out of 309)\n",
      "\tModeling junction J151: (146 out of 309)\n",
      "\tModeling junction J153: (147 out of 309)\n",
      "\tModeling junction J154: (148 out of 309)\n",
      "\tModeling junction J155: (149 out of 309)\n",
      "\tModeling junction J156: (150 out of 309)\n",
      "\tModeling junction J157: (151 out of 309)\n",
      "\tModeling junction J158: (152 out of 309)\n",
      "\tModeling junction J159: (153 out of 309)\n",
      "\tModeling junction J160: (154 out of 309)\n",
      "\tModeling junction J161: (155 out of 309)\n",
      "\tModeling junction J162: (156 out of 309)\n",
      "\tModeling junction J163: (157 out of 309)\n",
      "\tModeling junction J164: (158 out of 309)\n",
      "\tModeling junction J165: (159 out of 309)\n",
      "\tModeling junction J166: (160 out of 309)\n",
      "\tModeling junction J167: (161 out of 309)\n",
      "\tModeling junction J168: (162 out of 309)\n",
      "\tModeling junction J169: (163 out of 309)\n",
      "\tModeling junction J170: (164 out of 309)\n",
      "\tModeling junction J171: (165 out of 309)\n",
      "\tModeling junction J172: (166 out of 309)\n",
      "\tModeling junction J173: (167 out of 309)\n",
      "\tModeling junction J175: (168 out of 309)\n",
      "\tModeling junction J176: (169 out of 309)\n",
      "\tModeling junction J177: (170 out of 309)\n",
      "\tModeling junction J179: (171 out of 309)\n",
      "\tModeling junction J180: (172 out of 309)\n",
      "\tModeling junction J181: (173 out of 309)\n",
      "\tModeling junction J182: (174 out of 309)\n",
      "\tModeling junction J183: (175 out of 309)\n",
      "\tModeling junction J184: (176 out of 309)\n",
      "\tModeling junction J185: (177 out of 309)\n",
      "\tModeling junction J186: (178 out of 309)\n",
      "\tModeling junction J187: (179 out of 309)\n",
      "\tModeling junction J188: (180 out of 309)\n",
      "\tModeling junction J189: (181 out of 309)\n",
      "\tModeling junction J190: (182 out of 309)\n",
      "\tModeling junction J191: (183 out of 309)\n",
      "\tModeling junction J192: (184 out of 309)\n",
      "\tModeling junction J193: (185 out of 309)\n",
      "\tModeling junction J194: (186 out of 309)\n",
      "\tModeling junction J195: (187 out of 309)\n",
      "\tModeling junction J196: (188 out of 309)\n",
      "\tModeling junction J197: (189 out of 309)\n",
      "\tModeling junction J198: (190 out of 309)\n",
      "\tModeling junction J199: (191 out of 309)\n",
      "\tModeling junction J200: (192 out of 309)\n",
      "\tModeling junction J202: (193 out of 309)\n",
      "\tModeling junction J203: (194 out of 309)\n",
      "\tModeling junction J204: (195 out of 309)\n",
      "\tModeling junction J205: (196 out of 309)\n",
      "\tModeling junction J206: (197 out of 309)\n",
      "\tModeling junction J207: (198 out of 309)\n",
      "\tModeling junction J208: (199 out of 309)\n",
      "\tModeling junction J209: (200 out of 309)\n",
      "\tModeling junction J210: (201 out of 309)\n",
      "\tModeling junction J211: (202 out of 309)\n",
      "\tModeling junction J212: (203 out of 309)\n",
      "\tModeling junction J213: (204 out of 309)\n",
      "\tModeling junction J214: (205 out of 309)\n",
      "\tModeling junction J216: (206 out of 309)\n",
      "\tModeling junction J218: (207 out of 309)\n",
      "\tModeling junction J219: (208 out of 309)\n",
      "\tModeling junction J220: (209 out of 309)\n",
      "\tModeling junction J223: (210 out of 309)\n",
      "\tModeling junction J224: (211 out of 309)\n",
      "\tModeling junction J225: (212 out of 309)\n",
      "\tModeling junction J226: (213 out of 309)\n",
      "\tModeling junction J227: (214 out of 309)\n",
      "\tModeling junction J228: (215 out of 309)\n",
      "\tModeling junction J229: (216 out of 309)\n",
      "\tModeling junction J230: (217 out of 309)\n",
      "\tModeling junction J231: (218 out of 309)\n",
      "\tModeling junction J232: (219 out of 309)\n",
      "\tModeling junction J233: (220 out of 309)\n",
      "\tModeling junction J235: (221 out of 309)\n",
      "\tModeling junction J236: (222 out of 309)\n",
      "\tModeling junction J239: (223 out of 309)\n",
      "\tModeling junction J240: (224 out of 309)\n",
      "\tModeling junction J241: (225 out of 309)\n",
      "\tModeling junction J242: (226 out of 309)\n",
      "\tModeling junction J243: (227 out of 309)\n",
      "\tModeling junction J244: (228 out of 309)\n",
      "\tModeling junction J245: (229 out of 309)\n",
      "\tModeling junction J246: (230 out of 309)\n",
      "\tModeling junction J247: (231 out of 309)\n",
      "\tModeling junction J248: (232 out of 309)\n",
      "\tModeling junction J249: (233 out of 309)\n",
      "\tModeling junction J250: (234 out of 309)\n",
      "\tModeling junction J251: (235 out of 309)\n",
      "\tModeling junction J252: (236 out of 309)\n",
      "\tModeling junction J253: (237 out of 309)\n",
      "\tModeling junction J254: (238 out of 309)\n",
      "\tModeling junction J256: (239 out of 309)\n",
      "\tModeling junction J257: (240 out of 309)\n",
      "\tModeling junction J258: (241 out of 309)\n",
      "\tModeling junction J259: (242 out of 309)\n",
      "\tModeling junction J260: (243 out of 309)\n",
      "\tModeling junction J261: (244 out of 309)\n",
      "\tModeling junction J262: (245 out of 309)\n",
      "\tModeling junction J263: (246 out of 309)\n",
      "\tModeling junction J264: (247 out of 309)\n",
      "\tModeling junction J265: (248 out of 309)\n",
      "\tModeling junction J267: (249 out of 309)\n",
      "\tModeling junction J269: (250 out of 309)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\tModeling junction J273: (251 out of 309)\n",
      "\tModeling junction J274: (252 out of 309)\n",
      "\tModeling junction J275: (253 out of 309)\n",
      "\tModeling junction J276: (254 out of 309)\n",
      "\tModeling junction J277: (255 out of 309)\n",
      "\tModeling junction J278: (256 out of 309)\n",
      "\tModeling junction J279: (257 out of 309)\n",
      "\tModeling junction J280: (258 out of 309)\n",
      "\tModeling junction J281: (259 out of 309)\n",
      "\tModeling junction J282: (260 out of 309)\n",
      "\tModeling junction J283: (261 out of 309)\n",
      "\tModeling junction J284: (262 out of 309)\n",
      "\tModeling junction J285: (263 out of 309)\n",
      "\tModeling junction J286: (264 out of 309)\n",
      "\tModeling junction J287: (265 out of 309)\n",
      "\tModeling junction J288: (266 out of 309)\n",
      "\tModeling junction J289: (267 out of 309)\n",
      "\tModeling junction J290: (268 out of 309)\n",
      "\tModeling junction J291: (269 out of 309)\n",
      "\tModeling junction J292: (270 out of 309)\n",
      "\tModeling junction J293: (271 out of 309)\n",
      "\tModeling junction J294: (272 out of 309)\n",
      "\tModeling junction J295: (273 out of 309)\n",
      "\tModeling junction J296: (274 out of 309)\n",
      "\tModeling junction J297: (275 out of 309)\n",
      "\tModeling junction J298: (276 out of 309)\n",
      "\tModeling junction J299: (277 out of 309)\n",
      "\tModeling junction J300: (278 out of 309)\n",
      "\tModeling junction J301: (279 out of 309)\n",
      "\tModeling junction J302: (280 out of 309)\n",
      "\tModeling junction J303: (281 out of 309)\n",
      "\tModeling junction J304: (282 out of 309)\n",
      "\tModeling junction J305: (283 out of 309)\n",
      "\tModeling junction J306: (284 out of 309)\n",
      "\tModeling junction J307: (285 out of 309)\n",
      "\tModeling junction J308: (286 out of 309)\n",
      "\tModeling junction J309: (287 out of 309)\n",
      "\tModeling junction J310: (288 out of 309)\n",
      "\tModeling junction J311: (289 out of 309)\n",
      "\tModeling junction J312: (290 out of 309)\n",
      "\tModeling junction J313: (291 out of 309)\n",
      "\tModeling junction J314: (292 out of 309)\n",
      "\tModeling junction J315: (293 out of 309)\n",
      "\tModeling junction J316: (294 out of 309)\n",
      "\tModeling junction J317: (295 out of 309)\n",
      "\tModeling junction J318: (296 out of 309)\n",
      "\tModeling junction J319: (297 out of 309)\n",
      "\tModeling junction J320: (298 out of 309)\n",
      "\tModeling junction J321: (299 out of 309)\n",
      "\tModeling junction J322: (300 out of 309)\n",
      "\tModeling junction J323: (301 out of 309)\n",
      "\tModeling junction J324: (302 out of 309)\n",
      "\tModeling junction J326: (303 out of 309)\n",
      "\tModeling junction J327: (304 out of 309)\n",
      "\tModeling junction J539: (305 out of 309)\n",
      "\tModeling junction J544: (306 out of 309)\n",
      "\tModeling junction J546: (307 out of 309)\n",
      "\tModeling junction T1: (308 out of 309)\n",
      "\tModeling junction T2: (309 out of 309)\n"
     ]
    }
   ],
   "source": [
    "linreg_models, normalization_params = lr.model_network_with_linreg(n=3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "surface-taiwan",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[{'node': 'J1',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0003094505413836043,\n",
       "  'r2': 0.999779023346383,\n",
       "  'sensors': ['J284', 'P9', 'P460']},\n",
       " {'node': 'J2',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.2718571746158296e-05,\n",
       "  'r2': 0.9999860527644338,\n",
       "  'sensors': ['J284', 'P460', 'P9']},\n",
       " {'node': 'J3',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.4049235727063265e-05,\n",
       "  'r2': 0.9999835140029287,\n",
       "  'sensors': ['J284', 'P460', 'P9']},\n",
       " {'node': 'J4',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.979415354764518e-05,\n",
       "  'r2': 0.9999437589457921,\n",
       "  'sensors': ['J284', 'P460', 'P9']},\n",
       " {'node': 'J5',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.921094605668927e-05,\n",
       "  'r2': 0.9999575212346724,\n",
       "  'sensors': ['J284', 'P460', 'P9']},\n",
       " {'node': 'J6',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.001315916827095419,\n",
       "  'r2': 0.9993899974078964,\n",
       "  'sensors': ['P9', 'J79', 'P500']},\n",
       " {'node': 'J7',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00027063805115237587,\n",
       "  'r2': 0.9998813887377228,\n",
       "  'sensors': ['J284', 'P9', 'P460']},\n",
       " {'node': 'J8',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00025383092036565904,\n",
       "  'r2': 0.9998484476306594,\n",
       "  'sensors': ['J284', 'P9', 'P460']},\n",
       " {'node': 'J9',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00031276677981475684,\n",
       "  'r2': 0.9998625916900282,\n",
       "  'sensors': ['J284', 'P9', 'P460']},\n",
       " {'node': 'J10',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.000967601264886546,\n",
       "  'r2': 0.9994893211964133,\n",
       "  'sensors': ['P9', 'J79', 'P500']},\n",
       " {'node': 'J11',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0005484068847349475,\n",
       "  'r2': 0.9996674336502284,\n",
       "  'sensors': ['P9', 'J12', 'J79']},\n",
       " {'node': 'J12',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.048709793414476e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J12', 'P9', 'P13']},\n",
       " {'node': 'J13',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 9.996649036817023e-06,\n",
       "  'r2': 0.9999950194358924,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J14',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.483175393701755e-05,\n",
       "  'r2': 0.9999786597753552,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J15',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.0926223192945634e-06,\n",
       "  'r2': 0.999999399201519,\n",
       "  'sensors': ['P13', 'J17', 'J12']},\n",
       " {'node': 'J16',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.200021502730398e-07,\n",
       "  'r2': 0.9999997178383436,\n",
       "  'sensors': ['P13', 'J17', 'J12']},\n",
       " {'node': 'J17',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.048709793414476e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J17', 'P13', 'J12']},\n",
       " {'node': 'J18',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.215422198335973e-06,\n",
       "  'r2': 0.9999952246654162,\n",
       "  'sensors': ['J17', 'P13', 'P536']},\n",
       " {'node': 'J19',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.1712982609723712e-05,\n",
       "  'r2': 0.9999887174591027,\n",
       "  'sensors': ['J17', 'P536', 'P13']},\n",
       " {'node': 'J20',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.380061830000289e-05,\n",
       "  'r2': 0.9999860933823945,\n",
       "  'sensors': ['J17', 'P536', 'P13']},\n",
       " {'node': 'J21',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.474717566267128e-06,\n",
       "  'r2': 0.9999983454320741,\n",
       "  'sensors': ['P536', 'J163', 'J17']},\n",
       " {'node': 'J22',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.047445567390268e-07,\n",
       "  'r2': 0.999999529794353,\n",
       "  'sensors': ['P536', 'J163', 'J17']},\n",
       " {'node': 'J24',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.001344239201029277,\n",
       "  'r2': 0.9993863050755517,\n",
       "  'sensors': ['P539', 'J539', 'P460']},\n",
       " {'node': 'J25',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0004703120680675365,\n",
       "  'r2': 0.9998767758046666,\n",
       "  'sensors': ['P22', 'J316', 'J28']},\n",
       " {'node': 'J26',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.455982471323399e-05,\n",
       "  'r2': 0.9999773284191099,\n",
       "  'sensors': ['P22', 'J28', 'J316']},\n",
       " {'node': 'J27',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 9.450851838671285e-05,\n",
       "  'r2': 0.9999639369722124,\n",
       "  'sensors': ['P22', 'J28', 'J316']},\n",
       " {'node': 'J28',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.048709793414476e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J28', 'P22', 'J316']},\n",
       " {'node': 'J29',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.242931142569591e-05,\n",
       "  'r2': 0.999972856764367,\n",
       "  'sensors': ['J28', 'P22', 'P27']},\n",
       " {'node': 'J30',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0002447535290838555,\n",
       "  'r2': 0.9999190695352701,\n",
       "  'sensors': ['J28', 'P27', 'P22']},\n",
       " {'node': 'J31',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.5537138877471384e-06,\n",
       "  'r2': 0.9999990302548533,\n",
       "  'sensors': ['P27', 'J33', 'J28']},\n",
       " {'node': 'J32',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.049688047313418e-07,\n",
       "  'r2': 0.9999997140523526,\n",
       "  'sensors': ['P27', 'J33', 'J28']},\n",
       " {'node': 'J33',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J33', 'P27', 'J151']},\n",
       " {'node': 'J34',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.5314859182154094e-05,\n",
       "  'r2': 0.9999900855362918,\n",
       "  'sensors': ['J316', 'P22', 'P520']},\n",
       " {'node': 'J35',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 9.274934037450648e-07,\n",
       "  'r2': 0.999999724147211,\n",
       "  'sensors': ['J316', 'P520', 'P22']},\n",
       " {'node': 'J36',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.150771469279833e-06,\n",
       "  'r2': 0.9999969550047124,\n",
       "  'sensors': ['P520', 'J316', 'J39']},\n",
       " {'node': 'J37',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.1017803536520094e-05,\n",
       "  'r2': 0.9999965245450793,\n",
       "  'sensors': ['P520', 'J316', 'J39']},\n",
       " {'node': 'J38',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.038761426132093e-05,\n",
       "  'r2': 0.9999640228569393,\n",
       "  'sensors': ['P520', 'J39', 'P478']},\n",
       " {'node': 'J39',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J40',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.341309032537044e-05,\n",
       "  'r2': 0.999990051810908,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J41',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.0386293560323352e-05,\n",
       "  'r2': 0.9999952096313989,\n",
       "  'sensors': ['P174', 'J172', 'J39']},\n",
       " {'node': 'J42',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.597408250328944e-05,\n",
       "  'r2': 0.9999876819423837,\n",
       "  'sensors': ['J43', 'P174', 'P511']},\n",
       " {'node': 'J43',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J43', 'P511', 'P46']},\n",
       " {'node': 'J44',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.6557949198848845e-05,\n",
       "  'r2': 0.9999930122698921,\n",
       "  'sensors': ['J43', 'P511', 'P46']},\n",
       " {'node': 'J45',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.5749172626485926e-08,\n",
       "  'r2': 0.9999999772741769,\n",
       "  'sensors': ['P46', 'J62', 'P55']},\n",
       " {'node': 'J46',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.969693376383129e-05,\n",
       "  'r2': 0.9999909803158696,\n",
       "  'sensors': ['J43', 'P511', 'P174']},\n",
       " {'node': 'J47',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.4760331993319954e-05,\n",
       "  'r2': 0.9999873990027696,\n",
       "  'sensors': ['J172', 'J43', 'P511']},\n",
       " {'node': 'J48',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.9146040974762316e-05,\n",
       "  'r2': 0.9999894274155003,\n",
       "  'sensors': ['J172', 'P174', 'J43']},\n",
       " {'node': 'J49',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.2123438565381882e-05,\n",
       "  'r2': 0.9999957688989652,\n",
       "  'sensors': ['J137', 'P43', 'J33']},\n",
       " {'node': 'J50',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00018843046168956666,\n",
       "  'r2': 0.9999252650610226,\n",
       "  'sensors': ['J137', 'P43', 'J33']},\n",
       " {'node': 'J52',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.048709793414476e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J53',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.7661034194382057e-14,\n",
       "  'r2': 0.9999999999999892,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J54',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.7480319128711703e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J54', 'P46', 'J62']},\n",
       " {'node': 'J56',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.0359445957261123e-09,\n",
       "  'r2': 0.9999999994157706,\n",
       "  'sensors': ['P52', 'J60', 'J54']},\n",
       " {'node': 'J57',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.661074497767976e-10,\n",
       "  'r2': 0.9999999997620511,\n",
       "  'sensors': ['P52', 'J60', 'J54']},\n",
       " {'node': 'J58',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.58141664203778e-11,\n",
       "  'r2': 0.9999999999694438,\n",
       "  'sensors': ['P52', 'J60', 'J54']},\n",
       " {'node': 'J59',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.6219603706453107e-11,\n",
       "  'r2': 0.9999999999750315,\n",
       "  'sensors': ['P52', 'J60', 'J54']},\n",
       " {'node': 'J60',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.048709793414476e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J60', 'P52', 'P520']},\n",
       " {'node': 'J61',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.6526170816113934e-11,\n",
       "  'r2': 0.9999999999887017,\n",
       "  'sensors': ['J60', 'P52', 'P520']},\n",
       " {'node': 'J62',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J62', 'P55', 'P204']},\n",
       " {'node': 'J63',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.738329496805738e-08,\n",
       "  'r2': 0.9999999867357642,\n",
       "  'sensors': ['P56', 'J64', 'J62']},\n",
       " {'node': 'J64',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.048709793414476e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J64', 'P56', 'J129']},\n",
       " {'node': 'J65',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.306926291360694e-07,\n",
       "  'r2': 0.9999999427983561,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J66',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.986539301629791e-05,\n",
       "  'r2': 0.999972402593206,\n",
       "  'sensors': ['J52', 'P45', 'J67']},\n",
       " {'node': 'J67',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J67', 'P459', 'J79']},\n",
       " {'node': 'J68',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.63270435436336e-08,\n",
       "  'r2': 0.9999999936950821,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J69',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.542261073017471e-07,\n",
       "  'r2': 0.9999999351382833,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J70',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.652817355353378e-05,\n",
       "  'r2': 0.999972532717188,\n",
       "  'sensors': ['J33', 'J151', 'P213']},\n",
       " {'node': 'J71',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.074294186869159e-05,\n",
       "  'r2': 0.9999768587382079,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J72',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.3936046208056114e-05,\n",
       "  'r2': 0.9999843631617031,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J73',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.761208637188742e-05,\n",
       "  'r2': 0.9999868950443915,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J75',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0004920381264396507,\n",
       "  'r2': 0.9997624474931932,\n",
       "  'sensors': ['J79', 'P500', 'P9']},\n",
       " {'node': 'J76',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00022123445651693212,\n",
       "  'r2': 0.9999016598670031,\n",
       "  'sensors': ['J79', 'P500', 'P9']},\n",
       " {'node': 'J77',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.807566774882057e-05,\n",
       "  'r2': 0.9999645909462678,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J78',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.000400787448541e-06,\n",
       "  'r2': 0.9999966849337295,\n",
       "  'sensors': ['J67', 'P459', 'J79']},\n",
       " {'node': 'J79',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J80',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.6389481979971035e-07,\n",
       "  'r2': 0.9999998641336237,\n",
       "  'sensors': ['J79', 'P500', 'J67']},\n",
       " {'node': 'J81',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0004853892451705881,\n",
       "  'r2': 0.9997681738830559,\n",
       "  'sensors': ['J79', 'P500', 'P9']},\n",
       " {'node': 'J82',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00038742771930648506,\n",
       "  'r2': 0.9998615227023512,\n",
       "  'sensors': ['J79', 'P500', 'P9']},\n",
       " {'node': 'J83',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00048677361871871273,\n",
       "  'r2': 0.9997450444132178,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J84',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0007524665196843946,\n",
       "  'r2': 0.9996546519505071,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J85',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0005603066815644632,\n",
       "  'r2': 0.9996672123097808,\n",
       "  'sensors': ['J79', 'P500', 'P9']},\n",
       " {'node': 'J86',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0004615575009067705,\n",
       "  'r2': 0.9997128219543098,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J87',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.128228276735326e-05,\n",
       "  'r2': 0.999989725619763,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J88',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.4586056543202752e-05,\n",
       "  'r2': 0.9999838142540389,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J89',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.0211744221270188e-05,\n",
       "  'r2': 0.9999851337384364,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J90',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.2153940956599702e-05,\n",
       "  'r2': 0.9999902980172007,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J91',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.653365071394493e-05,\n",
       "  'r2': 0.9999860680913212,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J92',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.73626627887566e-05,\n",
       "  'r2': 0.999989718326059,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J93',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.3593782483951216e-05,\n",
       "  'r2': 0.9999886983816298,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J94',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.7355506703232603e-05,\n",
       "  'r2': 0.9999868966962034,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J95',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.9331022199334736e-08,\n",
       "  'r2': 0.9999999724915982,\n",
       "  'sensors': ['J12', 'P9', 'P13']},\n",
       " {'node': 'J96',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.309668704688895e-06,\n",
       "  'r2': 0.9999980420153904,\n",
       "  'sensors': ['P13', 'J17', 'J67']},\n",
       " {'node': 'J97',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00035466851793021417,\n",
       "  'r2': 0.9998272250419779,\n",
       "  'sensors': ['P13', 'J67', 'P459']},\n",
       " {'node': 'J98',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00017560875197262645,\n",
       "  'r2': 0.9999006489881112,\n",
       "  'sensors': ['P13', 'J67', 'P459']},\n",
       " {'node': 'J99',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00022238984588480244,\n",
       "  'r2': 0.9998791159674647,\n",
       "  'sensors': ['J67', 'P459', 'P13']},\n",
       " {'node': 'J100',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00016517599050054366,\n",
       "  'r2': 0.9999284166886067,\n",
       "  'sensors': ['J67', 'P459', 'P13']},\n",
       " {'node': 'J102',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0001785776938072101,\n",
       "  'r2': 0.9999041680781081,\n",
       "  'sensors': ['J67', 'P459', 'P536']},\n",
       " {'node': 'J103',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.021286957296424e-05,\n",
       "  'r2': 0.9999860682382572,\n",
       "  'sensors': ['J67', 'P459', 'J163']},\n",
       " {'node': 'J104',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.9080180027920025e-05,\n",
       "  'r2': 0.9999899143967654,\n",
       "  'sensors': ['J163', 'P536', 'J172']},\n",
       " {'node': 'J105',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 9.697300368099563e-06,\n",
       "  'r2': 0.9999952349201447,\n",
       "  'sensors': ['J172', 'J163', 'P174']},\n",
       " {'node': 'J106',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.3157605087875871e-05,\n",
       "  'r2': 0.9999946465665996,\n",
       "  'sensors': ['J172', 'P174', 'J163']},\n",
       " {'node': 'J107',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.5411947363498066e-05,\n",
       "  'r2': 0.999991522314046,\n",
       "  'sensors': ['J172', 'P174', 'J43']},\n",
       " {'node': 'J108',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.571135967340238e-05,\n",
       "  'r2': 0.9999751829441477,\n",
       "  'sensors': ['J67', 'P459', 'J151']},\n",
       " {'node': 'J109',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 9.149234541168969e-05,\n",
       "  'r2': 0.9999654218219229,\n",
       "  'sensors': ['J67', 'P459', 'J151']},\n",
       " {'node': 'J110',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0008402155446009716,\n",
       "  'r2': 0.999644897630201,\n",
       "  'sensors': ['J151', 'P213', 'J67']},\n",
       " {'node': 'J111',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0017742277831137498,\n",
       "  'r2': 0.9991663296527359,\n",
       "  'sensors': ['J151', 'P213', 'J113']},\n",
       " {'node': 'J112',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.5801186667772643e-05,\n",
       "  'r2': 0.9999938076937248,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J113',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J114',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.7184292322373657e-05,\n",
       "  'r2': 0.9999926430013888,\n",
       "  'sensors': ['J113', 'P110', 'J39']},\n",
       " {'node': 'J115',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00013176046973010404,\n",
       "  'r2': 0.9999270764941358,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J116',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.3386681369535347e-05,\n",
       "  'r2': 0.9999841491751876,\n",
       "  'sensors': ['J113', 'P110', 'J146']},\n",
       " {'node': 'J117',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.228497121274496e-05,\n",
       "  'r2': 0.9999815515183552,\n",
       "  'sensors': ['J113', 'P110', 'J146']},\n",
       " {'node': 'J118',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00010806236239460373,\n",
       "  'r2': 0.999949176696505,\n",
       "  'sensors': ['J172', 'P174', 'J43']},\n",
       " {'node': 'J119',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.195731706830594e-06,\n",
       "  'r2': 0.9999977189424392,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J120',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.4829818524125973e-05,\n",
       "  'r2': 0.9999883396906001,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J121',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.127777714013021e-06,\n",
       "  'r2': 0.9999991729743282,\n",
       "  'sensors': ['P174', 'J172', 'J39']},\n",
       " {'node': 'J122',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.704397754330465e-06,\n",
       "  'r2': 0.9999986305045072,\n",
       "  'sensors': ['P174', 'J172', 'J39']},\n",
       " {'node': 'J123',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.371526028702894e-08,\n",
       "  'r2': 0.9999999754153427,\n",
       "  'sensors': ['P56', 'J54', 'J64']},\n",
       " {'node': 'J124',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.731518000862999e-08,\n",
       "  'r2': 0.9999999668702495,\n",
       "  'sensors': ['P56', 'J54', 'J64']},\n",
       " {'node': 'J125',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.5162991585660086e-08,\n",
       "  'r2': 0.9999999773996789,\n",
       "  'sensors': ['P56', 'J54', 'J64']},\n",
       " {'node': 'J126',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.405321986057694e-08,\n",
       "  'r2': 0.9999999606283742,\n",
       "  'sensors': ['P56', 'J54', 'J64']},\n",
       " {'node': 'J127',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.615100456573543e-10,\n",
       "  'r2': 0.999999999779467,\n",
       "  'sensors': ['P244', 'J129', 'J207']},\n",
       " {'node': 'J128',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.568907803918457e-10,\n",
       "  'r2': 0.9999999995740728,\n",
       "  'sensors': ['P244', 'J129', 'J207']},\n",
       " {'node': 'J129',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J129', 'P244', 'J64']},\n",
       " {'node': 'J130',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.401532292806226e-10,\n",
       "  'r2': 0.9999999997907855,\n",
       "  'sensors': ['J129', 'P244', 'J64']},\n",
       " {'node': 'J131',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.70244082005889e-13,\n",
       "  'r2': 0.9999999999997008,\n",
       "  'sensors': ['J60', 'P52', 'P520']},\n",
       " {'node': 'J132',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.934000124674556e-10,\n",
       "  'r2': 0.9999999998428051,\n",
       "  'sensors': ['J28', 'P22', 'J316']},\n",
       " {'node': 'J133',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.832850957115586e-05,\n",
       "  'r2': 0.9999750255374282,\n",
       "  'sensors': ['J28', 'P22', 'P27']},\n",
       " {'node': 'J134',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0003095097491859559,\n",
       "  'r2': 0.9998765205687736,\n",
       "  'sensors': ['J28', 'P27', 'P22']},\n",
       " {'node': 'J135',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.29672905747933e-06,\n",
       "  'r2': 0.9999994844645388,\n",
       "  'sensors': ['P27', 'J33', 'J28']},\n",
       " {'node': 'J136',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.4313702376931591e-06,\n",
       "  'r2': 0.9999994591146651,\n",
       "  'sensors': ['P27', 'J33', 'J28']},\n",
       " {'node': 'J137',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J137', 'P43', 'J52']},\n",
       " {'node': 'J138',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.779358598676585e-06,\n",
       "  'r2': 0.9999979474896441,\n",
       "  'sensors': ['J33', 'J137', 'P43']},\n",
       " {'node': 'J139',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0002814195437604446,\n",
       "  'r2': 0.9999117348951604,\n",
       "  'sensors': ['J33', 'J151', 'P213']},\n",
       " {'node': 'J140',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0006023113726631208,\n",
       "  'r2': 0.9997443515301468,\n",
       "  'sensors': ['J33', 'J151', 'P213']},\n",
       " {'node': 'J141',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.505944606223703e-06,\n",
       "  'r2': 0.9999965284489745,\n",
       "  'sensors': ['J33', 'J146', 'P280']},\n",
       " {'node': 'J142',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.219478244903574e-06,\n",
       "  'r2': 0.9999988624380457,\n",
       "  'sensors': ['J316', 'P520', 'J146']},\n",
       " {'node': 'J143',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.4099884089692123e-06,\n",
       "  'r2': 0.9999994339067554,\n",
       "  'sensors': ['J146', 'P280', 'J316']},\n",
       " {'node': 'J144',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.588691815950113e-07,\n",
       "  'r2': 0.9999996700958405,\n",
       "  'sensors': ['J146', 'P280', 'J316']},\n",
       " {'node': 'J145',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.980850457932984e-07,\n",
       "  'r2': 0.9999997055866862,\n",
       "  'sensors': ['J146', 'P280', 'J33']},\n",
       " {'node': 'J146',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J146', 'P280', 'J316']},\n",
       " {'node': 'J147',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.8659697681588794e-08,\n",
       "  'r2': 0.9999999866413876,\n",
       "  'sensors': ['J146', 'P280', 'J316']},\n",
       " {'node': 'J148',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.772601938369866e-05,\n",
       "  'r2': 0.9999649176170667,\n",
       "  'sensors': ['J67', 'P459', 'J151']},\n",
       " {'node': 'J149',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.560906729011211e-05,\n",
       "  'r2': 0.9999665462536836,\n",
       "  'sensors': ['J67', 'P459', 'J151']},\n",
       " {'node': 'J150',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.930047713829576e-05,\n",
       "  'r2': 0.9999822943507852,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J151',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J153',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.861848159880964e-06,\n",
       "  'r2': 0.9999968425626624,\n",
       "  'sensors': ['J33', 'J146', 'P280']},\n",
       " {'node': 'J154',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 9.516490232640812e-07,\n",
       "  'r2': 0.9999996650201399,\n",
       "  'sensors': ['J146', 'P280', 'J33']},\n",
       " {'node': 'J155',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.1379810757523192e-06,\n",
       "  'r2': 0.9999996245799212,\n",
       "  'sensors': ['J146', 'P280', 'J316']},\n",
       " {'node': 'J156',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.596539176642965e-05,\n",
       "  'r2': 0.9999845002434241,\n",
       "  'sensors': ['J146', 'P280', 'J316']},\n",
       " {'node': 'J157',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.931569736042464e-05,\n",
       "  'r2': 0.9999691309446287,\n",
       "  'sensors': ['J146', 'P280', 'J113']},\n",
       " {'node': 'J158',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.273222453410115e-05,\n",
       "  'r2': 0.999971143908286,\n",
       "  'sensors': ['J146', 'P280', 'J113']},\n",
       " {'node': 'J159',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.517671150765822e-05,\n",
       "  'r2': 0.9999703268295954,\n",
       "  'sensors': ['J146', 'P280', 'J113']},\n",
       " {'node': 'J160',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.9135292115944567e-05,\n",
       "  'r2': 0.9999830902079154,\n",
       "  'sensors': ['J43', 'P174', 'P511']},\n",
       " {'node': 'J161',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.9137965542129806e-05,\n",
       "  'r2': 0.9999858024164527,\n",
       "  'sensors': ['J43', 'P174', 'P511']},\n",
       " {'node': 'J162',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.949742161449265e-06,\n",
       "  'r2': 0.999995904995016,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J163',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.048709793414476e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J164',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.0752485448719343e-05,\n",
       "  'r2': 0.9999904560508865,\n",
       "  'sensors': ['J163', 'P536', 'J172']},\n",
       " {'node': 'J165',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00010247466161095954,\n",
       "  'r2': 0.9999362377585138,\n",
       "  'sensors': ['J163', 'P536', 'J67']},\n",
       " {'node': 'J166',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.998982083636653e-06,\n",
       "  'r2': 0.9999965862781034,\n",
       "  'sensors': ['P536', 'J163', 'J17']},\n",
       " {'node': 'J167',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.790363734959225e-06,\n",
       "  'r2': 0.9999965189638003,\n",
       "  'sensors': ['P536', 'J163', 'J17']},\n",
       " {'node': 'J168',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.3715671560067987e-06,\n",
       "  'r2': 0.9999994443905164,\n",
       "  'sensors': ['J67', 'P459', 'J79']},\n",
       " {'node': 'J169',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0002619620388547669,\n",
       "  'r2': 0.999865879757053,\n",
       "  'sensors': ['J67', 'P459', 'J172']},\n",
       " {'node': 'J170',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00023881378808216763,\n",
       "  'r2': 0.9998895160412041,\n",
       "  'sensors': ['J67', 'P459', 'J172']},\n",
       " {'node': 'J171',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.2758648020466248e-05,\n",
       "  'r2': 0.9999899252742401,\n",
       "  'sensors': ['J172', 'P174', 'J67']},\n",
       " {'node': 'J172',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J172', 'P174', 'J39']},\n",
       " {'node': 'J173',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.000153571583355292,\n",
       "  'r2': 0.9999453392583657,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J175',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.3136624406882335e-07,\n",
       "  'r2': 0.9999998501098465,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J176',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.5098451365345562e-08,\n",
       "  'r2': 0.9999999924958849,\n",
       "  'sensors': ['J137', 'P43', 'J52']},\n",
       " {'node': 'J177',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.728520579272681e-09,\n",
       "  'r2': 0.9999999992615488,\n",
       "  'sensors': ['J137', 'P43', 'J52']},\n",
       " {'node': 'J179',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.06124940272460691,\n",
       "  'r2': 0.9830569244654552,\n",
       "  'sensors': ['J316', 'P22', 'P520']},\n",
       " {'node': 'J180',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0010979294788283807,\n",
       "  'r2': 0.9996328381221977,\n",
       "  'sensors': ['J146', 'P280', 'J316']},\n",
       " {'node': 'J181',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0011399183717909897,\n",
       "  'r2': 0.9996770013306597,\n",
       "  'sensors': ['J316', 'J146', 'P280']},\n",
       " {'node': 'J182',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0015869070006646416,\n",
       "  'r2': 0.9994950089417007,\n",
       "  'sensors': ['P520', 'J316', 'J146']},\n",
       " {'node': 'J183',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.1804041710245325e-05,\n",
       "  'r2': 0.9999958536613037,\n",
       "  'sensors': ['P520', 'J316', 'J39']},\n",
       " {'node': 'J184',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.60901409023244e-05,\n",
       "  'r2': 0.9999871009408673,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J185',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0018663617144833477,\n",
       "  'r2': 0.9992876735153008,\n",
       "  'sensors': ['J113', 'P110', 'J146']},\n",
       " {'node': 'J186',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.014230436836224904,\n",
       "  'r2': 0.9947814050466948,\n",
       "  'sensors': ['J113', 'P110', 'J146']},\n",
       " {'node': 'J187',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0008745789892500552,\n",
       "  'r2': 0.9997251208184567,\n",
       "  'sensors': ['J113', 'P110', 'J146']},\n",
       " {'node': 'J188',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.420157175220146e-06,\n",
       "  'r2': 0.9999971171248861,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J189',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.8815971007243241,\n",
       "  'r2': 0.8119807751479047,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J190',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.8545628474459361,\n",
       "  'r2': 0.770963527923422,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J191',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.07722700440898325,\n",
       "  'r2': 0.9654358994682762,\n",
       "  'sensors': ['J163', 'P536', 'J67']},\n",
       " {'node': 'J192',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.10758643756603131,\n",
       "  'r2': 0.9519645902305385,\n",
       "  'sensors': ['J163', 'P536', 'J67']},\n",
       " {'node': 'J193',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.001472801604997264,\n",
       "  'r2': 0.9993130829600435,\n",
       "  'sensors': ['J67', 'P459', 'P536']},\n",
       " {'node': 'J194',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0004745917343606279,\n",
       "  'r2': 0.9997485539793385,\n",
       "  'sensors': ['J151', 'P213', 'J67']},\n",
       " {'node': 'J195',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.008828452697630695,\n",
       "  'r2': 0.9958881574815306,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J196',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0010926221192865196,\n",
       "  'r2': 0.9995902125929177,\n",
       "  'sensors': ['J151', 'P213', 'J113']},\n",
       " {'node': 'J197',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00012704452800859386,\n",
       "  'r2': 0.9999461212349353,\n",
       "  'sensors': ['P520', 'J60', 'J316']},\n",
       " {'node': 'J198',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0018063947779459063,\n",
       "  'r2': 0.9992514736972985,\n",
       "  'sensors': ['P520', 'J60', 'J316']},\n",
       " {'node': 'J199',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.722586347139571e-05,\n",
       "  'r2': 0.999972420433668,\n",
       "  'sensors': ['J113', 'P110', 'J146']},\n",
       " {'node': 'J200',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.7465239889924205e-05,\n",
       "  'r2': 0.9999887251798776,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J202',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0002511741596746751,\n",
       "  'r2': 0.9999069914343518,\n",
       "  'sensors': ['P520', 'J316', 'J60']},\n",
       " {'node': 'J203',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.057260481991854e-09,\n",
       "  'r2': 0.9999999981719188,\n",
       "  'sensors': ['J54', 'P46', 'J62']},\n",
       " {'node': 'J204',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.321806765492906e-09,\n",
       "  'r2': 0.9999999945534861,\n",
       "  'sensors': ['J54', 'P46', 'J62']},\n",
       " {'node': 'J205',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.473871703450945e-09,\n",
       "  'r2': 0.9999999982000339,\n",
       "  'sensors': ['P244', 'J129', 'J207']},\n",
       " {'node': 'J206',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.189236846384284e-08,\n",
       "  'r2': 0.9999999788433623,\n",
       "  'sensors': ['P244', 'J129', 'J207']},\n",
       " {'node': 'J207',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J207', 'J62', 'P55']},\n",
       " {'node': 'J208',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.4025205443798806e-07,\n",
       "  'r2': 0.99999980441786,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J209',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.7029900095278145e-05,\n",
       "  'r2': 0.9999885648241367,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J210',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.004841150926946122,\n",
       "  'r2': 0.9979098103029428,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J211',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.335715471806339e-05,\n",
       "  'r2': 0.9999713588550363,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J212',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0005960824174549384,\n",
       "  'r2': 0.9997746031953848,\n",
       "  'sensors': ['J33', 'P27', 'J151']},\n",
       " {'node': 'J213',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.3100866382410244e-05,\n",
       "  'r2': 0.9999891477806669,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J214',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.020592192783040395,\n",
       "  'r2': 0.9913849489440388,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J216',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0005622003664559073,\n",
       "  'r2': 0.9996657675534891,\n",
       "  'sensors': ['P536', 'J67', 'P459']},\n",
       " {'node': 'J218',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0005628769287314014,\n",
       "  'r2': 0.9997185452220778,\n",
       "  'sensors': ['P536', 'J163', 'J17']},\n",
       " {'node': 'J219',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.39303460696077974,\n",
       "  'r2': 0.8703839927858528,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J220',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.214878434841415e-06,\n",
       "  'r2': 0.9999965603539998,\n",
       "  'sensors': ['J33', 'J137', 'P43']},\n",
       " {'node': 'J223',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.007433477468291703,\n",
       "  'r2': 0.9966960164920549,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J224',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00046607351107260184,\n",
       "  'r2': 0.9997828979067922,\n",
       "  'sensors': ['J163', 'P536', 'J172']},\n",
       " {'node': 'J225',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.0780424081708886,\n",
       "  'r2': 0.7691903903750713,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J226',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0014417211529271154,\n",
       "  'r2': 0.9993803805964531,\n",
       "  'sensors': ['J43', 'P174', 'P511']},\n",
       " {'node': 'J227',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.004942665850217382,\n",
       "  'r2': 0.9977549026265989,\n",
       "  'sensors': ['J172', 'P174', 'J39']},\n",
       " {'node': 'J228',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0014971720855011135,\n",
       "  'r2': 0.9993218970471414,\n",
       "  'sensors': ['J113', 'P110', 'J151']},\n",
       " {'node': 'J229',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.24467790275339513,\n",
       "  'r2': 0.9054179303699389,\n",
       "  'sensors': ['J172', 'P174', 'J67']},\n",
       " {'node': 'J230',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.045357810409903e-08,\n",
       "  'r2': 0.9999999629587504,\n",
       "  'sensors': ['P56', 'J54', 'J64']},\n",
       " {'node': 'J231',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.2806630704856338e-14,\n",
       "  'r2': 0.9999999999999936,\n",
       "  'sensors': ['J64', 'P56', 'J129']},\n",
       " {'node': 'J232',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.000980818808977671,\n",
       "  'r2': 0.9994938452093411,\n",
       "  'sensors': ['P56', 'J64', 'J62']},\n",
       " {'node': 'J233',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.535423311737599e-08,\n",
       "  'r2': 0.9999999707390954,\n",
       "  'sensors': ['J137', 'P43', 'J33']},\n",
       " {'node': 'J235',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.306961155055793e-08,\n",
       "  'r2': 0.9999999898448159,\n",
       "  'sensors': ['J137', 'P43', 'J52']},\n",
       " {'node': 'J236',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.884677797814177e-08,\n",
       "  'r2': 0.9999999777112225,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J239',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.916001989582635e-06,\n",
       "  'r2': 0.9999976964120979,\n",
       "  'sensors': ['J316', 'J146', 'P280']},\n",
       " {'node': 'J240',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.005442614465622082,\n",
       "  'r2': 0.9978449983696953,\n",
       "  'sensors': ['P520', 'J316', 'J146']},\n",
       " {'node': 'J241',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.974354965979175e-07,\n",
       "  'r2': 0.9999997822495357,\n",
       "  'sensors': ['J146', 'P280', 'J33']},\n",
       " {'node': 'J242',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.0521437677197579e-06,\n",
       "  'r2': 0.9999995975926069,\n",
       "  'sensors': ['J146', 'P280', 'J33']},\n",
       " {'node': 'J243',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.006432081234945e-08,\n",
       "  'r2': 0.9999999782984543,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J244',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 9.6021347653231e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J244', 'P266', 'P22']},\n",
       " {'node': 'J245',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.013305083970289048,\n",
       "  'r2': 0.9999856278933872,\n",
       "  'sensors': ['P544', 'J244', 'P266']},\n",
       " {'node': 'J246',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.354843640764391e-28,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J244', 'P266', 'P22']},\n",
       " {'node': 'J247',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0006324032658338907,\n",
       "  'r2': 0.9997752346889472,\n",
       "  'sensors': ['J113', 'P110', 'J146']},\n",
       " {'node': 'J248',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0001266228785666625,\n",
       "  'r2': 0.999952799204551,\n",
       "  'sensors': ['P520', 'J39', 'P478']},\n",
       " {'node': 'J249',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.608012621954379e-05,\n",
       "  'r2': 0.9999843447452066,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J250',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.6042682327569933e-05,\n",
       "  'r2': 0.999992681877482,\n",
       "  'sensors': ['P13', 'J17', 'J67']},\n",
       " {'node': 'J251',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.139151332013642e-05,\n",
       "  'r2': 0.9999886260716152,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J252',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0002169644614793642,\n",
       "  'r2': 0.9998852076275105,\n",
       "  'sensors': ['J67', 'P459', 'P13']},\n",
       " {'node': 'J253',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.0710514544838576e-05,\n",
       "  'r2': 0.9999799235928402,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J254',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.4014635972564627e-05,\n",
       "  'r2': 0.9999897323374105,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J256',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.3703049532132415e-06,\n",
       "  'r2': 0.9999995025293622,\n",
       "  'sensors': ['P27', 'J33', 'J28']},\n",
       " {'node': 'J257',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.057790454316844e-11,\n",
       "  'r2': 0.9999999999768736,\n",
       "  'sensors': ['J54', 'P56', 'P46']},\n",
       " {'node': 'J258',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.567696647759022e-12,\n",
       "  'r2': 0.9999999999981268,\n",
       "  'sensors': ['J54', 'P46', 'J62']},\n",
       " {'node': 'J259',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.001737804569586838,\n",
       "  'r2': 0.9992252768982857,\n",
       "  'sensors': ['P539', 'J539', 'P460']},\n",
       " {'node': 'J260',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0020766921859079693,\n",
       "  'r2': 0.999173223781608,\n",
       "  'sensors': ['P539', 'J539', 'J137']},\n",
       " {'node': 'J261',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.877322022330363e-05,\n",
       "  'r2': 0.999994956688614,\n",
       "  'sensors': ['P539', 'J539', 'J137']},\n",
       " {'node': 'J262',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0022284259578712203,\n",
       "  'r2': 0.9989324368849812,\n",
       "  'sensors': ['P539', 'J539', 'J137']},\n",
       " {'node': 'J263',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.04151594406227607,\n",
       "  'r2': 0.9999518821801977,\n",
       "  'sensors': ['P544', 'J244', 'P266']},\n",
       " {'node': 'J264',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.015828136498021766,\n",
       "  'r2': 0.9999825421993567,\n",
       "  'sensors': ['P544', 'J244', 'P266']},\n",
       " {'node': 'J265',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.013196502456835494,\n",
       "  'r2': 0.9999858987146477,\n",
       "  'sensors': ['P544', 'J244', 'P266']},\n",
       " {'node': 'J267',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.89507737529789e-05,\n",
       "  'r2': 0.9999493458473329,\n",
       "  'sensors': ['P460', 'J284', 'P539']},\n",
       " {'node': 'J269',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.38600051306001e-08,\n",
       "  'r2': 0.9999999735105968,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J273',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.879059780504579e-07,\n",
       "  'r2': 0.9999997602850533,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J274',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0001911090423841159,\n",
       "  'r2': 0.9999266164357943,\n",
       "  'sensors': ['P520', 'J316', 'J60']},\n",
       " {'node': 'J275',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00020738825160168933,\n",
       "  'r2': 0.999920008695059,\n",
       "  'sensors': ['P520', 'J316', 'J60']},\n",
       " {'node': 'J276',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 8.686865361382852e-06,\n",
       "  'r2': 0.9999967403645169,\n",
       "  'sensors': ['J316', 'J146', 'P280']},\n",
       " {'node': 'J277',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.023576323804381e-05,\n",
       "  'r2': 0.9999920199339306,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J278',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.003329075837260173,\n",
       "  'r2': 0.9988842932592445,\n",
       "  'sensors': ['J151', 'P213', 'J33']},\n",
       " {'node': 'J279',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.6538459739937424e-10,\n",
       "  'r2': 0.9999999997828037,\n",
       "  'sensors': ['P52', 'J60', 'J54']},\n",
       " {'node': 'J280',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.0256013222868118e-06,\n",
       "  'r2': 0.9999991835262719,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J281',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.9056377860552365e-06,\n",
       "  'r2': 0.9999984771658815,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J282',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.7961784003576465e-05,\n",
       "  'r2': 0.9999708936806921,\n",
       "  'sensors': ['J67', 'P459', 'J52']},\n",
       " {'node': 'J283',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.309978943734572e-05,\n",
       "  'r2': 0.9999777085939653,\n",
       "  'sensors': ['J67', 'P459', 'J52']},\n",
       " {'node': 'J284',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.262177448353619e-29,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J284', 'P460', 'P539']},\n",
       " {'node': 'J285',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.5629513452694785e-10,\n",
       "  'r2': 0.9999999999330126,\n",
       "  'sensors': ['J284', 'P460', 'P539']},\n",
       " {'node': 'J286',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J284', 'P460', 'P539']},\n",
       " {'node': 'J287',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.027963667757043092,\n",
       "  'r2': 0.990972307297612,\n",
       "  'sensors': ['J316', 'P22', 'P520']},\n",
       " {'node': 'J288',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.7068458063885802e-05,\n",
       "  'r2': 0.999991757674896,\n",
       "  'sensors': ['J316', 'P22', 'P520']},\n",
       " {'node': 'J289',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.2414983312292027,\n",
       "  'r2': 0.9279999452865866,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J290',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0001523214208764417,\n",
       "  'r2': 0.9999361672149325,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J291',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0001492344840589687,\n",
       "  'r2': 0.9999417519497584,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J292',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00010018646002718767,\n",
       "  'r2': 0.9999474957649701,\n",
       "  'sensors': ['J172', 'P174', 'J43']},\n",
       " {'node': 'J293',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.09946389035007e-06,\n",
       "  'r2': 0.9999971066415507,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J294',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.492590061363636e-06,\n",
       "  'r2': 0.9999976692736156,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J295',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 9.71189216657201e-06,\n",
       "  'r2': 0.9999970873401406,\n",
       "  'sensors': ['J39', 'P478', 'P174']},\n",
       " {'node': 'J296',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.152322426113252e-06,\n",
       "  'r2': 0.9999982616293513,\n",
       "  'sensors': ['P174', 'J172', 'J39']},\n",
       " {'node': 'J297',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.121659238068613e-06,\n",
       "  'r2': 0.9999987343779857,\n",
       "  'sensors': ['P174', 'J172', 'J39']},\n",
       " {'node': 'J298',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.6336017517682262e-05,\n",
       "  'r2': 0.9999945552831017,\n",
       "  'sensors': ['P520', 'J316', 'J39']},\n",
       " {'node': 'J299',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00028061214713008965,\n",
       "  'r2': 0.9998832978480167,\n",
       "  'sensors': ['P520', 'J39', 'P478']},\n",
       " {'node': 'J300',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0003412332416223126,\n",
       "  'r2': 0.9998481268328747,\n",
       "  'sensors': ['J284', 'P9', 'P460']},\n",
       " {'node': 'J301',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0003815344567717908,\n",
       "  'r2': 0.9997907621948937,\n",
       "  'sensors': ['J284', 'P9', 'P460']},\n",
       " {'node': 'J302',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0004698792648810825,\n",
       "  'r2': 0.9997477661902137,\n",
       "  'sensors': ['J79', 'P500', 'P9']},\n",
       " {'node': 'J303',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0008163928659866165,\n",
       "  'r2': 0.9996064870830957,\n",
       "  'sensors': ['J79', 'P500', 'P9']},\n",
       " {'node': 'J304',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0013009699825104475,\n",
       "  'r2': 0.9993346432284048,\n",
       "  'sensors': ['J163', 'P536', 'J43']},\n",
       " {'node': 'J305',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.043594893260729e-06,\n",
       "  'r2': 0.9999974044697066,\n",
       "  'sensors': ['P13', 'J17', 'J12']},\n",
       " {'node': 'J306',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.0862424626674731e-05,\n",
       "  'r2': 0.9999941419918389,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J307',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.5802402903412713e-07,\n",
       "  'r2': 0.9999999014100589,\n",
       "  'sensors': ['J79', 'P500', 'J12']},\n",
       " {'node': 'J308',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 7.20758123314891e-06,\n",
       "  'r2': 0.9999967025085064,\n",
       "  'sensors': ['J79', 'P500', 'J67']},\n",
       " {'node': 'J309',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 6.682200864289395e-08,\n",
       "  'r2': 0.999999969542588,\n",
       "  'sensors': ['J67', 'P459', 'J79']},\n",
       " {'node': 'J310',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.1212546993924613e-06,\n",
       "  'r2': 0.9999994083080762,\n",
       "  'sensors': ['J67', 'P459', 'J79']},\n",
       " {'node': 'J311',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.845273636033316e-05,\n",
       "  'r2': 0.9999853810519428,\n",
       "  'sensors': ['J172', 'P174', 'J67']},\n",
       " {'node': 'J312',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.3681818104890295e-06,\n",
       "  'r2': 0.9999986551593889,\n",
       "  'sensors': ['J43', 'P511', 'P46']},\n",
       " {'node': 'J313',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.0119804076813552e-05,\n",
       "  'r2': 0.9999957020662392,\n",
       "  'sensors': ['J43', 'P511', 'J163']},\n",
       " {'node': 'J314',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.709377446375991e-05,\n",
       "  'r2': 0.999991037861208,\n",
       "  'sensors': ['J113', 'P110', 'J39']},\n",
       " {'node': 'J315',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 4.1945876935802504e-05,\n",
       "  'r2': 0.9999814169879125,\n",
       "  'sensors': ['J113', 'P110', 'J39']},\n",
       " {'node': 'J316',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J316', 'P520', 'P22']},\n",
       " {'node': 'J317',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.603457345474872e-06,\n",
       "  'r2': 0.9999985843435296,\n",
       "  'sensors': ['J316', 'P520', 'P22']},\n",
       " {'node': 'J318',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0003852508437197719,\n",
       "  'r2': 0.9998208185628795,\n",
       "  'sensors': ['J79', 'P500', 'P9']},\n",
       " {'node': 'J319',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 3.124605161938362e-06,\n",
       "  'r2': 0.9999989775288554,\n",
       "  'sensors': ['J316', 'P520', 'P22']},\n",
       " {'node': 'J320',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.8937223726027567e-07,\n",
       "  'r2': 0.9999999379538008,\n",
       "  'sensors': ['J316', 'P520', 'P22']},\n",
       " {'node': 'J321',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.33423790469084563,\n",
       "  'r2': 0.9852682630998942,\n",
       "  'sensors': ['P22', 'J316', 'J28']},\n",
       " {'node': 'J322',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 1.794568154239284e-09,\n",
       "  'r2': 0.999999999094777,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J323',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.00011271672461465454,\n",
       "  'r2': 0.9999573619156374,\n",
       "  'sensors': ['P520', 'J39', 'P478']},\n",
       " {'node': 'J324',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 5.156740820181028e-05,\n",
       "  'r2': 0.9999741134589979,\n",
       "  'sensors': ['J67', 'P459', 'J52']},\n",
       " {'node': 'J326',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.0077959826727576e-05,\n",
       "  'r2': 0.9999901342433322,\n",
       "  'sensors': ['J12', 'P13', 'J79']},\n",
       " {'node': 'J327',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 2.5013486680919158e-14,\n",
       "  'r2': 0.9999999999999889,\n",
       "  'sensors': ['J52', 'P45', 'J137']},\n",
       " {'node': 'J539',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.0,\n",
       "  'r2': 1.0,\n",
       "  'sensors': ['J539', 'P539', 'J137']},\n",
       " {'node': 'J544',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.013403462278741582,\n",
       "  'r2': 0.9999845713667007,\n",
       "  'sensors': ['P544', 'J244', 'P266']},\n",
       " {'node': 'J546',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.03500491732055697,\n",
       "  'r2': 0.9999638438725439,\n",
       "  'sensors': ['P544', 'J244', 'P266']},\n",
       " {'node': 'T1',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.1277801020426148,\n",
       "  'r2': 0.9160019175985321,\n",
       "  'sensors': ['P460', 'J284', 'P539']},\n",
       " {'node': 'T2',\n",
       "  'linreg': LinearRegression(),\n",
       "  'msq': 0.003172214927225728,\n",
       "  'r2': 0.9981125395936037,\n",
       "  'sensors': ['J54', 'P52', 'P46']}]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "linreg_models"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "higher-worth",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\users\\pc\\anaconda3\\envs\\wodociagi\\lib\\site-packages\\scipy\\optimize\\minpack.py:828: OptimizeWarning: Covariance of the parameters could not be estimated\n",
      "  warnings.warn('Covariance of the parameters could not be estimated',\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Creating a NetworkX graph from .inp file (networkx_graph.py - create_graph())\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>J1</th>\n",
       "      <th>J2</th>\n",
       "      <th>J3</th>\n",
       "      <th>J4</th>\n",
       "      <th>J5</th>\n",
       "      <th>J6</th>\n",
       "      <th>J7</th>\n",
       "      <th>J8</th>\n",
       "      <th>J9</th>\n",
       "      <th>J10</th>\n",
       "      <th>...</th>\n",
       "      <th>J322</th>\n",
       "      <th>J323</th>\n",
       "      <th>J324</th>\n",
       "      <th>J326</th>\n",
       "      <th>J327</th>\n",
       "      <th>J539</th>\n",
       "      <th>J544</th>\n",
       "      <th>J546</th>\n",
       "      <th>T1</th>\n",
       "      <th>T2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>36.007728</td>\n",
       "      <td>30.043946</td>\n",
       "      <td>29.924401</td>\n",
       "      <td>31.958709</td>\n",
       "      <td>32.148558</td>\n",
       "      <td>45.656096</td>\n",
       "      <td>31.842500</td>\n",
       "      <td>28.841451</td>\n",
       "      <td>33.860558</td>\n",
       "      <td>43.419479</td>\n",
       "      <td>...</td>\n",
       "      <td>39.770512</td>\n",
       "      <td>70.297167</td>\n",
       "      <td>43.725289</td>\n",
       "      <td>49.239213</td>\n",
       "      <td>27.804030</td>\n",
       "      <td>30.955127</td>\n",
       "      <td>83.333916</td>\n",
       "      <td>79.320369</td>\n",
       "      <td>25.459113</td>\n",
       "      <td>19.839908</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>600</th>\n",
       "      <td>38.995160</td>\n",
       "      <td>33.155931</td>\n",
       "      <td>33.036006</td>\n",
       "      <td>35.025148</td>\n",
       "      <td>35.215815</td>\n",
       "      <td>48.320933</td>\n",
       "      <td>34.824782</td>\n",
       "      <td>31.824736</td>\n",
       "      <td>36.845347</td>\n",
       "      <td>46.002609</td>\n",
       "      <td>...</td>\n",
       "      <td>42.109475</td>\n",
       "      <td>72.072027</td>\n",
       "      <td>45.880900</td>\n",
       "      <td>51.501378</td>\n",
       "      <td>30.138199</td>\n",
       "      <td>30.583525</td>\n",
       "      <td>85.060845</td>\n",
       "      <td>81.056026</td>\n",
       "      <td>28.321491</td>\n",
       "      <td>20.400044</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1200</th>\n",
       "      <td>39.116052</td>\n",
       "      <td>33.277189</td>\n",
       "      <td>33.157258</td>\n",
       "      <td>35.146255</td>\n",
       "      <td>35.336938</td>\n",
       "      <td>48.442713</td>\n",
       "      <td>34.945636</td>\n",
       "      <td>31.945608</td>\n",
       "      <td>36.966216</td>\n",
       "      <td>46.124199</td>\n",
       "      <td>...</td>\n",
       "      <td>42.230851</td>\n",
       "      <td>72.193733</td>\n",
       "      <td>46.002615</td>\n",
       "      <td>51.622721</td>\n",
       "      <td>30.259576</td>\n",
       "      <td>30.704902</td>\n",
       "      <td>85.172661</td>\n",
       "      <td>81.168453</td>\n",
       "      <td>28.419035</td>\n",
       "      <td>20.520383</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1800</th>\n",
       "      <td>36.748178</td>\n",
       "      <td>30.906448</td>\n",
       "      <td>30.786441</td>\n",
       "      <td>32.776993</td>\n",
       "      <td>32.967176</td>\n",
       "      <td>46.008250</td>\n",
       "      <td>32.578082</td>\n",
       "      <td>29.578119</td>\n",
       "      <td>34.598261</td>\n",
       "      <td>43.710317</td>\n",
       "      <td>...</td>\n",
       "      <td>39.060846</td>\n",
       "      <td>68.986034</td>\n",
       "      <td>43.241612</td>\n",
       "      <td>49.267862</td>\n",
       "      <td>27.100760</td>\n",
       "      <td>27.511339</td>\n",
       "      <td>24.344780</td>\n",
       "      <td>21.084591</td>\n",
       "      <td>26.512908</td>\n",
       "      <td>20.167304</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2400</th>\n",
       "      <td>36.714150</td>\n",
       "      <td>30.872316</td>\n",
       "      <td>30.752312</td>\n",
       "      <td>32.742904</td>\n",
       "      <td>32.933083</td>\n",
       "      <td>45.973971</td>\n",
       "      <td>32.544065</td>\n",
       "      <td>29.544096</td>\n",
       "      <td>34.564239</td>\n",
       "      <td>43.676092</td>\n",
       "      <td>...</td>\n",
       "      <td>39.026682</td>\n",
       "      <td>68.951776</td>\n",
       "      <td>43.207352</td>\n",
       "      <td>49.233707</td>\n",
       "      <td>27.066595</td>\n",
       "      <td>27.477174</td>\n",
       "      <td>24.313286</td>\n",
       "      <td>21.052974</td>\n",
       "      <td>26.485451</td>\n",
       "      <td>20.133432</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>465600</th>\n",
       "      <td>32.274378</td>\n",
       "      <td>26.460614</td>\n",
       "      <td>26.340631</td>\n",
       "      <td>28.321744</td>\n",
       "      <td>28.511380</td>\n",
       "      <td>41.543831</td>\n",
       "      <td>28.104025</td>\n",
       "      <td>25.103830</td>\n",
       "      <td>30.124328</td>\n",
       "      <td>39.232900</td>\n",
       "      <td>...</td>\n",
       "      <td>34.504603</td>\n",
       "      <td>64.433670</td>\n",
       "      <td>38.671457</td>\n",
       "      <td>44.749682</td>\n",
       "      <td>22.543039</td>\n",
       "      <td>23.069615</td>\n",
       "      <td>20.160332</td>\n",
       "      <td>16.878780</td>\n",
       "      <td>23.059765</td>\n",
       "      <td>16.318970</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>466200</th>\n",
       "      <td>32.185867</td>\n",
       "      <td>26.373761</td>\n",
       "      <td>26.253760</td>\n",
       "      <td>28.234313</td>\n",
       "      <td>28.423923</td>\n",
       "      <td>41.455189</td>\n",
       "      <td>28.015443</td>\n",
       "      <td>25.015286</td>\n",
       "      <td>30.035751</td>\n",
       "      <td>39.144698</td>\n",
       "      <td>...</td>\n",
       "      <td>34.419044</td>\n",
       "      <td>64.348648</td>\n",
       "      <td>38.584743</td>\n",
       "      <td>44.663115</td>\n",
       "      <td>22.457443</td>\n",
       "      <td>22.984958</td>\n",
       "      <td>20.082207</td>\n",
       "      <td>16.800167</td>\n",
       "      <td>22.994623</td>\n",
       "      <td>16.240604</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>466800</th>\n",
       "      <td>32.100813</td>\n",
       "      <td>26.288083</td>\n",
       "      <td>26.168086</td>\n",
       "      <td>28.148874</td>\n",
       "      <td>28.338471</td>\n",
       "      <td>41.369494</td>\n",
       "      <td>27.930430</td>\n",
       "      <td>24.930259</td>\n",
       "      <td>29.950720</td>\n",
       "      <td>39.059392</td>\n",
       "      <td>...</td>\n",
       "      <td>34.334990</td>\n",
       "      <td>64.265235</td>\n",
       "      <td>38.500540</td>\n",
       "      <td>44.579043</td>\n",
       "      <td>22.373392</td>\n",
       "      <td>22.899159</td>\n",
       "      <td>20.005425</td>\n",
       "      <td>16.723079</td>\n",
       "      <td>22.924540</td>\n",
       "      <td>16.162398</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>467400</th>\n",
       "      <td>32.016400</td>\n",
       "      <td>26.203089</td>\n",
       "      <td>26.083097</td>\n",
       "      <td>28.064103</td>\n",
       "      <td>28.253688</td>\n",
       "      <td>41.284434</td>\n",
       "      <td>27.846055</td>\n",
       "      <td>24.845870</td>\n",
       "      <td>29.866329</td>\n",
       "      <td>38.974692</td>\n",
       "      <td>...</td>\n",
       "      <td>34.251417</td>\n",
       "      <td>64.182202</td>\n",
       "      <td>38.416806</td>\n",
       "      <td>44.495449</td>\n",
       "      <td>22.289821</td>\n",
       "      <td>22.814050</td>\n",
       "      <td>19.929131</td>\n",
       "      <td>16.646297</td>\n",
       "      <td>22.855141</td>\n",
       "      <td>16.084085</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>468000</th>\n",
       "      <td>31.932555</td>\n",
       "      <td>26.118703</td>\n",
       "      <td>25.998716</td>\n",
       "      <td>27.979926</td>\n",
       "      <td>28.169498</td>\n",
       "      <td>41.199940</td>\n",
       "      <td>27.762246</td>\n",
       "      <td>24.762049</td>\n",
       "      <td>29.782505</td>\n",
       "      <td>38.890530</td>\n",
       "      <td>...</td>\n",
       "      <td>34.168272</td>\n",
       "      <td>64.099508</td>\n",
       "      <td>38.333491</td>\n",
       "      <td>44.412284</td>\n",
       "      <td>22.206678</td>\n",
       "      <td>22.729553</td>\n",
       "      <td>19.853081</td>\n",
       "      <td>16.569881</td>\n",
       "      <td>22.786349</td>\n",
       "      <td>16.005680</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>781 rows × 309 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "               J1         J2         J3         J4         J5         J6  \\\n",
       "0       36.007728  30.043946  29.924401  31.958709  32.148558  45.656096   \n",
       "600     38.995160  33.155931  33.036006  35.025148  35.215815  48.320933   \n",
       "1200    39.116052  33.277189  33.157258  35.146255  35.336938  48.442713   \n",
       "1800    36.748178  30.906448  30.786441  32.776993  32.967176  46.008250   \n",
       "2400    36.714150  30.872316  30.752312  32.742904  32.933083  45.973971   \n",
       "...           ...        ...        ...        ...        ...        ...   \n",
       "465600  32.274378  26.460614  26.340631  28.321744  28.511380  41.543831   \n",
       "466200  32.185867  26.373761  26.253760  28.234313  28.423923  41.455189   \n",
       "466800  32.100813  26.288083  26.168086  28.148874  28.338471  41.369494   \n",
       "467400  32.016400  26.203089  26.083097  28.064103  28.253688  41.284434   \n",
       "468000  31.932555  26.118703  25.998716  27.979926  28.169498  41.199940   \n",
       "\n",
       "               J7         J8         J9        J10  ...       J322       J323  \\\n",
       "0       31.842500  28.841451  33.860558  43.419479  ...  39.770512  70.297167   \n",
       "600     34.824782  31.824736  36.845347  46.002609  ...  42.109475  72.072027   \n",
       "1200    34.945636  31.945608  36.966216  46.124199  ...  42.230851  72.193733   \n",
       "1800    32.578082  29.578119  34.598261  43.710317  ...  39.060846  68.986034   \n",
       "2400    32.544065  29.544096  34.564239  43.676092  ...  39.026682  68.951776   \n",
       "...           ...        ...        ...        ...  ...        ...        ...   \n",
       "465600  28.104025  25.103830  30.124328  39.232900  ...  34.504603  64.433670   \n",
       "466200  28.015443  25.015286  30.035751  39.144698  ...  34.419044  64.348648   \n",
       "466800  27.930430  24.930259  29.950720  39.059392  ...  34.334990  64.265235   \n",
       "467400  27.846055  24.845870  29.866329  38.974692  ...  34.251417  64.182202   \n",
       "468000  27.762246  24.762049  29.782505  38.890530  ...  34.168272  64.099508   \n",
       "\n",
       "             J324       J326       J327       J539       J544       J546  \\\n",
       "0       43.725289  49.239213  27.804030  30.955127  83.333916  79.320369   \n",
       "600     45.880900  51.501378  30.138199  30.583525  85.060845  81.056026   \n",
       "1200    46.002615  51.622721  30.259576  30.704902  85.172661  81.168453   \n",
       "1800    43.241612  49.267862  27.100760  27.511339  24.344780  21.084591   \n",
       "2400    43.207352  49.233707  27.066595  27.477174  24.313286  21.052974   \n",
       "...           ...        ...        ...        ...        ...        ...   \n",
       "465600  38.671457  44.749682  22.543039  23.069615  20.160332  16.878780   \n",
       "466200  38.584743  44.663115  22.457443  22.984958  20.082207  16.800167   \n",
       "466800  38.500540  44.579043  22.373392  22.899159  20.005425  16.723079   \n",
       "467400  38.416806  44.495449  22.289821  22.814050  19.929131  16.646297   \n",
       "468000  38.333491  44.412284  22.206678  22.729553  19.853081  16.569881   \n",
       "\n",
       "               T1         T2  \n",
       "0       25.459113  19.839908  \n",
       "600     28.321491  20.400044  \n",
       "1200    28.419035  20.520383  \n",
       "1800    26.512908  20.167304  \n",
       "2400    26.485451  20.133432  \n",
       "...           ...        ...  \n",
       "465600  23.059765  16.318970  \n",
       "466200  22.994623  16.240604  \n",
       "466800  22.924540  16.162398  \n",
       "467400  22.855141  16.084085  \n",
       "468000  22.786349  16.005680  \n",
       "\n",
       "[781 rows x 309 columns]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# simulate a leak scenario\n",
    "sim_LEAK = wntr_WSN.get_sim_results_LEAK(leak_node, leak_area, leak_start, leak_end)\n",
    "\n",
    "# join pressures and flowrates\n",
    "sim_LEAK_raw = pd.concat([sim_LEAK.node['pressure'], sim_LEAK.link['flowrate']], axis=1)\n",
    "predictions = {}\n",
    "\n",
    "# for each node in graph\n",
    "G = ng.create_graph()\n",
    "for node in G.nodes():\n",
    "\n",
    "    # find linear regression model for this particular node\n",
    "    for model in linreg_models:\n",
    "        if model['node'] == node:\n",
    "            linreg_dict = model\n",
    "\n",
    "    # choose only the available sensors\n",
    "    sim_LEAK_scope = sim_LEAK_raw[linreg_dict['sensors']]\n",
    "\n",
    "    # normalise the input data with respective parameters\n",
    "    sim_LEAK_scope = (sim_LEAK_scope - normalization_params['mean'].loc[linreg_dict['sensors']]) / normalization_params['std'].loc[linreg_dict['sensors']]\n",
    "\n",
    "    # predict pressure for each node\n",
    "    predictions[node] = linreg_dict['linreg'].predict(sim_LEAK_scope)\n",
    "\n",
    "# \"sim_duration*6\" because (60 minutes = 10 minutes * 6)\n",
    "# it's related to simulation step time \n",
    "index = [600*i for i in range(0, (sim_duration*6)+1)]\n",
    "\n",
    "pred_LEAK_results = pd.DataFrame(data=predictions, index=index)\n",
    "pred_LEAK_results"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "enormous-installation",
   "metadata": {},
   "source": [
    "## Get results of leak scenario"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "frank-fountain",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>J1</th>\n",
       "      <th>J2</th>\n",
       "      <th>J3</th>\n",
       "      <th>J4</th>\n",
       "      <th>J5</th>\n",
       "      <th>J6</th>\n",
       "      <th>J7</th>\n",
       "      <th>J8</th>\n",
       "      <th>J9</th>\n",
       "      <th>J10</th>\n",
       "      <th>...</th>\n",
       "      <th>P548</th>\n",
       "      <th>P549</th>\n",
       "      <th>P551</th>\n",
       "      <th>P553</th>\n",
       "      <th>PMP51</th>\n",
       "      <th>PMP63</th>\n",
       "      <th>PMP72</th>\n",
       "      <th>FCV51</th>\n",
       "      <th>FCV63</th>\n",
       "      <th>FCV72</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>36.011416</td>\n",
       "      <td>30.044852</td>\n",
       "      <td>29.924852</td>\n",
       "      <td>31.958688</td>\n",
       "      <td>32.148688</td>\n",
       "      <td>45.531102</td>\n",
       "      <td>31.841401</td>\n",
       "      <td>28.841400</td>\n",
       "      <td>33.861399</td>\n",
       "      <td>43.310270</td>\n",
       "      <td>...</td>\n",
       "      <td>0.005835</td>\n",
       "      <td>-0.001642</td>\n",
       "      <td>0.001183</td>\n",
       "      <td>-0.000037</td>\n",
       "      <td>0.020833</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0434</td>\n",
       "      <td>0.020833</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0434</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>600</th>\n",
       "      <td>38.996412</td>\n",
       "      <td>33.156409</td>\n",
       "      <td>33.036409</td>\n",
       "      <td>35.026409</td>\n",
       "      <td>35.216409</td>\n",
       "      <td>48.256448</td>\n",
       "      <td>34.826398</td>\n",
       "      <td>31.826397</td>\n",
       "      <td>36.846395</td>\n",
       "      <td>45.955310</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.003068</td>\n",
       "      <td>0.000596</td>\n",
       "      <td>-0.000168</td>\n",
       "      <td>-0.000035</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0434</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0434</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1200</th>\n",
       "      <td>39.117789</td>\n",
       "      <td>33.277785</td>\n",
       "      <td>33.157785</td>\n",
       "      <td>35.147786</td>\n",
       "      <td>35.337786</td>\n",
       "      <td>48.377824</td>\n",
       "      <td>34.947774</td>\n",
       "      <td>31.947773</td>\n",
       "      <td>36.967772</td>\n",
       "      <td>46.076686</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.003068</td>\n",
       "      <td>0.000596</td>\n",
       "      <td>-0.000168</td>\n",
       "      <td>-0.000035</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0434</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0434</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1800</th>\n",
       "      <td>36.745965</td>\n",
       "      <td>30.905961</td>\n",
       "      <td>30.785961</td>\n",
       "      <td>32.775962</td>\n",
       "      <td>32.965962</td>\n",
       "      <td>46.006000</td>\n",
       "      <td>32.575951</td>\n",
       "      <td>29.575950</td>\n",
       "      <td>34.595948</td>\n",
       "      <td>43.706562</td>\n",
       "      <td>...</td>\n",
       "      <td>0.001599</td>\n",
       "      <td>-0.000348</td>\n",
       "      <td>0.000072</td>\n",
       "      <td>-0.000008</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2400</th>\n",
       "      <td>36.711800</td>\n",
       "      <td>30.871797</td>\n",
       "      <td>30.751797</td>\n",
       "      <td>32.741797</td>\n",
       "      <td>32.931797</td>\n",
       "      <td>45.971836</td>\n",
       "      <td>32.541786</td>\n",
       "      <td>29.541785</td>\n",
       "      <td>34.561783</td>\n",
       "      <td>43.672397</td>\n",
       "      <td>...</td>\n",
       "      <td>0.001599</td>\n",
       "      <td>-0.000348</td>\n",
       "      <td>0.000072</td>\n",
       "      <td>-0.000008</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>465600</th>\n",
       "      <td>32.287546</td>\n",
       "      <td>26.463700</td>\n",
       "      <td>26.343700</td>\n",
       "      <td>28.328067</td>\n",
       "      <td>28.518067</td>\n",
       "      <td>41.516770</td>\n",
       "      <td>28.117472</td>\n",
       "      <td>25.117466</td>\n",
       "      <td>30.137446</td>\n",
       "      <td>39.212196</td>\n",
       "      <td>...</td>\n",
       "      <td>0.002238</td>\n",
       "      <td>-0.000472</td>\n",
       "      <td>-0.000261</td>\n",
       "      <td>-0.000033</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>466200</th>\n",
       "      <td>32.199657</td>\n",
       "      <td>26.376986</td>\n",
       "      <td>26.256986</td>\n",
       "      <td>28.240945</td>\n",
       "      <td>28.430945</td>\n",
       "      <td>41.429355</td>\n",
       "      <td>28.029243</td>\n",
       "      <td>25.029191</td>\n",
       "      <td>30.049191</td>\n",
       "      <td>39.124934</td>\n",
       "      <td>...</td>\n",
       "      <td>0.002277</td>\n",
       "      <td>-0.000481</td>\n",
       "      <td>-0.000258</td>\n",
       "      <td>-0.000033</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>466800</th>\n",
       "      <td>32.114172</td>\n",
       "      <td>26.291202</td>\n",
       "      <td>26.171202</td>\n",
       "      <td>28.155264</td>\n",
       "      <td>28.345264</td>\n",
       "      <td>41.344436</td>\n",
       "      <td>27.943758</td>\n",
       "      <td>24.943705</td>\n",
       "      <td>29.963705</td>\n",
       "      <td>39.040154</td>\n",
       "      <td>...</td>\n",
       "      <td>0.002303</td>\n",
       "      <td>-0.000487</td>\n",
       "      <td>-0.000255</td>\n",
       "      <td>-0.000033</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>467400</th>\n",
       "      <td>32.029338</td>\n",
       "      <td>26.206106</td>\n",
       "      <td>26.086106</td>\n",
       "      <td>28.070259</td>\n",
       "      <td>28.260259</td>\n",
       "      <td>41.260099</td>\n",
       "      <td>27.858923</td>\n",
       "      <td>24.858871</td>\n",
       "      <td>29.878871</td>\n",
       "      <td>38.955939</td>\n",
       "      <td>...</td>\n",
       "      <td>0.002325</td>\n",
       "      <td>-0.000493</td>\n",
       "      <td>-0.000252</td>\n",
       "      <td>-0.000033</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>468000</th>\n",
       "      <td>31.945081</td>\n",
       "      <td>26.121619</td>\n",
       "      <td>26.001619</td>\n",
       "      <td>27.985852</td>\n",
       "      <td>28.175852</td>\n",
       "      <td>41.176279</td>\n",
       "      <td>27.774667</td>\n",
       "      <td>24.774614</td>\n",
       "      <td>29.794615</td>\n",
       "      <td>38.872225</td>\n",
       "      <td>...</td>\n",
       "      <td>0.002345</td>\n",
       "      <td>-0.000497</td>\n",
       "      <td>-0.000250</td>\n",
       "      <td>-0.000033</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>781 rows × 679 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "               J1         J2         J3         J4         J5         J6  \\\n",
       "0       36.011416  30.044852  29.924852  31.958688  32.148688  45.531102   \n",
       "600     38.996412  33.156409  33.036409  35.026409  35.216409  48.256448   \n",
       "1200    39.117789  33.277785  33.157785  35.147786  35.337786  48.377824   \n",
       "1800    36.745965  30.905961  30.785961  32.775962  32.965962  46.006000   \n",
       "2400    36.711800  30.871797  30.751797  32.741797  32.931797  45.971836   \n",
       "...           ...        ...        ...        ...        ...        ...   \n",
       "465600  32.287546  26.463700  26.343700  28.328067  28.518067  41.516770   \n",
       "466200  32.199657  26.376986  26.256986  28.240945  28.430945  41.429355   \n",
       "466800  32.114172  26.291202  26.171202  28.155264  28.345264  41.344436   \n",
       "467400  32.029338  26.206106  26.086106  28.070259  28.260259  41.260099   \n",
       "468000  31.945081  26.121619  26.001619  27.985852  28.175852  41.176279   \n",
       "\n",
       "               J7         J8         J9        J10  ...      P548      P549  \\\n",
       "0       31.841401  28.841400  33.861399  43.310270  ...  0.005835 -0.001642   \n",
       "600     34.826398  31.826397  36.846395  45.955310  ... -0.003068  0.000596   \n",
       "1200    34.947774  31.947773  36.967772  46.076686  ... -0.003068  0.000596   \n",
       "1800    32.575951  29.575950  34.595948  43.706562  ...  0.001599 -0.000348   \n",
       "2400    32.541786  29.541785  34.561783  43.672397  ...  0.001599 -0.000348   \n",
       "...           ...        ...        ...        ...  ...       ...       ...   \n",
       "465600  28.117472  25.117466  30.137446  39.212196  ...  0.002238 -0.000472   \n",
       "466200  28.029243  25.029191  30.049191  39.124934  ...  0.002277 -0.000481   \n",
       "466800  27.943758  24.943705  29.963705  39.040154  ...  0.002303 -0.000487   \n",
       "467400  27.858923  24.858871  29.878871  38.955939  ...  0.002325 -0.000493   \n",
       "468000  27.774667  24.774614  29.794615  38.872225  ...  0.002345 -0.000497   \n",
       "\n",
       "            P551      P553     PMP51  PMP63   PMP72     FCV51  FCV63   FCV72  \n",
       "0       0.001183 -0.000037  0.020833    0.0  0.0434  0.020833    0.0  0.0434  \n",
       "600    -0.000168 -0.000035  0.000000    0.0  0.0434  0.000000    0.0  0.0434  \n",
       "1200   -0.000168 -0.000035  0.000000    0.0  0.0434  0.000000    0.0  0.0434  \n",
       "1800    0.000072 -0.000008  0.000000    0.0  0.0000  0.000000    0.0  0.0000  \n",
       "2400    0.000072 -0.000008  0.000000    0.0  0.0000  0.000000    0.0  0.0000  \n",
       "...          ...       ...       ...    ...     ...       ...    ...     ...  \n",
       "465600 -0.000261 -0.000033  0.000000    0.0  0.0000  0.000000    0.0  0.0000  \n",
       "466200 -0.000258 -0.000033  0.000000    0.0  0.0000  0.000000    0.0  0.0000  \n",
       "466800 -0.000255 -0.000033  0.000000    0.0  0.0000  0.000000    0.0  0.0000  \n",
       "467400 -0.000252 -0.000033  0.000000    0.0  0.0000  0.000000    0.0  0.0000  \n",
       "468000 -0.000250 -0.000033  0.000000    0.0  0.0000  0.000000    0.0  0.0000  \n",
       "\n",
       "[781 rows x 679 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sim_LEAK_results = pd.concat([sim_LEAK.node['pressure'], sim_LEAK.link['flowrate']], axis=1)\n",
    "sim_LEAK_results"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "powerful-session",
   "metadata": {},
   "source": [
    "## Calculate residuum"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "nearby-standard",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>J1</th>\n",
       "      <th>J10</th>\n",
       "      <th>J100</th>\n",
       "      <th>J102</th>\n",
       "      <th>J103</th>\n",
       "      <th>J104</th>\n",
       "      <th>J105</th>\n",
       "      <th>J106</th>\n",
       "      <th>J107</th>\n",
       "      <th>J108</th>\n",
       "      <th>...</th>\n",
       "      <th>J98</th>\n",
       "      <th>J99</th>\n",
       "      <th>JPMP51</th>\n",
       "      <th>JPMP63</th>\n",
       "      <th>JPMP72</th>\n",
       "      <th>R5</th>\n",
       "      <th>R6</th>\n",
       "      <th>R7</th>\n",
       "      <th>T1</th>\n",
       "      <th>T2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0.003687</td>\n",
       "      <td>-0.109209</td>\n",
       "      <td>0.011408</td>\n",
       "      <td>0.027907</td>\n",
       "      <td>0.013708</td>\n",
       "      <td>0.012505</td>\n",
       "      <td>-0.005332</td>\n",
       "      <td>-0.006691</td>\n",
       "      <td>-0.017601</td>\n",
       "      <td>0.035667</td>\n",
       "      <td>...</td>\n",
       "      <td>0.008589</td>\n",
       "      <td>0.027087</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.540887</td>\n",
       "      <td>0.160092</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>600</th>\n",
       "      <td>0.001252</td>\n",
       "      <td>-0.047299</td>\n",
       "      <td>-0.039084</td>\n",
       "      <td>-0.027194</td>\n",
       "      <td>-0.009116</td>\n",
       "      <td>0.003952</td>\n",
       "      <td>-0.004576</td>\n",
       "      <td>-0.005254</td>\n",
       "      <td>-0.007135</td>\n",
       "      <td>0.002183</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.040772</td>\n",
       "      <td>-0.033877</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-2.101202</td>\n",
       "      <td>-0.315379</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1200</th>\n",
       "      <td>0.001736</td>\n",
       "      <td>-0.047513</td>\n",
       "      <td>-0.038637</td>\n",
       "      <td>-0.026882</td>\n",
       "      <td>-0.009232</td>\n",
       "      <td>0.003735</td>\n",
       "      <td>-0.004569</td>\n",
       "      <td>-0.005303</td>\n",
       "      <td>-0.007103</td>\n",
       "      <td>0.002171</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.040328</td>\n",
       "      <td>-0.033429</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-2.198747</td>\n",
       "      <td>-0.314342</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1800</th>\n",
       "      <td>-0.002213</td>\n",
       "      <td>-0.003755</td>\n",
       "      <td>-0.003068</td>\n",
       "      <td>-0.002196</td>\n",
       "      <td>0.000417</td>\n",
       "      <td>0.002022</td>\n",
       "      <td>0.003481</td>\n",
       "      <td>0.002464</td>\n",
       "      <td>0.003436</td>\n",
       "      <td>0.001874</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.003166</td>\n",
       "      <td>-0.001375</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.292619</td>\n",
       "      <td>0.009757</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2400</th>\n",
       "      <td>-0.002349</td>\n",
       "      <td>-0.003694</td>\n",
       "      <td>-0.003194</td>\n",
       "      <td>-0.002284</td>\n",
       "      <td>0.000450</td>\n",
       "      <td>0.002083</td>\n",
       "      <td>0.003479</td>\n",
       "      <td>0.002477</td>\n",
       "      <td>0.003427</td>\n",
       "      <td>0.001877</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.003291</td>\n",
       "      <td>-0.001501</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.265163</td>\n",
       "      <td>0.009464</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>465600</th>\n",
       "      <td>0.013168</td>\n",
       "      <td>-0.020704</td>\n",
       "      <td>-0.009884</td>\n",
       "      <td>-0.011965</td>\n",
       "      <td>-0.008830</td>\n",
       "      <td>-0.002837</td>\n",
       "      <td>-0.001542</td>\n",
       "      <td>-0.000557</td>\n",
       "      <td>-0.001875</td>\n",
       "      <td>-0.004576</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.009853</td>\n",
       "      <td>-0.009901</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.363185</td>\n",
       "      <td>-0.027191</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>466200</th>\n",
       "      <td>0.013791</td>\n",
       "      <td>-0.019764</td>\n",
       "      <td>-0.009016</td>\n",
       "      <td>-0.011202</td>\n",
       "      <td>-0.008611</td>\n",
       "      <td>-0.002678</td>\n",
       "      <td>-0.001509</td>\n",
       "      <td>-0.000502</td>\n",
       "      <td>-0.001896</td>\n",
       "      <td>-0.004208</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.008908</td>\n",
       "      <td>-0.009230</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.382425</td>\n",
       "      <td>-0.025979</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>466800</th>\n",
       "      <td>0.013359</td>\n",
       "      <td>-0.019238</td>\n",
       "      <td>-0.009300</td>\n",
       "      <td>-0.011429</td>\n",
       "      <td>-0.008618</td>\n",
       "      <td>-0.002499</td>\n",
       "      <td>-0.001507</td>\n",
       "      <td>-0.000471</td>\n",
       "      <td>-0.001955</td>\n",
       "      <td>-0.004175</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.009180</td>\n",
       "      <td>-0.009518</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.398718</td>\n",
       "      <td>-0.025284</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>467400</th>\n",
       "      <td>0.012937</td>\n",
       "      <td>-0.018753</td>\n",
       "      <td>-0.009584</td>\n",
       "      <td>-0.011652</td>\n",
       "      <td>-0.008616</td>\n",
       "      <td>-0.002324</td>\n",
       "      <td>-0.001506</td>\n",
       "      <td>-0.000441</td>\n",
       "      <td>-0.002010</td>\n",
       "      <td>-0.004145</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.009452</td>\n",
       "      <td>-0.009806</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.414936</td>\n",
       "      <td>-0.024740</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>468000</th>\n",
       "      <td>0.012526</td>\n",
       "      <td>-0.018305</td>\n",
       "      <td>-0.009867</td>\n",
       "      <td>-0.011870</td>\n",
       "      <td>-0.008605</td>\n",
       "      <td>-0.002154</td>\n",
       "      <td>-0.001505</td>\n",
       "      <td>-0.000411</td>\n",
       "      <td>-0.002062</td>\n",
       "      <td>-0.004117</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.009726</td>\n",
       "      <td>-0.010093</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-0.431089</td>\n",
       "      <td>-0.024331</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>781 rows × 315 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "              J1       J10      J100      J102      J103      J104      J105  \\\n",
       "0       0.003687 -0.109209  0.011408  0.027907  0.013708  0.012505 -0.005332   \n",
       "600     0.001252 -0.047299 -0.039084 -0.027194 -0.009116  0.003952 -0.004576   \n",
       "1200    0.001736 -0.047513 -0.038637 -0.026882 -0.009232  0.003735 -0.004569   \n",
       "1800   -0.002213 -0.003755 -0.003068 -0.002196  0.000417  0.002022  0.003481   \n",
       "2400   -0.002349 -0.003694 -0.003194 -0.002284  0.000450  0.002083  0.003479   \n",
       "...          ...       ...       ...       ...       ...       ...       ...   \n",
       "465600  0.013168 -0.020704 -0.009884 -0.011965 -0.008830 -0.002837 -0.001542   \n",
       "466200  0.013791 -0.019764 -0.009016 -0.011202 -0.008611 -0.002678 -0.001509   \n",
       "466800  0.013359 -0.019238 -0.009300 -0.011429 -0.008618 -0.002499 -0.001507   \n",
       "467400  0.012937 -0.018753 -0.009584 -0.011652 -0.008616 -0.002324 -0.001506   \n",
       "468000  0.012526 -0.018305 -0.009867 -0.011870 -0.008605 -0.002154 -0.001505   \n",
       "\n",
       "            J106      J107      J108  ...       J98       J99  JPMP51  JPMP63  \\\n",
       "0      -0.006691 -0.017601  0.035667  ...  0.008589  0.027087     NaN     NaN   \n",
       "600    -0.005254 -0.007135  0.002183  ... -0.040772 -0.033877     NaN     NaN   \n",
       "1200   -0.005303 -0.007103  0.002171  ... -0.040328 -0.033429     NaN     NaN   \n",
       "1800    0.002464  0.003436  0.001874  ... -0.003166 -0.001375     NaN     NaN   \n",
       "2400    0.002477  0.003427  0.001877  ... -0.003291 -0.001501     NaN     NaN   \n",
       "...          ...       ...       ...  ...       ...       ...     ...     ...   \n",
       "465600 -0.000557 -0.001875 -0.004576  ... -0.009853 -0.009901     NaN     NaN   \n",
       "466200 -0.000502 -0.001896 -0.004208  ... -0.008908 -0.009230     NaN     NaN   \n",
       "466800 -0.000471 -0.001955 -0.004175  ... -0.009180 -0.009518     NaN     NaN   \n",
       "467400 -0.000441 -0.002010 -0.004145  ... -0.009452 -0.009806     NaN     NaN   \n",
       "468000 -0.000411 -0.002062 -0.004117  ... -0.009726 -0.010093     NaN     NaN   \n",
       "\n",
       "        JPMP72  R5  R6  R7        T1        T2  \n",
       "0          NaN NaN NaN NaN  0.540887  0.160092  \n",
       "600        NaN NaN NaN NaN -2.101202 -0.315379  \n",
       "1200       NaN NaN NaN NaN -2.198747 -0.314342  \n",
       "1800       NaN NaN NaN NaN -0.292619  0.009757  \n",
       "2400       NaN NaN NaN NaN -0.265163  0.009464  \n",
       "...        ...  ..  ..  ..       ...       ...  \n",
       "465600     NaN NaN NaN NaN -0.363185 -0.027191  \n",
       "466200     NaN NaN NaN NaN -0.382425 -0.025979  \n",
       "466800     NaN NaN NaN NaN -0.398718 -0.025284  \n",
       "467400     NaN NaN NaN NaN -0.414936 -0.024740  \n",
       "468000     NaN NaN NaN NaN -0.431089 -0.024331  \n",
       "\n",
       "[781 rows x 315 columns]"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "residuals = sim_LEAK.node['pressure'] - pred_LEAK_results\n",
    "residuals"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "deluxe-revelation",
   "metadata": {},
   "source": [
    "## Take care of column naming convention"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "interested-bangladesh",
   "metadata": {},
   "outputs": [],
   "source": [
    "residuals = residuals.add_suffix('_residuum')\n",
    "\n",
    "sim_LEAK.node['pressure'] = sim_LEAK.node['pressure'].add_suffix('_pressure')\n",
    "sim_LEAK.link['flowrate'] = sim_LEAK.link['flowrate'].add_suffix('_flow')\n",
    "sim_LEAK = pd.concat([sim_LEAK.node['pressure'], sim_LEAK.link['flowrate']], axis=1)\n",
    "\n",
    "pred_LEAK_results = pred_LEAK_results.add_suffix('_model_output')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "worth-yugoslavia",
   "metadata": {},
   "source": [
    "## Sanity check"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "early-prisoner",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of residuals columns: 315\n",
      "Number of sim_LEAK columns: 679\n",
      "Number of pred_LEAK_results columns: 309\n",
      "Sum: 1303\n"
     ]
    }
   ],
   "source": [
    "print(f'Number of residuals columns: {len(residuals.columns)}')\n",
    "print(f'Number of sim_LEAK columns: {len(sim_LEAK.columns)}')\n",
    "print(f'Number of pred_LEAK_results columns: {len(pred_LEAK_results.columns)}')\n",
    "print(f'Sum: {len(residuals.columns) + len(sim_LEAK.columns) + len(pred_LEAK_results.columns)}')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "institutional-chest",
   "metadata": {},
   "source": [
    "## Form final dataframe"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "noted-renewal",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>J1_residuum</th>\n",
       "      <th>J10_residuum</th>\n",
       "      <th>J100_residuum</th>\n",
       "      <th>J102_residuum</th>\n",
       "      <th>J103_residuum</th>\n",
       "      <th>J104_residuum</th>\n",
       "      <th>J105_residuum</th>\n",
       "      <th>J106_residuum</th>\n",
       "      <th>J107_residuum</th>\n",
       "      <th>J108_residuum</th>\n",
       "      <th>...</th>\n",
       "      <th>J322_model_output</th>\n",
       "      <th>J323_model_output</th>\n",
       "      <th>J324_model_output</th>\n",
       "      <th>J326_model_output</th>\n",
       "      <th>J327_model_output</th>\n",
       "      <th>J539_model_output</th>\n",
       "      <th>J544_model_output</th>\n",
       "      <th>J546_model_output</th>\n",
       "      <th>T1_model_output</th>\n",
       "      <th>T2_model_output</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0.003687</td>\n",
       "      <td>-0.109209</td>\n",
       "      <td>0.011408</td>\n",
       "      <td>0.027907</td>\n",
       "      <td>0.013708</td>\n",
       "      <td>0.012505</td>\n",
       "      <td>-0.005332</td>\n",
       "      <td>-0.006691</td>\n",
       "      <td>-0.017601</td>\n",
       "      <td>0.035667</td>\n",
       "      <td>...</td>\n",
       "      <td>39.770512</td>\n",
       "      <td>70.297167</td>\n",
       "      <td>43.725289</td>\n",
       "      <td>49.239213</td>\n",
       "      <td>27.804030</td>\n",
       "      <td>30.955127</td>\n",
       "      <td>83.333916</td>\n",
       "      <td>79.320369</td>\n",
       "      <td>25.459113</td>\n",
       "      <td>19.839908</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>600</th>\n",
       "      <td>0.001252</td>\n",
       "      <td>-0.047299</td>\n",
       "      <td>-0.039084</td>\n",
       "      <td>-0.027194</td>\n",
       "      <td>-0.009116</td>\n",
       "      <td>0.003952</td>\n",
       "      <td>-0.004576</td>\n",
       "      <td>-0.005254</td>\n",
       "      <td>-0.007135</td>\n",
       "      <td>0.002183</td>\n",
       "      <td>...</td>\n",
       "      <td>42.109475</td>\n",
       "      <td>72.072027</td>\n",
       "      <td>45.880900</td>\n",
       "      <td>51.501378</td>\n",
       "      <td>30.138199</td>\n",
       "      <td>30.583525</td>\n",
       "      <td>85.060845</td>\n",
       "      <td>81.056026</td>\n",
       "      <td>28.321491</td>\n",
       "      <td>20.400044</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1200</th>\n",
       "      <td>0.001736</td>\n",
       "      <td>-0.047513</td>\n",
       "      <td>-0.038637</td>\n",
       "      <td>-0.026882</td>\n",
       "      <td>-0.009232</td>\n",
       "      <td>0.003735</td>\n",
       "      <td>-0.004569</td>\n",
       "      <td>-0.005303</td>\n",
       "      <td>-0.007103</td>\n",
       "      <td>0.002171</td>\n",
       "      <td>...</td>\n",
       "      <td>42.230851</td>\n",
       "      <td>72.193733</td>\n",
       "      <td>46.002615</td>\n",
       "      <td>51.622721</td>\n",
       "      <td>30.259576</td>\n",
       "      <td>30.704902</td>\n",
       "      <td>85.172661</td>\n",
       "      <td>81.168453</td>\n",
       "      <td>28.419035</td>\n",
       "      <td>20.520383</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1800</th>\n",
       "      <td>-0.002213</td>\n",
       "      <td>-0.003755</td>\n",
       "      <td>-0.003068</td>\n",
       "      <td>-0.002196</td>\n",
       "      <td>0.000417</td>\n",
       "      <td>0.002022</td>\n",
       "      <td>0.003481</td>\n",
       "      <td>0.002464</td>\n",
       "      <td>0.003436</td>\n",
       "      <td>0.001874</td>\n",
       "      <td>...</td>\n",
       "      <td>39.060846</td>\n",
       "      <td>68.986034</td>\n",
       "      <td>43.241612</td>\n",
       "      <td>49.267862</td>\n",
       "      <td>27.100760</td>\n",
       "      <td>27.511339</td>\n",
       "      <td>24.344780</td>\n",
       "      <td>21.084591</td>\n",
       "      <td>26.512908</td>\n",
       "      <td>20.167304</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2400</th>\n",
       "      <td>-0.002349</td>\n",
       "      <td>-0.003694</td>\n",
       "      <td>-0.003194</td>\n",
       "      <td>-0.002284</td>\n",
       "      <td>0.000450</td>\n",
       "      <td>0.002083</td>\n",
       "      <td>0.003479</td>\n",
       "      <td>0.002477</td>\n",
       "      <td>0.003427</td>\n",
       "      <td>0.001877</td>\n",
       "      <td>...</td>\n",
       "      <td>39.026682</td>\n",
       "      <td>68.951776</td>\n",
       "      <td>43.207352</td>\n",
       "      <td>49.233707</td>\n",
       "      <td>27.066595</td>\n",
       "      <td>27.477174</td>\n",
       "      <td>24.313286</td>\n",
       "      <td>21.052974</td>\n",
       "      <td>26.485451</td>\n",
       "      <td>20.133432</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>465600</th>\n",
       "      <td>0.013168</td>\n",
       "      <td>-0.020704</td>\n",
       "      <td>-0.009884</td>\n",
       "      <td>-0.011965</td>\n",
       "      <td>-0.008830</td>\n",
       "      <td>-0.002837</td>\n",
       "      <td>-0.001542</td>\n",
       "      <td>-0.000557</td>\n",
       "      <td>-0.001875</td>\n",
       "      <td>-0.004576</td>\n",
       "      <td>...</td>\n",
       "      <td>34.504603</td>\n",
       "      <td>64.433670</td>\n",
       "      <td>38.671457</td>\n",
       "      <td>44.749682</td>\n",
       "      <td>22.543039</td>\n",
       "      <td>23.069615</td>\n",
       "      <td>20.160332</td>\n",
       "      <td>16.878780</td>\n",
       "      <td>23.059765</td>\n",
       "      <td>16.318970</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>466200</th>\n",
       "      <td>0.013791</td>\n",
       "      <td>-0.019764</td>\n",
       "      <td>-0.009016</td>\n",
       "      <td>-0.011202</td>\n",
       "      <td>-0.008611</td>\n",
       "      <td>-0.002678</td>\n",
       "      <td>-0.001509</td>\n",
       "      <td>-0.000502</td>\n",
       "      <td>-0.001896</td>\n",
       "      <td>-0.004208</td>\n",
       "      <td>...</td>\n",
       "      <td>34.419044</td>\n",
       "      <td>64.348648</td>\n",
       "      <td>38.584743</td>\n",
       "      <td>44.663115</td>\n",
       "      <td>22.457443</td>\n",
       "      <td>22.984958</td>\n",
       "      <td>20.082207</td>\n",
       "      <td>16.800167</td>\n",
       "      <td>22.994623</td>\n",
       "      <td>16.240604</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>466800</th>\n",
       "      <td>0.013359</td>\n",
       "      <td>-0.019238</td>\n",
       "      <td>-0.009300</td>\n",
       "      <td>-0.011429</td>\n",
       "      <td>-0.008618</td>\n",
       "      <td>-0.002499</td>\n",
       "      <td>-0.001507</td>\n",
       "      <td>-0.000471</td>\n",
       "      <td>-0.001955</td>\n",
       "      <td>-0.004175</td>\n",
       "      <td>...</td>\n",
       "      <td>34.334990</td>\n",
       "      <td>64.265235</td>\n",
       "      <td>38.500540</td>\n",
       "      <td>44.579043</td>\n",
       "      <td>22.373392</td>\n",
       "      <td>22.899159</td>\n",
       "      <td>20.005425</td>\n",
       "      <td>16.723079</td>\n",
       "      <td>22.924540</td>\n",
       "      <td>16.162398</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>467400</th>\n",
       "      <td>0.012937</td>\n",
       "      <td>-0.018753</td>\n",
       "      <td>-0.009584</td>\n",
       "      <td>-0.011652</td>\n",
       "      <td>-0.008616</td>\n",
       "      <td>-0.002324</td>\n",
       "      <td>-0.001506</td>\n",
       "      <td>-0.000441</td>\n",
       "      <td>-0.002010</td>\n",
       "      <td>-0.004145</td>\n",
       "      <td>...</td>\n",
       "      <td>34.251417</td>\n",
       "      <td>64.182202</td>\n",
       "      <td>38.416806</td>\n",
       "      <td>44.495449</td>\n",
       "      <td>22.289821</td>\n",
       "      <td>22.814050</td>\n",
       "      <td>19.929131</td>\n",
       "      <td>16.646297</td>\n",
       "      <td>22.855141</td>\n",
       "      <td>16.084085</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>468000</th>\n",
       "      <td>0.012526</td>\n",
       "      <td>-0.018305</td>\n",
       "      <td>-0.009867</td>\n",
       "      <td>-0.011870</td>\n",
       "      <td>-0.008605</td>\n",
       "      <td>-0.002154</td>\n",
       "      <td>-0.001505</td>\n",
       "      <td>-0.000411</td>\n",
       "      <td>-0.002062</td>\n",
       "      <td>-0.004117</td>\n",
       "      <td>...</td>\n",
       "      <td>34.168272</td>\n",
       "      <td>64.099508</td>\n",
       "      <td>38.333491</td>\n",
       "      <td>44.412284</td>\n",
       "      <td>22.206678</td>\n",
       "      <td>22.729553</td>\n",
       "      <td>19.853081</td>\n",
       "      <td>16.569881</td>\n",
       "      <td>22.786349</td>\n",
       "      <td>16.005680</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>781 rows × 1303 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        J1_residuum  J10_residuum  J100_residuum  J102_residuum  \\\n",
       "0          0.003687     -0.109209       0.011408       0.027907   \n",
       "600        0.001252     -0.047299      -0.039084      -0.027194   \n",
       "1200       0.001736     -0.047513      -0.038637      -0.026882   \n",
       "1800      -0.002213     -0.003755      -0.003068      -0.002196   \n",
       "2400      -0.002349     -0.003694      -0.003194      -0.002284   \n",
       "...             ...           ...            ...            ...   \n",
       "465600     0.013168     -0.020704      -0.009884      -0.011965   \n",
       "466200     0.013791     -0.019764      -0.009016      -0.011202   \n",
       "466800     0.013359     -0.019238      -0.009300      -0.011429   \n",
       "467400     0.012937     -0.018753      -0.009584      -0.011652   \n",
       "468000     0.012526     -0.018305      -0.009867      -0.011870   \n",
       "\n",
       "        J103_residuum  J104_residuum  J105_residuum  J106_residuum  \\\n",
       "0            0.013708       0.012505      -0.005332      -0.006691   \n",
       "600         -0.009116       0.003952      -0.004576      -0.005254   \n",
       "1200        -0.009232       0.003735      -0.004569      -0.005303   \n",
       "1800         0.000417       0.002022       0.003481       0.002464   \n",
       "2400         0.000450       0.002083       0.003479       0.002477   \n",
       "...               ...            ...            ...            ...   \n",
       "465600      -0.008830      -0.002837      -0.001542      -0.000557   \n",
       "466200      -0.008611      -0.002678      -0.001509      -0.000502   \n",
       "466800      -0.008618      -0.002499      -0.001507      -0.000471   \n",
       "467400      -0.008616      -0.002324      -0.001506      -0.000441   \n",
       "468000      -0.008605      -0.002154      -0.001505      -0.000411   \n",
       "\n",
       "        J107_residuum  J108_residuum  ...  J322_model_output  \\\n",
       "0           -0.017601       0.035667  ...          39.770512   \n",
       "600         -0.007135       0.002183  ...          42.109475   \n",
       "1200        -0.007103       0.002171  ...          42.230851   \n",
       "1800         0.003436       0.001874  ...          39.060846   \n",
       "2400         0.003427       0.001877  ...          39.026682   \n",
       "...               ...            ...  ...                ...   \n",
       "465600      -0.001875      -0.004576  ...          34.504603   \n",
       "466200      -0.001896      -0.004208  ...          34.419044   \n",
       "466800      -0.001955      -0.004175  ...          34.334990   \n",
       "467400      -0.002010      -0.004145  ...          34.251417   \n",
       "468000      -0.002062      -0.004117  ...          34.168272   \n",
       "\n",
       "        J323_model_output  J324_model_output  J326_model_output  \\\n",
       "0               70.297167          43.725289          49.239213   \n",
       "600             72.072027          45.880900          51.501378   \n",
       "1200            72.193733          46.002615          51.622721   \n",
       "1800            68.986034          43.241612          49.267862   \n",
       "2400            68.951776          43.207352          49.233707   \n",
       "...                   ...                ...                ...   \n",
       "465600          64.433670          38.671457          44.749682   \n",
       "466200          64.348648          38.584743          44.663115   \n",
       "466800          64.265235          38.500540          44.579043   \n",
       "467400          64.182202          38.416806          44.495449   \n",
       "468000          64.099508          38.333491          44.412284   \n",
       "\n",
       "        J327_model_output  J539_model_output  J544_model_output  \\\n",
       "0               27.804030          30.955127          83.333916   \n",
       "600             30.138199          30.583525          85.060845   \n",
       "1200            30.259576          30.704902          85.172661   \n",
       "1800            27.100760          27.511339          24.344780   \n",
       "2400            27.066595          27.477174          24.313286   \n",
       "...                   ...                ...                ...   \n",
       "465600          22.543039          23.069615          20.160332   \n",
       "466200          22.457443          22.984958          20.082207   \n",
       "466800          22.373392          22.899159          20.005425   \n",
       "467400          22.289821          22.814050          19.929131   \n",
       "468000          22.206678          22.729553          19.853081   \n",
       "\n",
       "        J546_model_output  T1_model_output  T2_model_output  \n",
       "0               79.320369        25.459113        19.839908  \n",
       "600             81.056026        28.321491        20.400044  \n",
       "1200            81.168453        28.419035        20.520383  \n",
       "1800            21.084591        26.512908        20.167304  \n",
       "2400            21.052974        26.485451        20.133432  \n",
       "...                   ...              ...              ...  \n",
       "465600          16.878780        23.059765        16.318970  \n",
       "466200          16.800167        22.994623        16.240604  \n",
       "466800          16.723079        22.924540        16.162398  \n",
       "467400          16.646297        22.855141        16.084085  \n",
       "468000          16.569881        22.786349        16.005680  \n",
       "\n",
       "[781 rows x 1303 columns]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "output = pd.concat([residuals, sim_LEAK, pred_LEAK_results], axis=1)\n",
    "output"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "prompt-delivery",
   "metadata": {},
   "source": [
    "## Export the final dataframe to .xlsx"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "impressed-joint",
   "metadata": {},
   "outputs": [],
   "source": [
    "if comment:\n",
    "    output.to_excel(os.path.join(os.getcwd(), f'{method}_n{leak_node}_la{leak_area}_lst{leak_start}_let{leak_end}_{comment}.xlsx'))\n",
    "else:\n",
    "    output.to_excel(os.path.join(os.getcwd(), f'{method}_n{leak_node}_la{leak_area}_lst{leak_start}_let{leak_end}.xlsx'))"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}