Dear Sayan,
I've been unable to replicate the error you have below. When I do
demo(robust), I get:
summary(z.out3)
Call: zelig(formula =
Surv(duration, ciep12) ~ polar + numst2 + crisis,
model = "weibull", data = coalition, cluster = "invest",
robust = TRUE)
Value Std. Err (Naive SE) z p
(Intercept) 3.26773 0.099103 0.10924 32.97 1.98e-238
polar -0.04197 0.000937 0.00431 -44.78 0.00e+00
numst2 0.40133 0.046546 0.11095 8.62 6.57e-18
crisis 0.00482 0.000350 0.00188 13.78 3.26e-43
Log(scale) -0.12891 0.032200 0.04997 -4.00 6.24e-05
Scale= 0.879
Weibull distribution
Loglik(model)= -1052 Loglik(intercept only)= -1101
Chisq= 97.95 on 3 degrees of freedom, p= 0
(Loglikelihood assumes independent observations)
Number of Newton-Raphson Iterations: 5
n= 314
Can you please make sure that you are running the latest version of Zelig on
R 2.9.2?
Yours,
Olivia
On Tue, Nov 3, 2009 at 5:54 AM, sayan dasgupta <kittudg(a)gmail.com> wrote:
Dear Zelig developers and users,
I want to find robust standard errors after fitting a model
So I checked out package Zelig
Here are the following commands
R >library(Zelig)
data(coalition)
z.out3 <- zelig(Surv(duration, ciep12) ~ polar + numst2 + crisis, model
=
"weibull", data = coalition, cluster = "invest", robust = TRUE)
It gives the following error
Error in rowsum.default(resid(fit, "dfbeta"), cluster) :
'x' must be numeric
Please Help
Thanking You
Sayan Dasgupta