similar to: motif search

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]]

Dear all, do you know if there is any particular version of R to implement with windows 64 bit, in such a way to increase the amount of memory it can use? How should I increase the memory, and more importantly to set a higher max vector size? It still stops me saying "Could not allocate vector of size 145" thanks to all alessia

The myoglobin sequence, with reference number NM_005368 in Gen bank, has the following frequencies of DNA nucleotides: A C G T 237 278 309 242 Do these data provide sufficient evidence, at the 1% level of significance, that the DNA nucleotides have an unequal distribution, that is the DNA nucleotides are not evenly utilised? Clearly state your hypothesis, test statistic and conclusion.

Hi Function sequence() repeatedly concatenates its output, and this is slow. It is possible to improve on the performance of sequence by defining myseq <- function(x){unlist(sapply(x,function(i){1:i}))} The following session compares the performance of myseq(), and sequence(), at least on my G5: > identical(sequence(1:50),myseq(1:50)) [1] TRUE > system.time(ignore <-

Hi Function sequence() repeatedly concatenates its output, and this is slow. It is possible to improve on the performance of sequence by defining myseq <- function(x){unlist(sapply(x,function(i){1:i}))} The following session compares the performance of myseq(), and sequence(), at least on my G5: > identical(sequence(1:50),myseq(1:50)) [1] TRUE > system.time(ignore <-

Hi, I have encountered this problem quite a few times and thought I would ask. Let's say that I have two endpoints, a and b, which are integers. If a <= b, I would like to get a:b, but if a > b, then numeric(0), for example: myseq(3,5) => 3:5 myseq(3,3) => 3 myseq(3,2) => numeric(0) The operator : just gives decreasing sequences in the latter case, and I could not coax seq

Dear all, I just installed the new version of R, 2.10.1, and I am currently using the package sparseM. (I also use a 64 bit windows version) I got a problem that I never had: when I try to multiply with a kronecker product (%x%) two sparse matrixes I get the following message: Error in dim(x) <- length(x) : invalid first argument I never had this problem with previous versions of R. May

Hi. Suppose I have a vector that I partition into disjoint, contiguous subvectors. For example, let v = c(1,4,2,6,7,5), partition it into three subvectors, v1 = v[1:3], v2 = v[4], v3 = v[5:6]. I want to find the maximum element of each subvector. In this example, max(v1) is 4, max(v2) is 6, max(v3) is 7. If I knew that the successive subvector maxima would never decrease, as in the example,

Last Friday, I noticed that it is difficult to work with regression models in which there are factors. It is easier to do the old fashioned thing of coding up "dummy" variables with 0-1 values. The predict function's newdata argument is not suited to insertion of hypothetical values for the factor, whereas it has no trouble with numeric variables. For example, if one uses a

Dear all, I 'm implementing the koenker procedure for quantile fixed effects. I would like also to apply weights to the procedure, so that to give more weight to the observation that better represent my original sample (much larger than it is possible to use in R). Do you know if it is possible? How could I solve this problem? Thank you alessia

Dear All, I am using the plotrix library to plot some matrices. I have a problem: some of my data are outliers, hence using a linear color scale does not work very well (you would see too many cells having a similar, indistinguishable color). See the code snipped at the end of the email. Plotting the logarithm of the data gets the job done, but my problem is that I would like to write in every

Hi, It might be a trivial question but I just wonder if you could advise good theoretical tutorials, reviews on NGS (different platforms) quality control like nucleotides quality by cycle, nucleotides frequency by cycle, GC content and distribution, K-mer frequency by cycle as well as quality control after alignment and mapping. Look forward to your reply, Carol [[alternative HTML version

Dear R users: > append(1:5, 0:1, after=2) [1] 1 2 0 1 3 4 5 If I want to repeat the appended value every 2 like the following: [1] 1 2 0 1 3 4 0 1 5 How should I modify? Thank you for any help.

I am a new R user. As a test, I want to write a simple code that does the following simulation: 1. Randomly generate a number from a distribution, say, Poisson. Let's say that number is 3. 2. Randomly generate 3 numbers from another distribution, say, Normal. 3. Compute the sum of the numbers generated in step 2 and read it into a vector, V. 4. Repeat steps 1 through 3 for 100,000 times. 5.

Dear all, I am running a quantile estimation with Sparse matrix and when I run the procedure rq.fit.sfn I receive the following warning: tiny diagonals replaced with Inf when calling blkfct. Does anyone knows exactly what does it mean? What's this kind of error? Should I get worried about this message? The matrix I use is a full rank matrix (and indeed I do not have any message about

Apparently there is one or more concepts that I do not fully understand from the descriptions of a function and the apply material. I have been reading the mail from this forum and have learned much but, in this case, what I have been reading here and from the manual isn't enough. The following code produces what I want with the for loop. From what I have read from this forum, a for

All Is there a package or library that will, given a nucleotide sequence 1. calculate the extinction coefficient at 260 nm for (Beer-Lambert's law) 2. calculate molecular weight 3. return it's complementary sequence I was able to find several packages that can do similar calculations for an amino acid sequence for proteins but none for nucleic acids. Any pointers, etc. would be

library(sos) (mp <- findFn('{molecular properties}')) ????? ** found 7 matches in 4 packages and opened two web pages in my default browser with (a) the 7 matches and (b) the 4 packages. The first function was something for amino acids, like you suggested.? Two others returned compound and substance information from PubChem. ????? Does this help? ????? Spencer On

In my database the id for a field is 810. RoR thinks that it is 809. My log file (below) shows that just before the insert it grabs the next sequence value. I dont think it should do that. This causes problems further down the line when trying to access the Person object. Development Log file : "SQL (0.016000) select people_seq.nextval id from dual Person Create (0.015000) INSERT

Hi, I have a piece of code originally written for s-plus - I am trying to run it in R now. The code was obtained from someone who is now not available to give any pointers and I am a beginner in R. Here is where it is getting stuck: > +names(good.motifs[,1]) Error in +names(good.motifs[, 1]) : Invalid argument to unary operator here is now names(good.motifs,1]) looks: >