Evaluating shoreline retreat rate (SRR) on different spatial‐temporal scales is critical for effective coastal management. Large‐scale evaluations typically rely on data‐driven methods such as Discrete Bayesian networks (BNs). However, these BNs require discretization of continuous variables which can lead to information loss. Here, we propose a new method, the Hybrid BN to incorporate continuous variables without discretization. Both Discrete and Hybrid BNs were developed and compared to evaluate large‐scale (continental scale) SRR in Australia, using Digital Earth Australia data set. These BNs used forcing parameters (e.g., waves, tide, sediment sink/source, and sea level rise [SLR]) and geomorphic settings (e.g., geomorphology, backshore profile, and surfzone slope) to predict SRR. Validation of the BNs showed that Hybrid BNs, which provide a more realistic assessment of the range of SRR, outperform in predicting continuous variables, when compared with Discrete BNs. However, Discrete and Hybrid BNs provide consistent qualitative findings for the SRR of Australia. Among forcing parameters, the sediment sink/source was found to be the most informative variable to indicate the shoreline retreat, followed by tide, SLR rate, and wave processes. In the scenario of an increased SLR rate, tropical tidal flats were predicted as the most at risk coasts in Australia. We found that BNs can reflect the impact of different factors on coastal evolution, and predict future shoreline change by exploring historical data. The performance of these models can be further improved when more data sets become available.