The first-order least-squares meshfree method for linear elasticity is presented. The conventional and the compatibility-imposed least-squares formulations are studied on the convergence behavior of the solution and the robustness to integration error. Since the least-squares formulation is a type of mixed formulation and induces positive-definite system matrix, by using shape functions of same order for both primal and dual variables, higher rate of convergence is obtained for dual variables than Galerkin formulation. Numerical examples also show that the presented formulations do not exhibit any volumetric locking for the incompressible materials.