selectBonds: Function to select bonds based on colors and parameters

In ick003/SpTMixture: Mixture models with spatio-temporal dependencies

Description

Function to select bonds based on colors and parameters

Usage

selectBonds(Bds, Cols, Betas)

Arguments

Bds	matrix of neighbouring bonds.
Cols	vector of colors per vertex.
Betas	vector of parameters for the Potts model.

Examples

Obs.loc = expand.grid(1:10,1:10)
Bds = getBonds(Obs.loc)
Bds = Bds[which(Bds[,1]>0),]
Betas = c(0.8, 0.8, 0.2)
Cols = sample(1:3, 100, replace = TRUE)
Bd = selectBonds(Bds, Cols, Betas)
par(mfrow = c(1,1))
col = grey.colors(3)
image(matrix(Cols, ncol = 10), col = col)
segments((Obs.loc[Bd[,1],1]-1)/9, (Obs.loc[Bd[,1],2]-1)/9, (Obs.loc[Bd[,2],1]-1)/9, (Obs.loc[Bd[,2],2]-1)/9)
segments((Obs.loc[,1]-1)/9, (Obs.loc[,2]-1)/9, (Obs.loc[,1]-1)/9, (Obs.loc[,2]-1)/9, lwd = 4)

ick003/SpTMixture documentation built on May 18, 2019, 2:32 a.m.