Displaying 2 results from an estimated 2 matches for "samp_i".
2006 Dec 07
2
Simulation in R
Hello!
I have the following problem.
My code:
--snip--
ergebnisse <- rep(0, each=2)
stichproben <- rep(0, each=2)
for (i in seq(1:2)) {
n <- dim(daten)[1]
ix <- sample(n,200) # producing samples
samp_i <- daten[ix,]
stichproben[i] <- samp_i # doesn???t works
# Calculation of the model:
posterior_i <- MCMClogit(y ~ fbl.ind + fekq3 + febitda4 + fuvs + fkru +
fzd + fur3, data=samp_i, b0=prior, B0=precision, tune=0.5) # calculation
ergebnisse[i] <- summary(posterior_i) # saving...
2006 Dec 07
1
Simulation in R - Part 2
...awing 100 samples and then calculate the
model for each sample). Here is the code:
--snip--
# sample size n=200
ergebnisse200 <- rep(0, each=100)
stichproben200 <- vector(?list?, 100)
default200 <- rep(0, each=100)
for (i in seq(1:100)) {
n <- dim(daten)[1]
ix <- sample(n,200)
samp_i <- daten[ix,] # draw samples
y <- sum(samp_i$y) # number of defaults
stichproben200[[i]] <- samp_i # saving the samples
default200[i] <- y # saving the number of defaults
# Modell berechnen:
posterior_i <- MCMClogit(y ~ fbl.ind + fekq3 + febitda4 + fuvs + fkru +
fzd + fur3,...