Blog Archives

When can we expect the last damn microarray paper?

January 18, 2012
By
When can we expect the last damn microarray paper?

With bonus R codeIt came as a shock to learn from PubMed that almost 900 papers were published with the word "microarray" in their titles last year alone, just 12 shy of the 2010 count. More alarming, many of these papers were not of the innocuous "Microarray study of gene expression in dog scrotal tissue" variety, but dry...

Read more »

Making R’s paste act more like CONCAT

October 31, 2011
By

While vector-friendly, R's paste function has a few behaviors I don't particularly like. One is using a space as the default separator:> adjectives> paste(adjectives,"er")> paste(adjectives,"er") "lean er" "fast er" "strong er" #d'oh> paste(adjectives,"er",sep="") "leaner" "faster" "stronger"Empty vectors get an undeserved first class treatment: > paste(indelPositions,"i",sep="") "i"> indelPositions> paste(indelPositions,"i",sep="") "5i"...

Read more »

Installing RStudio Server on Scientific Linux 6: My bash notebook

August 18, 2011
By
Installing RStudio Server on Scientific Linux 6: My bash notebook

Granted, not a brilliant sysadmin mind at work here, but this might help someone someday.Scientific Linux (SL) is built from Red Hat Enterprise LinuxSee installation instructions here:http://rstudio.org/download/server$ sudo rpm -Uvhhttp://download.fedoraproject.org/pub/epel/6/x86_64/epel-release-6-5.noarch.rpm password for leipzig: Retrievinghttp://download.fedoraproject.org/pub/epel/6/x86_64/epel-release-6-5.noarch.rpmwarning: /var/tmp/rpm-tmp.S2RQAH: Header V3 RSA/SHA256 Signature, key ID0608b895: NOKEYPreparing... ...

Read more »

RStudio: My thoughts

March 15, 2011
By
RStudio: My thoughts

Let me get this out of the way: I just love RStudio.Created by a team lead by JJ Allaire, a name that should ring a bell if you were involved in web development during the Clinton administration, RStudio is an R IDE that is actually designed for R from...

Read more »

Chromosome bias in R, my notebook

December 23, 2010
By
Chromosome bias in R, my notebook

My goal is to develop a means of detecting chromosome bias from a human BAM file.Because I've been working with proprietary and novel plant genomes for the last three years, I haven't had the chance to use any of the awesome UCSC-based annotational features that have been introduced and refined in Bioconductor until now. I've returned to biomedical research...

Read more »

Getting the basics from readAligned

March 9, 2010
By

The UCR guide is a little sparse with regard to getting basic information from readAligned.I'd like to add to the general cookbook. If some bioc people out there can contribute some alignment recipes can fill me in on some more basics please comment:alignedReads #how many reads did I attempt to align#i don't think you can't get this from...

Read more »

Quality trimming in R using ShortRead and Biostrings

March 3, 2010
By

I wrote an R function to do soft-trimming, right clipping FastQ reads based on quality.This function has the option of leaving out sequences trimmed to extinction and will do left-side fixed trimming as well.#softTrim#trim first position lower than minQuality and all subsequent positions#omit sequences that after trimming are shorter than minLength#left trim to firstBase, (1 implies no left trim)#input:...

Read more »

R’s xtabs for total weighted read coverage

November 4, 2009
By

Samtools and its BioPerl wrapper Bio::DB:Sam prefer to give read coverage on a depth per base pair basis. This is typically an array of depths, one for every position that has at least one read aligned. OK, works for me. But how can we quickly see which targets (in my case transcripts) have the greatest total weighted read coverage...

Read more »

Standardized Velvet Assembly Report

August 25, 2009
By
Standardized Velvet Assembly Report

http://code.google.com/p/standardized-velvet-assembly-report/I finally got my Velvet Assembler report script up on google code. This "program" consists of some short scripts and a Sweave report designed to help Velvet users identify the optimal kmer an...

Read more »

Calculating an N50 from Velvet output

November 25, 2008
By

In sequencing circles the N50 length is a useful heuristic for judging the quality of an assembly. Here is my definition of N50 length, which you may or may not find intuitive: N50 length is the length of the shortest contig such that the sum of contigs of equal length or longer is at least 50% of the total length...

Read more »

Sponsors

Mango solutions



RStudio homepage



Zero Inflated Models and Generalized Linear Mixed Models with R

Quantide: statistical consulting and training

datasociety

http://www.eoda.de





ODSC

ODSC

CRC R books series





Six Sigma Online Training









Contact us if you wish to help support R-bloggers, and place your banner here.

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)