Mercurial > repos > stevecassidy > wrassp

annotate util.R @ 1:dbcf9bc275e3 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
author stevecassidy
date Wed, 07 Dec 2016 19:10:51 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
1 assp_to_tsv <- function(assp, outfile) {
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
2
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
3 # convert to a dataframe, need to add one column at a time
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
4 result <- data.frame(assp[1])
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
5 if (length(attr(assp,'names'))>1) {
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
6 for(idx in seq(2,length(attr(assp,'names')))) {
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
7 result <- cbind(result, assp[idx])
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
8 }
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
9 }
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
10 # add a column of timestamps
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
11 start <- attr(assp, 'startTime')
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
12 rate <- attr(assp, 'sampleRate')
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
13 time <- seq(start, by=1/rate, length.out=nrow(result))
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
14 result <- cbind(time, result)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
15
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
16 write.table(result, file=outfile, sep="\t", quote=F, row.names=F, col.names=T, append=F)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
17 }