Riesgo de ganglios linfáticos comprometidos en cáncer de próstata

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

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")

intercept = -4.71127143

psa_m1 = 0.04427068

g1_3_m1 = 1.90797075

g2_3_m1 = 0.89152421

s2a_m1 = 0.54386908

s2b_m1 = 0.92753512

s2c_m1 = 1.24526645

s3_m1 = 1.63966106

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>> ")

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 = intercept + float(psa) * psa_m1 + float(g1) * g1_3_m1 + float(g2) * g2_3_m1 + float(s2a) * s2a_m1 + float(s2b) * s2b_m1 + float(s2c) * s2c_m1 + float(s3) * s3_m1

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

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