Question: use MotifDb package with rGADEM?
0
gravatar for Paul Shannon
13 months ago by
Paul Shannon370
Paul Shannon370 wrote:

Reading through the rGADEM man pages, I have yet to figure out how to supply PWMs of my own choosing for rGADEM to match against.  MotifDb provides a method ("query") for extracting a list of matrices by various criteria, for instance:

as.list(query(MotifDb, "MA0079.3"))[[1]]

           1          2          3 4           5         6 7 8         9        10         11
A 0.09812228 0.00000000 0.01133501 0 0.000000000 0.1555988 0 0 0.0000000 0.1206778 0.07465079
C 0.20448821 0.71158690 0.76746050 1 0.991641859 0.0000000 1 1 0.7647126 0.7278452 0.56892604
G 0.48900847 0.07350584 0.00000000 0 0.000000000 0.5294252 0 0 0.0000000 0.0000000 0.08403939
T 0.20838104 0.21490726 0.22120449 0 0.008358141 0.3149760 0 0 0.2352874 0.1514770 0.27238379

It would be very handy to pass such a matrix to the GADEM function, along with the DNAStringSet to match against.  Is this possible?  Sorry if it is there, documented, and I missed it.

 - Paul

matchpwm • 162 views
ADD COMMENTlink written 13 months ago by Paul Shannon370
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 16.09
Traffic: 262 users visited in the last hour