Fix for JuliaLang/julia#9098

REQUIRE

@@ -1,3 +1,3 @@
 julia 0.3.0-
 StatsBase 0.3
-Compat 0.2.6
+Compat 0.2.11-

src/indexing.jl

@@ -107,6 +107,13 @@ getindex(t::AbstractDataArray, i::Real) =
 ## getindex: DataArray
 
 # Scalar case
+function getindex(da::DataArray, I::Real)
+    if getindex(da.na, I)
+        return NA
+    else
+        return getindex(da.data, I)
+    end
+end
 @nsplat N function getindex(da::DataArray, I::NTuple{N,Real}...)
     if getindex(da.na, I...)
         return NA
@@ -125,7 +132,7 @@ end
     k = 1
     srcextr = daextract(src)
     destextr = daextract(dest)
-    @nloops N i dest d->(@inbounds offset_{d-1} = offset_d + (Base.unsafe_getindex(I_d, i_d)-1)*stride_d) begin
+    @nloops N i dest d->(@inbounds offset_{d-1} = offset_d + (Base.unsafe_getindex(I[d], i_d)-1)*stride_d) begin
         @inbounds if unsafe_isna(src, srcextr, offset_0)
             unsafe_dasetindex!(dest, destextr, NA, k)
         else
@@ -159,6 +166,13 @@ end
 ## getindex: PooledDataArray
 
 # Scalar case
+function getindex(pda::PooledDataArray, I::Real)
+    if getindex(pda.refs, I) == 0
+        return NA
+    else
+        return pda.pool[getindex(pda.refs, I)]
+    end
+end
 @nsplat N function getindex(pda::PooledDataArray, I::NTuple{N,Real}...)
     if getindex(pda.refs, I...) == 0
         return NA
@@ -215,8 +229,8 @@ end
     end
 
     Aextr = daextract(A)
-    @ncall N checkbounds A J
-    @nexprs N d->(I_d = Base.to_index(J_d))
+    checkbounds(A, J...)
+    @nexprs N d->(I_d = Base.to_index(J[d]))
     stride_1 = 1
     @nexprs N d->(stride_{d+1} = stride_d*size(A,d))
     @nexprs N d->(offset_d = 1)  # really only need offset_$N = 1