2020-01-01から1年間の記事一覧

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 = "…

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(…

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,…

日米のGDP統計(SNAとNIPA)の細かな違いについて

財政政策(政府支出)のデータ分析に関する自分用のメモ: 日本のGDP統計は正確には国民経済計算体系(System of National Accounts)といい、国際的な基準(国連が定めるSNA)に準拠している。それに対してアメリカのGDP統計はNIPA (National Income and Pr…

dplyrのmutateでlagと差分(自分の使うところだけ)

color <- c("R","B","B","R","R") weight <- c(51,55,72,57,64) dfrm <- data.frame(COL=color, WEI=weight) dfrm %>% mutate(WEILag = lag(WEI)) メモ:lagを使うと、データは下に1段下がる(1行下がる)。 差分を計算するとき、diff()は使えず、G_delta = …

embedの使い方(自分の使うところだけ)

embedの使い方(自分用のメモ) Rにはembedという関数があり、ベクトル・時系列を指定期間ずらす作業に使える。言葉で聞くとなんだかわからないと思うので、たとえば以下のような時系列データがあったとする。 デモ用データ 使い方は以下のとおり demoTs <- …

成長率の計算(改訂版)

以下の記事の続編: 成長率の計算 - ハリ・セルダンになりたくて(はてダから移行中) 少し改訂版を作ったので記録のために以下に(前期比・前年同期比も計算できるように)。 growthRate <- function(tsData, lagLen=1, percent=FALSE) { numSize <- length…

Rでダミー変数を作る(fastDummies)

そもそもlm()(とglm())には「data.frameの列がfactorである場合、自動的にダミー変数を作って、回帰してくれる」という機能がある(manにきちんと書かれてないので知らない人が多いかも)。 stats.idre.ucla.edu が、lmやglmの「自動でダミー変数を作る機…

openxlsxの使い方(自分が使うところだけ)

require(openxlsx) tmpdat<- read.xlsx("./data.xlsx", sheet = 2, startRow = 9, detectDates=TRUE)head(tmpdat) wb <- createWorkbook()addWorksheet(wb, 'Sheet1')addWorksheet(wb, 'Sheet2')writeData(wb, sheet = 'Sheet1', x = lm.res1, withFilter=F)…