Skip to contents

This function is a constructor for S3 class gaussmodel, which represents Gaussian-based model. It usually takes data and optionally freq as arguments and also optionally stepsize. Members of interest in practice are result and coeffs, which maintain the information of estimates and coefficients of polynomials, respectively.

Usage

gaussmodel(data = data, freq = NULL, stepsize = c(0.5, 0.3))

Arguments

data

A numeric vector of a data set to be estimated.

freq

A frequency vector corresponding to the data vector. The default value is NULL, which means all frequencies are one. If supplied, the length of a vector should be same as data and each element should be a nonnegative integer.

stepsize

A numeric vector whose element is larger than 0 and smaller than 1, and decreasing order. The default value is c(0.5, 0.3). If you encounter numerical difficulties, decreasing its values, for example, to c(0.4, 0.2), might help to estimate a model.

Value

An object of Gaussian-based model gaussmodel.

Examples

## Create `gaussmodel` object from a data set `mix2gauss$n200`.
 gmodel <- gaussmodel(mix2gauss$n200)
## Create `gaussmodel` object from a data set `mix2gaussHist$n200p` and
## its frequencies `mix2gaussHist$n200f`.
 gmodel <- gaussmodel(mix2gaussHist$n200p, mix2gaussHist$n200f)