R拟合面板logit模型时报错

panel是数据框，其数据形式如下：

1. > dput(head(panel,n = 2))
   
2. structure(list(time = 10:11, code = c("000002.SZ", "000002.SZ"
   
3. ), Y = c(0, 0), jgcg = c(40.0698, 43.7155), mv = c(90380627991.96,
   
4. 82134220328.46), eps = c(0.66, 0.88), qr = c(0.5568, 0.3702),
   
5.     TAI = c(56.7036, 37.364), bps = c(4.02, 4.82), a2d = c(74.6861,
   
6.     77.0997), FAT = c(39.3809, 50.9921), TAT = c(0.2871, 0.2805
   
7.     ), SAT = c(0.2692, 0.253)), .Names = c("time", "code", "Y",
   
8. "jgcg", "mv", "eps", "qr", "TAI", "bps", "a2d", "FAT", "TAT",
   
9. "SAT"), row.names = 1:2, class = "data.frame")

复制代码

1. > head(panel)
   
2.   time      code Y    jgcg           mv  eps     qr     TAI  bps     a2d     FAT    TAT    SAT
   
3. 1   10 000002.SZ 0 40.0698  90380627992 0.66 0.5568 56.7036 4.02 74.6861 39.3809 0.2871 0.2692
   
4. 2   11 000002.SZ 0 43.7155  82134220328 0.88 0.3702 37.3640 4.82 77.0997 50.9921 0.2805 0.2530
   
5. 3   12 000002.SZ 0 49.9849 111274997554 1.14 0.4141 27.8835 5.80 78.3163 64.2845 0.3055 0.2823
   
6. 4   13 000002.SZ 0 38.3900  88450200420 1.37 0.3372 26.5056 6.98 77.9970 72.3773 0.3157 0.3166
   
7. 5   14 000002.SZ 0 42.6325 153421350664 1.43 0.4255  6.0941 7.99 77.2046 65.9685 0.2964 0.3161
   
8. 6   15 000002.SZ 0 53.0831 269990888489 1.64 0.4259 20.2370 9.08 77.7015 54.1250 0.3493 0.4029

复制代码

以下是我的代码：

1. form <- formula(panel$Y~panel$jgcg+panel$mv+panel$eps+panel$qr+panel$TAI+panel$bps
   
2.                 +panel$a2d+panel$FAT+panel$TAT+panel$SAT)
   
3. pool <- pglm(form,data = panel,model = 'pooling',family = binomial())
   
4. rand <- pglm(form,data = panel,model = 'random',family = binomial())
   
5. fixed <- pglm(form,data = panel,model = 'within',family = binomial(link = 'logit'))

复制代码

混合效应和随机效应都没有报错，固定效应出现如下报错：

1. > fixed <- pglm(form,data = panel,model = 'within',family = binomial(link = 'logit'))
   
2. Error in maxRoutine(fn = logLik, grad = grad, hess = hess, start = start,  :
   
3.   argument "start" is missing, with no default

复制代码

不明白这里的start是怎么回事，要如何修改代码来正确拟合？
谢谢咯~

自己顶
求大神解答

唤月 发表于 2017-8-17 08:26
自己顶
求大神解答

请问这个问题解决了吗