fitpoly1 {tcplfit2}R Documentation

Polynomial 1 (Linear) Model Fit

Description

Function that fits to f(x) = a*x and returns generic model outputs.

Usage

fitpoly1(conc, resp, bidirectional = TRUE, verbose = FALSE, nofit = FALSE)

Arguments

conc

Vector of concentration values NOT in log units.

resp

Vector of corresponding responses.

bidirectional

If TRUE, model can be positive or negative; if FALSE, it will be positive only.

verbose

If TRUE, gives optimization and hessian inversion details.

nofit

If nofit = TRUE, returns formatted output filled with missing values.

Details

Zero background and increasing absolute response are assumed. Parameters are "a" (y scale) and error term "er". success = 1 for a successful fit, 0 if optimization failed, and NA if nofit = TRUE. cov = 1 for a successful hessian inversion, 0 if it fails, and NA if nofit = TRUE. aic, rme, modl, parameters, and parameter sds are set to NA in case of nofit or failure.

Value

Named list containing: success, aic (Akaike Information Criteria), cov (success of covariance calculation), rme (root mean square error), modl (vector of model values at given concentrations), parameters values, parameter sd (standard deviation) estimates, pars (vector of parameter names), sds (vector of parameter sd names).

Examples

fitpoly1(c(.03,.1,.3,1,3,10,30,100), c(0,.01,.1, .1, .2, .5, 2, 5))

[Package tcplfit2 version 0.1.6 Index]