search for: genmod

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 bin_data; inpu...

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

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 in R which is equivalent to the SAS code below? PROC GENMOD DATA=X; CLASS FLH; MODEL BS/OCCUPANCY = distcrop distfor flh distcrop*flh /D=B LINK=LOGIT TYPE3; RUN; Thanks for answers! Nick [[...

...e repeated() library by Lindsey seems to handle just one nested random effect, or so I believe I read while scanning backlogs of the R-Help list. glmmPQL() is in the ballpark of what I need, but once again seems to lack the "outer" variable specification that nlme() has, and which PROC GENMOD also has---and which I need. I read someplace of yags() that apparently uses GEE to estimate parameters of nonlinear models including GLIMs/mixed models, just the way PROC GENMOD (and many another program) does. But on trying to install it (either v4.0-1.zip or v4.0-2.tar.gz from Carey's...

...p/chenwj/install/powerpc64-unknown-linux-gnu/bin/ -g -O2 -DIN_GCC -W -Wall -Wwrite-strings -Wstrict-prototypes -Wmissing-prototypes -pedantic -Wno-long-long -Wno-variadic-macros -Wno-overlength-strings -Wold-style-definition -Wmissing-format-attribute -DHAVE_CONFIG_H -DGENERATOR_FILE -o build/genmodes \ build/genmodes.o build/errors.o .././libiberty/libiberty.a build/genmodes -h > tmp-modes.h /bin/sh: line 1: 9623 Segmentation fault build/genmodes -h > tmp-modes.h make[3]: *** [s-modes-h] Error 139 make[3]: Leaving directory `/tmp/chenwj/build/gcc' make[2]: *** [all...

Hi, folks I am trying to build LLVM-GCC 4.2 on a Linux/PowerPC machine. Actually, the PowerPC is a PS3. The source was downloaded from http://llvm.org/releases/2.8/llvm-gcc-4.2-2.8.source.tgz Here is my system information and configuration options: ---------------------------------------------------------------------- $ uname -a Linux ps3 2.6.32-rc2-00995-g96ebbe6-dirty #2 SMP Fri Oct 2

...g 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. The length of weights should be the same as the number of observations. This weights is not (yet) the weight as in sas proc genmod, and hence is not recommended to use.” Now my question is, is there any other package you know that might allow me to fit GEE with weight option and the results would be same as SAS PROC GENMOD with weights? Any help would be sincerely appreciated. Thank you very much. Sincerely, Abdus Sat...

...egypt-linux/build/llvm-initial/src/include -DENABLE_LLVM -I/tmp/ptmp/dag/buildbot/slaves/llvm-validate/llvm-validate-release-egypt-linux/build/llvm-initial/obj/../src/include -DNDEBUG -D_GNU_SOURCE -D__STDC_LIMIT_MACROS -D__STDC_CONSTANT_MACROS -o build/errors.o ../../src/gcc/errors.c build/genmodes -h > tmp-modes.h /bin/sh: line 1: build/genmodes: No such file or directory make[4]: *** [s-modes-h] Error 127 Google turned up a gcc bug in which a maintainer explained that doing a ./configure in the source directory could cause this. Except the buildbot doesn't do that. Anyone seen...

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 each subject, I need the REPEATED statement in or...

...itial/src/include > -DENABLE_LLVM > -I/tmp/ptmp/dag/buildbot/slaves/llvm-validate/llvm-validate-release-egypt-linux/build/llvm-initial/obj/../src/include > -DNDEBUG -D_GNU_SOURCE -D__STDC_LIMIT_MACROS -D__STDC_CONSTANT_MACROS -o > build/errors.o ../../src/gcc/errors.c > build/genmodes -h > tmp-modes.h > /bin/sh: line 1: build/genmodes: No such file or directory > make[4]: *** [s-modes-h] Error 127 > > Google turned up a gcc bug in which a maintainer explained that doing > a ./configure in the source directory could cause this. Except the buildbot > does...

Greetings. I'm trying to recreate in R some regression models I've done in SAS, but I'm not getting the same results. My advisor suspects this may be due to differences in precision between R and SAS. Does anyone know where I can find specifications for R's type double? (It doesn't seem to be in the R Language Definition.) Thanks in advance for any help anyone can

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 [[alternative HTML version deleted]]

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

...ility... yes gcc -static doesn't work? libtool doesn't support shared libraries? complex.h? Wow. I did find this: http://www.mail-archive.com/autoconf at gnu.org/msg09268.html , but I don't know what to do about it. Eventually the script crashes, complaining about the absence of build/genmodes .h . >echo timestamp > s-options >gawk -f ../../../source/gcc-4.5.1/gcc/opt-functions.awk -f ../../../source/gcc-4.5.1/gcc/opth-gen.awk \ > < optionlist > tmp-options.h >/bin/sh ../../../source/gcc-4.5.1/gcc/../move-if-change tmp-options.h options.h >echo timestamp...

...n), these intervals are very approximate indeed. Has anyone done any work on calculating more precise likelihood ratio - based confidence intervals for parameter estimates from generalised linear models? These are also known as "profile likelihood confidence intervals", I believe. PROC GENMOD in SAS can calculate them - I am happy to email the relevant page from the SAS online documentation which describes the way in which they are calculated in more detail to anyone who is interested - but we would much rather use R for this work... Regards, Tim Churches Epidemiology and Surveillance...

...;***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 (Dispersion parameter for Gamma family taken to be 0.2745989) Null deviance: 63.635 on 161 degrees of freedom Residual deviance: 43.214 on 153 degrees of freedom AIC: 1527.6 Results From SAS *proc* *genmod* data=a order=data; class grass altitude; model snail = grass gheight humidity altitude soiltemr airtemr / dist=gamma type3; *run*; Analysis Of Parameter Estimates Standard Wald 95% Confi...

...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 using lrm from the Hmisc library without restructuring the dataset so that for each hospital there is one row with y = 1 and one row with y = 0 and then using the weight option in lrm to weight these two responses by...

    Hello Everyone,   I'm trying to use SAS to replicate some results obtained in R. I was wondering if anyone call tell me the SAS equivalent of the code that appears below.   fm.glm.x <- glm(resp ~ . - 1, data = as.data.frame(mm.x),  na.action = na.exclude, family = binomial(link = "probit")) summary(fm.glm.x) Thanks,   Paul [[alternative HTML version deleted]]

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!) Questions: 1. How...

Hi, we try to do a logistic regression with the function glm. But we notice that this function don't give the same results as the SAS proc catmod (differents estimate given). We try to change the contrast on R system with: > options(contrasts=c(unordered="contr.SAS",ordered="contr.poly")) We also try with brlr and logistf functions. Unfortunately, the estimate