Displaying 1 result from an estimated 1 matches for "cus1".
Did you mean:
cus
2004 Sep 04
0
Non-Markovian Behaviour of a Cusum?
...xTime
# Q[i] <- max(0, Q[i-1]+z[i]), z[i] ~ N(offset, 1),
# Q[0] = FIR
# Store only the mean of Q[i] for each i
Qmean <- rep(NA, maxTime)
Q <- rep(FIR, nSims)
for(i in 1:maxTime){
Q <- pmax(0, Q+rnorm(nSims, mean=offset))
Qmean[i] <- mean(Q)
}
Qmean
}
set.seed(321)
Cus1 <- CusumSim()
plot(Cus1)
# Different simulation, essentially the same behavior
Cus2 <- CusumSim()
plot(Cus1, ylim=range(Cus1, Cus2))
lines(Cus2)
# Different random number generator, same behavior
RNGkind("Wichmann-Hill")
CusWH <- CusumSim()
plot(Cus1, ylim=range(Cus1, Cus2, CusW...