View source: R/lib_ErdosRenyExponential_MCMC.R
steps_ERE | R Documentation |
Runs a Gibbs sampler in the Erdos Reny model with Exponential weights (ERE model) and fixed marginals. The algorithm starts from a given matrix.
steps_ERE(L, p, lambda, nsamples = 10000, thin = 1000, burnin = 10000)
L |
Starting matrix for the Gibbs sampler. Implicitly defines the fixed marginals. |
p |
A matrix with entries in [0,1] |
lambda |
A matrix with nonnegative entries |
nsamples |
Number of samples to return. |
thin |
Frequency at which samples should be generated (default=1, every step) |
burnin |
Number of initial steps to discard. |
List of simulation results
sample_ERE
L <- matrix(rexp(4*4),nrow=4,ncol=4); diag(L)=0;
p <- matrix(0.5,nrow=4,ncol=4); diag(p) <-0;
lambda <- matrix(1,nrow=4,ncol=4); diag(lambda)<-0;
L <- steps_ERE(L=L,p=p,lambda=lambda,nsamples=5,thin=50,burnin=20)
L
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.