成長率の計算

Rを使って成長率(%表示)を計算する方法についての備忘録。 前提:tsDataをベクトルもしくはtsオブジェクトとする。 その1:対数差分 diff(log(tsData))*100 メリット:成長率の近似であるが、プログラムが簡単 デメリット:tsDataの要素が正である必要が…

RでX-13ARIMA-SEATSを使う(2019年6月版)

自分用の覚書CRANにあるx13binaryが使えるようになったので、自分でX-13ARIMA-SEATSをコンパイルする必要がなくなり、MacとWinの違いもなくなりました。 [インストール]GNU Rで二つのパッケージをインストール install.packages("seasonal"); install.packag…

GNU R (Windows10)にR Interface to 'Keras'をインストールする

GNU R (実際の作業はRstudio)でKerasをインストールしたときの覚書。基本的には(1)Anacondaインストール、(2)Anacondaでtensorflowインストール、(3)R上でdevtoolsインストール、(4)devtoolsでkerasインストール、(5)Windows上でRtoolsインス…

Rstudio (Windows10)でRのでパッケージをインストールする

Rstudio (Windows10)でRのでパッケージ(devtoolsなど)をインストールするとき、パスに日本語がある場合、Rstudioが認識できない(日本語パスの部分が????になる)。 libpathsでパスを指定する方法もあるが、一番簡単なのはRstudioを管理者権限で起動してイ…

matlab 警告: 名前が存在しないか、ディレクトリではありません

Windows10でDynareをインストールしたあとに「matlab 警告: 名前が存在しないか、ディレクトリではありません」という警告がでる問題の解決方法。 cd C:\dynare\4.5.7\contrib\ms-sbvar\TZcode\MatlabFiles\ edit startup.m ms-sbvarを使わない場合は、以下…

シンボリック関数を関数ハンドルに直して、arrayfunを適用する

arrayfunの使い方 - ハリ・セルダンになりたくて(はてダから移行中) と Matlabのシンボリック式を無名関数(関数ハンドル)に変換する関数 - ハリ・セルダンになりたくて(はてダから移行中) の合わせ技。 sf(x,y)というシンボリック関数があった場合、そ…

ある数字の繰り返しである行列を作成するのに効率的方法は?

たとえば、7が単に並んだ列ベクトルをつくるとき、Matlabで一番効率的な方法はなにか? tic; repmat(7, 1000000,1); toc;経過時間は 0.006551 秒です。tic; ones(1000000,1) * 7; toc;経過時間は 0.008067 秒です。 repmatの方が微妙に早い?のか?

arrayfunの使い方

Matlabに詳しい人には当たり前の話だと思うがarrayfunの使い方を自分用にメモ(Rではapply系の関数に該当するのがarrayfun、cellfun、bsxfunだと思えばいいのかな・・・): plusFun = @(x, y) x + 2 * yplusFun(1,2)arrayfun(plusFun, [1 2], [3 4])arrayfu…

Symbolic mathで微分演算子っぽいものを作る

Symbolic mathで微分演算子っぽいものを作ることができる。なんの役に立つか不明だが、面白い機能なので覚書。 [問題の定式化] 2回微分可能な関数について以下の簡単な常微分方程式を考える。 微分演算子を以下のように定義する。 微分演算子を使えば、常微…

Symbolic mathではシンボリック関数の配列が作成できる。微分もできる。

使い方次第だと思うが、面白い機能なので覚書。 シンボリック関数はsymfunで作ることができる(フォーマルなやり方)。簡易なやり方としてsym f(x,y)のようにすることもできる(xxとyyは自動生成される)。 syms phi1(xx) phi2(xx); phi1(xx) = xx; phi2(xx)…

シンボリック式・シンボリック関数はstringで文字列に変換できる

何に役立つのか分からないが、面白い機能なので覚書(この機能はMatlabドキュメントに載っていないようなので)。 syms xx yy; % シンボリック式はstringで文字列に変換できるstring(xx + yy) % シンボリック関数もstringで文字列に変換できるtestFun = @(xx…

Matlabのシンボリック式を無名関数(関数ハンドル)に変換する関数

Matlabのsymbolic math toolboxを使わない人にとってはどうでもいい話題ですが、シンボリック式を最適化しようとして少し困ったので、覚書。 Matlabのシンボリック式またはシンボリック関数を無名関数(関数ハンドル)に変換する関数が存在する。使い方はこ…

Miranda-Fackler (2002) p. 132のテンソル積による計算結果がMatlabのkronの結果と同一であることの確認。

% Miranda-Fackler (2002) p. 132のテンル積による計算結果がMatlabのkronの結果と同一であることの確認。 syms x1 x2; aa = [ 1, x1, x1^2] bb = [ 1, x2] kron(bb, aa) 注:ベクトルの場合、テンソル積とクロネッカー積の結果は同一となる(だったはず)。…

MacOSでAnaconda (python)のバージョンアップ

condaを使う場合、 $ conda update conda $ conda update anaconda $ conda update --all Anaconda Navigatorを起動して、バージョンが最新版になっていることを確認する。TensorFlowとKerasのインストール $ conda install -c jjhelmus tensorflow $ conda …

RでX-13ARIMA-SEATSを使う(2019年6月版)

自分用の備忘録です。CRANにあるx13binaryが使えるようになったので、自分でX-13ARIMA-SEATSをコンパイルする必要がなくなりました。[インストール]GNU Rで二つのパッケージをインストール install.packages("seasonal"); install.packages("x13binary") [GN…

macOS Sierraでgfortranを使う

Mac

方法(その1):homebrewでbrew install gccを実行する(かなり前にgfortranはgccに統合された)方法(その2):CRANに行き、"Download R for (Mac) OS X"をクリック、"tools"をクリック、そこにある"gfortran-4.2.3.pkg"をインストールする

Rでxlsxを読み書きする

以下のサイトが非常に分かりやすかったです。 R で openxlsx パッケージを用いて xlsx ファイルを扱う - Qiita http://qiita.com/zero310/items/71b6e5bc0731680d19dfRで解析:セル体裁!大きなデータも大丈夫!エクセル操作の「openxlsx」パッケージ http:/…

マイナス金利に関するまとめ(2016/09/23初稿、09/26修正)

以下、自分用の覚書: 概況: 2016年9月23日現在、「超過準備に対するマイナス金利(negative interest rate on excess reserves)」を実施しているのは日銀(BOJ)、デンマーク国立銀行(DNK)、欧州中央銀行(ECB)、スウェーデン中央銀行(Riksbank)、ス…

GNU Rで月次データを四半期に集約する

時系列解析で月次データを四半期に集約する(3ヶ月の合計や平均を計算する)ことはよくある。以下はstats, zoo, xtsを使った場合。まず準備: # 2010年1月から12月までの月次データを生成する rm(list=ls()) dataLength set.seed(1) # to make it reproducib…

Emacsでjulia-mode

ESS (Emacs Speaks Statistics) にはすでにjulia-modeが入っているはずなのですが、OS Xでうまく動かないので、しょうがなくJulia-lang公式のjulia-mode.elを入れてみました。 Julia-lang公式のgithubから julia-mode.el をダウンロード(ページの「RAW」を…

ESS (Emacs Speaks Statistics) をgitからインストール

cd ~/bin/ git clone https://github.com/emacs-ess/ESS.git cd ess make ~/.emacsに以下のように記述*1 (add-to-list 'load-path "~/bin/ess/lisp") (load "ess-site") 以上でemacs上でM-RでRを起動できるようになる*2。 *1:~/.emacs.d/init.elに書いておく…

OS X 10.9とMatlab R2016aでgfortranを使う

(1)OS X 10.9 (Mavericks)にgfortranをインストール Xcodeをインストール:App Storeの検索窓に"xcode"と入力して、xcodeを探して、インストールする Xcode 6.1以降ではcommand-line toolsがプレインストールされるようになった模様 [参考] Mac での Xcod…

Stateでダミー変数の作り方(の一例)

たとえばこんな感じ(詳しいことは出典参照) gen year1=1 if year==2011 replace year1=0 if year1 == . [出典] http://www.stata.com/support/faqs/data-management/creating-dummy-variables/

varsでのベクトル自己回帰の結果をxtableでLaTeXに出力する

変数それぞれに対して係数の推定結果を出力することは可能*1。たとえばこんな感じ: require(vars) data(Canada) var.canada summary(var.canada) ここまではvarsの機能。Canadaはe, prod, rw, Uの4変数からなるデータセット。たとえばUに関する回帰式の推定…

install.packagesでCRANを指定する

optionsコマンドでreposを指定する。 options(repos="http://cran.ism.ac.jp/")

配列の1列目をインデックスとしてソートする

Matlabで配列の1列目をインデックスとしてソート(並べ替え)するたとえば A = [2 5 8 ; 3 4 7 ; 1 6 9] という配列の1列目(1列目の[2;3;1])をインデックスとして昇順でソートする。つまり以下のような配列Cを得たい。 C = [1 6 9 ; 2 5 8 ; 3 4 7] この場…

OS X (10.9.5)にHomebrewでEmacs 24.5をインストール

$ brew install --japanese --cocoa -v emacs なお、Emacs 24.4は/usr/local/Cellar/emacs/24.5配下にインストールされる。ちなみに/usrは不可視になっているので、termianl上で以下のように入力すればフォルダが開ける: open /usr/local/Cellar/emacs/24.5…

OS Xでアプリ内のウィンドウ切り替え

IT

通常は「Command + F1」もしくは「Command + Fn + F1」だが、面倒(F1がキーボードの上の方にあるので届きにくい&矢野は手が不器用)なので、キーアサインを入れ替えることにした。とりあえず「Command + E」に設定してみた。後日、再設定する可能性あり。…

家計・所得・格差に関する統計

国内 家計調査(総務省統計局) http://www.stat.go.jp/data/kakei/ 国民生活基礎調査(厚生労働省) http://www.mhlw.go.jp/toukei/list/20-21.html 統計情報(国税庁)https://www.nta.go.jp/kohyo/tokei/kokuzeicho/tokei.htm 国民経済計算(内閣府) htt…

GNU Rのパッケージvarsのvarresultから係数を取り出す

varresultから係数を取り出す。tmp tmp["enlargePrisonTime"]