Using the Scripts
In this section you will learn how to run an analysis using the scripts.
1) load your signal in the workspace
FORMAT of input signal: the input signal need to be in a double of <nx1> ionic current values in nA.
Your input signal MUST contain baseline information, at least 100 points of baseline before and after each event.
2) Check your signal by plotting it
Note that the baseline shouldn't contain any sudden jump...
If you have a problem loading the signal please email us! See address on website: http://lben.epfl.ch/page-79460-en.html
2) Open main.m and fill in the input parameters
line 5: give a name to the folder in which your results will be stored
line 6: input the name of the signal you just loaded in the Workspace
line 7: input the value of delta: most likely current drop in nA that you want to detect in an efficient manner (in the demo case 0.2 nA is a typical current blockage value for λ DNA translocation)
line 8: input the value of sigma: RMS current noise value, can be calculated using std() function of MATLAB on 1000 points of baseline or more
line 9: input the sampling frequency used when recording the signal in Hz
3) run the main
You can either press the button shown above or write main in the Command Window.
4) look at the results
A .mat file called ExperimentData is created containing all results of an analysis:
EventDatabase, ConcatenatedEvents and ConcatenatedFits and the settings of the analysis.
EventDatabase is a structure array with level and dwell time information for each event (see Event Database for more information) ConcatenatedEvents is a concatenation of all the detected events with 100 points of baseline between each event, ConcatenatedFits are its corresponding fits.