Nothing
"pcnm" <-
function(matdist,thresh=give.thresh(as.dist(matdist)))
{
matdist <- as.matrix(matdist)
mattrunc <- ifelse(matdist >thresh, 4*thresh,matdist)
wa.old <- options()$warn
options(warn = -1)
mypcnm <- cmdscale(mattrunc,k=min(dim(matdist)) -1 ,eig=TRUE)
## now cmdscale returns only positive vectors but all eigenvalues...
## but it could returns eigenvector corresponding to null eigenvalue
mypcnm$eig <- mypcnm$eig[1:ncol(mypcnm$points)]
eq0 <- apply(as.matrix(mypcnm$eig/max((mypcnm$eig))),1,function(x) identical(all.equal(x, 0), TRUE))
res <- list()
res$values <- mypcnm$eig[!eq0]
res$vectors <- mypcnm$points[,!eq0]
res$vectors <- sweep(res$vectors,2,sqrt(res$values),"/")
options(warn = wa.old)
return(res)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.