is it possible to take the statistics for extracted region using R? -


i have binary file size of (360 720 )for globe.i wrote code given below read , extract area file. when use summary whole file got.

summary(a, na.rm=false) min. 1st qu.  median    mean 3rd qu.    max.    na's  0.00    1.00    3.00    4.15    7.00   20.00  200083 . 

but when used summary region(b) extracted, got many v1,v2. not right should have got 1 line (as a)not many v1,v2.

here code:

x <- c(200:300) y <- c(150:190) conne <- file("c:\\initial-wtd.bin", "rb") a=readbin(conne, numeric(), size=4,  n=360*720, signed=true) a[a == -9999] <- na y <- matrix(data=a,ncol=360,nrow=720) image(t(t(y[x,y])),ylim=c(1,0)) b = y[x,y] summary(b,na.rm=false)     v1            v2            v3            v4            v5           v6            v7         min.   : na   min.   : na   min.   : na   min.   : na   min.   : 8   min.   : na   min.   :     1st qu.: na   1st qu.: na   1st qu.: na   1st qu.: na   1st qu.:11   1st qu.: na   1st qu.:    median : na   median : na   median : na   median : na   median :14   median : na   median    mean   :nan   mean   :nan   mean   :nan   mean   :nan   mean   :14   mean   :nan   mean       3rd qu.: na   3rd qu.: na   3rd qu.: na   3rd qu.: na   3rd qu.:17   3rd qu.: na   3rd     max.   : na   max.   : na   max.   : na   max.   : na   max.   :20   max.   : na   max.       na's   :101   na's   :101   na's   :101   na's   :101   na's   :99   na's   :101   na's   : 

the problem not in indexing of matrix, place prior accessing it:

a <- matrix(1:100, 10, 10) summary( a[1:3,1:3] )        v1            v2             v3        min.   :1.0   min.   :11.0   min.   :21.0    1st qu.:1.5   1st qu.:11.5   1st qu.:21.5    median :2.0   median :12.0   median :22.0    mean   :2.0   mean   :12.0   mean   :22.0    3rd qu.:2.5   3rd qu.:12.5   3rd qu.:22.5    max.   :3.0   max.   :13.0   max.   :23.0   

you managed hit few non-na values (apparently 2) why doing such sparse data? scaled 100 columns (out of 1000) , still got expected results.


Comments

Popular posts from this blog

django - How can I change user group without delete record -

java - Need to add SOAP security token -

java - EclipseLink JPA Object is not a known entity type -