PurposeThe purpose of this paper is to attempt the buckling analysis of a laminated composite skew plate using the C0 finite element (FE) model based on higher-order shear deformation theory (HSDT) in conjunction with minimax probability machine regression (MPMR) and multivariate adaptive regression spline (MARS).Design/methodology/approachHSDT considers the third-order variation of in-plane displacements which eliminates the use of shear correction factor owing to realistic parabolic transverse shear stresses across the thickness coordinate. At the top and bottom of the plate, zero transverse shear stress condition is imposed. C0 FE model based on HSDT is developed and coded in formula translation (FORTRAN). FE model is validated and found efficient to create new results. MPMR and MARS models are coded in MATLAB. Using skew angle (α), stacking sequence (Ai) and buckling strength (Y) as input parameters, a regression problem is formulated using MPMR and MARS to predict the buckling strength of laminated composite skew plates.FindingsThe results of the MPMR and MARS models are in good agreement with the FE model result. MPMR is a better tool than MARS to analyze the buckling problem.Research limitations/implicationsThe present work considers the linear behavior of the laminated composite skew plate.Originality/valueTo the authors’ best of knowledge, there is no work in the literature on the buckling analysis of a laminated composite skew plate using C0 FE formulation based on third-order shear deformation theory in conjunction with MPMR and MARS. These machine-learning techniques increase efficiency, reduce the computational time and reduce the cost of analysis. Further, an equation is generated with the MARS model via which the buckling strength of the laminated composite skew plate can be predicted with ease and simplicity.