ERE_step_cycle: Does one Gibbs Step on a cycle
| ERE_step_cycle | R Documentation |
Does one Gibbs Step on a cycle
Description
Execute one Gibbs step on a cycle keeping row and column sums fixed
Usage
ERE_step_cycle(r, c, L, lambda, p, eps = 1e-10)
Arguments
r |
Row indies of cycle, starting at 0 (vector of length k) |
c |
Column indices of cycle, starting at 0 (vector of length k) |
L |
nxn matrix with nonnegative values (will be modified) |
lambda |
nxn matrix of intensities |
p |
nxn matrix of probabilities (must be in [0,1] and 0 on diagonal) |
eps |
Threshold for values to be interpreted as equal to 0 (default = 1e-10) |
Value
no return value
Examples
L=matrix(rexp(9),nrow=3)
lambda <- matrix(0.5,nrow=3,ncol=3)
p <- matrix(0.7, nrow=3,ncol=3)
ERE_step_cycle(r=c(0,1),c=c(1,2),L=L,lambda=lambda,p=p)
ERE_step_cycle(r=c(0,1,2),c=c(0,1,2),L=L,lambda=lambda,p=p)
ERE_step_cycle(r=c(0,1,2),c=c(2,1,0),L=L,lambda=lambda,p=p)
