Reliability analysis of stress–strength models usually assumes that the stress and strength variables are independent. However, in numerous real-world scenarios, stress and strength variables exhibit dependence. This paper investigates the reliability estimation in a multicomponent stress–strength model for parallel-series system assuming that the dependence between stress and strength is based on the Clayton copula. The estimators for the unknown parameters and system reliability are derived using the two-step maximum likelihood estimation and the maximum product spacing methods. Additionally, confidence intervals are constructed by utilizing asymptotically normal distribution theory and bootstrap method. Furthermore, Monte Carlo simulations are conducted to compare the effectiveness of the proposed inference methods. Finally, a real dataset is analyzed for illustrative purposes.