2013年6月より統計解析向けプログラミング言語「R」を使い始めました． きっかけは，仕事上のニーズです．しばらく前に仕事で簡単な統計解析のレポートを書いたのですが，そのときはエクセルとFortranプログラムで対応しました．しかしもう少し複雑な解析になってくると自分のプログラムでは心もとないし，エクセルを使って複雑なことをやる気はない，ということで，前から名前だけ知っていた「R」を試してみようと思った次第です．
使用環境はインストールしたままですが，自分が使う中での注意事項は今のところ以下のとおりです．
| Contents |
|---|
| 重回帰分析 |
| 確率楕円 |
| 主成分分析 |
| 多次元尺度法 |
| クラスター分析 |
| 判別分析 |
| 水文統計解析 |
| FFT を用いた地震応答スペクトル計算 |
| Filename | Description |
|---|---|
| R_MRA.txt | R-script (using 'lm') |
| R_MRAorg.txt | R-script (using 'solve') |
| inp_MRA.txt | Input data sample (csv) |
| fig_R_MRA.pdf | Image sample |
| Filename | Description |
|---|---|
| R_ELLIPSE.txt | R-script (using 'ellipse') |
| inp_reg1.txt | Input data sample |
| fig_R_ellipse.pdf | Image sample |
| Filename | Description |
|---|---|
| R_PCA.txt | R-script (using 'princomp') |
| R_PCAorg.txt | R-script (using 'eigen') |
| inp_iris.txt | Input data sample |
| fig_R_PCAorg.pdf | Image sample |
| Filename | Description |
|---|---|
| R_MDS.txt | R-script (using 'cmdscale') |
| R_MDSorg.txt | R-script (using 'eigen') |
| inp_iris.txt | Input data sample |
| fig_R_MDSorg.pdf | Image sample |
| Filename | Description |
|---|---|
| R_CLUSTER.txt | R-script (using 'kmeans') |
| R_CLUSTERorg.txt | R-script (using 'eigen') |
| inp_iris.txt | Input data sample |
| fig_R_CLUSPLOT.pdf | Image sample |
| fig_R_CLUSorg.pdf | Image sample |
| Filename | Description |
|---|---|
| R_DCA.txt | R-script (using 'lda') |
| R_DCA_MAHA.txt | R-script (using 'mahalanobis') |
| inp_iris.txt | Input data sample |
| fig_R_DCA.pdf | Image sample (by 'lda') |
| fig_R_MAHA.pdf | Image sample (by 'mahalanobis') |
| Filename | Description |
|---|---|
| inp_RF_S.txt | 札幌雨量データ |
| inp_RF_M.txt | 前橋雨量データ |
| inp_RF_K.txt | 京都雨量データ |
| Filename | Description |
|---|---|
| R_HF_qq_LN3_IWAI.txt | R-script (LN3, IWAI) |
| R_HF_qq_LN3_mom.txt | R-script (LN3, moment) |
| R_HF_qq_LP3_mom.txt | R-script (LP3, moment) |
| R_HF_qq_GEV.txt | R-script (GEV) |
| R_HF_qq_Gum.txt | R-script (Gumbel) |
| R_HF_qq_SQRT.txt | R-script (SQRT-ET) |
| fig_R_HF_qq_LN3_IWAI.png | Image sample |
| fig_R_HF_qq_LN3_mom.png | Image sample |
| fig_R_HF_qq_LP3_mom.png | Image sample |
| fig_R_HF_qq_GEV.png | Image sample |
| fig_R_HF_qq_Gum.png | Image sample |
| fig_R_HF_qq_SQRT.png | Image sample |
| Filename | Description |
|---|---|
| R_HF_pdf_GEV.txt | R-script (Histogram & GEV pdf) |
| fig_R_HF_hist_GEV.png | Image sample |
| Filename | Description |
|---|---|
| R_HF_paper_GUM.txt | R-script (Gumbel probability paper) |
| fig_R_HF_paper_GUM.png | Image sample |
| Filename | Description |
|---|---|
| R_HF_RT_LN3_IWAI.txt | R-script (Rejection test, LN3, IWAI) |
| R_HF_RT_LP3.txt | R-script (Rejection test, LP3) |
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
| dat_acc_TCGH16_EW2.txt | Input data sample (1) |
| dat_acc_TCGH16_NS2.txt | Input data sample (2) |
| dat_acc_TCGH16_UD2.txt | Input data sample (3) |
| R_EQdat.txt | R-script for drawing (oroginal acc.wave) |
| fig_R_acc_org.pdf | Image sample (original acc.wave) |
| R_FSP.txt | R-script for calculation of Fourier spectrum |
| R_NL_graph.txt | R-script for drawing (Fourier spectrum) |
| dat_spc_1.txt | Output data sample (1) |
| dat_spc_2.txt | Output data sample (2) |
| dat_spc_3.txt | Output data sample (3) |
| fig_R_fsp_1.png | Image sample (original acc.wave 1) |
| fig_R_fsp_2.png | Image sample (original acc.wave 2) |
| fig_R_fsp_3.png | Image sample (original acc.wave 3) |
| Filename | Description |
|---|---|
| R_RSA.txt | R-script for calculation of spectrum |
| R_LL_graph.txt | R-script for drawing (response spectrum) |
| dat_spc_1.txt | Output data sample (1) |
| dat_spc_2.txt | Output data sample (2) |
| dat_spc_3.txt | Output data sample (3) |
| fig_R_spc.pdf | Image sample (spectrum) |
