changePointAnalysis: Get the offset for specific RiboSeq read width

Description Usage Arguments Details Value See Also

View source: R/shift_footprints_helpers.R

Description

Creates sliding windows of transcript normalized counts per position and check which window has most in upstream window vs downstream window. Pick the position with highest absolute value maximum of the window difference. Checks windows with split sites between positions -17 to -7, where 0 is TIS. Normally you expect the shift around -12.

Usage

1
2
3
4
5
6
changePointAnalysis(
  x,
  feature = "start",
  max.pos = 40L,
  interval = seq.int(14L, 24L)
)

Arguments

x

a vector with count per position to analyse, assumes the zero position (TIS) is in the middle + 1 (position 0). Default it is size 60, from -30 to 29 in p-shifting

feature

(character) either "start" or "stop"

max.pos

integer, default 40L, subset x to go from index 1 to max.pos, if tail is not relevant.

interval

integer vector , default seq.int(14L, 24L). Seperation points for upstream and downstream windows. That is (+/- 5 from -12) position.

Details

Transcript normalized means per CDS TIS region, count reads per position, divide that number per position by the total of that transcript, then sum up these numbers per position for all transcripts.

Value

a single numeric offset, -12 would mean p-site is 12 bases upstream

See Also

Other pshifting: detectRibosomeShifts(), shiftFootprintsByExperiment(), shiftFootprints()


ORFik documentation built on March 27, 2021, 6 p.m.