Re: [R] SVYPLOT

2014-11-20 Thread Raphael Fraser
at<-svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, > fpc=~fpc) > svyplot(api00~api99, design=dstrat, style="bubble") > svyplot(api00~api99, design=dstrat, style="bubble", xlim=c(500, 700)) > > > On Thu, Nov 20, 2014 at 12:54 AM, Raphae

[R] SVYPLOT

2014-11-19 Thread Raphael Fraser
How do I set the limits of my x and y axis in svyplot? xlim and ylim does not work. Regards, Raphael [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posti

Re: [R] Re : Elegant Code

2012-03-16 Thread Raphael Fraser
gt; > On Fri, 16 Mar 2012 10:02:48 -0400 Raphael Fraser > wrote: > >> Thank you all. Your responses were very helpful. >> >> Raphael >> >> On Fri, Mar 16, 2012 at 5:04 AM, Daniel Nordlund >> wrote: >> > >> >> -Original Message-

Re: [R] Re : Elegant Code

2012-03-16 Thread Raphael Fraser
ay, March 16, 2012 12:36 AM >> To: Raphael Fraser >> Cc: r-help@r-project.org >> Subject: [R] Re : Elegant Code >> >> Hi Raphael, >> >> Something like that? >> >> require(pscl)sim <- numeric() >> for(i in 1:5){ >>   eval(parse(text

[R] Elegant Code

2012-03-16 Thread Raphael Fraser
Hi, Can anyone help to write a more elegant version of my code? I am sure this can be put into a loop but I am having trouble creating the objects b1,b2,b3,...,etc. b1 <- rigamma(50,1,1) theta1 <- rgamma(50,0.5,(1/b1)) sim1 <- rpois(50,theta1) b2 <- rigamma(50,1,1) theta2 <- rgamma(50,0.5,(1/b2)

Re: [R] tukey.1

2010-11-09 Thread Raphael Fraser
x27;s in it? > What's aov.pen? > > On Tue, Nov 9, 2010 at 10:13 AM, Raphael Fraser > wrote: >> I have been trying to do tukey's test to no avail. This is what I have >> and the error produced: >> >> pen.df = data.frame(blend, treatment, y) >> source

[R] tukey.1

2010-11-09 Thread Raphael Fraser
I have been trying to do tukey's test to no avail. This is what I have and the error produced: pen.df = data.frame(blend, treatment, y) source("tukey.1.r") tukey.1(aov.pen, pen.df) Error in tukey.1(aov.pen, pen.df) : the model must be two-way This is a two-way design. Therefore I am confused. Can

[R] Tinn R

2010-10-02 Thread Raphael Fraser
Dear ALL, I am trying to run Tinn R for the first time with R. I have already selected the SDI option. When ever I run the code I get the following error: > source(.trPaths[4]) Error in source(.trPaths[4]) : object '.trPaths' not found Can any one help? Raphael ___

Re: [R] Rounding Elements of a Matrix

2010-09-20 Thread Raphael Fraser
a MD > > Work: 301-451-8575 > Mobile: 301-204-5642 > Email: zoppo...@mail.nih.gov > ________ > From: Raphael Fraser [raphael.fra...@gmail.com] > Sent: Monday, September 20, 2010 12:21 PM > To: r-help@r-project.org > Subject: [R] Rounding

[R] Rounding Elements of a Matrix

2010-09-20 Thread Raphael Fraser
I have a 2x50 matrix and would like to round all the elements to 2 decimal places. Can any one help? Thanks, Raphael __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.

Re: [R] Mantel-Byar test

2010-05-11 Thread Raphael Fraser
Thank you Terry. On Tue, May 11, 2010 at 8:28 AM, Terry Therneau wrote: > >> The Mantel-Byar test is a simple modification to the logrank test. Can >> R perform this test or any other similar test? > > The logrank test = "score test" from coxph > The Mantel-Byar is the score test from a Cox mod

[R] Comparing Survival: Mantel-Byar Test

2010-05-10 Thread Raphael Fraser
The Mantel-Byar test is a simple modification to the logrank test. Can R perform this test or any other similar test? Raphael __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] Skew-Normal Linear Mixed Model (SNLMM)

2009-07-28 Thread Raphael Fraser
Dear ALL, In the Normal Linear Mixed Model (NLMM) we make the assumption that the random effects are normally distributed. Is it possible to fit a linear mixed model in R where the random effects are assumed to have a Skew-Normal distribution? I am trying to reproduce the results in the paper refe

[R] Skew-Normal Linear Mixed Model (SNLMM)

2009-07-28 Thread Raphael Fraser
Dear ALL, In the Normal Linear Mixed Model (NLMM) we make the assumption that the random effects are normally distributed. Is it possible to fit a linear mixed model in R where the random effects are assumed to have a Skew-Normal distribution? I am trying to reproduce the results in the paper refe

[R] 3D Graphs

2008-09-02 Thread Raphael Fraser
Can R plot 3D graphs? If so, how would you plot x^2 + y^2 + z^2 = 1? Thanks Raphael __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provid