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
Post a Comment