Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • lagrace
    Junior Member
    • May 2017
    • 4

    #1

    NA's in data (using TopHat-HTSeq-DESeq)

    Hello all,

    I have some data that I have retrieved from GEO. I am using the TopHat-HTSeq-DESeq pipeline to analyze this data. Within linux, I am using R and ran into an issue with NA's in the data. I am trying to determine why there are so many NA's. I was told to view the annotations. Is that possible in R? If so, how?
    Not sure if there is some other reason. Below are the codes used:

    samTab=read.table(file="samTab.txt",sep="\t")
    samTab
    df=newCountDataSetFromHTSeqCount(samTab)
    df

    cds<-estimateSizeFactors(df)
    sizeFactors(cds)
    head(counts(cds,normalized=TRUE))

    cds=estimateDispersions(cds)
    head(cds)
    str(fitInfo(cds))
    head(fData(cds))

    res=nbinomTest(cds,"Nor","Can")
    head(res)
    id baseMean baseMeanA baseMeanB foldChange log2FoldChange pval padj
    A1 138.33103 129.97950 146.68255 1.1285052 0.17441313 0.9057131 1


    resSig = res[res$padj < 0.1, ]
    head(resSig[order(resSig$pval), ])
    id baseMean baseMeanA baseMeanB foldChange log2FoldChange pval padj
    NA <NA> NA NA NA NA NA NA NA
    NA.1 <NA> NA NA NA NA NA NA NA
    NA.2 <NA> NA NA NA NA NA NA NA
    NA.3 <NA> NA NA NA NA NA NA NA
    NA.4 <NA> NA NA NA NA NA NA NA
    NA.5 <NA> NA NA NA NA NA NA NA

    I don't understand why the padj is 1. I suppoose that is the reason for the NA's. Is there a way to view the annotations while in R? Any suggestions? Please assist.

Latest Articles

Collapse

  • SEQadmin2
    Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
    by SEQadmin2



    CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

    Despite this, “CRISPR helped turn genome editing from a specialized technique into
    ...
    07-31-2026, 11:01 AM
  • SEQadmin2
    Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
    by SEQadmin2


    Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

    The systematic characterization of the human proteome has
    ...
    07-20-2026, 11:48 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, Today, 12:22 PM
0 responses
12 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-11-2026, 10:35 AM
0 responses
11 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-06-2026, 07:41 AM
0 responses
30 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-03-2026, 10:13 AM
0 responses
48 views
0 reactions
Last Post SEQadmin2  
Working...