and effective distance are The redshift z$_{d}$ at the baryon drag epoch is fitted with the formula [19] where Ω$_{m}$$_{0}$ is the current value of the dimensionless matter energy density, Ω$_{b}$ is the di- mensionless baryon matter energy density, and h = H$_{0}$/ 100. The sound speed c$_{s}$ ( z ) = 1 / √ 3[1 + ¯ R$_{b}$/ (1 + z )], where ¯ R$_{b}$ = 315000Ω$_{b}$ h $^{2}$( T$_{cmb}$/ 2 . 7K) − 4 and T$_{cmb}$ = 2 . 726K. In [8], two points d$_{0}$$_{.}$$_{2}$ and d$_{0}$$_{.}$$_{35}$ and their covariance matrix are given. We will use d$_{0}$$_{.}$$_{2}$ = 0 . 1905 ± 0 . 0061, d$_{0}$$_{.}$$_{35}$ = 0 . 1097 ± 0 . 0036 and their covariance matrix (hereafter Bao2) as the second sample of BAO data in the data fitting. When we use these data, we add two more parameters Ω$_{b}$ h 2 and h . Based on these two BAO distance measurements, we can derive the BAO distance ratio D$_{V}$ (0 . 35) /D$_{V}$ (0 . 2) = 1 . 736 ± 0 . 065 (hereafter BaoR), which is relatively model inde- pendent quantity. This BAO distance ratio was employed in the analysis in [11]. In our case we will further use another two radial BAO data at redshifts z = 0 . 24 and z = 0 . 43 [20] by using Δ z ( z ) = H ( z ) r$_{s}$ ( z$_{d}$ ) /c (hereafter BaoZ). Combining BAO data sets Bao2 and BaoZ, we have the data set called Bao4. Therefore, for different BAO data sets, we can perform χ 2 statistics for the model parameter p as follows and where i, j = 0 . 2 , 0 . 35, Δ d$_{0}$$_{.}$$_{2}$ = d$_{0}$$_{.}$$_{2}$ − 0 . 1905, Δ d$_{0}$$_{.}$$_{35}$ = d$_{0}$$_{.}$$_{35}$ − 0 . 1097. Since the SNIa and BAO data contain information about the universe at relatively low redshifts, we will include the CMB information by implementing the Wilkinson microwave anisotropy probe 5 year (WMAP5) data to probe the entire expansion history up to the last 5 χ 2 Bao $_{2}$( p , h, Ω$_{b}$ h $^{2}$) = ∑ ij Δ d$_{i}$ Cov − $^{1}$( d$_{i}$, d$_{j}$ )Δ d$_{j}$, (10) χ 2 $_{BaoZ}$( p , h, Ω$_{b}$ h $^{2}$) = [Δ z (0 . 24) − 0 . 0407] 2 0 . 0011 2 + [Δ z (0 . 43) − 0 . 0442] 2 0 . 0015 2 , (9) χ 2 $_{Baor}$( p ) = [ D$_{V}$ (0 . 35) /D$_{V}$ (0 . 2) − 1 . 736] 2 0 . 065 2 , (8) b$_{1}$ = 0 . 313(Ω$_{m}$$_{0}$ h $^{2}$) − 0 . $^{419}$[1 + 0 . 607(Ω$_{m}$$_{0}$ h $^{2}$)0 . $^{674}$] , b$_{2}$ = 0 . 238(Ω$_{m}$$_{0}$ h $^{2}$)0 . $^{223}$, (7) z$_{d}$ = 1291(Ω$_{m}$$_{0}$ h $^{2}$)0 . 251 1 + 0 . 659(Ω$_{m}$$_{0}$ h $^{2}$)0 . $_{828}$[1 + b$_{1}$ (Ω$_{b}$ h $^{2}$) b 2 ] , (6) r$_{s}$ ( z, p ) = ∫ ∞ z dx c$_{s}$ ( x ) E ( x ) , (4) D$_{V}$ ( z, p , H$_{0}$ ) = [ d 2 $_{L}$( z ) (1 + z ) 2 z H ( z ) ] 1 / 3 = H − 1 0 [ z E ( z ) (∫ z 0 dx E ( x ) ) 2 ] 1 / 3 . (5)