changePointAnalysis {ORFik} | R Documentation

## Get the offset for specific RiboSeq read width

### 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

    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()`

[Package ORFik version 1.8.6 Index]