Take the 2-minute tour ×
Cross Validated is a question and answer site for people interested in statistics, machine learning, data analysis, data mining, and data visualization. It's 100% free, no registration required.

How to transform a normal random variable such that I can simulate normal samples between the range of 1 and 45?

Do I need to do a jacobian transformation ?

share|improve this question

2 Answers 2

up vote 12 down vote accepted

Any normal distribution is defined on the entire real line $(-\infty,\infty)$.

This means that while we can choose some probability of being between 1 and 45, you can't have it both be actually normal and definitely in that range.

You have to decide which of the two of those can be compromised.

However, it also sounds like perhaps you're after an integer result; if that's so, that, again, means you don't have normality, and have to compromise one or the other.

Are you after a distribution with mean 23 ($\frac{_{45+1}}{^2}$)?

Can you say any more about the situation?


Whoever asked you to do this seems to be operating under a misunderstanding because the conditions are inconsistent.

I have no idea how to choose the variance parameter.

The smaller the variance, the less chance of going outside the range. But the smaller the variance, the less chance of taking up most of the range.

For example, if you choose $\sigma=5$, you'll have about 11 of your 1000 samples with at least one value outside the bounds, but on the other hand, most of the time you won't have any values very close to 1 or 45.

You can use simulation to judge the tradeoff.

Here's a simulated distribution of the largest observation for n=1000 with mean 23 and sd 5:

enter image description here

As you see almost all of the time, it's less than 45 but then often it's less than 40.


One possibility is to truncate the normal; e.g. try mean 45 sd around 6 (say) and then regenerate values that fall outside the range (or generate slightly too many and then omit what you don't need).

Here's an example with sd=6 and another with sd=9:

enter image description here

enter image description here

(This is not a single sample of 1000, but many samples of 1000 so you can see the distribution)

The narrower one looks more normal but you're less likely to get values near 1 or 45; the wider one is a little less normal looking but generally gives values right out near the ends.

You may want to choose something in between.


Another possibility: the person may actually want you to fit it to 1 to 45 post-hoc (that is to make the smallest value 1 and the largest 45). This is fine, but you don't actually have normality any more (though with n=1000 it will look very close).

For example, here's what happens when you use this strategy with n=6 instead of n=1000:

enter image description here

(1000 samples of n=6)

There are spikes at 1 and 45 because those values are certain to occur, but leaving that aside, the distribution of values in between 1 and 45 is fairly flat, only gently bulging in the center. [Macro makes a similar point about the n=4 case here; it's even more nearly flat, though even the n=3 case still shows a slight bulge in the center.]

For small samples this isn't much use. As sample sizes increase, this improves. Here it is for samples of size 100:

enter image description here

(1000 samples of n=100)

It looks like a Winsorized normal (or ignoring the endpoints, it's like a truncated normal). This will still be having an impact at n=1000:

enter image description here

Note that this isn't a single sample, but 1000 samples of 1000. I did that to show you what the actual distribution you're ending up with is like.

I think you're better off with a truncated normal than doing this. Nevertheless, if this is what they want, it is very easy to do.

share|improve this answer
    
That's what I thought as well. A normal distribution is defined on the entire real line. But I was asked to draw normal samples between the range of 1 and 45. And I am not after an integer result. So I think I should just set the mean to 23. But I have no idea how to choose the variance parameter. –  user1769197 Aug 1 at 9:50
    
Do you have a specific sample size or range of sample sizes? –  Glen_b Aug 1 at 9:54
    
Sure, 1000 samples but that's the number of simulations --- how many observations in a sample? –  Glen_b Aug 1 at 10:00
    
1000 samples of size 1000? Okay. –  Glen_b Aug 1 at 10:02
    
Sorry, I meant 1000 samples of size 1. Sorry for the confusion. –  user1769197 Aug 1 at 10:12

How about just sampling 1000 from the standard normal distribution and then scaling the outcome to [1,45] by 44*(x - min(x))/(max(x)-min(x))+1.

share|improve this answer
    
This is a useful shortcut, even though the conflict @Glen_b describes still must be reckoned with. –  rolando2 Aug 1 at 10:45
    
The result, unfortunately, won't be normal any longer, since the min and max are random variables. It will, however, be fairly close with a large sample. –  Glen_b Aug 1 at 10:56
    
I've done some simulations to show what happens with this approach in my answer. In some cases, it's a reasonable option (in particular, at n=1000, many people may be happy with the result). –  Glen_b Aug 1 at 11:55

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.