similar to: How to get the "Naive SE" of coefficients from the zelig output

Dear list, I am trying to do a repeated analysis using lme in R and a little bit unsure if I have set up the right statement. The problem is the IL6 (interleukin 6) was measured 5 times on each individual in each of 6 companies. The hypotheses are to see whether there is a relationship between IL6 and the total dust in each of the companies and if there is any change in IL6 across time

Dear R-Users: I am trying to find the robust (or sandwich) estimates of the standard error of fixed effects parameter estimates using the package "lmer2". In model-1, I used "robust=TRUE" on the other, in model-2, I used "robust=FALSE". Both models giving me the same estimates. So my question is, does the robust option works in lmer2 to get the robust estimates of

Hello R-Users: I am want to use tobit regression for left censored panel/longitudinal data. Could you please provide me the name of "library" and/or "package" that will give me option of fitting tobit regression model for longitudinal data? Thank you. Sattar __________________________________________________ [[alternative HTML version deleted]]

Hello Everybody: I am searching for an R package for fitting Generalized Estimating Equations (GEE) with weights (i.e. Weighted Estimating Equations). From the R documentation I found "geese(geepack)" for fitting Generalized Estimating Equations. In this documentation, under the paragraph “weights” it has been written, “an optional vector of weights to be used in the fitting process.

Hello EVERYONE, I need an URGENT help from you please! How can I see the "estfun" (empirical estimating function) and "df" (degree of freedom) from the following mixed-model please? (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) Many thanks in advance for your kind help. Sattar

Dear R-user: I am using survreg(Surv()) for fitting a Tobit model of left-censored longitudinal data. For logarithmic transformation of y data, I am trying use survreg.distributions in the following way: tfit=survreg(Surv(y, y>=-5, type="left")~x + cluster(id), dist="gaussian", data=y.data, scale=0, weights=w) my.gaussian<-survreg.distributions$gaussian

how can i put the row labels on the left hand of the heatmap in heatmap.2? abnr<-structure(c(1, 0.678622097406395, 0.670294749652918, -0.0016314464654279, -0.000519068106572792, 0.199581999119988, -0.0106623494189115, 0.0840111691399559, -0.0461494399639137, 0.249279171677728, NA, 1, 0.757114062773504, 0.0352642759270137, -0.0255518450373996, 0.0943268190664674, -0.0536269679247722,

Hello R Gurus, I have a matrix for which I am doing a heatmap using heatmap.2. I want to put the rownames on the lefthand side instead of the right side of the heatmap. how can i put the rownames on left hand side: I have already tried axis but could not make it work ccc<-structure(c(1, 0.283300333755851, 0.237863231117007, 0.0148696794159253, -0.0780756406815149, -0.106161465097659,

Dear R-user: I have a left censored longitudinally measured data set with 4 variables such as sub (which is id), x (only covariate), y (repeatedly measured response) and w (weights) (note, ?-5? indicates the left censored value in the attached data set). I am using following R codes (?survival? library and ?survreg? package) for fitting a random effect tobit model for the left censored

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]]

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

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 anything, minus a question by me on this same

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

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,

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)

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

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)

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"),

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

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 +