-1*f()より-f()

R言語

Rで「-1*f()が気持ち悪い」と言う指摘を受けて、-f()だとSyntax errorになりそうで逆に気持ち悪いんじゃないかと思っていたら、年寄りの戯言でした。どういう事かと言うと、この二つがRではS式が異なるわけです。 f <- function(n){ n + 1 } as.list(quote(…

2012-02-11

加算後の除算≿除算後の加算

R言語

かなり大量に演算しないと目立ちませんが、500万サンプルで比較すると、0.10秒と0.16秒で差が計測できました。 obs <- 5000000 x <- runif(obs) y <- runif(obs) ### まとめて割り算 ### gc();gc() system.time({ n <- x + y m <- n / 2 print(sum(m)) }) ##…

2012-02-11

Rは行列演算≿apply関数≿ループ演算

R言語

Rのプロビット回帰を用いたマイクロ・ベンチマークを書いたので、ソースコードを公開しておきます。 1. セットアップ推定結果が常に同じになるように、セットアップします。 ### セットアップ ### set.seed(20120209) obs crt x1 x2 y frml # 説明変数を行…

2012-02-11

ifelse()と論理値の掛け算

R言語

Rでは論理値TRUEとFALSEは、それぞれ1と0として演算する事ができます。ゆえにyとzがxの値に関わらず計算可能な値になるのであれば、ifelse()の代わりに論理値で演算を行う事ができます。 ifelse()はベクター演算にならないので、論理値を使う方が速いです。…

2012-02-11

配列の拡張とパフォーマンス

R言語

動的に配列サイズが変更されるほど、パフォーマンスは劣化します。 num = 50000### 配列を拡張していく ### gc(); gc() system.time({ x <- numeric(0) i <- 1 while(i <= num){ x[i] <- i i <- i + 1 } })### 配列を一度に初期化 ### gc(); gc() system.tim…

2012-02-10

Rでトービット・モデルによる打ち切りデータの推定

R言語

金額で測った不動産の保有量を考えましょう。大半の貧乏人はゼロでありますが、一定以上の所得がある人から保有量が増加してきます。集合住宅、持ち家、そして別荘や投資物件など増加していくわけです。こういう説明変数の値に対して、一定まではゼロ（もし…

2012-02-04

地獄のミサワ的Rユーザーのソースコード

R言語

別のブログでアイナちゃんの課題をミサワ君とサエコ先輩が助ける話を書いたのですが、検証用にその物語中のソースコードを公開しておきます。 1. データセットの作成データセットは課題で出されたモノになりますが、ファイルを用意するのも面倒なので、今回…

2012-01-28

Rでマルコフ連鎖モンテカルロ法を試す

R言語

地味にここ5年間ぐらい、マルコフ連鎖モンテカルロ法（MCMC）が流行っているようです。汎用的な分布でベイズ推定を行う時に有用な数値解析アルゴリズムの総称で、Metropolis-Hastings algorithm（M-Hアルゴリズム）などが主要なメソッドとして使われています…

2012-01-21

Rによる数値微分

R言語

Rはderiv()関数で記号/微分ができるものの、和や積、つまりsum()やprod()があると微分できません。以下のようなコードを実行すると、『関数 'sum' は導関数の表中にありません』のようにエラーメッセージが出ます。 # ベクトルxを代入 x # 変数bを2とする b …

2012-01-10

Rと手作業で覚える最尤法

R言語

OLSより進んだ統計手法で最初に覚えるのは最尤法だと思います。大半の人はツールとして知っていて、あまり中身を意識していない気がするのですが、「尤度」の説明無しで『尤度が最大になるパラメーターを求める方法』と言う説明が横行しているのは、問題があ…

2011-12-27

Rでニュートン・ラフソン法

R言語

文系の学生などで特に目的無くプログラミングを覚えたい人には、Rが最適です。統計解析に使わなくても、ちょっと計算した気になるプログラミングの練習をしてみましょう。 1. よくある近似解を求める関数f(x)=0があるとして、テーラー展開をしてみます。f(x…

2011-12-08

Rで年次ダミーを作成

R言語

実際の計量分析ではダミー変数は無くてはならないテクニックで、技術進歩の計測でもしない限りは、複数年にまたがるデータには年次ダミーが含まれている事が通例です。Rでは簡単にダミー変数を作る事ができます。例えば以下のようなデータフレームdfがある…

2011-11-20

Rでデータフレームに、t-1期のデータを追加する（ミクロデータ向き）

R言語

被説明変数が説明変数に影響を与えている事を同時性と言います。アイドルのCMやドラマ露出時間と、アイドルの人気度の関係などが例にあげられるでしょう。単回帰では問題ないのですが、重回帰分析では分析に深刻な影響を及ぼす可能性があります。これを解決…

2011-11-19

Rでバランスド・データの作り方

R言語

複数年に渡るミクロデータでは、個体が入ったり抜けたりします。例えば地域住人のデータでは、転居や死亡などで個体が減りますし、逆に転入や誕生で個体が増えたりします。サンプル数が膨大な場合は誤差になりますが、限定されたサンプルでは推定結果に影響…

2011-11-18

Rでミクロデータにマクロデータを紐付ける

R言語

ミクロデータへのマクロデータへの紐付けはデータ分析では良くある作業です。例えば企業の財務データを、その年のインフレ率で割り引く作業があります。Rの演算は各年各企業のデータにインフレ率が無いと演算しにくいので、Excelのvlookup関数で紐付けを済ま…

2011-06-08

PidginでYahoo!Japanに接続不能のときの対処方法

Windows

PidginでYahoo!Japanに接続できなくなりました。原因は、認証サーバーのIPアドレスをもどすマスター・サーバーが停止したためのようです。Pidginはマスター・サーバーしか問い合わせに行かないようで、サーバー応答が無いとメッセージが出て接続できなくな…

2011-05-24

今日は第何週？

R言語

R言語は、他の人気スクリプト言語に比較すると、ベクトル演算以外は記述が煩雑な気がします。特に日付関係は特殊なオブジェクトが乱立しており、少なくとも他の言語とは扱い方がちょっと違います。例えば、今日が第何週かを調べるには以下のようにすれば良…

2011-05-03

64bits版Linuxでは、64bits版Flashを使う

Linux

64bits版Ubuntu Linuxでは、64bits版Flashを使う方が良いようです。 Synapticパッケージマネージャーでは32bits版をインストールするのですが、ia32-libsをインストールするためか、C++/Qtでアプリケーションを作成すると、ELF class: ELFCLASS64エラーが出…