A note to my own post above:
> psCsv <- read.table(HuEx.small.csv,skip=1,sep=" ",comment.char="")
> colnames(psCsv) <- c("probeset_id","start", "stop", "probe_count")
If you are reading more columns in (like "strand") which are
characters, you should keep in mind, that R usually treats them as
Hi Frank.
Comments below.
On 05/03/2009, at 12:34 PM, Frank Ziegner wrote:
> Hi Mark,
>
> I have some questions regarding you plot and some more general
> questions regarding probesets/exons.
>
> I assume your plot shows the WNK1 (ENSG0060237) gene.
>
> When doing:
> u<-which(getUnitNa
Hi everybody,
based on Mark's hints I wrote some code to plot FirmaScores with
GenomeGraphs. Mark's way of extracting the data (d <-log2(extractMatrix
(cs,cells=ind,verbose=verbose)) ) is very elegant and effective. But
somehow I didn't manage to pull the firma scores data off in the same
way.
Th
Hi Mark,
thanks for your help. The awk-hint is wonderful. The annotation file
is MUCH smaller and better manageable in R.
I want to give some advice for Windows-Users:
- get the UnxUtils (update your PATH variable to use the commands
system-wide)
- the awk command is called gawk in UnxUtils
- on
Hi Frank.
A very good question and unfortunately, I haven't had time to put a
vignette together. I have given a working example below that will at
least to get you started.
On 03/03/2009, at 9:34 AM, Mark Robinson wrote:
> [this was sent to me offline. i've created a new thread and will t