Probability: Normal Distribution

 Normal Distribution

Generate a Normal Distribution

We can use the rnorm() function in R to generate random values from a normal distribution with a specified mean and standard deviation. The syntax of the function is:


rnorm(n, mean = 0, sd = 1)

where n is the number of values to generate, mean is the mean of the distribution, and sd is the standard deviation of the distribution.

For example, to generate 1000 random values from a normal distribution with a mean of 0 and a standard deviation of 1, we can use the following code:


set.seed(123) # set a seed for reproducibility
x <- rnorm(1000, mean = 0, sd = 1)

The set.seed() function is used to set a seed value for the random number generator, which ensures that the same set of random numbers will be generated each time the code is run. This can be useful for reproducibility.

Plot a Histogram of the Distribution

We can use the hist() function to plot a histogram of the generated values. The syntax of the function is:


hist(x, breaks = ..., freq = TRUE, main = ..., xlab = ..., ylab = ...)

where x is the vector of values to plot, breaks is the number of intervals in the histogram, freq is a logical value indicating whether to plot frequencies or densities, main is the main title of the plot, xlab is the x-axis label, and ylab is the y-axis label.

For example, to plot a histogram of the 1000 random values we generated earlier, we can use the following code:


hist(x, breaks = 30, freq = TRUE, main = "Normal Distribution", xlab = "Value", ylab = "Frequency")

This will plot a histogram with 30 bars, showing the frequency of each value in the generated distribution. The main title of the plot is "Normal Distribution", and the x-axis label is "Value", and the y-axis label is "Frequency".

Calculate Summary Statistics of the Distribution

We can use the mean() and sd() functions to calculate the mean and standard deviation of the generated values, respectively.

For example, to calculate the mean and standard deviation of the 1000 random values we generated earlier, we can use the following code:


mean(x) # calculate the mean 
sd(x) # calculate the standard deviation

This will output the mean and standard deviation of the generated values.

Generate a Normal Probability Density Function

We can use the dnorm() function to generate a probability density function (PDF) of a normal distribution. The syntax of the function is:


dnorm(x, mean = 0, sd = 1)

where x is the vector of values to calculate the PDF for, mean is the mean of the distribution, and sd is the standard deviation of the distribution.

For example, to generate a PDF of a normal distribution with a mean of 0 and a standard deviation of 1, we can use the following code:


x <- seq(-4, 4, length.out = 1000) # create a sequence of values 
y <- dnorm(x, mean = 0, sd = 1) # calculate the PDF 
plot(x, y, type = "l", main = "Normal Probability Density Function", xlab = "Value", ylab = "Density") 

pnorm Function

The pnorm() function is used to calculate the cumulative probability density function (CDF) of a normal distribution. The CDF represents the probability that a random variable from the distribution takes on a value less than or equal to a specified value. The syntax of the pnorm() function is:

pnorm(x, mean = 0, sd = 1, lower.tail = TRUE)

where x is the vector of values to calculate the CDF for, mean is the mean of the distribution, sd is the standard deviation of the distribution, and lower.tail is a logical value indicating whether to calculate the probability of values less than or equal to x (lower.tail = TRUE) or greater than x (lower.tail = FALSE).

For example, if we want to calculate the probability that a random variable from a normal distribution with a mean of 0 and a standard deviation of 1 is less than or equal to 1, we can use the following code:

pnorm(1, mean = 0, sd = 1, lower.tail = TRUE)

This will return the value of the CDF at x = 1, which is approximately 0.8413.

qnorm Function

The qnorm() function is used to calculate the quantile function of a normal distribution. The quantile function represents the value of the distribution for a given probability. In other words, it gives us the inverse of the CDF. The syntax of the qnorm() function is:

qnorm(p, mean = 0, sd = 1, lower.tail = TRUE)

where p is the probability value for which to calculate the quantile, mean is the mean of the distribution, sd is the standard deviation of the distribution, and lower.tail is a logical value indicating whether the quantile should be calculated for probabilities less than or equal to p (lower.tail = TRUE) or greater than p (lower.tail = FALSE).

For example, if we want to find the value of a normal distribution with a mean of 0 and a standard deviation of 1 that has a cumulative probability of 0.8413 (as we found earlier), we can use the following code:

qnorm(0.8413, mean = 0, sd = 1, lower.tail = TRUE)

This will return the value of the distribution that has a cumulative probability of 0.8413, which is approximately 1.


https://rstudio-pubs-static.s3.amazonaws.com/533598_9c53bf874d9b41ea8390d0c9046bea71.html

Comments

Post a Comment

Popular posts from this blog

Probability: Binomial Distribution

 Binomial Distribution The binomial distribution is a discrete probability distribution that describes the number of successes in a fixed number of independent trials, where each trial has a constant probability of success. The binomial distribution has two parameters: n, the number of trials, and p, the probability of success in each trial. The probability of getting exactly k successes in n trials is given by the formula: P(X=k) = (n choose k) * p^k * (1-p)^(n-k) where (n choose k) is the binomial coefficient, which is equal to n! / (k! * (n-k)!) and represents the number of ways to choose k successes out of n trials. The mean or expected value of a binomial distribution is given by E(X) = n * p, and the variance is given by Var(X) = n * p * (1-p). The binomial distribution is widely used in statistics and probability theory, particularly in fields such as genetics, engineering, and finance, where it is used to model the probability of success or failure in a series of trials. Co...
Read more

Hypothesis testing

Drive:  https://drive.google.com/drive/folders/1XSJoE-3rCDDJ6S4mRb8rYCgw5CQ_PYxP?usp=share_link Problem Statement Paired dataset: Suppose we have a dataset of 20 students who took a math test before and after attending a tutoring program. We want to test if the tutoring program had a significant effect on the students' math scores. Our null hypothesis is that there is no significant difference between the mean scores before and after the tutoring program, and our alternative hypothesis is that there is a significant difference. Data is stored in a file called "math_scores.csv" with two columns: "pre_tutoring_scores" and "post_tutoring_scores".  Load the data: math_scores <- read.csv("math_scores.csv") Calculate the sample mean and standard deviation for both pre-tutoring and post-tutoring scores: pre_tutoring_mean <- mean(math_scores$pre_tutoring_scores) pre_tutoring_sd <- sd(math_scores$pre_tutoring_scores) post_tutoring_mean <- ...
Read more