search for: zelig

Displaying 20 results from an estimated 34 matches for "zelig".

Did you mean: zelis
2007 May 26
How to get the "Naive SE" of coefficients from the zelig output
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...
2009 Feb 19
Zelig method setx()
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...
2011 Jul 06
matching, treatment effect-ATT and Zelig package
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 +...
2011 Jan 25
Problem with matchit() and zelig()
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,, "control"), model="ls" , weight...
2009 Dec 07
zelig logit survey
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...
2005 Jan 06
package Zelig problem with setx
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...
2006 Apr 04
simulation with Zelig
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 =, weights = my.weights) Since the regression is weighted before running the regression I used the procedure: zelig2ls <- function(formula, model, data, M, ...)...
2009 May 04
Zelig, oprobit error
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",...
2010 Jan 10
lme4 and function 'cholmod_start' not provided by package 'Matrix' / Ubuntu
Hello all, Using Ubuntu 9.04 and R 2.8.1. For a project I need to use the Zelig package, which in turn wants to use the lme4 package. When trying to use Zelig and it tries to its required packages I get the following error message. Error in dyn.load(file, DLLpath = DLLpath, ...) : function 'cholmod_start' not provided by package 'Matrix' Error in loadModelD...
2009 Feb 12
gamma regression (Zelig package)
Hey, I want to estimate to regressions. First I want to estimate simple OLS. Since my dependent variable looks like a gamma distribution I want to estimate a gamme regression with the "zelig" package. But, I have some problems with the interpretation of the estimated coefficient. In the OLS the estimated coefficient is negative. In the gamma regression the estimated coefficient is positive. What is the interpretation of the estimated parameters within a gamma regression? Rega...
2011 Dec 16
Incorrect Number of Dimensions in Zelig with setx()
I'm running an ordered logit in R with the Zelig package and am trying to calculate some predicted probabilities. However, I get the following error message. > x.low <- setx(mod, cars=1)Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] : incorrect number of dimensions I googled this problem and couldn't find anyt...
2012 Aug 22
Clustered standard errors in Relogit (Zelig)
How do I apply clustered standard errors (by country) when using "relogit" in Zelig? I can not find a function covering this. [[alternative HTML version deleted]]
2011 Oct 18
getting basic descriptive stats off multiple imputation data
Hi, all, I'm running multiple imputation to handle missing data and I'm running into a problem. I can generate the MI data sets in both amelia and the mi package (they look fine), but I can't figure out how to get pooled results. The examples from the mi package, zelig, etc., all seem to go right to something like a regression, though all I want are the mean and SE for all the variables in the data set, and the cov and for among them. How would I request those basic descriptives (using mi, zelig, or any other package). Thanks! Mike
2020 Oct 05
Simultaneous Equation Model with Dichotomous Dependent Variables
...ooks like: glm(DV1~ IV1+ Control1+ Control2+ DV4+ DV5+ DV6, family="binomial", data= data) which I repeat for each of the 6 DVs. I realize that I need to run a series of simultaneous equation models, but I am not having much luck figuring out how to run an SEM with binary data in R. The zelig package used to have a bivariate logit command, but zelig does not work with the most current version of R. Additionally, I would run into issues later on using that command since 2 of my variables are not binary. Any ideas? Thanks for your help! Casey --- Ph.D student, Political Science Unive...
2005 Oct 26
Pb encountered with demos
...demo(lattice) ---- ~~~~~~~ Type <Return> to start : > require(grid) [1] TRUE > old.prompt <- grid.prompt(TRUE) > old.settings <- trellis.par.get() Erreur dans eval.with.vis(expr, envir, enclos) : impossible de trouver la fonction "trellis.par.get" - Zelig > demo(poisson,package='Zelig') demo(poisson) ---- ~~~~~~~ Type <Return> to start : > data(sanction) > user.prompt() Erreur dans eval.with.vis(expr, envir, enclos) : impossible de trouver la fonction "user.prompt" De plus : Warning message: data set...
2011 Dec 16
Zellig Error Message
I'm trying to calculate predicted probabilities in R with Zelig and keep getting the following error. Can anyone help? > x.low <- setx(mod, type=1)Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] : incorrect number of dimensions When I ran the model, I ran everything but the explanatory variable as a numeric variable. Now, I...
2011 Dec 22
Finding predicted probabilities
I ran three logit models in R with the Zelig package and I'm trying to compute the predicted probabilities for a number of different values on the independent variable. My dep variable was accepted or decline and my indep variable is bid amount, and varies. So for a bid amount of 3, what's the expected probability of winning. For a...
2009 Jul 26
Sweave, cacheSweave, and data frame
...the whoel thing every time I made changes. Then I discover the "cacheSweave" package, which seems the right solution to my problem. I only have on problem. Here is what I did: ---------------------------------------------- <<results=hide,echo=false>>= library(foreign) library(Zelig) library(memisc) options(digits=4) @ <<echo=false>>= d <- read.dta("~/project/abortion/data/data_transition_wide.dta") @ ... ---------------------------------------------- It can be compiled using " foo.Rnw", but when I tried "Swe...
2008 Jul 23
Time series reliability questions
...correct) results. R on the other hand, has the exogenous regressor with a negative coefficient. If I use other data I encounter the same problem - agreement between EViews and Gretl, disagreement with R (for identical data sets). Are there any known bugs with arima estimation in R? If I use the Zelig package, I get the same results as the arima{stats} function call. If I remove the exogenous regressor from the estimations then I have agreement between R, Gretl and EViews, but with the exogenous regressor (basically an ARMAX model) the estimation results are substantially different.
2009 Apr 22
Gee with nested desgin
...I'm wondering if this is correct. library(geepack) geeglm(formula = Y ~ Year, id = PlotID, waves = TreeID, family = binomial, corstr = "exchangeable") #using a unstructured correlation on the plot level. geeglm with unstructured correlation resulted in an out of memory error. library(Zelig) zelig(formula = Y ~ Year, model = "logit.gee", id = "PlotID", corstr = "unstructured") #Ideally I think I need a correlation matrix structured like below (given for a plot with 3 trees). Here a1 is the correlation within a tree, a2 the correlation between trees from...