# R/wygalg.R In robeth: R Functions for Robust Statistics

```"wygalg" <-
function(x,a,exu=ucv,exup=upcv,nobs=nrow(x),maxit=.dFvGet()\$mxg,
nitmon=.dFvGet()\$ntm,icnv=.dFvGet()\$icv,tol=.dFvGet()\$tlo,xfud=.dFvGet()\$xfd) {
if (missing(x)) messagena("x")
nvar <- ncol(x)
ncov <- length(a)
mdx <- nrow(x)
mdz <- nobs
if (missing(a)) a <- double(ncov)
nit <- integer(1)
dist <- single(nobs)
sa <- double(ncov)
ss <- double(ncov)
sz <- matrix(single(1),mdz,nvar)
su <- double(nobs)
sup <- double(nobs)
sy1 <- double(nvar)
sy2 <- double(nvar)
st <- double(ncov)
st1 <- double(ncov)
f.res <- .Fortran("int55",
x=to.single(x),
a=as.double(a),
as.integer(exu()),
as.integer(exup()),
nobs=to.integer(nobs),
nvar=to.integer(nvar),
ncov=to.integer(ncov),
mdx=to.integer(mdx),
mdz=to.integer(mdz),
maxit=to.integer(maxit),
nitmon=to.integer(nitmon),
icnv=to.integer(icnv),
tol=to.single(tol),
xfud=to.single(xfud),
nit=to.integer(nit),
dist=to.single(dist),
sa=as.double(sa),
ss=as.double(ss),
sz=to.single(sz),
su=as.double(su),
sup=as.double(sup),
sy1=as.double(sy1),
sy2=as.double(sy2),
st=as.double(st),
st1=as.double(st1))
list(a=f.res\$a,nit=f.res\$nit,dist=f.res\$dist,su=f.res\$su)
}
```

## Try the robeth package in your browser

Any scripts or data that you put into this service are public.

robeth documentation built on June 13, 2018, 5:04 p.m.