[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