mlogit Wald检验，R

Question

我已经使用mlogit软件包估算了R中的多项式logit回归（mlog）。代码

coef(mlog)

给出值，

-1:(intercept)  1:(intercept)     -1:Coefficient1      1:Coefficient2 
   -0.09696444     0.00569650          1.17099302          1.04593858 

因此，该模型具有两个截距和两个系数。现在，我想使用Wald检验来检验两个系数同时为1和两个截距为0的零假设。我已经尝试了以下代码，

vec0col <- matrix(rbind(0,0,1,1),nrow=4)
Rmat <- cbind(vec0col,diag(4))
wald.test(b=coef(mlog), Sigma=vcov(mlog), L=Rmat)

但是我收到错误消息，

Error in dimnames(L) <- list(paste("L", as.character(seq(NROW(L))), sep = ""),  : length of 'dimnames' [2] not equal to array extent

我的问题是：我该如何解决这个问题或如何以更合适的方式进行这项测试？

谢谢！