Zheng Zheng (Yale), Renyue Cen (Princeton), Hy Trac (CfA), and Jordi Miralda-Escude (ICREA/UB) We combine a cosmological reionization simulation with box size of 100Mpc/h on a side and a Monte Carlo Lyman-alpha radiative transfer code to model Lyman Alpha Emitters (LAEs) at z~5.7. The model essentially adopts a constant halo mass to (UV or intrinsic Lya) light ratio and introduces Lya radiative transfer as the single factor for transforming the intrinsic Lya emission properties into the observed ones. The spatial and frequency diffusions of Lya photons prove to be crucial in determining the observed properties of high-z LAEs. Radiative transfer depends on the environment of intergalactic medium surrounding LAEs and on the viewing angle, which leads to a broad distribution of apparent (observed) Lya luminosity for a given intrinsic Lya luminosity. This environment and viewing angle dependent appearance of LAEs introduced by Lya radiative transfer is able to explain an array of observed properties of LAEs, including their morphology, Lya spectra, Lya luminosity function, UV luminosity function, and distribution of Lya equivalent width. The model also predicts new features in the clustering of LAEs. I will discuss implications of the results for our understanding of LAEs and for studying reionization and cosmology with LAEs.