similar to: Dotmatrix Plots

Hi, I am very new to R and wanted to know if there is a package that, given very long nucleotide sequences, searches and identifies short (7-10nt) motifs.. I would like to look for enrichment of certain motifs in genomic sequences. I tried using MEME (not an R package, I know), but the online version only allows sequences up to MAX 60000 nucleotides, and that's too short for my needs..

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

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, I have a data which is in text file and i would like to import the data to R. From the manual, i?ve found the read.table command function is the most appropriate but when i wrote the command an error had occur. It say ?Error in read.table"C:/Users/user/Documents/cfa-1.txt", header = T, sep = "\t",skip=10) :more columns than column names?. Please help me with this as

Hi - When I'm trying to read in a text file into a labeled character array, the memory stamp/footprint of R will exceed 4 gigs or more. I've seen this behavior on Mac OS X, Linux for AMD_64 and X86_64., and the R versions are 2.4, 2.4 and 2.2, respectively. So, it would seem that this is platform and R version independant. The file that I'm reading contains the upstream regions

Dear R-help, I am using Sweave and pdflatex to generate a large report from data contained in my database (Postgres via RODBC). Currently, I work with a single R/Sweave file, containing several "chapter" indications for the Latex engine. My master tex file sets the document class, and includes the introduction, the main Sweave file, and a conclusions and reference file. I use a makefile

A couple days ago, Mark Leeds asked about a solution that would basically stagger two lists, a and b, to return a list in the form of a[1], b[1], a[2], b[2], a[3].... In particular, the summary of his question was in reference to lists defined by x <- 5 tempin <- seq(1,1411, by=30) a <- tempin b <- tempin + x I offered the following function everyOther <-

Hi all, I am learning R by "doing". And this is my first post. I want to use R: 1- to fetch a DNA sequence from a databank (see bellow) and 2- store it as FASTA file. The problem: neither an error is prompted nor the fasta file is created. Testing the code (see bellow), I notice that everything works until the *"write.dna" *command - which is not creating the fasta file.

All, Recently I ran across a URL documenting published research using R: http://www.cgd.ucar.edu/ccr/ammann/millennium/CODES_MBH.html A note on the site indicates that the code is being revised. The code and data are provided, so that one could reproduce the results without having to buy a proprietary software program. In poking around the R website it is clear that a lot of thought has gone

Dear Mag. Ferri Leberl, I'm using something like: ----------------------- tex.tex --------------------------- \documentclass{article} \usepackage{graphicx} \usepackage{fancyvrb} \newcommand{\Rlogo}{\protect\includegraphics[height=1.8ex,keepaspectratio]{Rlogo.pdf}} \newcommand{\myinput}[1] {\begin{scriptsize} \VerbatimInput[frame=single,label=#1]{#1} \end{scriptsize}} \title{The R logo,

Dear Mag. Ferri Leberl, I'm using something like: ----------------------- tex.tex --------------------------- \documentclass{article} \usepackage{graphicx} \usepackage{fancyvrb} \newcommand{\Rlogo}{\protect\includegraphics[height=1.8ex,keepaspectratio]{Rlogo.pdf}} \newcommand{\myinput}[1] {\begin{scriptsize} \VerbatimInput[frame=single,label=#1]{#1} \end{scriptsize}} \title{The R logo,

Dear all, I'm trying to produce two dotcharts side-by-side within a Sweave document. When I'm compiling this example: \documentclass{article} \begin{document} <<fig=T,width=8,height=4>>= par(mfrow = c(1, 2), cex = 0.7) for(i in 1:2) dotchart(1:10) @ <<fig=T,width=8,height=4>>= par(mfrow = c(1, 2), cex = 0.7) for(i in 1:2) hist(1:10) @

Dear R-help, After making a typo (reminiscent of FORTRAN 77, I guess) I found the following: > identical(all.equal(2^(-10:10), 2**(-10:10)), TRUE) [1] TRUE I have tried to find the documentation about the ** operator but I was unsuccesful this way: > sessionInfo() Version 2.3.1 (2006-06-01) powerpc-apple-darwin8.6.0 attached base packages: [1] "methods" "stats"

> Dear useRs, > > by a circular diagram representation I mean what you will get by entering > this at your R promt: > > pie(1:5) > > Nice to have R as a lingua franca :-) > > The folowing quote is from page 360 in this very interesting paper: > > @article{SpenceI2005, > title = {No Humble Pie: The Origins and Usage of a Statistical Chart}, >

Hallo everyone, excuse me if this is not a genuine R question but I do not know where to ask else. Referring to e.g. https://stat.ethz.ch/pipermail/r-help/2004-December/062114.html I wonder if these measurements of 3000 criminals (raw data) are available anywhere. At least I didn't find them in the R datasets package or by means of Google. What I did find was a table of frequencies of

Dear R-help, I would like to include the documentation of an R function in an *.rnw document processed by Sweave. Because I'm sharing my *.rnw files with colleagues under Linux and Windows (I'm on Mac OS X), I would like a pure R solution. The naive approach doesn't work, because Sweaving this *.rnw file: -------- tmp.rnw -------- \documentclass{article} \begin{document}

I suspect you meant WD <- "~/Documents/Scripting/R_Studio/Sequences/" but I am entirely unfamiliar with the packages you are using, and know nothing about what is on your hard drive. For future reference: A) Read the Posting Guide. This is a plain text email list, and your html formatting gets removed leaving a mess that is not always readable. B) Most frequent users of R

Dear all, I'm confused by the time serie parameters in data(lh) : sueoka:~ lobry$ R --vanilla --quiet > tsp(lh) [1] 1 48 1 because documentation says: QUOTE A regular time series giving the luteinizing hormone in blood samples at 10 mins intervals from a human female, 48 samples. UNQUOTE So that I would expect the time serie to end at 480 minutes or 8 hours. Shouldn't we have

Hi, I want to analyse DNA sequence data (mtDNA) in R as in calculate Fst, Heterozygosity and such summary statistics. Package Adagenet converts the DNA sequence into retaining only retaining the polymorphic sites and then calcuates Fst.. but is there any other way to do this? I mean analyse the DNA sequence as it is.. and calculate the statistics? Thanks! [[alternative HTML version deleted]]

Dear All, I'am looking for examples showing that correlation does not imply causality, the targeted audience consists of undergraduate students (their first year at the university but in the BioMathStat track). All practicals are under R. I was able to extract this from R datasets: ### begin data(sunspots) data(lynx) spots <- window(sunspots, freq = 1, start = 1880, end = 1900) lnx <-