WANtaroHP (R tips)

 

Contents

2013年6月より統計解析向けプログラミング言語「R」を使い始めました. きっかけは,仕事上のニーズです.しばらく前に仕事で簡単な統計解析のレポートを書いたのですが,そのときはエクセルとFortranプログラムで対応しました.しかしもう少し複雑な解析になってくると自分のプログラムでは心もとないし,エクセルを使って複雑なことをやる気はない,ということで,前から名前だけ知っていた「R」を試してみようと思った次第です.

使用環境はインストールしたままですが,自分が使う中での注意事項は今のところ以下のとおりです.

  • 言語は英語のほうが良い.コマンドプロンプトで使うとき文字化けして何が書かれているかわからなくなるため.
  • 64bit版と32bit版の両方をインストールする.コマンドプロンプトからrscriptを使う場合に64bitのみでは動作しない.
  • 楕円の描画のため,package:ellipseを追加インストールする.

紹介されている内容は以下のとおりです.















水文統計解析

入力サンプルデータ

FilenameDescription
inp_RF_S.txt札幌雨量データ
inp_RF_M.txt前橋雨量データ
inp_RF_K.txt京都雨量データ

Q-Qプロット

FilenameDescription
R_HF_qq_LN3_IWAI.txtR-script (LN3, IWAI)
R_HF_qq_LN3_mom.txtR-script (LN3, moment)
R_HF_qq_LP3_mom.txtR-script (LP3, moment)
R_HF_qq_GEV.txtR-script (GEV)
R_HF_qq_Gum.txtR-script (Gumbel)
R_HF_qq_SQRT.txtR-script (SQRT-ET)
fig_R_HF_qq_LN3_IWAI.pngImage sample
fig_R_HF_qq_LN3_mom.pngImage sample
fig_R_HF_qq_LP3_mom.pngImage sample
fig_R_HF_qq_GEV.pngImage sample
fig_R_HF_qq_Gum.pngImage sample
fig_R_HF_qq_SQRT.pngImage sample

ヒストグラムとGEV確率密度関数

FilenameDescription
R_HF_pdf_GEV.txtR-script (Histogram & GEV pdf)
fig_R_HF_hist_GEV.pngImage sample

Gumbel 確率紙プロット

FilenameDescription
R_HF_paper_GUM.txtR-script (Gumbel probability paper)
fig_R_HF_paper_GUM.pngImage sample

棄却検定

FilenameDescription
R_HF_RT_LN3_IWAI.txtR-script (Rejection test, LN3, IWAI)
R_HF_RT_LP3.txtR-script (Rejection test, LP3)


FFT を用いた地震応答スペクトル計算

FFTを使うための基本的なRスクリプト

Rスクリプト事例

d  <- scan("inp_FFT_1.csv",skip=3)
nn <- length(d)
z  <- fft(d)
z1 <- Re(z)/nn
z2 <- Im(z)/nn
cc <- sqrt(z1*z1+z2*z2)
zz <- complex(re=z1*nn,im=z2*nn)
y  <- Re(fft(zz,inverse=TRUE))/nn
out<- cbind(d,z1,z2,cc,y)
write.table(out,"out_FFT_1.txt",quote=F,append=F,row.names=F)

サンプル入力データ

Example
dt,1.0
ndata,16
5
32
38
-33
-19
-10
1
-8
-20
10
-1
4
11
-1
-7
-2

サンプル出力データ

d  z1                z2               cc                 y
5  0                 0                0                  5
32  3.87961256380478  2.07134405386938 4.3979381344821   32
38  2.74448051533946 -4.19010191002141 5.00890480199445  38
-33  2.47896579998472 -5.9762205702964  6.46996783163006 -33
-19 -3.375            -4.375            5.52550902632509 -19
-10 -2.09400457708259  1.92803176917456 2.84642963584018 -10
1 -3.61948051533946  1.18489808997859 3.80849349278122   0.999999999999996
-8  1.98542621329309  2.47559639334034 3.17340428422712  -8
-20  1                 0                1                -20
10  1.98542621329309 -2.47559639334034 3.17340428422712  10
-1 -3.61948051533946 -1.18489808997859 3.80849349278123  -1
4 -2.0940045770826  -1.92803176917456 2.84642963584018   4
11 -3.375             4.375            5.52550902632509  11
-1  2.47896579998472  5.9762205702964  6.46996783163006  -1
-7  2.74448051533946  4.19010191002141 5.00890480199445  -7
-2  3.87961256380478 -2.07134405386938 4.3979381344821   -2

地震応答スペクトル計算のためのコードと文書

入力地震波データ

FilenameDescription
dat_acc_TCGH16_EW2.txtInput data sample (1)
dat_acc_TCGH16_NS2.txtInput data sample (2)
dat_acc_TCGH16_UD2.txtInput data sample (3)

入力地震波時刻歴

FilenameDescription
R_EQdat.txtR-script for drawing (oroginal acc.wave)
fig_R_acc_org.pdfImage sample (original acc.wave)

フーリエスペクトル

FilenameDescription
R_FSP.txtR-script for calculation of Fourier spectrum
R_NL_graph.txtR-script for drawing (Fourier spectrum)
dat_spc_1.txtOutput data sample (1)
dat_spc_2.txtOutput data sample (2)
dat_spc_3.txtOutput data sample (3)
fig_R_fsp_1.pngImage sample (original acc.wave 1)
fig_R_fsp_2.pngImage sample (original acc.wave 2)
fig_R_fsp_3.pngImage sample (original acc.wave 3)

加速度応答スペクトル

FilenameDescription
R_RSA.txtR-script for calculation of spectrum
R_LL_graph.txtR-script for drawing (response spectrum)
dat_spc_1.txtOutput data sample (1)
dat_spc_2.txtOutput data sample (2)
dat_spc_3.txtOutput data sample (3)
fig_R_spc.pdfImage sample (spectrum)


pic
inserted by FC2 system