[BioC] maSigPro Problems in see.genes

Marcelo Laia marcelolaia at gmail.com
Thu Feb 21 16:41:51 CET 2008


Dear,

We have a Single Series Time Course and we would like to analyse it on maSigPro.

But we got an error in see.genes:

> see.genes(sigs, design, main = "Exp Prof MI", show.fit = T,
+ dis = design$dis, cluster.method="kmeans", cluster.data = 1, k = 9)
Error in edesign[, time.col] : incorrect number of dimensions

sig.genes function is OK.

> sigs <- get.siggenes(tstep, rsq = 0.6, vars = "each")
> sigs$summary
    if (stringsAsFactors) factor(x) else x
1                            AD1X1FD8 (1)
2                            AD1X1FD8 (2)
(...)
136                          ST3X9BD3 (1)
137                          ST3X7DD0 (2)

my design file (degree = 4):

> design
$dis
        Time Time2 Time3 Time4
Array1     1     1     1     1
Array2     1     1     1     1
Array3     1     1     1     1
Array4     5    25   125   625
Array5     5    25   125   625
Array6     5    25   125   625
Array7     9    81   729  6561
Array8     9    81   729  6561
Array9     9    81   729  6561
Array10   13   169  2197 28561
Array11   13   169  2197 28561
Array12   13   169  2197 28561
Array13   17   289  4913 83521
Array14   17   289  4913 83521
Array15   17   289  4913 83521

$groups.vector
[1] "Group" "Group" "Group" "Group"

$edesign
        Time Replicates Group
Array1     1          1     1
Array2     1          1     1
Array3     1          1     1
Array4     5          2     1
Array5     5          2     1
Array6     5          2     1
Array7     9          3     1
Array8     9          3     1
Array9     9          3     1
Array10   13          4     1
Array11   13          4     1
Array12   13          4     1
Array13   17          5     1
Array14   17          5     1
Array15   17          5     1

>
> sessionInfo()
R version 2.6.2 (2008-02-08)
i386-pc-mingw32

locale:
LC_COLLATE=Portuguese_Brazil.1252;LC_CTYPE=Portuguese_Brazil.1252;LC_MONETARY=Portuguese_Brazil.1252;LC_NUMERIC=C;LC_TIME=Portuguese_Brazil.1252

attached base packages:
[1] tools     stats     graphics  grDevices utils     datasets
methods   base

other attached packages:
[1] maSigPro_1.10.0 Biobase_1.16.3
>

Could you help me? What I could do wrong?

Thank you very much.

-- 
Marcelo Luiz de Laia
Jaboticabal - SP - Brazil

"Você vê as coisas como elas são e pergunta: por quê? Mas eu sonho com
coisas que nunca foram e pergunto: por que não? " - Bernard Shaw



More information about the Bioconductor mailing list