The editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: Ballgown for differential expressed genes
gravatar for marasodi
8 months ago by
marasodi0 wrote:

 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




ballgown • 178 views
ADD COMMENTlink modified 8 months ago by James W. MacDonald49k • written 8 months ago by marasodi0
Answer: Ballgown for differential expressed genes
gravatar for James W. MacDonald
8 months ago by
United States
James W. MacDonald49k wrote:

The help page for stattest says this:

timecourse: if 'TRUE', tests whether or not the expression profiles of
          genomic features vary over time (or another continuous
          covariate) in the study.  Default 'FALSE'.  Natural splines
          are used to fit time profiles, so you must have more
          timepoints than degrees of freedom used to fit the splines.
          The default df is 4.

covariate: string representing the name of the covariate of interest
          for the differential expression tests.  Must correspond to
          the name of a column of 'pData(gown)'. If 'timecourse=TRUE',
          this should be the study's time variable.

If you are doing a timecourse study, you have to tell stattest that's what you are doing. And if there is no 'Time' variable in the pData object for your gown object, then stattest is going to tell you that it can't find it.


ADD COMMENTlink written 8 months ago by James W. MacDonald49k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 169 users visited in the last hour