[BioC] interaction effects for DEXSeq
Simon Anders
anders at embl.de
Wed Jun 12 17:10:10 CEST 2013
Just looking over my own plot and noticed a typo:
On 12/06/13 09:15, Simon Anders wrote:
> Not quite. This tests for main effects and interaction in one go, so it
> returns as hits all exons whose usage differs between colonies _and/or_
> between strains. You are looking for interactions, i.e., for exons whose
> usage is different between strains _and_ where this difference itself
> differs between colonies. So, you need
>
> formuladispersion<-count~sample+(colony*strain)+exon
^ should be "*exon"
>
> formula0<-count~sample+exon+(colony+strain)*exon
>
> formula1<-count~sample+exon+(colony+strain)*exon+(colony:strain)*I(exon==exonID)
>
>
> Things look a bit simpler if you use the new "TRT" functions:
>
> formuladispersion<-count~sample+(colony*strain)+exon
^ here too
>
> formula0<-count~sample+exon+(colony+strain)*exon
> formula1<-count~sample+exon+(colony*strain)*exon
More information about the Bioconductor
mailing list