Displaying 20 results from an estimated 900 matches similar to: "Zelig and the "blogit" model"
2010 Apr 11
1
Matrix is not symmetric under lme4
Dear all,
My code is presented as the following.
library(MASS)
library(rmutil)
library(repeated)
library(lme4)
library(arm)
#install.packages("Zelig", repos = "http://gking.harvard.edu")
library(Zelig)
rm(list = ls())
beta0<-2.5
beta1<--0.3
sigs2<-0.5
I<-4
#numberpatients<-c(40,100,160,200,400,600)
numberpatients<-c(1000)
#numberpatients<-3
times<-1
2013 Mar 05
2
Zelig package: Coxph model problems
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 available with the currently
loaded packages,
** and is not an
2009 Oct 30
1
Package zelig
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 <- zelig(Surv(duration, ciep12) ~ polar + numst2 +
2003 Aug 06
1
demo() function returns alphabetical matches, not unique matches?
Hi,
I'm working on a set of demo files for a package, and I'm having a problem because I have two demo files (one called "logit", the other called "blogit") and when I type demo(logit), it gives me demo(blogit) instead. (And calls it the the demo for logit.)
I've figured out that this is because the demo() function searches through all the demo directories (in
2010 Mar 29
1
Question about 'logit' and 'mlogit' in Zelig
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 + education962 + personal962 + economy962 +
2007 May 26
1
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.data$il6 ~ il6.data$apache, model =
2010 Dec 21
0
"variable lengths differ (found for '(weights)')" error in Zelig library
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
2009 Feb 19
0
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 parameters
2011 Dec 16
0
Error constructing probabilities in Zelig
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 generating
predicted probabilities for
2011 Jul 06
0
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 + hispan + married + nodegree + re74 + re75, data = lalonde)
2011 Jan 25
0
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,
data=match.data(m.out1, "control"),
2012 Oct 26
0
combined output with zelig is not working!?!
Hi everyone,
I have carried out a multiple imputation in R using Amelia II and have
created 5 multiply imputed datasets. The purpose of my research is to fit a
Poisson Model to the data to estimate numbers of hospital 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
2009 Dec 07
0
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 "quasibinomial" (see R-News
2003, Analyzing Survey Data in R, by Thomas Lumley). Is there
2004 Jun 25
2
R 1.9.1 package installation problems
Hello,
I am writing as an administrator, not as an R user, so forgive me if I
am not completely knowledgeable 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
2010 Feb 15
1
argh .. if/else .. why?
Hello, would someone please help explain the following inconsistency
in the if/else statement to me?
The format of the if/else #3 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
2005 Jan 06
0
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' must be specified
#??? I have no date
2009 May 04
0
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,
2006 Apr 04
0
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 =
my.data, weights = my.weights)
Since the regression is weighted before running the regression I used
the procedure:
zelig2ls <- function(formula, model, data, M, ...) {
mf <- match.call(expand.dots = TRUE)
2010 Mar 24
0
Zelig: Error message for 'mlogit'
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")$depFactorVar :
# $ operator is invalid
2009 Apr 13
2
joint estimation of two poisson equations
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