Riesgo de ganglios linfáticos positivos con parámetros clínicos y patológicos completos

https://repl.it/@mauriciolema/ProstateCaNM2

print("Riesgo de ganglios linfáticos positivos en cáncer de próstata no metastásico.\n")

print("Parámetros: PSA, Gleason primario, Gleason secundario, estadío, número de cores negativos, número de cores positivos.\n")

intercept2 = -3.9238186

psa_m2 = 0.04127647

g1_3_m2 = 1.69555954

g2_3_m2 = 0.73627978

s2a_m2 = 0.41819731

s2b_m2 = 0.56954741

s2c_m2 = 0.77240836

s3_m2 = 1.00598066

NegCores_m2 = -0.10567544

PositCores_m2 = 0.06398027

psa = input("Valor del PSA>> ")

g1 = input("Gleason primario >3, 0. No, 1. Sí>> ")

g2 = input("Gleason secundario > 3, 0. No, 1. Sí>> ")

stage = input("Estadío 1. I, 2. IIa, 3. IIb, 4. IIc, 5. III o mayor>> ")

NegCores = input("Número de cores negativos: ")

PositCores = input("Número de cores positivos: ")

if int(stage) == 1:

s2a = 0

s2b = 0

s2c = 0

s3 = 0

if int(stage) == 2:

s2a = 1

s2b = 0

s2c = 0

s3 = 0

if int(stage) == 3:

s2a = 0

s2b = 1

s2c = 0

s3 = 0

if int(stage) == 4:

s2a = 0

s2b = 0

s2c = 1

s3 = 0

if int(stage) == 5:

s2a = 0

s2b = 0

s2c = 0

s3 = 1

ln_risk = intercept2 + float(psa) * psa_m2 + float(g1) * g1_3_m2 + float(g2) * g2_3_m2 + float(s2a) * s2a_m2 + float(s2b) * s2b_m2 + float(s2c) * s2c_m2 + float(s3) * s3_m2 + NegCores_m2 * float(NegCores) + PositCores_m2 * float(PositCores)

predprob = ((2.7172**(ln_risk)/(1+2.7172**(ln_risk))) * 100)

print("Probabilidad de ganglios linfáticos regionales comprometidos: " + str(int(predprob)) + "%.")