Traffic Intensity Estimation in Finite Markovian Queueing Systems

Listing 6

Monte Carlo simulation.

MtCaRoMM1K<−function(ssize, rho, K, fEst,...)

rep<−10000

samp<−numeric(ssize)

est<−numeric(rep)

for (i in 1:rep)

samp<−rmm1k(ssize, rho, K)

oldseed<−.GlobalEnv$.Random.seed

est [i]<−fEst (K, samp,…)

.GlobalEnv$.Random.seed<−oldseed

return(c(mean(est), var(est)))