datatab<- read.table(system.file("doc", "exp2_2waycfa.txt", package="CFAssay"), header=TRUE, sep="\t")
names(datatab) # has columns "x5fuCis" and "siRNA"
fitcomp<- cfa2way(datatab, A="siRNA", B="x5fuCis", param="A/B")
print(fitcomp, labels=c(A="siRNA",B="x5fuCis"))
pdf("TwoWay_experimental_plots.pdf")
plotExp(fitcomp, labels=c(A="siRNA", B="x5fuCis"))
dev.off()
Run the code above in your browser using DataLab