The intrinsic radiosensitivity index (RSI) and genomic-adjusted radiation dose (GARD) were reported to be able to predict the surviving fraction at 2 Gy and therapeutic effect when delivering actual treatment doses using the gene expression profiles of clinical cases. Given the impact of p16 status, a surrogate marker of the human papillomavirus (HPV) infection, on radiosensitivity, we attempted to apply the RSI and GARD to estimate p16-associated radiosensitivity in head and neck squamous cell carcinoma (HNSC). For this purpose, The Cancer Genome Atlas (TCGA) dataset was employed. In the GARD calculation, we assumed that p16-positive patients received 60 Gy in 30 fractions, while p16-negative patients received 70 Gy in 35 fractions. p16 positivity was associated with favorable characteristics compared to negative patients. The RSI and GARD analyses demonstrated increased radiosensitivity and high therapeutic effect in p16-positive patients, compared to p16-negative patients. Additionally, tumor microenvironmental conditions predicted by other models were also significantly affected by p16 status. Collectively, the models used in this study could be a promising tool for estimating p16-associated radiosensitivity in HNSC.