data-raw/expset_translation_id.R

In Berghopper/GAPGOM: GAPGOM (novel Gene Annotation Prediction and other GO Metrics)

# this block serves as a test to see if the package's algorithm matches up with the original.
# For this it also uses the original lncRNA2Function dataset. While the original source is unavaible, there's still an intermediary version available.

## ORIGINAL DATA

# expression data containing fpkm expression values. fpkm is fragements per kilobase million. fragments means that is is for paired-end data.
options(stringsAsFactors = FALSE)
ExpressionData = read.table('/media/casper/USB_ccpeters/internship_thesis/data/rezvan_lncrna2function/lncRNA2function_data.txt', sep='\t', head=TRUE)

# Gene ID's and annotation term type.
EnsemblID2GOID = read.table('/media/casper/USB_ccpeters/internship_thesis/data/rezvan_lncrna2function/EG2GO.txt', sep='\t',head=TRUE)

# expression data selection <--- THIS IS WHERE EXAMPLE DATA IS MADE
# ENSG00000228630
# rowsample <- sample(1:nrow(ExpressionData),10000)
# ExpressionData <- unique(rbind(ExpressionData[rowsample,], 
#   ExpressionData[ExpressionData$GeneID=="ENSG00000228630",]))
# EnsemblID2GOID <- EnsemblID2GOID[EnsemblID2GOID$ensembl_gene_id %in% 
#   ExpressionData$GeneID,]

## Conversion to expressionset

expression_matrix <- as.matrix(ExpressionData[,4:ncol(ExpressionData)])
rownames(expression_matrix) <- ExpressionData$GeneID
featuredat <- as.data.frame(ExpressionData[,1:3])
rownames(featuredat) <- ExpressionData$GeneID
expset <- ExpressionSet(expression_matrix, 
                        featureData = new("AnnotatedDataFrame", 
                                          data=featuredat))

## selection of filter

# filter out everything that IS NOT a protein coding gene
filter_vector <- expset@featureData@data[(expset@featureData@data$
                                            GeneType=="protein_coding"),]$GeneID

## make translational data.frame (BP only)
bp_only <- EnsemblID2GOID[(EnsemblID2GOID[ ,3] == "biological_process"), ]
id_translation_df <- bp_only[,1:2]
colnames(id_translation_df) <- c("ORIGID","GO")

## algorithm run
gid <- "ENSG00000228630"

result <- GAPGOM::expression_prediction(gid, 
                                        expset, 
                                        filter_vector, 
                                        "human", 
                                        "BP",
                                        id_translation_df = id_translation_df,
                                        method = "combine", verbose = TRUE, 
                                        filter_pvals = TRUE
)
result