Dcoef {sojourner} | R Documentation |
Caclulate diffusion coefficient (Dcoef) for trajecotries.
Dcoef( MSD=NULL,trackll=NULL,dt=6,filter=c(min=7,max=Inf),rsquare=0.8, resolution=0.107,binwidth=NULL,method=c("static","percentage", "rolling.window"),plot=FALSE,output=FALSE,t.interval=0.01, profile=NULL)
MSD |
Mean Square Displacement calculated using msd() function. Either MSD or trackll can be passed into Dcoef for calculation of diffusion coefficient. |
trackll |
Track list output from readDiatrack(). |
dt |
Time intervals. Default 6. |
filter |
a vector specifies the minimum and max length of trajecotries to be analyzed. Take only trajectories that has number of frames greater than (>=) min and less than (<) max. |
rsquare |
rsquare filter on Dcoef results. Default to be 0.8. Set value to 0 if rsquare filter is not desired. |
resolution |
ratio of pixel to uM. |
binwidth |
binwidth used for histogram. Default NULL, automatically assign binwidth. |
method |
"static", uses time lags 2~5 to calculate diffusion coefficient; "percentage", uses (tierd) percentage (default 0.25) of time lags (see Details). "rolling.window", time lags uses for Dcoef follows a rolling window with specified window size (default 4). |
plot |
A parameter for plotting. Default FALSE, no plot; If TRUE, automatically plots 'histogram' with count information, binwidth can be set through parameter binwidth; as well as 'density' with density/frequency. |
output |
An Logical indicate if output should be generated. See Values for detail. |
t.interval |
time interval between frames, default 0.010 s (10ms). |
profile |
Location of preference file. By default (NULL), it is stored at : system.file("extdata","PREF","profile.csv",package="sojourner"). User can provide preference file by specifying the location of the file, e.g. profile="/Users/shengliu/Desktop/profile.csv". |
Generic parameters (parameter applied to all methods, such as resolution etc) are set in the function. Method dependent parameters (such as lag.start, lag.end for method = "static"), are stored in profile.csv in PREF folder under extdata. To change preference parameter, can either programably or manually go to folder system.file("extdata","PREF","profile.csv",package="sojourner"), and c hange the profile.csv.
lag.start: time lag used as start of dt for compute Dcoef. Default 2.
lag.end: Time lag used as end of dt for compute Dcoef. Default 2.
method for calculating Dcoef:
static stabilize the number of time lags used for fitting using time lag 2~ 5 despite the total time lags measured.
percentage "percentage", uses (tierd) percentage (default 0.25) of time lags.
[,1] | TrackLength | Percentage | TimeLagsForFitting |
[,2] | 31~ | 0.25 | 2~5-2~7 |
[,3] | 22~30 | 0.25 | 2~5-2~7 |
[,4] | 15~21 | 0.4 | 2~5-2~7 |
[,5] | 10~15 | 0.6 | 2~5-2~7 |
[,6] | 7~9 | 1 | 2~5-2~7 |
rolling.window time lags uses for Dcoef follows a rolling window with specified window size (default 4).
Dcoef A list of Dcoef for each file in trackll.
PDF Log.Dcoef histogram fitted with density curve, when plot = TRUE.
csv Dcoef output in csv format, when output = TRUE.
# compare files folder=system.file("extdata",'SWR1',package="sojourner") trackll = createTrackll(folder=folder, input=3) MSD=msd(trackll=trackll) Dcoef(MSD=MSD,method="static",plot=TRUE) # compare folders folder1=system.file("extdata",'SWR1',package="sojourner") folder2=system.file("extdata",'HTZ1',package="sojourner") trackll2=compareFolder(folders=c(folder1,folder2), input=3) Dcoef(trackll=trackll2,method="percentage",plot=TRUE) Dcoef(trackll=trackll2,method="rolling.window",plot=TRUE)