{ "cells": [ { "cell_type": "code", "execution_count": 12, "id": "6b84294b", "metadata": {}, "outputs": [], "source": [ "# import all packages\n", "import pandas as pd \n", "import numpy as np \n", "import textwrap \n", "\n", "import seaborn as sns \n", "import matplotlib.pyplot as plt \n", "\n", "import shap \n", "\n", "from sklearn.preprocessing import StandardScaler \n", "from sklearn import linear_model \n", "\n", "from sklearn.model_selection import KFold, cross_validate, cross_val_predict, RepeatedKFold, cross_val_score, GridSearchCV \n", "from sklearn.metrics import r2_score, mean_absolute_error \n", "\n", "import statsmodels.api as sm" ] }, { "cell_type": "code", "execution_count": 13, "id": "5662f441", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of Columns: 107\n" ] }, { "data": { "text/html": [ "
\n", " | Unnamed: 0 | \n", "Accretion Measurement 1 (mm) | \n", "Accretion Measurement 2 (mm) | \n", "Accretion Measurement 3 (mm) | \n", "Accretion Measurement 4 (mm) | \n", "Latitude | \n", "Longitude | \n", "Year (yyyy) | \n", "Month (mm) | \n", "Average Accretion (mm) | \n", "... | \n", "percent_waterlevel_complete | \n", "Unnamed: 6 | \n", "Surface Elevation Change Rate (cm/y) | \n", "Flood Freq (Floods/yr) | \n", "Avg. Flood Depth when Flooded (ft) | \n", "90th Percentile Flood Depth when Flooded (ft) | \n", "10th Percentile Flood Depth when Flooded (ft) | \n", "Std. Deviation Flood Depth when Flooded | \n", "Acc_rate_fullterm (cm/y) | \n", "Average_Marsh_Elevation (ft. NAVD88) | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
0 | \n", "CRMS0002 | \n", "32.545 | \n", "34.0 | \n", "34.0 | \n", "30.0 | \n", "30.100356 | \n", "-89.791784 | \n", "2016.0 | \n", "4.0 | \n", "28.75 | \n", "... | \n", "97.739726 | \n", "NaN | \n", "0.49 | \n", "163.104824 | \n", "0.442805 | \n", "1.104 | \n", "0.00 | \n", "0.506806 | \n", "1.05 | \n", "0.351500 | \n", "
1 | \n", "CRMS0003 | \n", "24.500 | \n", "23.0 | \n", "24.0 | \n", "11.0 | \n", "30.099396 | \n", "-89.252697 | \n", "2015.0 | \n", "4.0 | \n", "11.25 | \n", "... | \n", "99.965753 | \n", "NaN | \n", "0.55 | \n", "243.955101 | \n", "0.140492 | \n", "0.531 | \n", "0.00 | \n", "0.333084 | \n", "1.26 | \n", "0.853500 | \n", "
2 | \n", "CRMS0006 | \n", "28.000 | \n", "30.0 | \n", "27.0 | \n", "20.5 | \n", "30.263548 | \n", "-89.982233 | \n", "2016.0 | \n", "4.0 | \n", "18.50 | \n", "... | \n", "98.390411 | \n", "NaN | \n", "0.41 | \n", "57.171053 | \n", "0.334836 | \n", "0.898 | \n", "0.00 | \n", "0.445684 | \n", "1.18 | \n", "0.522500 | \n", "
3 | \n", "CRMS0008 | \n", "0.000 | \n", "0.0 | \n", "0.0 | \n", "0.0 | \n", "30.228348 | \n", "-90.685225 | \n", "2016.0 | \n", "5.0 | \n", "0.00 | \n", "... | \n", "0.000000 | \n", "NaN | \n", "0.42 | \n", "69.701166 | \n", "0.729878 | \n", "1.270 | \n", "0.27 | \n", "0.514597 | \n", "0.75 | \n", "0.441111 | \n", "
4 | \n", "CRMS0030 | \n", "44.000 | \n", "48.0 | \n", "49.0 | \n", "37.5 | \n", "30.303906 | \n", "-90.325364 | \n", "2016.0 | \n", "4.0 | \n", "37.75 | \n", "... | \n", "96.529680 | \n", "NaN | \n", "0.60 | \n", "142.737073 | \n", "0.271160 | \n", "0.850 | \n", "0.00 | \n", "0.446133 | \n", "1.73 | \n", "0.661000 | \n", "
5 rows × 107 columns
\n", "\n", " | Soil Porewater Salinity (ppt) | \n", "Average Height Dominant (cm) | \n", "NDVI | \n", "TSS (mg/l) | \n", "windspeed | \n", "avg_percentflooded (%) | \n", "Flood Freq (Floods/yr) | \n", "Distance to River (km) | \n", "Distance to Water (km) | \n", "Land Lost (km2) | \n", "Tidal Amplitude (cm) | \n", "90th Percentile Flood Depth (cm) | \n", "10th Percentile Flood Depth (cm) | \n", "Avg. Flood Depth (cm) | \n", "Std. Deviation Flood Depth (cm) | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
count | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "266.000000 | \n", "
mean | \n", "7.186861 | \n", "112.439643 | \n", "0.353358 | \n", "31.821080 | \n", "6.683651 | \n", "39.207975 | \n", "261.107686 | \n", "0.009934 | \n", "0.006082 | \n", "0.090947 | \n", "9.283287 | \n", "24.836617 | \n", "0.314712 | \n", "9.731371 | \n", "11.897189 | \n", "
std | \n", "6.039870 | \n", "28.938786 | \n", "0.110811 | \n", "13.033298 | \n", "0.272169 | \n", "28.512498 | \n", "549.612559 | \n", "0.001005 | \n", "0.002949 | \n", "0.198047 | \n", "9.204314 | \n", "9.858570 | \n", "1.560128 | \n", "5.785805 | \n", "3.026577 | \n", "
min | \n", "0.100000 | \n", "54.860000 | \n", "0.057249 | \n", "0.000000 | \n", "5.864558 | \n", "0.000000 | \n", "0.000000 | \n", "0.005079 | \n", "0.000000 | \n", "0.000000 | \n", "0.060442 | \n", "0.914400 | \n", "0.000000 | \n", "1.173515 | \n", "4.595155 | \n", "
25% | \n", "2.425000 | \n", "97.000000 | \n", "0.292554 | \n", "28.556742 | \n", "6.525110 | \n", "15.867526 | \n", "51.426390 | \n", "0.009427 | \n", "0.005693 | \n", "0.001554 | \n", "0.312294 | \n", "18.267426 | \n", "0.000000 | \n", "5.835615 | \n", "10.014167 | \n", "
50% | \n", "5.600000 | \n", "108.300000 | \n", "0.364691 | \n", "32.970487 | \n", "6.703245 | \n", "39.127680 | \n", "161.037234 | \n", "0.010165 | \n", "0.007177 | \n", "0.020488 | \n", "6.584424 | \n", "24.292560 | \n", "0.000000 | \n", "8.483845 | \n", "11.656643 | \n", "
75% | \n", "10.300000 | \n", "122.100000 | \n", "0.399708 | \n", "39.331860 | \n", "6.872427 | \n", "59.170743 | \n", "250.736063 | \n", "0.010722 | \n", "0.008095 | \n", "0.080221 | \n", "18.043368 | \n", "29.794200 | \n", "0.000000 | \n", "12.263635 | \n", "13.760751 | \n", "
max | \n", "27.500000 | \n", "315.800000 | \n", "0.627081 | \n", "65.011248 | \n", "7.315856 | \n", "99.999999 | \n", "4693.096330 | \n", "0.011222 | \n", "0.009412 | \n", "1.394380 | \n", "27.500165 | \n", "65.257680 | \n", "18.592800 | \n", "41.731574 | \n", "24.649485 | \n", "