use MotifDb package with rGADEM?
0
0
Entering edit mode
Paul Shannon ▴ 470
@paul-shannon-5944
Last seen 2.6 years ago
United States

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 • 750 views
ADD COMMENT

Login before adding your answer.

Traffic: 597 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6