Composites are known to display a considerable amount of scatter in their material properties due to large number of parameters associated with the manufacturing and fabrication processes. In the present work, the material properties have been taken as random variables for accurate prediction of the system behavior. Higher order shear theory including rotatory inertia e ects has been accounted for in the system dynamic equations. A ΓΏrst order perturbation technique has been employed to obtain the solution of the governing equations. An approach has been outlined for obtaining closed form expressions for the variances of eigen solutions. The e ects of side to thickness ratio and variation in standard deviation of the material properties have been investigated for cross-ply symmetric and anti-symmetric laminates. The mean and standard deviations of the ΓΏrst ΓΏve natural frequencies have been worked out for laminated rectangular plates with all edges simply supported. The higher order shear deformation theory results have been validated with Monte Carlo simulation results and compared with the results based on classical laminate and ΓΏrst order shear deformation theories.