Enter the body of text here

Code should be placed in three backticks as shown below

query.exp.hg38 <- GDCquery(project ="TCGA-BRCA",
                           data.category = "Transcriptome Profiling",
                           data.type = "Gene Expression Quantification",
                           workflow.type = "STAR - Counts",
                           legacy = FALSE)#legacy=FALSE means we are download the data against hg38, otherwise it will be hg19
GDCdownload(query.exp.hg38, method = "api")
# Step 3: prepare the downloaded data (by default reading data from local folder C:/Users/dhe2/Documents/GDCdata/): 
maf <- GDCprepare(query.exp.hg38)

STAR - Counts is the raw counts, which is the best to start from for any analysis. If you want TPM and there is the gene length in that object use DESeq2: Is it possible to convert read counts to expression values via TPM and return these values? or alternatively use any of the common normalization methods from packages such as DESeq2 or edgeR/limma, check their manuals. By the way, it is appreciated if a post is not just code but at least some words to describe the problem outside of the title.