Results
Rj Editor +
> #Jatin Goyal, 2021EEB1180, Set B
> #a)
> given = data$SetB
> u = (given-2)/3
> mean(u)
[1] 0.4979
> var(u)
[1] 0.08406
> Q1 = median(u) - IQR(u)/2
> Q2 = median(u)
> Q3 = median(u) + IQR(u)/2
> Q1
[1] 0.237
> Q2
[1] 0.4852
> Q3
[1] 0.7334
> #b)
> myexp = function(u)
{
x = log(1-u)*(-1/3)
return(x)
}
> x = myexp(u)
> mean(x)
[1] 0.3344
> var(x)
[1] 0.1141
> Q1 = median(x) - IQR(x)/2
> Q2 = median(x)
> Q3 = median(x) + IQR(x)/2
> Q1
[1] 0.04025
> Q2
[1] 0.2213
> Q3
[1] 0.4024
> #c)
> mygamma = function(x,alpha,n)
{
y = matrix(x, nrow = alpha, ncol = n)
colSums(y)
}
> y = mygamma(x,10,100)
> mean(y)
[1] 3.344
> var(y)
[1] 1.17
> Q1 = median(y) - IQR(y)/2
> Q2 = median(y)
> Q3 = median(y) + IQR(y)/2
> Q1
[1] 2.449
> Q2
[1] 3.224
> Q3
[1] 4
> #d)
> a = y[1:50]
> b = y[51:100]
> c = a+b
> ve = a/c
> mean(ve)
[1] 0.5042
> var(ve)
[1] 0.01118
> Q1 = median(ve) - IQR(ve)/2
> Q2 = median(ve)
> Q3 = median(ve) + IQR(ve)/2
> Q1
[1] 0.439
> Q2
[1] 0.5144
> Q3
[1] 0.5897
> hist(ve,col = "blue")
> # abline(v=c(0.5,mean(ve)), col = ("green","red"), lwd = 2, lty = c(1,2))
> plot(density(ve), col = "blue", lwd = 2)
> # abline(v=c(0.5,mean(ve)),lty = c(1,2), col = ("green","red"), lwd = 2)
References
[1] The jamovi project (2024). jamovi. (Version 2.5) [Computer Software]. Retrieved from https://www.jamovi.org.
[2] R Core Team (2023). R: A Language and environment for statistical computing. (Version 4.3) [Computer software]. Retrieved
from https://cran.r-project.org. (R packages retrieved from CRAN snapshot 2024-01-09).