The integrity of timber structures is mainly related to its capacity to resist crack propagation under various load conditions. However, this phenomenon is random by nature, and the need to incorporate statistical information is mandatory for practical use in structures. This paper aims at defining a probabilistic model in order to characterize the scatter of the toughness test results of timber. The instantaneous failure tests are performed using the mixed-mode crack growth specimen. The crack tip growth is recorded by a video camera for mixed-mode ratios of 15°, 30°and 60°, where the relative displacement of loading points is recorded by LVDT sensor. The experimental energy release rate is evaluated by the compliance method. As large scatter of the energy release rate is observed, the statistical analysis is performed by using the bootstrap simulation, in order to characterize the probabilistic models in the opening and shear crack modes. The reliability analysis is then performed in order to underline the impact of the statistical uncertainties on the rupture of wood material.