tests/testthat/test_fit_gaussian_optim.R

In pmbrophy/SummitR: Provides a set of tools for peak fitting

context("fit_gaussian_optim")

test_that("fit_gaussian_optim returns correct data types", {
  xVec <- seq(from = 1, to = 100, by = 0.1)
  gauss <- func_gaussian(x = xVec, mu = 10, sigma = 1, probDensity = FALSE, k = 10)

  gauss <- gauss + rnorm(n = length(gauss), mean = 0, sd = 1)

  fit <- fit_gaussian_optim(x = xVec, y = gauss, init_mu = 9, init_sigma = 0.5, init_k = 3, maxit = 1000)

  expect_type(fit, "list")
  expect_named(fit)
  expect_equal(length(fit$fitPeak), length(xVec))
})