Re: [R] How can I sample from a two-dimensional grid of points

From: Dimitris Rizopoulos <dimitris.rizopoulos_at_med.kuleuven.be>
Date: Mon, 10 Mar 2008 15:57:00 +0100

you can still use sample(), e.g.,

n.grid <- 500
muA.grid <- seq(-4, 4, length = n.grid)
muB.grid <- seq(-4, 4, length = n.grid)

vals <- data.matrix(expand.grid(muA.grid, muB.grid))
mu.p <- dnorm(vals[, 1]) * dnorm(vals[, 2], sd = 0.5)
mu.p <- mu.p / sum(mu.p)

ind <- 1:nrow(vals)
samp.ind <- sample(ind, 1000, TRUE, mu.p) mu <- vals[samp.ind, ]
head(mu, 10)

I hope it helps.

Best,
Dimitris



Dimitris Rizopoulos
Biostatistical Centre
School of Public Health
Catholic University of Leuven

Address: Kapucijnenvoer 35, Leuven, Belgium

Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
     http://www.student.kuleuven.be/~m0390867/dimitris.htm

> Hi everyone,
>
> My goal is to sample from a two-dimensional grid. Consider the
> following example of code:
>
> n.grid <- 500
> muA.grid <- seq(-4,4, length=n.grid)
> muB.grid <- seq(-4,4, length=n.grid)
> mu.p <- matrix(NA, nrow=n.grid, ncol=n.grid)
> for(i in 1:n.grid){
> for(j in 1:n.grid){
> mu.p[i,j] <- dnorm(muA.grid[i], 0, 1)*dnorm(muB.grid[j], 0, 0.5)
> }
> }
>
> mu.p <- mu.p/sum(mu.p)
>
> I would now like to sample the grid of points from the
> probabilities in mu.p. Im using the multivariate normal here for
> illustration as my real problem is a more complicated probability
> density. If this problem were only one-dimensional, this is easy:
>
> n.samples <- 1000
> # assuming mu.p and muA.grid are now the appropriate vectors
> mu <- sample(muA.grid, n.samples, replace=T, prob=mu.p)
>
> However, im not sure how to do this in two-dimensions in R.
>
> Thanks in advance for any help.
>
> All the best,
> Gregory Gentlemen
>
>
> ---------------------------------
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help_at_r-project.org 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.
>

Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm



R-help_at_r-project.org 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 Mon 10 Mar 2008 - 15:00:53 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 10 Mar 2008 - 15: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.

list of date sections of archive