######################################################################### ### PhD course, Advanced Computational Statistics ### ### Code for Problem 4.3 from lecture slides (Lecture 4) ### ### Frank Miller, 2023-04-18 ### ######################################################################### w <- seq(-1, 1, by=0.05) # design space lw <- length(w) des <- rep(0, lw) indices <- 1:lw for (i in 1:n){ ind <- sample(indices, size=1) des[ind] <- des[ind]+1 } irem <- sample(indices[des>0], size=1) iadd <- sample(indices, size=1) # this might be improved by excluding iadd=irem desnew <- des desnew[irem] <- desnew[irem]-1 desnew[iadd] <- desnew[iadd]+1