The purposes of this paper are to estimate the behavior of laterally loaded single piles by the back analysis and the effects of non-homogeneity consisted of two layers, constraint condition of pile head, relative density, and embedded length of pile on the behavior of single pile in sand. These effects can be quantified only by the results of model tests. Also, the model tests were compared with the results of the numerical methods (p-y method, modified Vlasov method; new γ parameter, Characteristic Load Method;CLM). In this study, the modulus of soil resistance and ultimate soil resistance are assumed to exponential function with depth for Nak-Dong River sand. Therefore, a new γ parameter equation based on the Vlasov method was derived to calculate the modulus of soil resistance (E h =n hi Z n ) of exponential function with depth. Thep-y method of nonlinear behavior analysis is relatively effective method of designing deep foundations subjected to lateral loads for only the reasonable of the modulus of soil resistance (E h ) and ultimate soil resistance (P u ). Therefore, the back analysis is used to determineE h andP u in this study. For this, theE h(DST) andP u(DST) are obtained from the stress-strain relationship of direct shear test (DST). TheE h andP u applicable to predict lateral behavior are obtained from back analysis usingp-y curve with hyperbolic function. The coefficients used to modify ofE h(DST) andP u(DST) are obtained with the range of 0.014–0.05 and 0.2–0.4, respectively. The predict of modified Vlasov method usedE h andP u suggested in this study is very good agreement with measured results blow y/D=0.06–0.14 and the predict ofp-y method used the recommend values of Terzaghi (1955) forE h and Reese et al.(1974) forP u is more under-estimated than the model test results.