为美好的一天开始解决错误

Question

让我们考虑以下数据

library(plm)
library(pglm)
data("EmplUK", package="plm")

我将添加随机放置 0 和 1 的新列。之后我想执行logit随机效应模型。

df1<-EmplUK
#adding 0's and 1's
df1<-cbind(df1,'binary'=sample(0:1,1031,replace=T))
#Performing logit regression 
pglm(binary~output+wage, data=df1, family=quasibinomial(link='logit'), start = NULL, model = 'random')

并且出现以下问题：

Error in maxRoutine(fn = logLik, grad = grad, hess = hess, start = start,  : 
  argument "start" is missing, with no default

我不确定到底是什么原因，我读过这个错误，似乎在你尝试估计 'within' 模型时出现了一些问题，但我对每个模型都遇到了这个错误类型。你能帮我指出这个错误的原因吗？

Answer 1

我认为这个函数中没有设置拟二项式族。在 pglm 中有一个函数 pglm:::starting.values 可以查找特定的家族：

"binomial"
"ordinal"
"poisson"
"negbin"
"gaussian"
"tobit"

负二项式允许对方差进行建模，以便可以满足您的需要，否则 binomial(link='logit') 如果没有过度分散的证据，则可以正常工作。

编辑：很高兴得到纠正，我以前没有用过这个包:)

为美好的一天开始解决错误

Solving error for the good beginning of the day

regression

r

panel