Search
Question: masigpro Error in terms.formula(formula, data = data) : '.' in formula and no 'data' argument
0
gravatar for john
2.5 years ago by
john40
Germany
john40 wrote:

Hello fellow researchers,

using masigpro i encountered an error which I cant seem to fix. Hopefully someone is working with masigpro as well or the developers are still around.

As you can see the p.vector function does what it should with "ncounts" beeing the normalized counts and d_1 the design matrix.

Calling T.fit() however gives me the following error:

 

> NBp_1 <- p.vector(ncounts, d_1, counts=TRUE,min.obs=3)
[1] "fitting gene 100 out of 165"
> NBt_1 <- T.fit(NBp_1)
[1] "fitting gene 100 out of 120"
Error in terms.formula(formula, data = data) : 
  '.' in formula and no 'data' argument
To show you the data used as input:

> d_1
$dis
           group.2vsgroup.1 time timexgroup.2 time2 time2xgroup.2
2h                        1    2            2     4             4
4h                        1    4            4    16            16
8h                        1    8            8    64            64
10h                       1   10           10   100           100
12h                       1   12           12   144           144
14h                       1   14           14   196           196
24h                       1   24           24   576           576
h2Control                 0    2            0     4             0
h4Control                 0    4            0    16             0
h8Control                 0    8            0    64             0
h10Control                0   10            0   100             0
h12Control                0   12            0   144             0
h14Control                0   14            0   196             0
h24Control                0   24            0   576             0

$groups.vector
[1] "group.2vsgroup.1" "group.1"          "group.2vsgroup.1" "group.1"          "group.2vsgroup.1"

$edesign
           time replicates group.1 group.2
2h            2          1       0       1
4h            4          2       0       1
8h            8          3       0       1
10h          10          4       0       1
12h          12          5       0       1
14h          14          6       0       1
24h          24          7       0       1
h2Control     2          8       1       0
h4Control     4          9       1       0
h8Control     8         10       1       0
h10Control   10         11       1       0
h12Control   12         12       1       0
h14Control   14         13       1       0
h24Control   24         14       1       0
> head(ncounts)
         2h 4h 8h 10h 12h 14h 24h h2Control h4Control h8Control h10Control h12Control h14Control h24Control
A1BG-AS1 14 17  8   4   4   7  15        17        17        17         17         17         17         17
A1CF      0  3  1   5   0   0   1         0         0         0          0          0          0          0
A2M      11 15 90  34  50  49  38         3         3         3          3          3          3          3
A2M-AS1   7  2 22   4   6  12   8         0         0         0          0          0          0          0
A2MP1    38 30 17  35  13  13   0         9         9         9          9          9          9          9
AAAS      6 59 83   4  67  97 150         1         1         1          1          1          1          1
> dim(ncounts)
[1] 165  14

 

 

Thanks a lot,

john

ADD COMMENTlink modified 18 months ago by paolo.uva0 • written 2.5 years ago by john40

Hi John,

I am wondering whether you found a solution to your problem in the end? I recently also started using maSigPro and have encountered the exact same problem. After a lot of trying and looking around I still haven't managed to fix it...

Kind regards,

Matthias

ADD REPLYlink written 11 days ago by naets_matthias0
0
gravatar for paolo.uva
18 months ago by
paolo.uva0
paolo.uva0 wrote:

Hi,

I had the same error with the T.fit function.
In my case the error was due the column names of the design matrix starting with numbers instead of letters (e.g. 10wt instead of wt10).
I simply renamed the columns and this time it worked.

Hope it helps,

Paolo

 

ADD COMMENTlink written 18 months ago by paolo.uva0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 146 users visited in the last hour