...
|
...
|
@@ -688,7 +688,6 @@ To create the GTF files, the algorithm uses the *EventPointer_RNASeq_TranRef_IGV
|
688
|
688
|
As EventDetection_transcriptome function returns the splicing graph information, this function does not need to create the splicing graph. (EventPointer_IGV is equivalent function for array platform).
|
689
|
689
|
|
690
|
690
|
```{r ep_tranref_igv, eval=TRUE, warning=FALSE, collapse=TRUE}
|
691
|
|
-
|
692
|
691
|
SG_List <- EventXtrans$SG_List
|
693
|
692
|
PathEventsTxt<-system.file('extdata',package='EventPointer')
|
694
|
693
|
PathEventsTxt <- paste0(PathEventsTxt,"/EventsFound_Gencode24_2genes.txt")
|
...
|
...
|
@@ -751,7 +750,7 @@ Internally, given the TxD matrix *Protein_Domain_Enrichment* relates protein dom
|
751
|
750
|
**NOTE**: Check that EventXtrans$transcritnames annotation type match to the rownames of the TxD matrix.
|
752
|
751
|
|
753
|
752
|
```{r pd_enrichmetn, eval=TRUE, warning=FALSE, collapse=TRUE}
|
754
|
|
-
|
|
753
|
+data(EventXtrans)
|
755
|
754
|
data("TxD")
|
756
|
755
|
|
757
|
756
|
#check same annotation for transcripts:
|
...
|
...
|
@@ -759,7 +758,7 @@ EventXtrans$transcritnames[1]
|
759
|
758
|
rownames(TxD)[1]
|
760
|
759
|
|
761
|
760
|
|
762
|
|
-## as si not the same, we change EventXtrans$transcritnames annotation
|
|
761
|
+## If they are different, we need to change EventXtrans$transcritnames annotation
|
763
|
762
|
transcriptnames <- EventXtrans$transcritnames
|
764
|
763
|
transcriptnames <- gsub("\\..*","",transcriptnames)
|
765
|
764
|
EventXtrans$transcritnames <- transcriptnames
|