similar to: Splitting vector

A version using regular expressions, regexpr() and substr() functions is attached. Finally everything is packed into splitSeq() function (chunk 14 in the attached file) Seq<- "GCCTCGATAGCTCAGTTGGGAGAGCGTACGACTGAAGATCGTAAGGtCACCAGTTCGATCCTGGTTCGGGGCA" Str<-

Hiya, I am trying to look at the similarities between a number of sequences, for example i am trying to see how similar "ababbbassdaa" is to "addffggssbbsbbs" I was wondering is the some way for me to see how similar they are in terms of, for example, number of a's, number of b's, how often a and ab are consecutive, how often abab is together etc. Any advice

Hello! I'm writing a script with a lot of loops and it executes really slowly over huge amounts of data. I assume it's because I don't know how to avoid using loops. Logical subscripts are more desirable, but I don't know how to implement them. One example of that issue: library(seqinr) GCsequence <- vector() for( i in 1:(length(data$sequence))) {

Hi, I am using ecdf-function and want to use the ecdf()-data-points for nls() as data-parameter. nls() expects 'list' or 'environment' as a data-type, knots(ecdf(mydata)) gives me 'numeric'. What should I do now? Thanks in advance - Jochen Here is the code: ################################################# # --- Fit --- # Gumbel-Dist-Function, cumulative,

A prebuilt binary. Needs non-free Watcom compiler to build. Signed-off-by: Lubomir Rintel <lkundrak at v3.sk> --- dosutil/mdiskchk.com | Bin 8106 -> 0 bytes 1 file changed, 0 insertions(+), 0 deletions(-) delete mode 100644 dosutil/mdiskchk.com diff --git a/dosutil/mdiskchk.com b/dosutil/mdiskchk.com deleted file mode 100644 index

Hi there, I am looking for R-packages that can help me visualize properties on nucleotide sequences. I want to display sequences in the 1-100K base range as lines and plot features above and below those lines. Any ideas would be welcome. Thanks, Bernd [[alternative HTML version deleted]]

I need to increment cells of a matrix (collusionM). The indexes to increment are in an index (matchIndex). This is some of the code for (j in 1:(rows-1)) matchIndex[[j]] <- which(mx[j]==mx) for (j in 1:(rows-1)) collisionM[j,matchIndex[[j]]] <- collisionM[j,matchIndex[[j]]] + 1 I could put them in the same loop but this is slower. The first for statement is fine. It finds the matches

Hi, I got two data point vectors. Now I want to make a ks.test(). I you print both vectors you will see, that they fit pretty fine. Here is a picture: http://www.jochen-bauer.net/downloads/kstest-r-help-list-plot.png As you can see there is one histogram and moreover there is the gumbel density function plotted. Now I took to bin-mids and the bin-height for vector1 and computed the

Hi R-Users, I read some texts related to KS-tests. Most of those authors stated, that KS-Tests are not suitable for binned data, but some of them refer to 'other' authors who are claiming that KS-Tests are okay for binned data. I searched for sources and can't find examples which approve that it is okay to use KS-Tests for binned data - do you have any links to articles or

Dear R-helpers: thank your for your attention. i am a newer to R and i am doing some protein category classification based on the amino acid sequence.while i have some questions urgently. 1. any packages for analysis amino acid sequence 2. given two sequences "AAA" and "BBB",how can i combine them into "AAABBB" 3. based on "AAABBB",how can i get some

Hello, I have a big number lets say of around hundred digits. I want to subset that big number into consecutive number of 5 digits and find the product of those 5 digits. For example my first 5 digit number would be 73167. I need to check the product of the individual numbers in 73167 and so on. The sample number is as follows:

Hello all, For some work I am doing on RNA, I want to use R to do string parsing that (I think) is like a simplistic HTML parsing. For example, let's say we have the following two variables: Seq <- "GCCTCGATAGCTCAGTTGGGAGAGCGTACGACTGAAGATCGTAAGGtCACCAGTTCGATCCTGGTTCGGGGCA" Str <-

Buenos dias a todos. Tengo este Barchart con dos tipos de datos. Desembarcado y Muestreado. Como los datos de Desembarcado son mucho mayores que los de Muestreado, estos ultimos datos apenas si se aprecian en la grafica. Lo que se me ocurre es cambiar la escala de este eje. ¿alguien sabria como hacerlo? O alguna otra sugerencia claro. Muchisimas gracias a todos los que me estan ayudando Jose Luis

When empty grids are shown (ex: no VMs in this pool) a specific graphic, message and improved action button are shown. For Task views, the user will see the current filter and an explanation that no tasks match that filter. Signed-off-by: Jeremy Perry <jeremy.perry at redhat.com> --- src/app/views/hardware/show_hosts.rhtml | 11 ++++------- src/app/views/hardware/show_storage.rhtml

Hey all, I''m not sure if this is a ZFS bug or a hardware issue I''m having - any pointers would be great! Following contents include: - high-level info about my system - my first thought to debugging this - stack trace - format output - zpool status output - dmesg output High-Level Info About My System --------------------------------------------- - fresh

> This message is in MIME format. Since your mail reader does not understand this format, some or all of this message may not be legible. --B_3294209241_521622 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit I have been unable to find any discussion regarding this issue, but I have been able to reproduce it on two different machines. Both computer are Macs running

Signed-off-by: Scott Seago <sseago at redhat.com> --- AUTHORS | 17 ++++++ README | 10 +++ conf/ovirt-agent | 12 ++++ conf/ovirt-db-omatic | 12 ++++ conf/ovirt-host-browser | 12 ++++