使用设置:
id <- c(1,1,2,2,3,4)
t <- c(1,2,1,2,1,1)
x <- c(1,2,2,1,2,1)
y <- c(1,0,0,0,1,0)
df <- data.frame(id, t, x, y)
从OLS回归中准备拟合值作为起始值:
tstart <- lm(y ~ x, data = df)
运行probit:
tfit <- pglm(y ~ x, data = df, index = c("id","t"),
model = "within", family = binomial('probit'), start = tstart$fitted.values)
返回错误
Error in lnl.binomial(param = start, y = y, X = X, id = id, model = model, :
object 'Li' not found
这个错误对我来说似乎非常缺乏信息。 我甚至没有在任何调用中寻找某个对象'Li',并且不知道该对象应该是什么。 回溯使它似乎出现在函数中:
9: lnl.binomial(param = start, y = y, X = X, id = id, model = model,
link = link, rn = rn) at <text>#1
但是试图查看错误发生时函数的代码会发现甚至没有函数
lnl.binomial()
我在哪里出错了?