R/SpecGap.R

Defines functions SpecGap

Documented in SpecGap

SpecGap <-
function(P)
{
	aux<-eigen(t(P),symmetric = FALSE)
	piE<-aux$vectors[,1]/sum(aux$vectors[,1])
	gap<-1-Mod(aux$values[2])
	return(list(SG=gap, pi=piE))
}

Try the IndTestPP package in your browser

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

IndTestPP documentation built on Aug. 29, 2020, 1:06 a.m.