[R] SAMR for paired samples

Haesook Kim kim at jimmy.harvard.edu
Tue Aug 24 18:16:57 CEST 2010


Hi R-help,

I am trying to use 'samr' for 10 pre and post paired samples to test 
whether
post is different from pre (i.e., the location shift for the delta of 
(post-pre)).
However, I got an error message saying

 > samr.obj<-samr(d, resp.type="Two class paired", nperms=100, 
random.seed=100)
perm= 1
Error in !logged2 : invalid argument type

Does anyone  know what this means? or how to  solve this error?

Also, the options 'One Sample paired' is for testing the differences 
directly (d1, d2,...), I suppose.

The source code is
*****************************************
function()
{
  library(samr)
  attach(proto_data)
  x <- cbind(pre1, pre2, pre3, pre4, pre5, ...
        post1, post2,post3,post4,post5, ....)
  y <- c(-(1:10),1:10)
 d <- list(x=x,y=y, geneid=array.id)
  samr.obj<-samr(d, resp.type="Two class paired", nperms=100)
  samr.plot(samr.obj)
}
*****************************************


Thank you

HS



More information about the R-help mailing list