View source: R/ModelGaussPeak.R
ModelGaussPeak | R Documentation |
ModelGaussPeak
will create a potentially overloaded Gaussian peak of requested width and height.
ModelGaussPeak(
height = 10^7,
width = 4,
scan_rate = 10,
e = 0,
ds = 10^7,
base_line = 10^2
)
height |
True peak height (=intensity counts). |
width |
Peak width in time units (preferably seconds). |
scan_rate |
Is determining the resolution of data points per time unit (preferably seconds). |
e |
Error term giving the percent amount of deviation from the ideal Gaussian curve for individual data points. |
ds |
Detector saturation. Intensity values will be cut off at this point if requested. |
base_line |
Defines if peak is supposed to have a higher base level. |
The main task of ModelGaussPeak
is to create peak data in Gaussian shape for testing.
Width is meant in the chromatographic sense, i.e. the time between peak front and tail hitting the baseline.
Dataframe with columns 'rt' and 'int'.
ylim <- c(0,10^7)
par(mfrow=c(1,5))
pk <- ModelGaussPeak(height=10^7, width=4, scan_rate=10, e=0, ds=10^7, base_line=10^2)
plot(pk,ylim=ylim,main="standard")
pk <- ModelGaussPeak(height=10^7, width=4, scan_rate=10, e=0, ds=8*10^6, base_line=10^2)
plot(pk,ylim=ylim,main="flat top")
pk <- ModelGaussPeak(height=10^7, width=4, scan_rate=10, e=0, ds=8*10^6, base_line=10^5)
plot(pk,ylim=ylim,main="high baseline")
pk <- ModelGaussPeak(height=10^7, width=4, scan_rate=10, e=0.05, ds=8*10^6, base_line=10^5)
plot(pk,ylim=ylim,main="e=5%")
pk <- ModelGaussPeak(height=10^7, width=4, scan_rate=5, e=0.05, ds=8*10^6, base_line=10^5)
plot(pk,ylim=ylim,main="sr=5")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.