Hi,

I have tried to find a beginners guide/tutorial on how to do a paired analysis in Ballgown. I have several samples pre and post treatment: ID<- seq(1,170, by=1) cond<-rep(c("post","pre"),85) pair<- rep(1:85, each = 2) test <- data.frame(id=ID,cond=cond, pair=factor(pair))

Am I right that I run the following? (after creating a ballgown object bg)

paired<- stattest(bg, feature='transcript', meas='FPKM', covariate='cond', adjustvars='pair')

Please help!