R、Pythonの勉強の過程をメモっていくノート。2020-01-23対数尤度が正で、AICが負になるケース初歩的な話ですが、今日学生から訊かれたので例を考えてみました。AICがマイナスの値になってもいいんだっけ?みたいな話です。統計モデルの最尤推定をする場合、尤度は1より小さい場合が多く、したがって対数尤度が負である場合が多いと...対数尤度が正で、AICが負になるケース2020-01-14researchmapに投入するデータをExcelからCSVに変換するスクリプト(Pythonとbash)以前、reserchmapにCSVでデータを投入するときの注意点を書きました。私は論文の業績はあまり無いので(笑)管理が楽なのですが、一般の雑誌に書いてる文章...researchmapに投入するデータをExcelからCSVに変換するスクリプト(Pythonとbash)2020-01-10Rでピボットテーブル学生の分析を手伝っていたところ、突然、Rでピボットテーブルみたいな集計をしたくなりました。度数分布表はtable関数でつくれますが、ピボットテーブルってどうやるのかな、と。最近、パネルデータを扱うことが増えてきたのだが、パネルデータで動学的な(つまりt-1期とかのラグ項が出てくる)分析をやろうとすると最小二乗法ではなくGMM推定量を用いる必要がある。備忘として、動学的パネルデータ...Rで動学的パネルデータ分析:plm、panelvarパッケージをつかったGMM推定2019-12-20Rで忘れがちな書き方(随時追記)よく忘れることのメモです。要素に含まれるかどうかの判定は、%in%かis.element()を使う。これはデータ全体の中から何かを抽出するときの条件を複数条件にしたい場面でも使えるときがある。Codeのアプデが必要だった)さっきハマったエラーを解決したので、備忘のためのメモです。RStudioとRの両方で試したのですが、{panelvar}というパッケージをインストールしようとしたら、依存パッケージのコンパイルのエラーがでて進めなくなりました。Codeのアプデが必要だった)2019-11-19Macで画像ファイルをまとめてPDF化する(ただし名前順にソート)Macで画像ファイルをまとめてPDFにしたいとき、従来はAutomatorを使ってワークフローを組んでいました。ところが、いつからか知らないですが、OSに標準でそういう機能が付けられていたようです。Q1"みたいな)を、文字列情報として取り出す方法が、ぱっとググって分からなかったのですが、とりあえず以下のようにしたらできました。.05」で「統計的に有意」なのか?――5%基準の由来について以前も書いたんですが、p値が0.05を下回るかどうかにとらわれる慣習を問題視する人が最近は増えてきていて、たしかにその理由はよく理解できる一方で、p値が過剰にバッシングされている気もします。.05」で「統計的に有意」なのか?――5%基準の由来について2019-08-17Macでひらがな・カタカナの濁点・半濁点が分離してしまったのを元に戻すスクリプト昨日のエントリで、テキストエディタの「CotEditor」のスクリプト機能の使い方を書きました。その機能をつかって、濁点・半濁点が分離してしまったテキストを元に戻すツールを作ります。いうやつをなんとかしたいわけです。英語のタブからダウンロードしたフォーマットを使っているなら、英語のタブからインポートしなければならない。説明書きに、情報がないセルにはnullを入れろと書いてあるが、べつに入れなく...researchmapにCSVで論文のデータを投入するときの注意点2019-05-15Rで二段階最小二乗法の演習(『gretlで計量経済分析』を参考に)『gretlで計量経済分析』という教科書があって、入門的な統計分析を非常にわかりやすく解説していると思いました。さっき就活エントリを書いたついでに、10年ぐらい前に別のブログに書いたネタを再掲しておこう。とある就活サイトの人から、リレーブログみたいなものへの参加を依頼されたので、今日は統計ではなく就活の話を書くことにします。この時は文中にも書いているとおり、差出人名を日本語表示するのがうまく出来ませんでした。ところがその問題は、下記の方法で解決しました。要するに、'...Python3でのメール送信時に日本語の差出人名を使う2018-05-24東京一極集中に関するデータ(ソースのメモ)以下は、分析とかではなくデータのソースについてのメモである。人口に関する「東京一極集中」がいかに凄いかは、世銀のサイトにある以下のデータをみるのが分かりやすい。主要国では、日...東京一極集中に関するデータ(ソースのメモ)2018-05-04RでVARモデル&インパルス応答関数を求めるとき、「ショックの値」をどうやって出すのか?RでVARモデルを推定してインパルス応答関数を出す時に、インパルス応答関数が対応しているところの「ショック」の大きさが幾らなのかを、どこから得たらいいのかという疑問がわきました。test)を実施してくれるadf.test()という関数があります。考慮するラグの次数を指定...Rでの単位根検定はadf.test()関数よりCADFtest()関数がいいのでは?2018-02-18Rで要素番号の指定の仕方をミスった考えてみればそりゃそうか、という感じではあるのですが、またいつかミスりそうなのでメモしておきます。以下のような感じで、startとendの値を変えて適切な期間を取りたいとします。メールを300人ぐらいに発信する必要がありまして、Toに全員入れるわけにはいかないし、BCCで送るのもダサいかなと思って、「1人1人を個別にToに指定して、同じ件名・同じ文面のメールを送る」ってのをPythonでやってみました。標準モ...Pythonの簡単なコードでメールを自動送信してみる2018-01-03Macでの年賀状作成環境について(2018年版)ブログの趣旨と全然違いますが、備忘のためにまとめておこうと思います。私は、サボる年もありますが、出す時は年賀状を200人以上に出すので、けっこうな大仕事になっています。Rのバージョンアップをしようとしたんですが、Macの場合はアップデート用のコマンドが使えないらしいので、アンインストールして新バージョンを再インストールをすることにしました。Pandas前回のエントリで、Pandasのデータフレームを文字列で検索すると遅かったと書いたんですが、実際に検索の方法を変えるとどれぐらい違うのかを計測してみました。前回のエントリで、著作権侵害にあたる違法アップロード動画を自分のTwitterで拡散してしまっている懸念を考えて、YouTube動画のリンクが貼ってあるツイートをまとめて削除しました。ツイッターで、YouTubeのリンクを貼り付けたツイートをしていますが、よく考えたら違法アップロードに該当するものを拡散してしまっている可能性があります。R、Pythonの勉強の過程をメモっていくノートです。初心者なのであまり信用なさらぬよう。間違いについてはコメント等でご指摘頂ければ幸いです。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。