# doublesigmoidalFitFunction: Double sigmoidal fit function. In sicegar: Analysis of Single-Cell Viral Growth Curves

## Description

The function fits a double sigmoidal curve to given data by using likelihood maximization (LM) algorithm and provides the parameters (maximum, final asymptote intensity, slope1Param, midpoint1Param, slope2Param, and midpoint distance) describing the double-sigmoidal fit as output. It also contains information about the goodness of fits such as AIC, BIC, residual sum of squares, and log likelihood.

## Usage

 ``` 1 2 3 4 5 6 7 8 9 10 11 12``` ```doublesigmoidalFitFunction( dataInput, tryCounter, startList = list(finalAsymptoteIntensityRatio = 0, maximum = 1, slope1Param = 1, midPoint1Param = 0.33, slope2Param = 1, midPointDistanceParam = 0.29), lowerBounds = c(finalAsymptoteIntensityRatio = 0, maximum = 0.3, slope1Param = 0.01, midPoint1Param = -0.52, slope2Param = 0.01, midPointDistanceParam = 0.04), upperBounds = c(finalAsymptoteIntensityRatio = 1, maximum = 1.5, slope1Param = 180, midPoint1Param = 1.15, slope2Param = 180, midPointDistanceParam = 0.63), min_Factor = 1/2^20, n_iterations = 1000 ) ```

## Arguments

 `dataInput` A data frame or a list contatining the dataframe. The data frame should be composed of at least two columns. One represents time, and the other represents intensity. The data should be normalized with the normalize data function sicegar::normalizeData() before imported into this function. `tryCounter` A counter that shows the number of times the data was fit via maximum likelihood function. `startList` The initial set of parameters vector that algorithm tries for the first fit attempt for the relevant parameters. The vector composes of six elements; 'finalAsymptoteIntensityRatio', 'maximum', 'slope1Param', 'midPoint1Param' , 'slope2Param', and 'midPointDistanceParam'. Detailed explanations of those parameters can be found in vignettes. Defaults are finalAsymptoteIntensityRatio = 0, maximum = 1, slope1Param = 1, midPoint1Param = 0.33, slope2Param = 1, and midPointDistanceParam=0.29. The numbers are in normalized time intensity scale. `lowerBounds` The lower bounds for the randomly generated start parameters. The vector composes of six elements; 'finalAsymptoteIntensityRatio', 'maximum', 'slope1Param', 'midPoint1Param' , 'slope2Param', and 'midPointDistanceParam'. Detailed explanations of those parameters can be found in vignettes. Defaults are finalAsymptoteIntensityRatio = 0, maximum = 0.3, slope1Param = .01, midPoint1Param = -0.52, slope2Param = .01, and midPointDistanceParam = 0.04. The numbers are in normalized time intensity scale. `upperBounds` The upper bounds for the randomly generated start parameters. The vector composes of six elements; 'finalAsymptoteIntensityRatio', 'maximum', 'slope1Param', 'midPoint1Param' , 'slope2Param', and 'midPointDistanceParam'. Detailed explanations of those parameters can be found in vignettes. Defaults are finalAsymptoteIntensityRatio = 1, maximum = 1.5, slope1Param = 180, midPoint1Param = 1.15, slope2Param = 180, and midPointDistanceParam = 0.63. The numbers are in normalized time intensity scale. `min_Factor` Defines the minimum step size used by the fitting algorithm. Default is 1/2^20. `n_iterations` Define maximum number of iterations used by the fitting algorithm. Default is 1000

## Value

Returns the fitted parameters and goodness of fit metrics.

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40``` ```time=seq(3, 24, 0.1) #simulate intensity data and add noise noise_parameter <- 0.2 intensity_noise <- stats::runif(n = length(time), min = 0, max = 1) * noise_parameter intensity <- doublesigmoidalFitFormula(time, finalAsymptoteIntensityRatio = .3, maximum = 4, slope1Param = 1, midPoint1Param = 7, slope2Param = 1, midPointDistanceParam = 8) intensity <- intensity + intensity_noise dataInput <- data.frame(intensity = intensity, time = time) normalizedInput <- normalizeData(dataInput) parameterVector <- doublesigmoidalFitFunction(normalizedInput, tryCounter = 2) #Check the results if(parameterVector\$isThisaFit){ intensityTheoretical <- doublesigmoidalFitFormula( time, finalAsymptoteIntensityRatio = parameterVector\$finalAsymptoteIntensityRatio_Estimate, maximum = parameterVector\$maximum_Estimate, slope1Param = parameterVector\$slope1Param_Estimate, midPoint1Param = parameterVector\$midPoint1Param_Estimate, slope2Param = parameterVector\$slope2Param_Estimate, midPointDistanceParam = parameterVector\$midPointDistanceParam_Estimate) comparisonData <- cbind(dataInput, intensityTheoretical) require(ggplot2) ggplot(comparisonData) + geom_point(aes(x = time, y = intensity)) + geom_line(aes(x = time, y = intensityTheoretical)) + expand_limits(x = 0, y = 0)} if(!parameterVector\$isThisaFit) {print(parameterVector)} ```

sicegar documentation built on May 8, 2021, 9:06 a.m.