regex - matching and counting strings (k-mer of DNA) in R -

-

I have a list of strings (DNA sequence including AN, T, C, G) I can find all the matches and In the table whose columns are all possible combinations of those DNA alphabets (4 ^ k; - Kashmir mer - "K" is the length of each match and should be specified by the user) and represent the number of rows Match the sequence in a list

Let's say 5 members in my list include: 

  DNAlst & LT; -List ("CAAACTGATTTT", "GATGAAAGTAAAATACCG", "ATTATGC", "TGGA", "CGCGCATCAA")

me k = 2 (2-mer) , So the 4 ^ 2 = 16 combinations are available including AA, AT, AC, AG, TA, TT, ...

Then my table will have 5 rows and 16 columns . I want to calculate the number of matches between my K-Mers and the list members.

My desired result: df: 

  AT AC AG TA TT TC lstMemb AA ... 1 2 1 1 0 0 3 0 2 ... 3 4 5

Can you help me implement it in R?

If you are looking for speed, then there is a clear solution package for calculation pattern stri_count_fixed function and now, code and benchmark check 

  DNAlst & LT; -List ("CAAACTGATTT", "GATGAAAGTAAAATACCG", "ATTATGC", "TGGA", "CGCGCATCAA") DNA & lt; -! Results (L (C, "A", "C", "G", "T"), Every = 4), C (Stri_paste ("A", "C", "G", "T")) "A", "L") [, 1] DNA result AA AC AG AT CA CC Coast Guard GTGC GG GT TA TC TG TT2 1 0 1 1 -; (Sapply (DNAlst, stri_count_fixed, pattern = DNA, overlap = TRUE)) colnames (results) <1 1 0 0 0 0 1 3 [2,] 5 1 1 2 0 1 1 0 2 0 0 1 2 0 1 0 [3,] 0 0 0 0 0 0 0 0 0 1 0 0 1 0 1 1 [4,] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 [5,] 0 0 1 2 0 2 0 0 0 0 0 0 0 0 fstri & lt; - Function (X) {T (sapply (x, stri_count_fixed, dna, T))} fbio & lt; - Function (x) {T (s, x, function (x) {x1 & lt; - DNAString (x); oligonucleotide frequency (x1,2)})) all (fstri (DNAlst) == fbio (DNAlst)) # Pills are the same [1] TRUE longDNA & lt; - Sample (DNAlst, 100, T) fstri (longDNA), fbio (longDNA) unit: Microsoft Expression Minimum LBA Median Average UK Maximum Naval FSTI (Long DNA) 68 9.378 738.184 825.3014 766.862 793.134 6027.039 100 FBOO (Long DNA ) 118371.825 125552.401 129543.6585 127245.489 129165.711 359335.294 100 127245.489 / 766.862 ## 165.9301

CA 165x times faster :)

< / Html>

Comments

Post a Comment

Popular posts from this blog

mysql - How to enter php data into a html multiple select box -

-
I am working on a page that will create a team. I just want to create a selection box which will list all the users of all the users. After that, the user can select more than one friend (i.e. 5) and can send the request to these 5 people to join the team but I am having problems I did this: Select $ sql = "Select i. *, Join * Connect with M * m.id = i.clique_id where adder_id =: id"; $ Stmt = $ db- & gt; Query ($ sql, array ('id' = & gt; $ _SESSION ['id'])); While ($ record = $ stmt-> fetch ()) {$ surname = $ record ['surname']; } and this: & lt; Select multiple = "true" name = "category" well id = "category" name = "category" category = "" & gt; & Lt; Option value = "& lt ;? php echo $ surname ;? & gt;" & Gt; & Lt ;? Php ek $ nickname; ? & Gt; & Lt; / Options & gt; & Lt; / Select & gt; But only a friend's nickname...
Read more

java - Can't add JTree to JPanel of a JInternalFrame -

-
Image
There are two panels in my JInternalframe. I want to add TopPanel to Jtree to nominate another. But I can not add Jaitari to the top of the panel. Please help me in this piece of code: DefaultMutableTreeNode root = new DefaultMutableTreeNode ("deck"); DefaultMutableTreeNode Item Clubs = New DefaultMutableTreeNode ("Club"); AddAllCard (itemClubs); Root.add (itemClubs); DefaultMutableTreeNode item = New DefaultMutableTreeNode ("Diamonds"); AddAllCard (itemDiamonds); Root.add (itemDiamonds); DefaultMutableTreeNode Itemspads = New DefaultMutableTreeNode ("Spades"); AddAllCard (itemSpades); Root.add (itemSpades); DefaultMutableTreeNode item increases = new DefaultMutableTreeNode ("heart"); AddAllCard (itemHearts); Root.add (itemHearts); Default Trimodel Tree Model = New DefaultTreeModel (Route); Tree = new jetty (tree modal); ScrollPane = new JScrollPane (tree); // scrollPane.setViewportView (tree); . ScrollPane.getViewport () add (tr...
Read more

c++ - Cassandra datastax cpp driver - avoiding unnecessary copies -

-
I use the following method to create my insert statement for Cassandra (using datastax cpp driver) : ssStmt1 & lt; & Lt; "Insert some table (...) value (....)" string sStatement1 = ssStmt1.str (); Castestring cass1 = cass_string_init (sStatement1.c_str ()); CassStatement * pStmt1 = cass_statement_new (cass1, 0); I have also tried the following: CassString cass1 = cass_string_init ("Enter some table (id, date, c) value (? ,?,?); "); CassStatement * pStmt1 = cass_statement_new (cass1, 3); Cass_statement_bind_uuid (pStmt1, 0, uuidKey); Cass_statement_bind_int64 (PSTMT1, 1, timestamp); Cass_statement_bind_string (pStmt1, 2, sSomething); The first pattern works with mixed keys, the second is not. The key has a timestamp and a UUID. Timestamp is one of two errors in the field (based on the used bond statement). If bound in this way: cass_statement_bind_int64_by_name (pStmt1, "date", timestamp); I get the following error: Inva...
Read more