The third assignment of Multivariate Statistics. The assignment is written in Rmarkdown, a smart syntax supported by RStudio helping with formula, plot visualization and plugin codes running.
most recommend: click here for html version of assignment, you can see codes as well as plots.
You may also find the PDF Version of this assignment from github. Or if you can cross the fire wall, just see below:
density function for the multivariate normal distribution:
the likelihood function for the two independent sample:
So the likelihood function can be defined as:
Using MLE, the maximum is:
2
(a)
It can be calculated that:
we can calculate S’s eigen value and its respective eigen vector:
The axes of the region are:
so the axes’ lengths are:
The directions of each aixs is determined by its corresponding eigen vector shown above.
(b)
1 |
# read the data |
1 |
plot(x1,x2) |
It seems that each variable’s normality is fine and they don’t have a significant relationship with each other, so the multivariate normal assumption seems justied.
3
the confidence region is defined as :
so the three regions are:
(b)
The Bonferroni region is defined as :
so the Bonferroni regions are:
which are smaller than $T^2$ confidence region because it focus on single confidence interval.
近期评论