[R] Simulation in R

From: Alexander Geisler <alexander.geisler_at_gmail.com>
Date: Thu 07 Dec 2006 - 00:20:49 GMT


I have the following problem.

My code:


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 the results (works) }


I have a data set called "daten". I produce samples of the size 200. The samples are saved in samp_i.

How is the easiest way to save this samples. My code doesn't work, only the first column of the sample samp_i is in stichproben[i]. I understand why.
My idea is to define a array stichproben and then save a matrix in the
particular fields of the array. But things like y <- matrix(c(samp_1),nrow=200, ncol=8)
and then saving y in the array stichproben doesn't work. How is it possible to define a matrix with the content of my samples (samp_i) and but them into an array stichproben? Or is there an easier way?

With kind regards


Alexander Geisler * Kaltenbach 151 * A-6272 Kaltenbach email: alexander.geisler@gmx.at | alexander.geisler@gmail.com phone: +43 650 / 811 61 90 | skpye: al1405ex

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu Dec 07 12:21:46 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 07 Dec 2006 - 18:30:21 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.