Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • DESeq/EdgeR on counts of reads supporting variant

    Is there precedent for and/or is it sensible to use DESeq or EdgeR to model counts of reads supporting variants in order to get significance?

    For example, suppose I have 10 cases and 10 controls in a 20 sample VCF, can I extract the counts of alternate alleles (after filtering) into a counts matrix to these tools without violating their assumptions?

    Some have recommended that I do a contingency table on number of cases / controls with / without the variant to get significance but that is discarding a lot of information.

    thanks,
    -Brent

    EDIT: I suppose the difficulty is adjusting for depth at each site.
    Last edited by brentp; 03-31-2014, 09:00 AM.

  • #2
    I would try a beta-binomial test (you can do this with the vglm function in the VGAM package). The beta-binomial distribution is similar to a regular binomial distribution, except the rates are derived from beta distributions (i.e, they're not a single number but a distribution that describes a population).

    The format would be something like:

    Code:
    library(VGAM)
    groups <- factor(c(rep("case",10),rep("control", 10)), level=c("control","case"))
    reference_count <- c(rbetabinom(10,seq(100,200,10),0.8), rbetabinom(10,seq(100,200,10),0.7)) #some example data
    alt_count <- rep(seq(100,190,10),2) - reference_count
    d <- data.frame(ref=reference_count, alt=alt_count, grp=groups)
    fit <- vglm(cbind(ref,alt)~grp, betabinomial, d)
    lrtest(fit)
    #One can get the coefficients with
    logit(coef(fit), inverse=T)
    #But keep in mind only the intercepts will be correct, so
    logit(sum(coef(fit)[c(1,3)]), inverse=T) #The actual "case" mu
    One can also fit the shape parameters, if you prefer. BTW, this is more or less how BiSeq works (the problem ends up being similar), though it uses a different package.

    Comment


    • #3
      Keep in mind that the distribution can be a bit tricky to fit, so you sometimes have to play around with parameters (don't blindly accept a fit without looking at the resulting deviations. You can also use the "binomial.ab" to fit the normal paramaterization of the distribution, though I suspect you're interested in changes in mu.

      Comment


      • #4
        Thanks Devon, I'll look into that more. I knew that BiSeq used the betaregression, but didn't think about the details.

        Some others have also suggested that I look into plink/seq or pseq.

        Comment


        • #5
          It's just a question of whether the data follows a binomial distribution or is over-dispersed. In the former case logistic regression (as done by Plink or in R) would be appropriate, but in the latter case not.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Choosing Between NGS and qPCR
            by seqadmin



            Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
            10-18-2024, 07:11 AM
          • seqadmin
            Non-Coding RNA Research and Technologies
            by seqadmin




            Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

            Nobel Prize for MicroRNA Discovery
            This week,...
            10-07-2024, 08:07 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 11-01-2024, 06:09 AM
          0 responses
          13 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-30-2024, 05:31 AM
          0 responses
          16 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-24-2024, 06:58 AM
          0 responses
          24 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-23-2024, 08:43 AM
          0 responses
          52 views
          0 likes
          Last Post seqadmin  
          Working...
          X