> source("ARMath.s")

> grep("\\.I",names(ARMath))
 [1]  1  3  6  8  9 10 11 12 13 14 15 16 18 20
> ARInit _ ARMath[,grep("\\.I",names(ARMath))]

> unlist(lapply(ARInit,mode))
     AR.I1    AR.I2C     AR.I3      AR.I3A      AR.I3B      AR.I3C      AR.I3D
 "numeric" "numeric" "numeric" "character" "character" "character" "character"

      AR.I3E     AR.I3ET     AR.I4       AR.I5     AR.I6     AR.I7     AR.I8
 "character" "character" "numeric" "character" "numeric" "numeric"
 "numeric"

> unlist(lapply(ARInit,mode))=="numeric"
 AR.I1 AR.I2C AR.I3 AR.I3A AR.I3B AR.I3C AR.I3D AR.I3E AR.I3ET AR.I4 AR.I5 
     T      T     T      F      F      F      F      F       F     T     F

 AR.I6 AR.I7 AR.I8 
     T     T     T

> ARInit _ ARInit[,unlist(lapply(ARInit,mode))=="numeric"]
>

> for (i in 1:7) hist(ARInit[,i])
> dev.print(postscript,command="cat > bj.eps",onefile=F)



N _ dim(ARInit)[1]
J _ dim(ARInit)[2]

for (i in 1:N) {
  for (j in 1:J) {
      if (ARInit[i,j]>10) {
        ARInit[i,j] <- NA
     }
  }
}