https://github.com/cran/smacof
Tip revision: 6e28d0110a861ff280afacd503d89827fcdc7ae3 authored by Jan de Leeuw on 29 November 2011, 16:49:27 UTC
version 1.2-1
version 1.2-1
Tip revision: 6e28d01
monregS.R
# secondary approach to ties
`monregS` <-
function(x, y, w = rep(1,length(x)), block = weighted.mean)
{
#x ... observed distance matrix
#y ... Guttman transformed distances
#w ... weights
wag <- tapply(w,x,sum)
yag <- tapply(y,x,mean)
xag <- tapply(x,x,mean)
o <- order(xag)
r <- order(o)
e <- pavasmacof(yag[o],wag[o])[r] #call pava
return(ifelse(outer(x,xag,"=="),1,0)%*%e)
}