Mercurial > repos > stevecassidy > wrassp

annotate outputAssp.R @ 3:f188eb0b526d draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 2d22bc12c66a3d1949d4bf4c84e7ce7c4610383c-dirty
author stevecassidy
date Wed, 14 Dec 2016 21:53:57 -0500
parents dbcf9bc275e3
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_est <- 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
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 for(idx in seq(2,length(attr(assp,'names')))) {
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
6 result <- cbind(result, assp[idx])
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
7 }
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
8 # add a column of timestamps
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
9 start <- attr(assp, 'startTime')
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
10 rate <- attr(assp, 'sampleRate')
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
11 sampletime <- seq(start, by=1/rate, length.out=nrow(result))
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
12 # and a column of 1s for some reason
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
13 #one <- rep(1, nrow(result))
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
14 result <- cbind(sampletime, result)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
15
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
16 # write header
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
17 #write("EST_File Track", file=outfile)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
18 #write("DataType ascii", file=outfile, append=T)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
19 #write("BreaksPresent true", file=outfile, append=T)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
20 #write(paste("NumFrames", nrow(result)), file=outfile, append=T)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
21
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
22 #count <- 0 # count overall number of columns
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
23 #vcount <- 1 # count variables in assp obj
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
24 #for(name in attr(assp,'names')) {
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
25 # for (coln in seq(ncol(data.frame(assp[vcount])))) {
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
26 # colname = paste(name, coln, sep="_")
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
27 # write(c(paste("Channel_", count, " ", colname, sep="")), file=outfile, append=T)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
28 # count <- count+1
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
29 # }
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
30 # vcount <- vcount+1
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
31 #}
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
32 #write(paste("NumChannels", count), file=outfile, append=T)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
33 #write("EST_Header_End", file=outfile, append=T)
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
34
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
35
dbcf9bc275e3 planemo upload commit 063e049d569aeb45e1008cbf044f2dad850eca3d-dirty
stevecassidy
parents:
diff changeset
36 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
37 }