Matlab program for calculating EUD-based NTCP and TCP in external beam radiotherapy

Function implementing Dr. Niemierko's EUD-based NTCP and TCP models.

Description

EUDMODEL(DVH), where DVH is a 2 column matrix corresponding to the cumulative, not differential, dose volume histogram. The 1st column corresponds to increasing absolute dose or percentage dose values, and the 2nd column to the corresponding absolute or relative volume value. The matrix must have a minimum of two rows, and both columns must be of equal length.
by Hiram A. Gay, MD
Revised July 8 2007

NTCP = normal tissue complication probability
TCP = tumor control probability
EUD = equivalent uniform dose

References:
Gay HA, Niemierko A. A free program for calculating EUD-based NTCP and TCP in external beam radiotherapy. Phys Med. 2007 Dec;23(3-4):115-25. Epub 2007 Sep 7.
PMID: 17825595 [PubMed - indexed for MEDLINE]

• MATLAB release: MATLAB 7 (R14)

[Heigrujam]

Dear Hiram, I just want to verify on line 161 and 163. According to LYMAN EUD model, inverse of 'a' should be in line 161 and normal 'a' should be in line 163. Can you just verify which formula for NTCP is used?

[Hiram]

Dear Heigrujam, the formula is for Niemierko's EUD based model, and not the Lyman model. It is described in the reference listed in the description which I can gladly provide.
Hiram

[Nicola Dinapoli]

Excellent paper abput TCP/NTCP topic