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
______________________________________________
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 provide commented, minimal, self-contained, reproducible code.
