VBLPCM (version 2.4.9)

predict.vblpcm: Find all link probabilities

Description

generate a matrix of link probabilities based on the fitted VB model.

Usage

# S3 method for vblpcm
predict(object, ...)

Value

The posterior predictive link probabilities given the fitted object

Arguments

object

The fitted values; output from vblpcmfit()

...

optional additional arguments.

Author

Michael Salter-Townshend

Examples

data(sampson)
v.fit<-vblpcmfit(vblpcmstart(samplike,G=3))
### create a matrix of link posterior probabilities given the fitted model
probs<-predict.vblpcm(v.fit)
# show this graphically; separation of the boxes implies a good fit to the data
boxplot(split(probs,v.fit$Y),
        ylab=expression(paste("P(",Y[i][j],"=1)")),xlab=expression(paste(Y[i][j])))