A coupled three-dimensional pile group analysis method (YSGroup) was developed considering nonlinear pile head stiffness matrices and compared with other analytical methods (elastic displacement method, Group 6.0 and FBPier 3.0). In this method, a pile cap was modelled by four-node flat shell element, a pier was modelled using 3 dimensional beam element, and individual piles were modelled as beam-column elements. Through the comparative studies on a piled pie. subjected to lateral loads in linear soil, it was found that present method (YSGroup), elastic displacement method and Group 6.0 gave similar results of lateral pile head displacement, but FBPier 3.0 was estimated to show somewhat larger displacements than those from the three methods. Displacements of superstructure (pier), including nonlinear soil behavior, could be estimated by present method (YSGroup) and FBPier 3.0 because these two methods modelled the superstructure directly by finite element techniques. It was found that pile groups in pinned pile head condition had a tendency to cause excessive rotation of the pile cap.