北極振動 Rで統計

北極振動　12月だけを抜き出しの続き

12月のデータだけ抜き出して、Rで統計処理する。 Rpy2を使う。
まずは以下のように準備。readAOとpickup_monthは以前作ったもの。

AO_rpy2.py 

from readAO import readAO # read AO data
from pick_month import pick_month # pick month
#
import scikits.timeseries as ts
#
import rpy2.robjects as robjects


# read AO data
AOindex_series=readAO()


# pickup December
AOindex_dec=pick_month(AOindex_series,month=12)


#remove missig value
data=AOindex_dec[AOindex_dec.mask==False].data


# for R
rdata=robjects.FloatVector(data)
robjects.globalEnv["rdatain"] = rdata
r=robjects.r

それで、

a=r('summary(rdatain)')

または

a=r.summary(rdata)

または

summary=r['summary']
a=summary(rdata)

いずれも同じ結果

print a
Min. 1st Qu. Median Mean 3rd Qu. Max.
-3.41300 -1.24200 -0.08762 -0.19740 0.82490 2.28200

print a[0]
-3.413

print a.names
[1] "Min." "1st Qu." "Median" "Mean" "3rd Qu." "Max."

b=a.r["Mean"]
print b
Mean
-0.1974

b[0]
-0.19739999999999999

print a.subset(1)
Min.
-3.413

参考
Rで統計: データ集合中の最大、最小、平均、中央値 - summary()関数http://www.yukun.info/blog/2008/09/r-summary-mean-median.html