Nothing
pois.rel.pca <-
function(x,lambda.min,lambda.max,len=10,plot=TRUE,seed=132)
{
if(!is.list(x)) {x <- lapply(seq_len(nrow(x)), function(i) x[i,])}
if(len > length(x)) stop("len exceeds row-length of data entered: For PCA number of columns cannot exceed number of rows.")
# Constructing the Distance matrix
mat.km <- matrix(NA,nrow = length(x),ncol = len)
rownames(mat.km) <- names(x)
colnames(mat.km) <- names(x)
for (i in 1:length(x))
{
y <- x[[i]]
mat.km[i,] <- wt.rel.pois(y,lambda.min,lambda.max,plot=FALSE,len=len)$Val
colnames(mat.km) <- wt.rel.pois(y,lambda.min,lambda.max,plot=FALSE,len=len)$Lambda
}
#---------------------------------------
# PCA
#---------------------------------------
lilkd.princomp <- princomp(mat.km)
if(plot==TRUE)biplot(lilkd.princomp,var.axes=FALSE,cex=0.7,col=c(1,0))
out <- list("PCA.Output" = lilkd.princomp)
return(out)
}
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.