similar to: Fitting Weighted Estimating Equations

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

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

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 =

I need capabilities, for my data analysis, like the Pinheiro & Bates S-Plus/R package nlme() but with binomial family and logit link. I need multiple crossed, possibly interacting fixed effects (age cohort of twin when entered study, sex of twin, sampling method used to acquire twin pair, and twin zygosity), a couple of random effects other than the cluster variable, and the ability to

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

Hello R-users: I am using "yags" for fitting GEE which is giving me the same result as "Proc GENMOD". Now I have couple of questions related to yags output. (By the way, someone told me to run the geeglm for the same analysis and I did run but did not get the same result as of genmod and don't know how to correct the geeglm codes so that all three will be same!)

Hello, I have different results from these two softwares for a simple binomial GLM problem. >From Genmod in SAS: LogLikelihood=-4.75, coeff(intercept)=-3.59, coeff(x)=0.95 >From glm in R: LogLikelihood=-0.94, coeff(intercept)=-3.99, coeff(x)=1.36 Is there anyone tell me what I did wrong? Here are the code and results, 1) SAS Genmod: % r: # of failure % k: size of a risk set data

List: I'm helping a colleague with some Poisson regression modeling. He uses SAS proc GENMOD and I'm using glm() in R. Note on the SAS and R output below that our estimates, standard errors, and deviances are identical but what we get for likelihoods differs considerably. I'm assuming that these must differ just by some constant but it would be nice to have some confirmation

I have included data at the bottom of this email. It can be read in by highlighting the data and then using this command: dat <- read.table("clipboard", header = TRUE,sep="\t") I can obtain solutions with both of these: library(gee) fit.gee<-gee(score ~ chem + time, id=id, family=gaussian,corstr="exchangeable",data=dat) and library(yags) fit.yags <-

Hello, I have to run a general linear mixed model which looks at 2 dependent variables at the same time (var1 divided by var2). I have tryed to search for such a kind of model structure but since I just started using R my search was not successful. Especielly since I only have an old SAS GENMOD code structure from my project supervisor as an indication. My question is no, does there exist a code

Hello. I am attempting to duplicate a negative binomial regression in R. SAS uses generalized estimating equations for model fitting in the GENMOD procedure. proc genmod data=mydata (where=(gender='F')); by agegroup; class id gender type; model count = var1 var2 var3 /dist=NB link=log offset=lregtm; repeated subject=id /type=exch; run; Since my dataset has several observations for

Hello, Is it possible to derive Pan's QIC (2001 Biometrics 57:120) from either a fitted gee() object in the gee package or from a geese() fit in the geepack package? If so, would anyone be kind enough to provide me with code to do so? I realize that QIC is part of the output from yags() but I would like to use one of the other functions. Thanks. Richard

Dear all, I'm working with count data following over-dispersed poisson distribution and have to work with mixed-models on them (like proc GENMOD on SAS sys.). I'm still not to sure about what function to use. It seems to me that a glmmPQL will do the job I want, but I'll be glad if people who worked on this type of data can share what they learned. Thanks for your time. simon

Hi All, I'm using glm() in R to perform Poisson regression, I'm wondering if its possible to get equivalent Type 1 / Type 3 Analysis (similar to one in PROC GENMOD). Thanks, Kim [[alternative HTML version deleted]]

Hi,?In your data file, the first row does not have an equal number of column like the rest of the row.Check your data file. Specially 1st row.?? Regards.............Tanvir Ahamed Stockholm, Sweden | mashranga at yahoo.com On Monday, 5 October 2020, 08:11:48 am GMT+2, Mir Md. Abdus Salam <mir.salam at uef.fi> wrote: Dear all, I need urgent help. I am a new user of R. I

i'm trying to learn R. i have no programing experience. i tried to add on a package and just picked "yags". i can't get it to work or get any help. --- > install.packages("yags") Warning in install.packages("yags") : argument 'lib' is missing: using 'C:\Documents and Settings\j\My Documents/R/win-library/2.6' trying URL

We are using glm() to models to counts of deaths due to rare causes using a log link and Poisson error distribution, with population as the offset. Approximate confidence intervals for the parameter estimates are easy to calculate using a standard normal deviate, but obviously when the counts of deaths are small (which is why we are using Poisson regression), these intervals are very approximate

I have a dataset at a hospital level (as opposed to the patient level) that contains number of patients experiencing events (call this number y), and the number of patients eligible for such events (call this number n). I am trying to model logit(y/n) = XBeta. In SAS this can be done in PROC LOGISTIC or GENMOD with a model statement such as: model y/n = <predictors>;. Can this be done

Here is an example given from ?yags library(methods) data(stackloss) Y1 <- yags(stack.loss~Air.Flow,id=1:21, data=stackloss) How can I access parts of the output. I tried: > str(Y1) Formal class 'yagsResult' [package "yags"] with 25 slots ..@ coefficients : num [1:2] -44.13 1.02 ..@ coefnames : chr(0) > Y1$coefnames Error in Y1$coefnames : $