Mercurial > repos > lecorguille > anova

comparison abims_anova.r @ 11:102049093b7d draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/workflow4metabolomics/anova commit 4922313a0e9569326b7723c41babb89f998dbfd9
author lecorguille
date Tue, 13 Mar 2018 09:47:21 -0400
parents b147b17759a6
children
comparison
equal deleted inserted replaced
10:b147b17759a6 11:102049093b7d
22 if (interaction) anova_formula_operator = "*" 22 if (interaction) anova_formula_operator = "*"
23 23
24 # -- import -- 24 # -- import --
25 data=read.table(file, header = TRUE, row.names=1, sep = sep, quote="\"", dec = dec, fill = TRUE, comment.char="",na.strings = "NA", check.names=FALSE) 25 data=read.table(file, header = TRUE, row.names=1, sep = sep, quote="\"", dec = dec, fill = TRUE, comment.char="",na.strings = "NA", check.names=FALSE)
26 26
27 if (mode == "row") data=t(data) 27 if (mode == "row") {data=t(data)} else {data=as.matrix(data)}
28 28
29 sampleinfoTab=read.table(sampleinfo, header = TRUE, row.names=1, sep = sep, quote="\"") 29 sampleinfoTab=read.table(sampleinfo, header = TRUE, row.names=1, sep = sep, quote="\"")
30 rownames(sampleinfoTab) = make.names(rownames(sampleinfoTab)) 30 rownames(sampleinfoTab) = make.names(rownames(sampleinfoTab))
31 31
32 varinfoTab=read.table(varinfo, header = TRUE, sep = sep, quote="\"") 32 varinfoTab=read.table(varinfo, header = TRUE, sep = sep, quote="\"")