REIDS (version 0.1.0)

REIDS_JunctionAssesment: REIDS_JunctionAssessment

Description

The REIDS_JunctionAssessment functions assess identified AS exons based on their 5'end and 3'end and exclusion junction support.

Usage

REIDS_JunctionAssesment(Indices, DataFile, ASProbeSets = c(),
  Juninfo = "User", JAnnotI = NULL, JAnnot = NULL, EandTrAnnotI = NULL,
  EandTrAnnot = NULL, PartiallyAnnotated = FALSE, positionData = NULL,
  transcriptData = NULL, Groups = list(c(3, 4, 5), c(6, 7, 8)),
  Low_AllSamples = c(), Low_GSamples = c(), Location = NULL,
  Name = "REIDS_Jun")

Arguments

Indices

The .csv file created by Line_Indexer.py which contains indices for every gene in geneIDs.

DataFile

The .csv file created by PivotTransformation.

ASProbeSets

The AS probe sets as identified by ASExons.

Juninfo

A parameter specifying wether the annotations are user of Ensembl defined. If JunInfo is "User" the annotations provided in EandTrAnnot are used. If JunInfo is "Ensembl" the annotations in EandTrAnnot are used to set up tje junction associations but the gene name and position in transcriptData and positionData are used to connect with the Ensembl data base and retrieve corresponding information.

JAnnotI

The file name with line indices for the junction associations.

JAnnot

The file name with the junction associations.

EandTrAnnotI

The file name with line indices for the exon and isoform annotations.

EandTrAnnot

The file name with the exon and isoform annotations.

PartiallyAnnotated

Logical. Should the exon annotations with partially annotated probe sets still be included? If FALSE, these are excluded. If TRUE, these are included. Default is FALSE.

positionData

The file with the chromosome start and ends for the probe sets. Only needed in JunInfo=Ensembl.

transcriptData

The file with gene name of the transcripts. Only needed in JunInfo=Ensembl.

Groups

A list with elements speficifing the columns of the data in each group.

Low_AllSamples

A character vector containing the probe sets which are not DABG in all samples.

Low_GSamples

A list with a character vector per group containing the probe sets which are not DABG in that group.

Location

A character string indication the place where the outputs are saved.

Name

A character string with the name of the ouput file. Defaults to "REIDS_Jun".

Value

The function returns four files. The first file has name "Name_ASInfo.txt" and contains a line per probe set. It shows the reached decisionregarding the probe set (Const/AS/not DABG),its linking and exclusion junctions, the fold change, the AS type and its annotated exons. The secondfile, "Name_Compositions.txt", is a list of all found transcripts for a particular TC ID. The third file,"Name_GroupTranscripts.txt" indicates whether a specific transcript is present or absent in a group. The fourth file "Name_NovelConnections.txt" contains junctions which are showing an undocumented connection between probe sets.

Examples

# NOT RUN {
data(TC1500264)
PivotTransformData(Data=TC1500264,GeneID=NULL,ExonID=NULL,
REMAPSplitFile="TC1500264_Gene_SplitFile.txt",Location="Output/",Name="TC1500264_Pivot")
REIDSFunction(ASPSR=c(), Indices="Output/TC1500264_LineIndex.csv",
DataFile="Output/TC1500264_Pivot.csv",nsim=50,informativeCalls=FALSE,Summarize=
c("WeightedAll","EqualAll"),rho=0.5,Low_AllSamples=c(),Groups=list(c(1:3),c(4:6)),
Location="Output",Name="TC1500264")

TC1500264_1vs2=ASExons(ExonScores="Output/TC1500264_ExonScores.txt",ArrayScores=
"Output/TC1500264_ArrayScores.txt",Exonthreshold=0.5,Groups=list(c(1:3),c(4:6)),paired=FALSE,
significancelevel=0.05)

ASPSR_PSR=TC1500264_1vs2[which(round(TC1500264_1vs2$ExonScore,2)>=0.50&
TC1500264_1vs2$adj.p.value<0.05),2]

REIDS_JunctionAssesment(Indices="Output/TC1500264_LineIndex.csv",DataFile=
"Output/TC1500264_Pivot.csv",ASProbeSets=ASPSR_PSR,Juninfo="User",JAnnotI=NULL,
JAnnot=NULL,EandTrAnnotI="Output/REMAP_Indices.txt",EandTrAnnot=
"Output/HJAY_REMAP.txt",positionData=NULL,transcriptData=NULL,Groups=
list(c(1:3),c(4:6)),Low_AllSamples=c(),Low_GSamples=c(),Location="Output",
Name="TC1500264")
# }