Learn R Programming

nlr (version 0.1-3)

loss.robchis: Hetroscedastic chi-square robust loss function.

Description

Compute the robustified chi-square pseoudo log likelihood,based on varaince model function. $$\ell (\boldsymbol{\theta},\tau )=\sum\{ w_i \log (H({x_i};\tau) + \rho \left[\sqrt {{z_i}/ H({x_i};{\sigma ^2},\lambda )} \right] \}$$

Usage

loss.robchis(formula, data, start, theta, varmodel, robfunc, ...)

Arguments

formula

nl.form object of nonlinear function \(f(x;\theta)\).

data

list of data include responce, predictor or possibly predictor of variance model function (t), if not represented then the predict of nonlinear model function will be replaced in predictor variable of nonlinear variance model function that is \(Var(\varepsilon)=\sigma^2 H(f(\theta),\tau)\)

start

list of parameter values of variance model function (\(\tau\) in \(H(t,\tau)\)), initial value or increament during optimization procedure.

theta

list of model function parameter (\(\theta\) in \(f(x,\theta)\)).

varmodel

nl.form object of heteroscedastic variance model function. varmodel must return variance function not standard deviation. Include gradient and hessian.

robfunc

nl.form of rho function. It must include tuning constants k0 and k1.

extra argument might pass to nonlinear regression, heteroscedastic functions, or tuning constant and other parameters to robust loss rho function.

Value

list od los function values:

value

value of minus loglikelihood of chi-square, include attribute "gradient"" and "hessian". These values use in optimization functions.

angvec

angular vector for checking the convergence.

angmat

angular matrix for checking convergence in optimization procedure.

refvar

refvar, sample variance \( \sum(wi * z) \over \sum(wi)\)

fmod

computed function model \(f(x,\theta)\), include response, predictor and their gradient and hessian depends on the defined form of nonlinear function model.

varcomp

computed variance function model \(H(t,\tau)\), include response or predictor and their gradient and hessian depends on the defined form of nonlinear function model.

vcmdata

list of data used in variance model function, that is varmodel$independent and varmodel$dependent typically is zi.

sourcefnc

source function from which this function is called. May be used in feature computing such as outlier detection measures.

rho

computed robust loss rho function. Include hessian and gradients.

zi

computed sample variance, which follows the chi-square distribution.

Details

For estimating variance model parameter \(\tau\), robustified form of chi-square pseudo likelihood is used as robust estimate. Based on calculating the sample variances.

References

Riazoshams H 2010 Outlier detection and robust estimation methods for nonlinear regression having autocorrelated and heteroscedastic errors.

See Also

nl.form, loss.chis

Examples

Run this code
# NOT RUN {
## The function is currently defined as
"loss.robchis"
# }

Run the code above in your browser using DataLab