I have differential expression RNA seq data, I need to run it on ball gown, the data is between two genotypes, at 4 time points each treated sample have its control and both the control and the treated plant are in triplicates, I used this
pData(bg_chrX_filt) = data.frame(id=sampleNames(bg_chrX_filt), group=rep(c(1,0), each=24)) when I run stattest
results_transcripts = stattest(bg_chrX_filt, feature="transcript", covariate="Time", getFC=TRUE, meas="FPKM")
I get this error
Error in stattest(bg_chrX_filt, feature = "transcript", covariate = "Time", :
invalid covariate name