Question: use MotifDb package with rGADEM?
0
gravatar for Paul Shannon
17 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 • 197 views
ADD COMMENTlink written 17 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: 212 users visited in the last hour