Description Slots Features Functions Author(s) See Also Examples
This S4 class represents a collection of Targeted Sequencing Experiments in R. All these experiments are characterized by a 'bed file' containing the specification of the explored 'features', as a 'feature panel'. These features could be amplicons, exons, transcripts, among others. In general each feature is associated to one gene but a gene could be related to many features. This class allows the representation and quality control of a set of Targeted Sequencing Experiment made over the same or different subjects but using always the same bed file'.
bedFile
GRanges object that models the bed file.
panels
GRanges object containing the feature/gene panels.
attribute
character indicates which attribute, 'coverage' or 'medianCounts' will be used to the analysis.
feature
character indicates the name of the analyzed features. E.g 'amplicon', 'exon', 'transcript', 'gene'.
Model sets of targeted sequencing experiments in R.
Evaluate the performance of the targeted sequencing technique across several experiments using coverage/read counts information.
Detect in early stage sequencing or library preparation errors.
Report quality control results.
TargetExperimentList S4 class includes the following functions:
constructor of TargetExperimentList to generate the feature panel starting from at least two TargetExperiment objects
return the respective TargetExperimentList slots
set the respective TargetExperimentList slot
generic output of the object
generic output of the object
print statistics summary for the set attribute
plot a summarized view of the attribute values achieved by each feature in each sample
plot the attribute distribution for each feature
plot the attribute distribution in each panel
plot the attribute distribution in each or pool
Gabriela A. Merino gmerino@bdmg.com.ar, Cristobal Fresno cfresno@bdmg.com.ar, Yanina Murua ymurua@leloir.org.ar, Andrea S. Llera allera@leloir.org.ar and Elmer A. Fernandez efernandez@bdmg.com.ar
Rsamtools
Other TargetExperimentList: TargetExperimentList
,
initialize,TargetExperimentList-method
,
object
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 | # Defining the set of TargetExperiment objects
data(ampliPanel, package="TarSeqQC")
data(ampliPanel2, package="TarSeqQC")
ampliList<-list(ampliPanel, ampliPanel2)
# Defining feature parameter
feature<-"amplicon"
# Defining attribute parameter
attribute<-"coverage"
##Calling the constructor
object<-TargetExperimentList(TEList=ampliList, attribute=attribute,
feature=feature)
setFeature(object)<-"amplicon"
## load the example dataset
data(TEList, package="TarSeqQC")
## Early exploration
# show/print
TEList
# summary
summary(TEList)
## Controlling low counts features
# Definition of the interval extreme values
attributeThres<-c(0,1,50,200,500, Inf)
# Do a frequency table for the attribute intervals
summaryIntervals(TEList, attributeThres)
# getting low counts features at gene level
getLowCtsFeatures(TEList, level="gene", threshold=50)
# exploring panel performance along several samples
g<-plot(TEList, attributeThres=attributeThres, featureLabs =TRUE)
if(interactive()){
g
}
g<-plotGlobalAttrExpl(TEList,log=FALSE)
# x11(type="cairo")
if(interactive()){
g
}
g<-plotPoolPerformance(TEList,log=FALSE)
if(interactive()){
g
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.