У меня была проблема при использовании кода pgmm в пакете plm.кажется, что он хорошо работает в чанках, но когда я вяжу его в файл pdf, он показывает мне сообщение об ошибке, показывающее
$Error in eval(predvars, data, env).
can not find the object "cigar_data63"
calls: <anonymous> ... model.frame.Formula ->
model.frame -> model.frame.default -> eval -> eval$
Данные, которые я использую, являются встроенным набором данных в пакете plm.это данные сигары.
---
title:
short:
journal: "AER" # AER, AEJ, PP, JEL
month: "`r format(Sys.Date(), '%m')`"
year: "`r format(Sys.Date(), '%Y')`"
author:
output: rticles::aea_article
---
library(plm)
data(Cigar)
cigar_data = log(Cigar)
cigar_data$state = Cigar$state
cigar_data$year = Cigar$year
cigar_data$lsales = lag(cigar_data$sales,1)
library(dummies)
cigar_data = cbind(cigar_data,dummy(cigar_data$year))
cigar_data = pdata.frame(cigar_data,index = c("state","year"))
q3_ols = plm(sales ~ lag(sales,1) + price + ndi + pimin,
data = cigar_data, effect = "time", model = "within")
q3_gmm = pgmm(sales ~ lag(sales,1) + price + ndi + pimin|
price + ndi + pimin + cigar_data63,
data = cigar_data, effect = "individual",
model = "twosteps")
summary(q3_ols)
summary(q3_gmm)
Результат должен дать оценку GMM.Заранее спасибо!