This paper presents the motion accuracy simulation considering loads such as workpiece weight, cutting force, cogging force of a linear motor, and force caused by misalignment and runout error of a ballscrew in linear motion units. The transfer function method is basically utilized to estimate 5-DOF motion errors, together with the equilibrium equations of force and moment on the table. The transfer function method is modified in order to consider clearance changed according to the loads in the double sided hydrostatic/aerostatic bearings. Then, the analytic model for predicting the 5-DOF motion errors is proposed with the modified transfer function method. Motion errors were simulated under different loading conditions in the linear motion units using hydrostatic, aerostatic, and linear motion bearings, respectively. And the proposed analytic model was verified by comparing the estimated and measured motion errors.