Chris,
Do summary(dat) and names(dat) -- is it as you expect? Do you have
repeating rows? Eliminate them using
dat1 <- unique(dat)
Then use dat1 (rather than dat) in the zelig() step.
Note, this is not a problem with zelig, this is a problem with your data.
Best,
Olivia
On Mon, Oct 22, 2012 at 8:41 PM, Chris Compton <CCompton(a)anexa.co.nz> wrote:
Dear Olivia and Kosuke,****
** **
Thank you for your prompt response. Unfortunately, the problem still
remained with the zelig functions after I made the suggested changes. Any
other suggestions or steps to diagnose the problem would be greatly
appreciated.****
** **
Thanks,****
** **
Chris****
** **
*From:* zelig-bounces(a)lists.gking.harvard.edu [mailto:
zelig-bounces(a)lists.gking.harvard.edu] *On Behalf Of *Olivia Lau
*Sent:* Saturday, 20 October 2012 11:08 a.m.
*To:* Kosuke Imai
*Cc:* Chris Compton; zelig(a)lists.gking.harvard.edu
*Subject:* Re: [zelig] "Error in data.frame" message****
** **
Hi, ****
** **
Try:****
** **
rownames(dat) <- NULL****
** **
then re-run the zelig step. ****
** **
Best, ****
Olivia****
On Fri, Oct 19, 2012 at 3:03 PM, Kosuke Imai <kimai(a)princeton.edu> wrote:*
***
Not sure what the problem is but I wonder if the variable name should not
contain underscores.
Kosuke
Department of Politics
Princeton University
http://imai.princeton.edu
On Oct 17, 2012, at 10:33 PM, Chris Compton <CCompton(a)anexa.co.nz> wrote:
Dear List Members,
I am running R 1.15.1 64 bit on Windows 7 64bit, and
'Zelig' version 4.0-11
'ZeligGEE' version 0.2
I wish to use Zelig to estimate mean effects from logistic regression
models with
confidence intervals, and differences between mean effects with
confidence intervals. Both standard and GEE logistic regression models are
operating in R packages base and geepack, respectively, but in Zelig are
returning the below messages:
z.out.1 <- zelig(cons1~rxgroup_id+herd_id,
model = "logit", data = dat)
Error in data.frame(visible = rep.int(FALSE,
n2), from = rep.int(msg, :
duplicate row.names: nobs.lm
library("ZeligGEE")
Loading
required package: gee
> z.out.1 <- zelig(cons1~rxgroup_id, model = "logit.gee", id =
"herd_id",
data = dat, corstr="exchangeable")
Beginning Cgee S-function, @(#) geeformula.q 4.13
98/01/27
running glm to get initial regression estimate
(Intercept) rxgroup_idGPG + P4 rxgroup_idGPG + P4
CoSynch
-0.08068891
-0.03956022
0.36057671
Error in data.frame(visible = rep.int(FALSE, n2),
from = rep.int(msg, :
duplicate row.names: nobs.default
I have run similar models in the past on previous versions of Zelig.
Look as I might, I can't find a problem with the dataframe. There are no
duplicate subject identifiers. I have tried using the subject identifier in
the first column in ascending order also. Any help on this problem would be
much appreciated.
Chris Compton BVSc MVS (Epi) (Hons)
Research Project Manager
Cognosco
Research division of Anexa Animal Health
PO Box 21 Morrinsville, New Zealand 3340
Ph +64 07 889 5159
Fax +64 07 889 3681
M +64 021 800 146
-
--
Zelig Mailing List, served by HUIT
Send messages: zelig(a)lists.gking.harvard.edu
[un]subscribe Options:
http://lists.gking.harvard.edu/mailman/listinfo/zelig
Zelig program information:
http://gking.harvard.edu/zelig/
Zelig mailing list
Zelig(a)lists.gking.harvard.edu
To unsubscribe from this list or get other information:
https://lists.gking.harvard.edu/mailman/listinfo/zelig
-
--
Zelig Mailing List, served by HUIT
Send messages: zelig(a)lists.gking.harvard.edu
[un]subscribe Options:
http://lists.gking.harvard.edu/mailman/listinfo/zelig
Zelig program information:
http://gking.harvard.edu/zelig/
Zelig mailing list
Zelig(a)lists.gking.harvard.edu
To unsubscribe from this list or get other information:
https://lists.gking.harvard.edu/mailman/listinfo/zelig****
** **