% Miranda-Fackler (2002) p. 132のテンル積による計算結果がMatlabのkronの結果と同一であることの確認。 syms x1 x2; aa = [ 1, x1, x1^2] bb = [ 1, x2] kron(bb, aa) 注：ベクトルの場合、テンソル積とクロネッカー積の結果は同一となる（だったはず）。…
condaを使う場合、 $ conda update conda $ conda update anaconda $ conda update --all Anaconda Navigatorを起動して、バージョンが最新版になっていることを確認する。TensorFlowとKerasのインストール $ conda install -c jjhelmus tensorflow $ conda …
自分用の備忘録です。CRANにあるx13binaryが使えるようになったので、自分でX-13ARIMA-SEATSをコンパイルする必要がなくなりました。[インストール]GNU Rで二つのパッケージをインストール install.packages("seasonal"); install.packages("x13binary") [GN…
方法（その１）：homebrewでbrew install gccを実行する（かなり前にgfortranはgccに統合された）方法（その２）：CRANに行き、"Download R for (Mac) OS X"をクリック、"tools"をクリック、そこにある"gfortran-4.2.3.pkg"をインストールする
以下のサイトが非常に分かりやすかったです。 R で openxlsx パッケージを用いて xlsx ファイルを扱う - Qiita http://qiita.com/zero310/items/71b6e5bc0731680d19dfRで解析：セル体裁！大きなデータも大丈夫！エクセル操作の「openxlsx」パッケージ http:/…
以下、自分用の覚書： 概況： 2016年9月23日現在、「超過準備に対するマイナス金利（negative interest rate on excess reserves）」を実施しているのは日銀（BOJ）、デンマーク国立銀行（DNK）、欧州中央銀行（ECB）、スウェーデン中央銀行（Riksbank）、ス…
時系列解析で月次データを四半期に集約する（3ヶ月の合計や平均を計算する）ことはよくある。以下はstats, zoo, xtsを使った場合。まず準備： # 2010年1月から12月までの月次データを生成する rm(list=ls()) dataLength set.seed(1) # to make it reproducib…
ESS (Emacs Speaks Statistics) にはすでにjulia-modeが入っているはずなのですが、OS Xでうまく動かないので、しょうがなくJulia-lang公式のjulia-mode.elを入れてみました。 Julia-lang公式のgithubから julia-mode.el をダウンロード（ページの「RAW」を…
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 (Mavericks)にgfortranをインストール Xcodeをインストール：App Storeの検索窓に"xcode"と入力して、xcodeを探して、インストールする Xcode 6.1以降ではcommand-line toolsがプレインストールされるようになった模様 [参考] Mac での Xcod…
たとえばこんな感じ（詳しいことは出典参照） gen year1=1 if year==2011 replace year1=0 if year1 == . [出典] http://www.stata.com/support/faqs/data-management/creating-dummy-variables/
変数それぞれに対して係数の推定結果を出力することは可能*1。たとえばこんな感じ： require(vars) data(Canada) var.canada summary(var.canada) ここまではvarsの機能。Canadaはe, prod, rw, Uの4変数からなるデータセット。たとえばUに関する回帰式の推定…
optionsコマンドでreposを指定する。 options(repos="http://cran.ism.ac.jp/")
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] この場…
$ brew install --japanese --cocoa -v emacs なお、Emacs 24.4は/usr/local/Cellar/emacs/24.5配下にインストールされる。ちなみに/usrは不可視になっているので、termianl上で以下のように入力すればフォルダが開ける： open /usr/local/Cellar/emacs/24.5…
通常は「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…
varresultから係数を取り出す。tmp tmp["enlargePrisonTime"]
ssc install outreg2
Emacsを起動すると「ESS（のlispが）が読み込めない」というエラーがでるようになった（最近EmacsとESSを使っていなかったので、いつからこのエラーが出ているのか不明）*1。原因を探るのが面倒なので、ホームディレクトリ配下にESSを保存・コンパイルするこ…
注意：内容が古いので、新しいエントリーをご覧ください（2017年5月19日記）。自分用の備忘録：[UNIX, OS Xで自分でコンパイルして使う場合] US CENSUS（米国の国勢調査局）からtar -zxvf x13assrc.tar.gzをダウンロードするhttps://www.census.gov/srd/www/…
自分用の備忘録：United Nations Statistics Division - National Accounts http://unstats.un.org/unsd/snaama/Introduction.asp総務省統計局、「国際連合（UN）National Accounts Main Aggregates Databaseの使い方」http://www.stat.go.jp/data/sekai/pdf…
ess-14.09.tgzをダウンロードし、以下の通り解凍する。 tar -zxvf ess-14.09.tgz HomebrewでインストールしたEmacs 24.4.1は/usr/local/Cellar/emacs/24.4配下にあるため、 ess-14.09をインストールするとき、MakeconfでPREFIXを以下の通りに書き換え、make;…
$ brew install --japanese --cocoa --with-gnutls -v emacs なお、Emacs 24.4は/usr/local/Cellar/emacs/24.4配下にインストールされる。[参考]http://www.sevencredit.com/2014/07/02/580/2015年7月20日追記：HomebrewでEmacs 24.5が使用可能です。http://…
毎度毎度忘れるので自分用にメモ： > nohup matlab -nodesktop -nodisplay output.txt 2>&1 &
CentOSに飽きてUbuntuに切り替えたため、CUDA toolkitをインストールする必要があり、その作業メモ。[インストール・ガイド]以下のLinux Getting Started Guideの通りにするだけ：http://developer.download.nvidia.com/compute/cuda/6_5/rel/docs/CUDA_Gett…
「CentOS 5が飽きたので、Ubuntuにしよう」と思ってUbuntu 14.04にしたらNVIDIA GeForce GTX 580のドライバでちょっと困ったので自分用のメモ：[問題点]NVIDIA GeForceのドライバがインストールされていない場合、しばらくするとハングアップする（Ubuntuに…
毎回、使い方を忘れるので、自分用の覚書。# クラスターでの並列プールの作成parpool('local',2);#これはローカル・プロファイルを使用して、2個のワーカーからなる並列プールを作る場合 # 走らせたい処理をparforで書く parfor i=1:1024 A(i) = sin(i*2*pi/…
Twitterでのshavetailさんの図に少し違和感があったので、自分でも似たような図を作ってみました。（Twitterのこの呟きの図を少し直したもの）
dvipdfmxオプションつければズレなかった・・・ \usepackage[dvipdfmx]{graphicx} [出典]http://heilman-jinroku.blogspot.jp/2013/10/matlabepsdvipdfmxpdf.html dvipdfmxオプションをつければ以下のコードは必要なし。 #!/bin/sh for f in `find . -name "…
