I was looking for the mutation data through TCGA portal using TCGAbiolinks and I have realized that sample size are not the same.
for instance TCGA-OV case TCGA data portal shows 419 cases, however TCGAbiolinks shows 462 samples. File counts are the same for both it is 482.
so why it is different?
this my query in TCGA data portal:
cases.project.project_id in ["TCGA-OV"] and files.analysis.workflow_type in ["Aliquot Ensemble Somatic Variant Merging and Masking"] and files.data_category in ["Simple Nucleotide Variation"] and files.data_type in ["Masked Somatic Mutation"]
this is same query in the TCGAbiolinks package:
#query
query <- GDCquery(
project = "TCGA-OV",
data.category = "Simple Nucleotide Variation",
access = "open",
data.type = "Masked Somatic Mutation",
workflow.type = "Aliquot Ensemble Somatic Variant Merging and Masking"
)
#download & read
GDCdownload(query)
maf <- GDCprepare(query)
mafr = maftools::read.maf(maf)
mutations = mafSummary(mafr)
print(as.numeric(mafr@summary[mafr@summary$ID=="Samples"]$summary))
1 answer
In the GDC query, you got 419 cases and 482 files (likely 482 aliquots). In the tcgabiolinks query, you got 462 samples. You are comparing apples to oranges.
Most of the cases in GDC have at least one tumor sample and one normal sample, and some could have more tumor samples such as metastasis and new primary, etc. So case count is not sample count.
In the link you have, there are only case tab and file tab. There are no summary tab for samples. If you really want to get samples, you can learn the GDC API, or add all files into cart, and download sample sheet from the cart.
Log in to answer this question.
You're comparing samples to cases. Can you check aliquot counts in both cases?
I thought 482 files = aliquots, isn't it like that? Or in another way to ask how can I find the sample number of given TCGA query in the portal? this is the query link TCGA-OV
I'm not entirely sure that num_files would equal num_aliquots. Please try and dig deeper to check if that's the case. I apologize, but I don't have the time to do a TCGA deep dive right now.