R/Ridi3.R

Ridi3 <-
function(X){
. <- "Shut up"
X=as.data.frame(X)
#X=scale(X,scale=F)
entree=apply(X,2,mean)

p=ncol(X)
#res=vector("list",p)
res=entree
res[1]=0.5*entree[1]



for(j in 2:p){
#if (j>1){ res[j]=entree[j]+entree[j-1]}
aux=entree[j-1]
res[j]=0.5*entree[j]+max(cumsum(entree[1:(j-1)]))
}

res=res/(sum(res))

return(list(p=p,entree=entree,Ridit=res))
}

Try the GraphPCA package in your browser

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

GraphPCA documentation built on May 2, 2019, 1:08 p.m.