[BioC] problem with pamr.listgenes in pamr package
Alice LE BARS
Alice.Le-Bars at univ-nantes.fr
Wed May 16 18:34:57 CEST 2007
Dear all,
I have trouble using the pamr.listgenes function in pamr package. When I
try to use it, I have the following error :
Error in dimnames(res) <- list(NULL, c("id", gnhdr, schdr)) :
length of 'dimnames' [2] not equal to array extent
I check if it was my data, but using the example for generating the data,
the problem remain.
Did someone have an idea?
I use a limited dataset of 12 samples and 14 genes
# With my originaly data
Val =
matrix(as.numeric(unlist(tableEntr[2:dim(tableEntr)[1],2:dim(tableEntr)[2]])),nrow=(dim(tableEntr)[1]-1))
LabelClasse = rep(c("CR","TOL"), each=6 )
ListeW = list(x=Val,y=LabelClasse)
entrainement = pamr.train(ListeW)
pamr.listgenes(entrainement,ListeW,1)
#with the example
set.seed(120)
x <- matrix(rnorm(1000*20),ncol=20)
y <- sample(c(1:4),size=20,replace=TRUE)
mydata <- list(x=x,y=factor(y))
train = pamr.train(mydata)
pamr.listgenes(train,mydata,1)
Thanks for help
Alice.
More information about the Bioconductor
mailing list