from numpy import array

v = array([0.51,0.4,0.54,0.54])
σ = array([10.62,5.07,13.07,16.23])
μ = array([42.90,19.30,66.43,86.14])

def CCS(X):
  return ((v/σ)*(X-μ)).sum()

#@title Prodromic average scores:

SDMT = 45 #@param {type:"integer"}
VFT = 25 #@param {type:"integer"}
SCR = 70 #@param {type:"integer"}
SWR = 90 #@param {type:"integer"}


X = array([SDMT, VFT, SCR, SWR], dtype=int)
if (X<0).any():
  print("Test scores can not be negative")
else:
  ccs = CCS(X)
  print("CCS = %.2f"%ccs)
  if ccs<-2.375:
    print("Group 1: Low CCS, Motor high risk")
  elif ccs<2.375:
    print("Group 2: medium CCS, Motor medium risk")
  else:
    print("Group 3: high CCS, Motor low risk")