shapiroAcceptor: shapiroAcceptor

Description Usage Arguments Value Author(s) See Also Examples

View source: R/main_splinter.R

Description

Shapiro's score of acceptor site (range is from -13 [intron] to +1 [exon]) is: 100 * ((t1 - l1)/(h1 - l1) + (t2 - l2)/(h2 - l2))/2, where t1 is the sum of the best 8 of 10 percentages at positions -13 to -4, l1 is the sum of the lowest 8 of 10 percentages at position -13 to -4, h1 is the sum of the highest 8 of 10 percentages at positions -13 to -4, t2 is the sum of percentages at positions -3 to +1, l2 is the sum of the lowest percentages at positions -3 to +1, and h2 is the sum of the highest percentages at positions -3 to +1

Usage

1
shapiroAcceptor(reference_fasta, target_fasta)

Arguments

reference_fasta

vector of strings or DNAStringSet of reference splice list

target_fasta

vector of strings or DNAStringSet of fasta to score

Value

data.frame with Shapiro scores

Author(s)

Diana Low

See Also

http://www.softberry.com/spldb/SpliceDB.html

Examples

1
2
3
4
5
6
7
library(BSgenome.Mmusculus.UCSC.mm9)
bsgenome <- BSgenome.Mmusculus.UCSC.mm9
data_path<-system.file("extdata",package="SPLINTER")
splice_data<-extractSpliceEvents(data=paste(data_path,"/skipped_exons.txt",sep=""))
splice_sites<-extractSpliceSites(splice_data,site="acceptor")
acceptor.ss<-getSeq(bsgenome,splice_sites)
##sacceptor<-shapiroAcceptor(acceptor.m,acceptor.ss)

SPLINTER documentation built on Nov. 8, 2020, 8:13 p.m.