Hi, I'm trying to fit a marginal model via GEE but I'm getting strange results and few problems. If I set the working correlation as exchangeable I'm getting the same fitting when I set as independent. Comparing to SAS results it shouldn't happen. If I try to use another working correlation (like AR-M or stat_M_dep), R just exits without giving any error message. Another doubt is how R estimates the scale parameter for Poisson? I get 1.36 in R, 1.25 in SAS (estimation by the square root of DEVIANCE/DOF) and 1.17 (estimation by the square root of Pearson's Chi-Square/DOF). I'm using R 1.4.1 but all the problems apply for S-Plus 2000 too. I hope anybody can help me. Sincerely, -- Frederico Zanqueta Poleto fred at poleto.com - ICQ# 4129787 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Dear all,
I want to chage the scale oder when I use the plot (x , y
,xlim=c(xmin,xmax), ylim=c(ymin,ymax) ) the default scale order was given like
x axis labels are 60,70,80,90,100 ,if I could , I hold axis labels are
60,65,70,75,80,85,90,95,100.
Thanks
Ken
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
ken_lee wrote:> > Dear all, > I want to chage the scale oder when I use the plot (x , y ,xlim=c(xmin,xmax), ylim=c(ymin,ymax) ) the default scale order was given like > x axis labels are 60,70,80,90,100 ,if I could , I hold axis labels are 60,65,70,75,80,85,90,95,100.Use axis(). Example: plot(60:100, 60:100, xaxt="n") axis(1, at=seq(60, 100, 5)) Uwe Ligges -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Dear all,
I just knew the prompt "#" can comment a line,but tow to comment
a paragraph more than one line?
Thanks
Ken
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
On Sun, 12 May 2002, Frederico Zanqueta Poleto wrote:> Hi, > > I'm trying to fit a marginal model via GEE but I'm getting strange > results and few problems. > If I set the working correlation as exchangeable I'm getting the same > fitting when I set as independent. Comparing to SAS results it shouldn't > happen. > If I try to use another working correlation (like AR-M or stat_M_dep), R > just exits without giving any error message.We probably need an example. The `OME' example in help(gee) gives different results using exchangeable and independence working correlations, so it isn't a universal problem.> Another doubt is how R estimates the scale parameter for Poisson? I get > 1.36 in R, 1.25 in SAS (estimation by the square root of DEVIANCE/DOF) > and 1.17 (estimation by the square root of Pearson's Chi-Square/DOF). > I'm using R 1.4.1 but all the problems apply for S-Plus 2000 too.This is because 1.36=1.17^2, I should think. The disperson parameter from gee() is a variance scale factor not a standard deviation scale factor. -thomas Thomas Lumley Asst. Professor, Biostatistics tlumley at u.washington.edu University of Washington, Seattle -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._