Проблема с использованием BLPestimatoR с демографическими данными
Я работаю с R BLPestimatoR
пакет. Мне удалось получить оценку пакета без использования демографических данных, но при добавлении демографии я получаю
"Ошибка в оценке BLP(Xlin = Xlin_example, Xrandom = Xrandom_example, Xexo = Xexo_example,: невозможно сопоставить демографические данные с рынками, поскольку по крайней мере одна демографическая переменная не имеет информации cdid ".
у меня есть cdid
переменная в обоих BLP_demo
а также BLP_data
как "двойной". Данные BLP_demo расположены следующим образом:
d1 d2 d3 cdid
4 2 8 1
3 7 4 1
5 6 3 2
8 2 4 2
. . . .
. . . .
Код R: этот код читает оба BLP_data
а также BLP_demo
(проверено с помощью View) и работает нормально, без демо-связанного кода.
library("openxlsx")
library("BLPestimatoR")
Xlin_example <- c("P1", "CP1", "SS", "WG", "TT")
Xexo_example <- c("SS", "WG", "TT")
Xrandom_example <- c("P1", "CP1")
instruments_example <- c("PC", "CL", "PA", "UE", "AGr", "VGr", "ApGr", "SGr")
demoghaphics_example <- c("Income", "Age", "HhSize")
BLP_data <- read.xlsx("ProdData.xlsx", 1)
BLP_demo <- read.xlsx("DemoData.xlsx", 1)
BLP_est<- estimateBLP(Xlin = Xlin_example,
Xrandom = Xrandom_example,
Xexo = Xexo_example,
instruments = instruments_example,
demographics = demoghaphics_example,
shares = "shares",
cdid = "cdid",
productData = BLP_data,
demographicData = BLP_demo,
starting.guesses.theta2 = c(.001, .001),
solver.control = list(maxeval = 5000),
solver.method = "BFGS_matlab",
starting.guesses.delta = rep(1, length(BLP_data$cdid)),
blp.control = list(inner.tol = 1e-6,
inner.maxit = 5000),
integration.control= list( method="MLHS",
amountNodes= 100,
seed= 3 ),
postEstimation.control= list(standardError = "robust",
extremumCheck = TRUE,
elasticities = "PhonePrice"),
printLevel = 2)
summary(BLP_est)