Displaying 1 result from an estimated 1 matches for "todos_v1_e_v2_back_c".
2007 Nov 13
0
Difficult to set a quiet formula in maanova
...ent time RNA was collected from treated and from not treated
plants for both varieties. So, I have:
Var: 2 (varieties)
Trat: 2 (treatment)
Time: 3
Sample: 3 (biological replicate)
Probes: 3575
Spots: 2 for each probe
I used the following code:
> library(vsn)
> todos <- read.matrix("todos_v1_e_v2_back_c.txt",sep="\t")
> todos.norm <- vsn2(todos)
> write.table(exprs(todos.norm),"todos.norm.txt",sep="\t")
> library(maanova)
> fabiana.raw <- read.madata("todos.norm.vsn2.maanova.txt", designfile="design.txt", header=TRUE, spotfla...