compare_RT_CDF-methods: Compare Residence time/Survival Curve (1-CDF)

Description Usage Arguments Details Value Examples

Description

Compare Residence time/Survival Curve of multiple trackll. Or simply plot the survival curve of one trackll.

Usage

1
2
compare_RT_CDF(trackll=NULL,x.max=30,
filter=c(min=3,max=Inf),t.interval=0.5,output=FALSE)

Arguments

trackll

trajectory list generated by createTrackll() and processing. if NULL, user will be prompted to enter the trackll name.

x.max

The maximum range of X axis, i.e. time, for the output plot. Default 30 sec.

filter

Filter the tracks by step/frame number (length). Only tracks pass through filter will be selected.

t.interval

time interval for image aquisition. Default 0.5 sec.

output

An Logical indicate if output should be generated. See Values for detail.

Details

Compare Residence time/Survival Curve of multiple track list (trackll). Or simply plot the survival curve of one trackll. The survival curve/probability is calculated as 1-CDF of the length of tracks/trajectories.

If the acquisition time interval of the tracklls are different, set argument trackll=NULL, users will be prompted to input the number of the track list (trackll) to compare/plot. Then users will be prompted to input the name and acquisition time interval of each trackll. The trackll should be masked and merged. The maximum time range to be plotted can be set using x.max, this will not change the curve/probabiltiy, which is determined by all track information in the trackll.

Value

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
# Generate trackll, and process, 
# e.g. mask region of interest, merge tracks from multiple files.
folder1=system.file('extdata','HSF',package='sojourner')
trackll1=createTrackll(folder1,input=3, cores = 2)
trackll1=maskTracks(folder1,trackll1)
trackll1=mergeTracks(folder1,trackll1)

folder2=system.file('extdata','HSF_2',package='sojourner')
trackll2=createTrackll(folder2,input=2, cores = 2)
trackll2=maskTracks(folder2,trackll2)
trackll2=mergeTracks(folder2,trackll2)

# Plot and output the survival curve,
compare_RT_CDF(trackll=c(trackll1,trackll2),x.max=30,
filter=c(min=3,max=Inf),t.interval=0.5,output=FALSE)

sheng-liu/sojourner documentation built on April 2, 2021, 2:30 a.m.