Hormone therapy in the form of androgen deprivation is a major treatment for advanced prostate cancer. However, if such therapy is overly prolonged, tumour cells may become resistant to this treatment and result in recurrent fatal disease. Long-term hormone deprivation also is associated with side effects poorly tolerated by patients. In contrast, intermittent hormone therapy with alternating on- and off-treatment periods is a possible clinical strategy to delay progression to hormone-refractory disease with the advantage of reduced side effects during the off-treatment periods. In this paper, we first overview previous studies on mathematical modelling of prostate tumour growth under intermittent hormone therapy. The model is categorized into a hybrid dynamical system because switching between on-treatment and off-treatment intervals is treated in addition to continuous dynamics of tumour growth. Next, we present an extended model of stochastic differential equations and examine how well the model is able to capture the characteristics of authentic serum prostate-specific antigen (PSA) data. We also highlight recent advances in time-series analysis and prediction of changes in serum PSA concentrations. Finally, we discuss practical issues to be considered towards establishment of mathematical model-based tailor-made medicine, which defines how to realize personalized hormone therapy for individual patients based on monitored serum PSA levels.