weights_input = None)
plt.plot(TCP_results['nom_doses'],TCP_results['TCP_pop'], color='darkgreen',lw=1,alpha=1)
plt.plot(TCP_results['dose_input'],TCP_results['TCP_input'],color='green',ls='',marker='o',alpha=0.5)
for i in range(TCP_results['n']):
plt.plot(TCP_results['nom_doses'],TCP_results['TCPs'][i],color='green',ls='-',alpha=0.1)
## should be able to get the list of doses etc from the TCP calc to
## then use within the NTCP calc.
if NTCP_calc == True:
NTCP_results = TCP_NTCP.complete_NTCP_calc(d_data=[45,70, 62, 67, 78, 65],
ntcp_data=[0.1,0.15,0.1,0.3,0.3, 0.19],
frac_doses = TCP_results['doses'],
initial_params_ntcp=None,
max_dose=100,
ntcp_params={'td50_1':(70,0),
'v':(1.0,0),
'm':(0.3,0),
'n':(0.3,0)},
fit_vals = True,
)
#print(NTCP_results['pop_fit'])
plt.plot(NTCP_results['fit_x'],NTCP_results['fit_y'],color='darkred',lw=1,alpha=1)
plt.plot(NTCP_results['d_data'],NTCP_results['ntcp_data'],color='red',ls='',marker='o',alpha=0.5)
for i in range(TCP_results['n']):
plt.plot(TCP_results['nom_doses'],NTCP_results['patient_ntcps'][i],color='red',ls='-',alpha=0.1)
plt.ylim(0,1)
plt.xlim(0,100)
plt.show()
