Question

In: Statistics and Probability

In the usual statistical softwares we have implemented several methods to generate (pseudo) random numbers according...

In the usual statistical softwares we have implemented several methods to generate (pseudo) random numbers according to a continuous uniform distribution in the interval [0, 1]. The method of the integral transformation of probability can be used to transform a random sample of such distribution into a random sample of essentially any other univariate distribution.

b) Lecture about some difficulties of the integral transformation of probability, in particular arguing that it is not applicable (nor can be easily adapted) to generate numbers according to a multivariate distribution.

Solutions

Expert Solution

This is a question related to understanding the difficulties /challenges of practicing /implementing the integral tranformation of probability which is used to transform a random sample of such distribution into a random sample of essentially any other univariate distribution.

We need to discuss the difficulties in its implementation particularly for multi variate distributions

But what is an integral probability transformation.

The probability integral transform states that if is a continuous random variable with cumulative distribution function , then the random variable has a uniform distribution on [0, 1].

The Probability integral transform (PIT) implies that the data values that are modelled as being random variables from any given continuous distribution can be converted to random variables having a standard uniform distribution

This holds exactly provided that the distribution being used is the true distribution of the random variables; if the distribution is one fitted to the data, the result will hold approximately in large samples.

So why is it important?

PIT helps to provide the basis for testing whether a set of observations can reasonably be modelled as arising from a specified distribution.

Let us take an example carry an exponential function .

If X has an exponential distribution with unit mean, then its CDF is

and the immediate result of the probability integral transform is that

has a uniform distribution.

The symmetry of the uniform distribution can then be used to show that

also has a uniform distribution.

Thus,

Suppose that a random variable X has a continuous distribution for which the cumulative distribution function (CDF) is FX. Then the random variable Y defined as

has a uniform distribution.

Now, how does one go about generating uniform random numbers? The short answer is that it’s impossible. This should not surprising as truly random uniform numbers from [0, 1] will take on irrational values with probability one and it is impossible to precisely represent irrational numbers.

However, one can approximate random numbers, or in other words, one can generate sets of rational numbers which in many ways appear to behave as if they were random uniform draws. (called as congruential pseudo-random number generation)

The intuition behind the congruential rule is that when one divides strings of very large numbers and then discards the integer part of the quotient, the result is a set of uniform-like numbers. More specifically, the rule generates a set of numbers

xi = (α + βxi−1) mod l, i = 1, 2, . . . ,

where x0, l, α and β < l are known beforehand, and mod denotes the modulus operator.

But this has a drawback .

One drawback is that it has a period of at most l, or in other words, it can generate at most l numbers before arriving back at x0.

The remedy to it can be to choosing l prime which helps to maximize the period of the algorithm.

A larger prime l also lengthens the period and generates a better approximation to true randomness

Another difficulty is that  The random uniform distribution created  is different for distinct values of the range, but the same range of random numbers  will repeatedly produce the same number. This effect multiplies as the application of the random distribution is increased in the process.


Related Solutions

Generate 1000 random numbers from ??3? starting with standard normal random numbers in R.
Generate 1000 random numbers from ??3? starting with standard normal random numbers in R.
Write a function that will generate an array of random numbers. It needs to:
DO IN C++Write a function that will generate an array of random numbers. It needs to:Take 3 integers as parameters-The first is the minimum value-the second is the maximum value-the third is the size of the new array-create a dynamically allocated array of the correct size-generate a random number (between min and max) for each element in the array-return a pointer to the arrayCreate a main() function that tests the above function and displays the values in the random array.
Generate 1000 random numbers from ??2, 5? starting with standard normal random numbers in R.
Generate 1000 random numbers from ??2, 5? starting with standard normal random numbers in R.
Write the R code to generate five independent uniform random numbers and use that to generate...
Write the R code to generate five independent uniform random numbers and use that to generate 5 independent Bernoulli random variables with probability of success p = 0.4. please use RStudio. please do not use lab nor rbern for calculations.
1) We have studied several types of common but faulty statistical reasoning (violation of the area...
1) We have studied several types of common but faulty statistical reasoning (violation of the area principle, supposing that correlation implies causation, etc.) and limitations of testing (sample bias, errors in hypothesis testing, etc.). Find TWO examples of these flaws/errors in advertisements, news articles, or research papers. Give the source (url or other citation). Explain the mistake in each, and explain what might be done differently to correct the mistake.
Write an Arduino code that does the following. Generate 50 random numbers between the numbers 100...
Write an Arduino code that does the following. Generate 50 random numbers between the numbers 100 and 300. Pick a number at random out of these 50 random variables. a. Determine the probability of the chosen number being greater than 200. This may be achieved by counting the numbers that are greater than 200 and dividing the count by 50. Make sure you, i.Formulate the appropriate if-conditions to check for a number being greater than 200 ii. Use a for-loop...
To generate 100 random numbers between 1-100 in a randomData.txt file To read the 100 random...
To generate 100 random numbers between 1-100 in a randomData.txt file To read the 100 random numbers from randomData.txt and store them in an array Print the data in the array Find the smallest and the largest of the random numbers and their array position Insert an element of value100 in the 51th position of the array Delete all the elements of the array having values between 50-80 and print the residual array Sort the data in the final array(residual)...
A computer was used to generate ten random numbers from a normal distribution with a set...
A computer was used to generate ten random numbers from a normal distribution with a set of unknown mean and variance: −1.1623, 0.2210, 1.6518, −1.1312, −0.2879, −1.0458, 1.3706, −0.7492, −0.1355, −1.2686. Eight more random normal numbers with the same variance perhaps a different mean were then generated (the mean may or may not actually be different): 0.3472, 2.2437, 1.0712, 2.5906, 0.5163, −1.1743, 0.0473, −0.8338. (a) What do you think the means of the random normal number generators were? What do...
Write a Java program to generate random numbers in the following range a. 1 <=n <=...
Write a Java program to generate random numbers in the following range a. 1 <=n <= 3 b. 1 <= n <= 200 c. 0 <= n <= 9 d. 1000 <= n <= 2112 e. -1 <= n <= 5 JAVA PROGRAMMING
Propose two of your own random number generation schemes. Please generate 100 random numbers ?? (?...
Propose two of your own random number generation schemes. Please generate 100 random numbers ?? (? = 1,2, … ,100) for each scheme and show the results on the same plot for comparison (i.e., x-axis of the plot will show the index ? and y-axis will show the generated random numbers ??. You can use different colors and/or symbols to distinguish one sequence from the other). Discuss which scheme will be preferred.
ADVERTISEMENT
ADVERTISEMENT
ADVERTISEMENT