search for: zelig

Hi, I'm having problems with the Zelig package - when using the below R displays the follwing message (I'm running R i386 2.15.3 for Windows and have updated all the Zelig packages): z.out<-zelig(Surv(psurv2, pcens2) ~ ren_sup3 + age, data=data_urgent, model="coxph") ** The model "coxph" is not availabl...

...about R. I have a user who is creating an R package for windows from a Linux environment using the crossbuild environment by Jun Yan and A.J. Rossini. The packages she generated worked fine until she tried to install in R 1.9.1 for Windows. Now when she installs with install.packages( "Zelig", CRAN="http://gking.harvard.edu" ) it results in two errors. First, it claims there are some missing files. When I look in the zip file, the files are there, but they coexist with other files that have the same name differing only in case. So there is both 'help/zelig' a...

Zelig: Everyone's Statistical Software Kosuke Imai, Gary King and Olivia Lau Version 1.0 (Available at http://gking.harvard.edu/zelig) A growing proportion of statisticians and methodologists from many disciplines are converging on R, a po...

hello all I am using the R package Zelig for some tobit regression with robust standard errors. I have got R version 2.9.2 (2009-08-24) and Zelig Version: 3.4-5 when i do demo(robust) It ends like this way data(coalition) > # Fit the model with robust standard error > user.prompt() Press <return> to continue: > z.out3...

Hi Folks, According to the PDF file blogit.pdf in the Zelig documentation: "Use the bivariate logistic regression model ["blogit"] if you have two binary dependent variables (Y1,Y2), and and wish to model them jointly as a function of some explanatory variables. Each pair of dependent variables (Yi1,Yi2) has four potential outc...

Dear R-user: After the fitting the Tobit model using zelig, if I use the following command then I can get the regression coefficents: beta=coefficients(il6.out) > beta (Intercept) apache 4.7826 0.9655 How may I extract the "Naive SE" from the following output please? > summary(il6w.out) Call: zelig(formula = il6.data$il6...

Dear R users, I am trying to estimate to estimate the average treatmen effect on the treated (ATT) using first the MatchIt software to weight the data set and, after this, the Zelig software as shown in Ho et al. (2007). See here for an explanation of how to apply this technique in R: http://imai.princeton.edu/research/files/matchit.pdf I encounter a slight problem when I apply the weights that are produced in the stage of preprocessing the data. The idea of this is to use t...

I'm running a multinomial logit in R using the Zelig packages. According to str(trade962a), my dependent variable is a factor with three levels. When I run the multinomial logit I get an error message. However, when I run 'model=logit' it works fine. any ideas on whats wrong? ## MULTINOMIAL LOGIT anes96two <- zelig(trade962a ~ age962 + ed...

Hello, I am attempting to "automate" a Bayesian normal linear regression using Zelig. Basically, I have a list containing several zelig() objects, each having a different formula, same data set, and same model (normal.bayes). My problem lies in the setx() method, where I am setting a numeric parameter to a value other than the mean. This is straightforward if I input the paramet...

I've run an ordered logistic regression model in R with Zelig and am looking to calculate predicted probabilities. Zelig has a series of simple one line commands to generate the information I want on first differences and so forth. Unfortunately, I keep getting an error when running the zelig function and was wondering if there was a quick alternative for gen...

Hi there, I'm wondering what Zelig in the following situation (code below) actually does. Is this considered as a so called regression adjustment after the propensity score matching? library(MatchIt) library(Zelig) data(lalonde) re78 represents the outcome variable 1. With Zelig m.out <- matchit(treat ~ age + educ + black +...

...below is ok, but if/else #1 is not? (I get an "unexpected else" type error.) In order for it to work I have to use if/else #2 Thanks .. maybe there is some reason for this, but this looks very inconsistent to me. R version 2.10.1 (2009-12-14) Ubuntu 9.04 Esmail --------------- #if (ZELIG) ######## if/else #1 # cat(sprintf("Zelig for %d runs - timeR\n", RUNS)) #else # cat(sprintf("lme for %d runs - timeR\n", RUNS)) if (ZELIG){ ####### if/else #2 cat(sprintf("Zeli...

Dear all, Does anybody know why the following code returns an error message? >library(MatchIt) >library(Zelig) >data(lalonde) > >m.out1<-matchit(treat~age+educ+black+hispan+nodegree+married +re74+re75, method="full", data=lalonde) > >z.out1<-zelig(re78~age+educ+black+hispan+nodegree+married+re74+re75, data=match.data(m.out1, "control"), model="ls" , weight...

...al admissions. Now that I have 5 completed datasets and I have to pool all the 5 datasets to get one combined output for a poisson model. I have checked previous queries about combining MI datasets and they are quite dated and none really related to my problem. I have been tirelessly working with zelig and unfortunately i am getting real stuck! here is my R code: ###these are my 5 imputed datasets#### d.1 <- read.table("imp1.csv", header=TRUE,sep=",") d.2 <- read.table("imp2.csv", header=TRUE,sep=",") d.3 <- read.table("imp3.csv", heade...

When I use zelig to compute parameter estimates with "model = logit.survey", I receive the following error: Nicholas Carnes. 2007. "logt.surveyWarning message: In eval(expr, envir, enclose) : non-integer #successes in a binomial glm! I believe this is because the model is not using "quasibinomi...

Dear list members, Is there a package somewhere for jointly estimating two poisson processes? I think the closest I've come is using the "SUR" option in the Zelig package (see below), but when I try the "poisson" option instead of the "SUR" optioin I get an error (error given below, and indeed, reading the documentation of the Zelig package, I get the impression "poisson" was not meant to handle a system of equations). I think...

Hi! Does somebody out there has experience with the Zelig package from Harvard uni? I have a problem when trying to set the explanatory variables with setx Polytomous logistic regression: >z.out <- zelig(OPARS ~ v1+v2+v3+...+vn, model = "mlogit", data=heb) that's OK >x.out<-set(z.out) Error in seq.Date(along = object) : `from...

Dear R-friends, I am trying to use the Zelig library to simulate the following model: z.out<-zelig(formula = y ~ x1*x2 + x3 + x4 + x5, model="ls", data = my.data, weights = my.weights) Since the regression is weighted before running the regression I used the procedure: zelig2ls <- function(formula, model, data, M, ...)...

Hello, I'm getting an error message when I use the ordered probit model "oprobit" in the zelig function. Using the same form as in the help file, we get an error message. It produces coefficients, but no standard errors. See results below. Any hints? Thanks! > > o.probit <- zelig(as.factor(checks.change) ~ Oda + Oil + sh_neg + sh_pos + checks1, model="oprobit",...

I'm running a multinomial logit in R using the Zelig package. However I get the following error. HELP anes96two <- zelig(trade962a ~ age962 + education962 + personal962 + economy962 + partisan962 + employment962 + union962 + home962 + market962 + race962 + income962, model="mlogit", data=data96) #Error in attr(tt, "depFactors&qu...