$1$ import arviz as az

$2$ from scipy import stats $3$ import os $4$ import pandas as pd $5$ import numpy as np $6$ import matplotlib. pyplot as plt

$7$ usecols $=$ I 'LOA', $\text{'}$LBP$\text{'},$ 'BEAM', 'DEPTH', 'DRAFT', $\text{'}$ PROP

POWER', 'AUX POWER', 'TEU', 'SPEED'$]$

$8$ cwd $=$

os$.$ path.abspath $(\text{'}/$content$/$gdrive$/$MyDrive$/$Container$\_$ShipData.xls

$\text{'})$

$9$ ships$\_$data$2=$ pd$.$read$\_$excel$($cwd$,$ usecols $=$ L$\text{'}$LBP$\text{'},$ 'TEU',

$\text{'}$ SPEED' $])$

$10$ ships$\_$data$2.$ head $()$

$11$ list$1=$ ships$\_$data$2[\text{'}$LBP$\text{'}]$

$12$ list$2=$ ships$\_$data$2[\text{'}$TEU$\text{'}]$

$13$ x $3=$ list$1.$ values

$14$ Y$\_3=$ listz. values

$15$ x$\_3=$x$\_3-$ x$\_3.$ mean $()16\_$ ax $=$ plt$.$ subplots $(1,2,$ Figsize$=(10,5))$

$17$ beta$\_$c$,$ alpha$\_$c $=$ stats.linregress $($x$\_3,$ Y$\_3)1$:$21$

$18$ import pymc$3$ as pm

$19$ with pm$.$ Model $()$ as ship $\_$model$\_$t:

$20$

x

$=$

pm$.$

Normal

$(\text{'}$a$\text{'},$

mu$=$y$\_$

$3.$

mean

B $=$ pm$.$HalfNormal $(\text{'}$B$\text{'},$ sd $=1)$

$=$ pm$.$ HalfNormal $(\text{'}\text{'},5)$

$(),$

sd$=2)$

$21$

$22$

$23$

y$\_$pred

$=$

pm$.$

Normall'y$\_$pred',

mu$=$a

$+$

B

$*$x$\_3,$

sd$=$e$,$ observed$=$y$\_3)$

$24$ trace$\_$t $=$ pm$.$ sample $(2000)$