Description Usage Arguments Value Examples
View source: R/PeakFunctions_multi_biGaussian.R
Calculate Bi-Gaussian Distribution(s) or Bi-Gaussian Peak(s)
1 2 3 4 5 6 7 8 9 | multi_biGaussian(
x,
mus,
sigma1s,
sigma2s,
probDensity,
returnComponentPks = FALSE,
ks
)
|
x |
A vector of x-coordinates from which the corrisponding y-coordinates are calculated |
mus |
A vector of mean values |
sigma1s |
A vector of standard deviations for the left side of the peaks |
sigma2s |
A vector of standard deviations for the right side of the peaks |
probDensity |
Should the function produce probability distributions 'TRUE' with integrated area from min(x) to max(x) equal to 1 or produce gaussian peaks 'FALSE' with amplitudes ks? default is 'TRUE'. |
returnComponentPks |
Should the function return a single vector containing the sum of each individual peak 'FALSE' or a data.frame containing the input vector 'x', each component peak 'peak_n', and the summed result of the component peaks 'peak_sum'. The default is FALSE. |
ks |
A vector of amplitudes for the gaussian peaks. Only used when 'probDensity == FALSE' |
returns either a single vector of y-coordinates the same length as x or a data.frame
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | #Exponentially Modified BiGaussian Probability Distribution from multiple peaks
xVec <- seq(from = 1, to = 100, by = 0.1)
pdist <- multi_biGaussian(x = xVec,
mus = c(5, 10, 15),
sigma1s = c(1, 2, 4),
sigma2s = c(2,2,2),
probDensity = TRUE,
returnComponentPks = FALSE)
pdist_components <- multi_biGaussian(x = xVec,
mus = c(5, 10, 15),
sigma1s = c(1, 2, 4),
sigma2s = c(2,2,2),
probDensity = TRUE,
returnComponentPks = TRUE)
plot(x = xVec, y = pdist)
plot(x = pdist_components$x, y = pdist_components$peak_sum)
points(x = pdist_components$x, y = pdist_components$peak_1, col = "red")
points(x = pdist_components$x, y = pdist_components$peak_2, col = "blue")
points(x = pdist_components$x, y = pdist_components$peak_3, col = "green")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.