weighted loss function is used to estimate Weighted M-estimate. It is a robustified form of Likelihood function for heteroscedastic variance case. loss.hetroWM will be used in optimization function for estimating parameters of nonlinear function model and variance function model simultaneously.
loss.hetroWM(formula, data, start, varmodel, robfunc, ...)nl.form object of nonlinear regression model.
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)\)
list of parameter values of nonlinear model function (\(\theta\) in \(f(x,\theta)\) and \(\tau\) in \(H(x,\tau)\)), initial values or increament during optimization procedure.
nl.form object of heteroscedastic variance model function. varmodel must return variance function not standard deviation. Include gradient and hessian.
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.
result <- list (value=value,angvec=angvec,angmat=angmat, refvar=refvar,sourcefnc= match.call(), rho=robvalue,fmod=fmod,varcomp=varcomp,correlation =nlrho,ri=ri)
Robustified form of log-likelihood function, use in optimize function. Include gradient and hessian attributes.
angular vector for checking the convergence.
angular matrix for checking convergence in optimization procedure.
refvar, sample variance \( \sum(wi * z) \over \sum(wi)\)
Object of class "callorNULL" source function called for fitt.
computed robust loss rho function. Include hessian and gradients.
computed function contains esponse and or its gradient and hessian predictor and or its gradient & hessian
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.
correlation of fitt.
residuals
Loss function in general form have robust rho function. l(tau,theta)=sum[log(h(mu;tau,sg))] + sum[rho(ri/h(mu;tau,sg))]
Lim, C., Sen, P. K., Peddada, S. D. (2010). Statistical inference in nonlinear regression under heteroscedasticity. Sankhya B 72:202-218.
# NOT RUN {
## The function is currently defined as
"loss.hetroWM"
# }
Run the code above in your browser using DataLab