We investigate the statistical anisotropy and gaussianity of temperature fluctuations of Cosmic Microwave Background (CMB) radiation data from the Wilkinson Microwave Anisotropy Probe survey, using the Multifractal Detrended Fluctuation Analysis, Rescaled Range, and Scaled Windowed Variance methods. Multifractal Detrended Fluctuation Analysis shows that CMB fluctuations has a long-range correlation function with a multifractal behavior. By comparing the shuffled and surrogate series of CMB data, we conclude that the multifractality nature of the temperature fluctuation of CMB radiation is mainly due to the long-range correlations, and the map is consistent with a gaussian distribution.