Rの質問です。
STATAではlogitやxtlogit（ランダム効果）でモデルを走らせたあと
mfx compute, eyex
などとすると各説明変数の平均でのelasticity（弾力性）が計算されます。
タバコの価格を説明変数と喫煙（0,1）を被説明変数としているので，
是非elasticity（価格が1%変わったとき，喫煙が何パーセントかわるか）
を計算したいのですが，R-projectのサイトで検索しても，
Rでどうやったらいいのかわかりません，ご教授お願いします。
STATAでやればいいですが，是非Rでやってみたいので，よろしくお願いします。
手で計算するにはすべて平均値を入れてlogit(p/1-p)を計算し，pを出して，
注目の変数の+1%と-1%を入れてp(+1)やp(-1)を出しなおして平均とれば
いいのでしょうか? もっとエレガントな方法ありますでしょうか? |