FKPnets.R
transMatrix.R
samp-matrices.R

# Playing around in general with matrices and their eigenspectrum:
> k <- play()$k
> j <- play()$j
> h <- play()$h
> eigen(h)
> egh <- eigen(h)$vec[,1]
> h %*% egh

# Plotting FKP nets and coloring nodes according to their positive or negative sign in eigenvector:

# vec2

> foo <- runFKP(N=100,alpha=1000000)
> mm <- transMatrix(outLink=foo$outLink)
> vec2 <- eigen(mm)$vec[,2]
> pls <- c(1:100)[vec2 > 0]
> mns <- c(1:100)[vec2 < 0]
> dev.set(2)
> points(foo$coords[pls,],pch="X",col="black")
> points(foo$coords[mns,],pch="M",col="green")

#vec4
> foo <- runFKP(N=100,alpha=1000000)
> mm <- transMatrix(outLink=foo$outLink)
> vec4 <- eigen(mm)$vec[,4]
> pls <- c(1:100)[vec4 > 0]
> mns <- c(1:100)[vec4 < 0]
> dev.set(2)
> points(foo$coords[pls,],pch="X",col="black")
> points(foo$coords[mns,],pch="M",col="green")