Skip to content
Prev 258350 / 398502 Next

how to generate a normal distribution with mean=1, min=0.2, max=0.8

That method (creating lots of samples and throwing most of them away) is 
usually frowned upon :-).

Try this:  (I haven't, so it may well have syntax errors)

% n28<- dnorm(seq(.2,.8,by=.001),mean=1,sd=1)

% x <- sample(seq(.2,.8,by=.001), size=500,replace=TRUE, prob=n28)

And I guess in retrospect this  will get really ugly if you want, say, a 
sampling grid resolution of 1e-6 or so.

Anyone know what's applicable from the "sampling" package?


Carl

-------<quote>__________________
From: David Winsemius <dwinsemius_at_comcast.net>
Date: Thu, 28 Apr 2011 13:06:21 -0400
On Apr 28, 2011, at 12:09 PM, Ravi Varadhan wrote:
> Surely you must be joking, Mr. Jianfeng.
 >

Perhaps not joking and perhaps not with correct statistical specification.

A truncated Normal could be simulated with:

set.seed(567)
x <- rnorm(n=50000, m=1, sd=1)
xtrunc <- x[x>=0.2 & x <=0.8]
require(logspline)
plot(logspline(xtrunc, lbound=0.2, ubound=0.8, nknots=7))