imputeTSを使って補完
na_interpolation(x = dataset$variable_name, option = "linear")
以下のサイトの情報から自分が使うところだけ抜き出した：
参考文献：
RでStataのdtaを読み込む
いくつか方法があるが、パッケージhavenを使う方法が比較的手軽に使えるようである。以下は覚書。
require(haven)
url <- "http://www.principlesofeconometrics.com/stata/broiler.dta"
data.df <- read_dta(url)
参考文献：
Rstudioでソースのあるディレクトリに移動する方法
Rstuidoを使っている人限定：
setwd(dirname(rstudioapi::getActiveDocumentContext()$path))
参考文献
dplyrのselectとlagが衝突する場合の対応
メモ：最初に書いた内容に間違いあったため、修正。
一年ほど前に書いたRのコード（dplyr）が動かないことがわかった。コードを追ったところ（１）dplyr::selectがMASS::selectと衝突、（２）dplyr::lagがstats::lagと衝突していることが分かった（1年前は問題なく動いたような気がするのだが・・・）。
やり方としては参考文献にあるように
select <- dplyr::select
と上書きする。
覚書１：「tidyverseを最後にロード（library(tidyverse) ）しても関数の衝突は解決しない」（場合がある？）ようだ。そのため上書きするのが確実な模様。
覚書２：Tidyverseのページに衝突の話はあった。
Conflicts between the tidyverse and other packages — tidyverse_conflicts • tidyverse
tidyverse_conflicts()
（これが分かるだけで30分以上費やし、なんだか少し悔しかったので、自戒のためにブログに書いておくことにする。将来、他の誰かの役に立ちますように。）
huxtableのhuxregを使ってlmの結果をワードに出力
自分用覚書：huxregと quick_docxは便利なのでオススメ
install.packages("huxtable")
require("huxtable")
lm.res1 <- lm(mpg ~ wt , data = mtcars)
lm.res2 <- lm(mpg ~ wt + qsec, data = mtcars)
lm.hux <- huxreg(lm.res1, lm.res2, number_format = "%.4f")
lm.hux
quick_docx(lm.hux, file="lmHuxSample.docx")
以下のマニュアルを読むと結構楽しい。
https://cran.r-project.org/web/packages/huxtable/vignettes/huxreg.html
jtoolsを使ってlmの結果をワードに出力
後日のメモ：huxtableのhuxregを使えばいいので、jtoolsは不要かもしれない。
install.packages("flextable")
install.packages("officer")
install.packages("huxtable")
require(“jtools”)
lm.res <- lm(mpg ~ wt + qsec, data = mtcars)
lm.summ <- summ(lm.res)
export_summs(lm.summ, to.file = "docx", file.name = "test.docx")
tidy()を使ってopenxlsxで結果を出力
library(ggplot2)
library(dplyr)
library(broom)
library(openxlsx)
lm.res <- lm(mpg ~ wt + qsec, data = mtcars)
lm.res %>% summary() %>% tidy() -> lm.tidy
lm.res %>% summary() %>% glance() -> lm.glance
wb <- createWorkbook()
addWorksheet(wb, 'Sheet1')
addWorksheet(wb, 'Sheet2')
writeData(wb, sheet = 'Sheet1', x = lm.tidy, withFilter=F)
writeData(wb, sheet = 'Sheet2', x = lm.glance, withFilter=F)
saveWorkbook(wb, "lm.res.xlsx", overwrite = T)
参考文献