R を使用して、回帰、シミュレーションを実行し、単純なデータ セット (Zelig では「承認」) で有意性をテストしています。
統計的有意性をテストするまでは、すべてうまくいっています。
これは私のコードです:
library(mvtnorm)
library(Zelig)
data(approval)
set.seed(12345)
xh1 <- setx(approval, avg.price = mean(approval$avg.price) + sd(approval$avg.price), sept.oct.2001=1)
xl1 <- setx(approval, avg.price= mean(approval$avg.price), sept.oct.2001=1)
xh0 <- setx(approval, avg.price = mean(approval$avg.price) + sd(approval$avg.price), sept.oct.2001 =0)
xl0 <- setx(approval, avg.price= mean(approval$avg.price), sept.oct.2001=0)
zh1 <- sim(approval, x=xh1)zl1 <- sim(approval, x=xl1)zh0 <-sim(approval, x=xh0)
zl0 <- sim(approval, x=xl0)
eff <- (zh1$qi$ev - zl1$qi$ev) -(zh0$qi$ev - zl0$qi$ev)
quantile(eff, c(.025, .975))
そして、私が得るエラーは
行 87 ～ 103 (Preview-151c5fba73e.Rmd) から終了します terms.default(object) のエラー: 用語コンポーネントも属性呼び出しもありません: ... eval -> setx -> setx.default -> terms -> terms.default 実行が停止しました
参考のため、
m1<- zelig(approve~avg.price, model="ls", data=approval)
m2<- zelig(approve~avg.price+sept.oct.2001+iraq.war, model="ls", data=approval)
m3<- zelig(approve~avg.price+sept.oct.2001+avg.price:sept.oct.2001, model="ls", data=approval)
そして私のシミュレーションのために
x1 <- setx(m2, sept.oct.2001= 1)
s1 <- sim(m2, x=x1)
summary(s1)
x1 <- setx(m2, sept.oct.2001= 0)
s1 <- sim(m2, x=x1)
summary(s1)
oilprice<-min(approval$avg.price):max(approval$avg.price)
x2 <- setx(m2, sept.oct.2001=0, avg.price=oilprice)
s2 <-sim (m2, x=x2)
plot.ci(s2)
oilprice<-min(approval$avg.price):max(approval$avg.price)
x2 <- setx(m2, sept.oct.2001=1, avg.price=oilprice)
s2 <-sim (m2, x=x2)
plot.ci(s2)
これはすべて完全にうまくいきました。