[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