はじめに

Rで開発を行っていると、途中までの実行結果を保存し後で呼び出ししたくなることがあった。その方法を備忘録として残しておく。

分析の例

データセットは「Rでの性別、身長、体重、血液型データの作成方法」で作成したものを「R shiny でデータの基礎分析用アプリの作成①」で csv に保存したので、それを用いて行う。

まずデータセットの読み込みと中身の確認を行う。

df <- read.csv("./data.csv")
str(df)

出力

'data.frame':	1000 obs. of  5 variables:
 $ X     : int  1 2 3 4 5 6 7 8 9 10 ...
 $ SEX   : chr  "female" "male" "male" "male" ...
 $ HEIGHT: num  154 178 163 176 156 ...
 $ WEIGHT: num  48 70.8 48.2 76.4 40.9 ...
 $ BLOOD : chr  "O" "A" "A" "O" ...

これから身長から体重を予測する線形回帰モデルを作成する。

result <- lm(WEIGHT~HEIGHT,data=df)
summary(result)

出力

Call:
lm(formula = WEIGHT ~ HEIGHT, data = df)

Residuals:
    Min      1Q  Median      3Q     Max 
-35.495  -7.378  -0.236   7.196  40.649 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) -85.32669    6.35537  -13.43   <2e-16 ***
HEIGHT        0.87842    0.03846   22.84   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 11.21 on 998 degrees of freedom
Multiple R-squared:  0.3433,	Adjusted R-squared:  0.3426 
F-statistic: 521.7 on 1 and 998 DF,  p-value: < 2.2e-16

この結果を保存する場合を考える。

変数（オブジェクト）を保存する場合

saveRDS関数を使って変数を保存する。

saveRDS(result, file = "./result.obj")

これで result.obj に変数の状態を保存できた。

保存した変数（オブジェクト）を呼び出す場合

呼び出すのも簡単で、readRDS関数を用いる。

result_new <- readRDS("./result.obj")
summary(result_new)

出力

Call:
lm(formula = WEIGHT ~ HEIGHT, data = df)

Residuals:
    Min      1Q  Median      3Q     Max 
-35.495  -7.378  -0.236   7.196  40.649 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) -85.32669    6.35537  -13.43   <2e-16 ***
HEIGHT        0.87842    0.03846   22.84   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 11.21 on 998 degrees of freedom
Multiple R-squared:  0.3433,	Adjusted R-squared:  0.3426 
F-statistic: 521.7 on 1 and 998 DF,  p-value: < 2.2e-16

オブジェクトごと読みだせているのが分かる。

おわりに

分析の結果を保存して途中で他の作業を行うことがあるので、saveRDS,readRDSは覚えておきたい関数の一つであった。今後も簡単な tips でも、何かあれば記事にしていきたいと思う。