# HỌ và TÊN: LÊ TUYẾT NHI
# LỚP 24NH1 MSSV: 2454020057
# BÀI KIỂM TRA CUỐI KỲ
#a) Thu thập - Nhập dữ liệu
PTDLCK<- file.choose()
options(width = 200)
options(scipen = 3)
tnhi<- read.csv(PTDLCK,header = TRUE)
save(tnhi,file="tnhi.rda")
attach(tnhi)
is.data.frame(tnhi)
## [1] TRUE
print(tnhi)
## STT Ngày VN.Index VN100 VN30
## 1 1 22/09/2025 1634.45 1780.83 1819.60
## 2 2 23/09/2025 1635.26 1781.08 1821.45
## 3 3 24/09/2025 1657.46 1814.83 1853.48
## 4 4 25/09/2025 1666.09 1821.95 1858.67
## 5 5 26/09/2025 1660.70 1812.71 1852.65
## 6 6 29/09/2025 1666.48 1820.38 1862.23
## 7 7 30/09/2025 1661.70 1817.70 1863.13
## 8 8 1/10/2025 1665.05 1826.30 1870.63
## 9 9 2/10/2025 1652.71 1810.69 1859.80
## 10 10 3/10/2025 1645.82 1803.15 1859.51
## 11 11 6/10/2025 1695.50 1864.34 1918.97
## 12 12 7/10/2025 1685.30 1850.55 1909.65
## 13 13 8/10/2025 1697.83 1861.87 1922.95
## 14 14 9/10/2025 1716.47 1879.62 1940.89
## 15 15 10/10/2025 1747.55 1911.55 1980.57
## 16 16 13/10/2025 1765.12 1935.65 2012.28
## 17 17 14/10/2025 1761.06 1935.25 2013.69
## 18 18 15/10/2025 1757.95 1937.17 2009.64
## 19 19 16/10/2025 1766.85 1950.03 2022.27
## 20 20 17/10/2025 1731.19 1908.40 1977.14
## 21 21 20/10/2025 1636.43 1800.51 1870.86
## 22 22 21/10/2025 1663.43 1832.98 1915.90
## 23 23 22/10/2025 1678.50 1847.12 1930.88
## 24 24 23/10/2025 1687.06 1854.20 1945.78
## 25 25 24/10/2025 1683.18 1848.01 1944.60
## 26 26 27/10/2025 1652.54 1810.65 1900.76
## 27 27 28/10/2025 1680.50 1853.69 1949.28
## 28 28 29/10/2025 1685.83 1859.11 1949.76
## 29 29 30/10/2025 1669.57 1835.82 1925.18
## 30 30 31/10/2025 1639.65 1798.62 1885.36
## 31 31 3/11/2025 1617.00 1765.36 1857.64
## 32 32 4/11/2025 1651.98 1812.10 1897.71
## 33 33 5/11/2025 1654.89 1802.59 1886.47
## 34 34 6/11/2025 1642.64 1785.91 1869.60
## 35 35 7/11/2025 1599.10 1735.99 1824.71
## 36 36 10/11/2025 1580.54 1715.56 1804.18
## 37 37 11/11/2025 1593.61 1732.96 1821.60
## 38 38 12/11/2025 1631.86 1779.78 1872.27
## 39 39 13/11/2025 1631.44 1773.67 1864.23
## 40 40 14/11/2025 1635.46 1780.35 1871.54
## 41 41 17/11/2025 1654.42 1804.18 1893.54
## 42 42 18/11/2025 1659.92 1810.72 1898.07
## 43 43 19/11/2025 1649.00 1796.70 1886.20
## 44 44 20/11/2025 1655.99 1804.39 1897.46
## 45 45 21/11/2025 1654.93 1802.47 1899.89
## 46 46 24/11/2025 1667.98 1813.34 1916.36
## 47 47 25/11/2025 1660.36 1803.41 1909.60
## 48 48 26/11/2025 1680.36 1826.35 1923.55
## 49 49 27/11/2025 1684.32 1827.32 1921.18
## 50 50 28/11/2025 1690.99 1829.88 1923.92
## 51 51 1/12/2025 1701.67 1834.56 1933.56
## 52 52 2/12/2025 1717.06 1847.40 1950.12
## 53 53 3/12/2025 1731.77 1868.56 1971.99
## 54 54 4/12/2025 1737.24 1878.72 1979.53
## 55 55 5/12/2025 1741.32 1876.20 1975.50
## 56 56 8/12/2025 1753.74 1880.12 1983.82
## 57 57 9/12/2025 1747.17 1869.86 1973.02
## 58 58 10/12/2025 1718.98 1844.94 1946.98
## 59 59 11/12/2025 1698.90 1827.05 1924.29
## 60 60 12/12/2025 1646.89 1767.08 1867.03
## 61 61 15/12/2025 1646.01 1765.08 1869.84
## 62 62 16/12/2025 1679.18 1807.62 1909.87
## 63 63 17/12/2025 1673.66 1798.31 1897.95
## 64 64 18/12/2025 1676.98 1804.06 1903.47
## 65 65 19/12/2025 1704.31 1831.78 1933.28
## 66 66 22/12/2025 1751.03 1885.64 1985.28
## 67 67 23/12/2025 1772.15 1907.33 2012.87
## 68 68 24/12/2025 1782.82 1922.80 2023.13
## 69 69 25/12/2025 1742.85 1873.79 1976.21
## 70 70 26/12/2025 1729.80 1861.57 1965.97
## 71 71 29/12/2025 1754.84 1884.24 1990.66
## 72 72 30/12/2025 1766.90 1900.27 2009.70
## 73 73 31/12/2025 1784.49 1918.00 2030.63
## 74 74 5/1/2026 1788.40 1912.69 2028.68
## 75 75 6/1/2026 1816.27 1935.93 2055.96
## 76 76 7/1/2026 1861.58 1973.69 2096.76
## 77 77 8/1/2026 1855.56 1956.73 2074.03
## 78 78 9/1/2026 1867.90 1944.22 2066.21
## 79 79 12/1/2026 1877.33 1961.77 2080.24
## 80 80 13/01/2026 1902.93 1976.38 2089.21
## 81 81 14/01/2026 1894.44 1950.48 2067.10
## 82 82 15/01/2026 1864.80 1935.22 2047.48
## 83 83 16/01/2026 1879.13 1956.26 2080.35
## 84 84 19/01/2026 1896.59 1971.81 2094.24
## 85 85 20/01/2026 1893.78 1964.23 2085.61
## 86 86 21/01/2026 1885.44 1954.03 2080.38
## 87 87 22/01/2026 1882.73 1961.84 2082.35
## 88 88 23/01/2026 1870.79 1956.64 2077.76
## 89 89 26/01/2026 1843.72 1909.38 2032.28
## 90 90 27/01/2026 1830.50 1893.38 2019.88
## 91 91 28/01/2026 1802.91 1869.14 1997.13
## 92 92 29/01/2026 1814.98 1881.82 2018.98
## 93 93 30/01/2026 1829.04 1895.40 2029.81
## 94 94 2/2/2026 1806.50 1869.17 2004.31
## 95 95 3/2/2026 1813.40 1872.59 1997.69
## 96 96 4/2/2026 1791.43 1861.24 1988.19
## 97 97 5/2/2026 1782.56 1847.81 1969.27
## 98 98 6/2/2026 1755.49 1817.49 1943.60
## 99 99 9/2/2026 1754.82 1819.63 1947.75
## 100 100 10/2/2026 1754.03 1826.92 1951.58
## 101 101 11/2/2026 1796.85 1879.24 2000.90
## 102 102 12/2/2026 1814.09 1895.81 2016.47
## 103 103 13/02/2026 1824.09 1904.16 2018.64
## 104 104 23/02/2026 1860.14 1927.92 2039.80
## 105 105 24/02/2026 1867.62 1936.82 2050.85
## 106 106 25/02/2026 1860.91 1932.71 2049.64
## 107 107 26/02/2026 1879.64 1952.69 2069.82
## 108 108 27/02/2026 1880.33 1948.72 2061.75
## 109 109 2/3/2026 1846.10 1903.50 2010.75
## 110 110 3/3/2026 1813.14 1864.78 1959.35
## 111 111 4/3/2026 1818.27 1856.66 1956.53
## 112 112 5/3/2026 1808.51 1850.66 1942.76
## 113 113 6/3/2026 1767.84 1817.73 1904.19
## 114 114 9/3/2026 1652.79 1699.60 1780.71
## 115 115 10/3/2026 1676.73 1740.71 1836.96
## 116 116 11/3/2026 1728.34 1796.23 1889.94
## 117 117 12/3/2026 1709.61 1776.44 1859.80
## 118 118 13/03/2026 1696.24 1772.56 1853.60
## 119 119 16/03/2026 1693.21 1772.60 1852.99
## 120 120 17/03/2026 1710.29 1788.73 1873.65
## 121 121 18/03/2026 1713.83 1783.56 1868.84
## 122 122 19/03/2026 1699.13 1771.74 1854.19
## 123 123 20/03/2026 1647.81 1726.81 1797.99
## 124 124 23/03/2026 1591.17 1664.96 1741.05
## 125 125 24/03/2026 1614.77 1696.93 1770.16
## 126 126 25/03/2026 1658.19 1744.40 1814.73
## 127 127 26/03/2026 1644.63 1728.57 1791.25
## 128 128 27/03/2026 1672.80 1759.10 1821.53
## 129 129 30/03/2026 1662.54 1749.23 1811.92
## 130 130 31/03/2026 1674.49 1767.05 1829.59
## 131 131 1/4/2026 1702.93 1792.16 1861.84
## 132 132 2/4/2026 1694.82 1781.33 1852.99
## 133 133 3/4/2026 1684.04 1766.97 1837.43
## 134 134 6/4/2026 1674.99 1760.97 1836.25
## 135 135 7/4/2026 1677.54 1769.19 1840.96
## 136 136 8/4/2026 1756.55 1856.15 1931.01
## 137 137 9/4/2026 1736.68 1843.49 1915.01
## 138 138 10/4/2026 1750.00 1856.76 1928.23
## 139 139 13/04/2026 1758.96 1860.33 1925.66
## 140 140 14/04/2026 1775.65 1878.19 1946.55
## 141 141 15/04/2026 1800.65 1892.92 1961.60
## 142 142 16/04/2026 1819.83 1908.08 1979.19
## 143 143 17/04/2026 1817.17 1907.34 1988.11
## 144 144 20/04/2026 1837.11 1927.96 2009.04
## 145 145 21/04/2026 1833.48 1925.79 2007.08
## 146 146 22/04/2026 1857.30 1943.91 2025.41
## 147 147 23/04/2026 1870.36 1944.05 2024.74
## 148 148 24/04/2026 1853.29 1930.67 2011.42
## 149 149 28/04/2026 1875.84 1952.67 2041.40
## 150 150 29/04/2026 1854.10 1936.31 2022.75
## 151 151 4/5/2026 1854.06 1928.89 2014.62
## 152 152 5/5/2026 1874.85 1944.47 2032.30
## 153 153 6/5/2026 1891.20 1967.64 2053.41
## 154 154 7/5/2026 1909.01 1992.48 2079.10
## 155 155 8/5/2026 1915.37 1997.02 2074.06
## 156 156 11/5/2026 1895.50 1978.37 2040.51
## 157 157 12/5/2026 1901.10 1988.39 2053.97
## 158 158 13/05/2026 1898.37 1976.98 2043.51
## 159 159 14/05/2026 1925.46 1996.22 2068.62
## 160 160 15/05/2026 1921.60 1982.89 2050.58
## 161 161 18/05/2026 1927.94 1978.45 2046.37
## 162 162 19/05/2026 1912.93 1965.64 2027.45
## 163 163 20/05/2026 1913.23 1958.35 2028.94
## 164 164 21/05/2026 1896.89 1952.00 2027.51
## 165 165 22/05/2026 1877.13 1937.14 2010.93
## 166 166 25/05/2026 1886.03 1950.39 2021.72
## 167 167 26/05/2026 1884.18 1956.87 2027.90
## 168 168 27/05/2026 1874.43 1951.75 2022.46
## 169 169 28/05/2026 1863.67 1933.02 1999.82
## 170 170 29/05/2026 1863.49 1925.55 1997.06
## 171 171 1/6/2026 1844.54 1915.27 1989.71
## 172 172 2/6/2026 1826.47 1894.18 1972.99
## 173 173 3/6/2026 1819.01 1893.09 1974.60
## 174 174 4/6/2026 1831.55 1902.15 1982.29
## 175 175 5/6/2026 1838.90 1904.58 1986.28
## 176 176 8/6/2026 1790.53 1859.67 1936.81
## 177 177 9/6/2026 1793.05 1871.55 1951.92
## 178 178 10/6/2026 1803.71 1881.42 1960.97
## 179 179 11/6/2026 1798.61 1871.39 1947.28
## 180 180 12/6/2026 1791.65 1865.93 1944.36
## 181 181 15/06/2026 1799.31 1882.63 1962.48
## 182 182 16/06/2026 1807.94 1887.98 1960.19
## 183 183 17/06/2026 1806.20 1887.29 1957.17
## 184 184 18/06/2026 1830.47 1899.67 1967.22
## 185 185 19/06/2026 1824.53 1894.86 1963.57
## 186 186 22/06/2026 1857.91 1915.25 1980.13
## 187 187 23/06/2026 1869.04 1928.14 1995.12
## 188 188 24/06/2026 1878.02 1940.44 2010.15
## 189 189 25/06/2026 1863.07 1932.68 2004.62
## 190 190 26/06/2026 1871.91 1937.79 2008.57
## 191 191 29/06/2026 1854.97 1933.17 2004.29
## 192 192 30/06/2026 1860.01 1931.52 1995.71
## 193 193 1/7/2026 1867.21 1942.19 2013.35
## 194 194 2/7/2026 1866.35 1938.99 2009.04
## 195 195 3/7/2026 1862.08 1932.28 2002.56
## 196 196 6/7/2026 1843.50 1914.29 1991.11
## 197 197 7/7/2026 1848.25 1922.44 1995.91
## 198 198 8/7/2026 1853.70 1926.12 1998.44
## 199 199 9/7/2026 1840.70 1908.47 1987.11
## 200 200 10/7/2026 1828.34 1891.94 1970.82
## 201 201 13/07/2026 1800.54 1861.08 1939.84
## 202 202 14/07/2026 1806.63 1868.47 1946.06
## 203 203 15/07/2026 1782.12 1839.44 1916.05
## 204 204 16/07/2026 1804.24 1864.31 1941.65
## 205 205 17/07/2026 1787.45 1852.65 1931.65
## 206 206 20/07/2026 1743.51 1805.16 1887.32
## 207 207 21/07/2026 1730.56 1797.80 1881.32
## 208 208 22/07/2026 1668.53 1743.88 1826.90
## 209 209 23/07/2026 1699.38 1767.70 1845.32
## 210 210 24/07/2026 1686.11 1751.49 1829.36
## 211 211 27/07/2026 1669.01 1729.22 1806.50
## 212 212 28/07/2026 1680.62 1745.18 1824.35
## 213 213 29/07/2026 1704.68 1762.60 1849.12
## 214 214 30/07/2026 1744.66 1800.95 1886.27
## 215 215 31/07/2026 1735.78 1786.90 1872.07
## 216 216 3/8/2026 1762.84 1825.30 1917.69
## 217 217 4/8/2026 1777.23 1841.63 1927.35
## 218 218 5/8/2026 1776.46 1836.22 1916.88
## 219 219 6/8/2026 1764.78 1822.07 1902.79
## 220 220 7/8/2026 1768.06 1823.70 1911.09
## 221 221 10/8/2026 1776.77 1838.02 1925.18
## 222 222 11/8/2026 1773.41 1837.69 1922.53
## 223 223 12/8/2026 1793.18 1850.94 1936.46
## 224 224 13/08/2026 1765.63 1824.51 1909.23
## 225 225 14/08/2026 1729.08 1795.69 1876.81
## 226 226 17/08/2026 1727.46 1793.86 1877.68
## 227 227 18/08/2026 1732.02 1791.35 1876.14
## 228 228 19/08/2026 1726.69 1788.39 1875.73
## 229 229 20/08/2026 1734.24 1795.57 1887.06
## 230 230 21/08/2026 1768.12 1837.77 1927.79
## 231 231 24/08/2026 1788.78 1853.99 1942.00
## 232 232 25/08/2026 1791.41 1851.76 1936.15
## 233 233 26/08/2026 1821.32 1879.52 1970.01
## 234 234 27/08/2026 1831.56 1888.06 1979.23
## 235 235 28/08/2026 1832.12 1891.44 1982.96
## 236 236 3/9/2026 1827.72 1874.20 1961.57
## 237 237 4/9/2026 1853.08 1894.93 1984.89
## 238 238 7/9/2026 1821.64 1870.98 1963.01
## 239 239 8/9/2026 1830.44 1873.51 1968.52
## 240 240 9/9/2026 1827.12 1869.84 1967.64
## 241 241 10/9/2026 1829.23 1876.33 1976.82
## 242 242 11/9/2026 1795.21 1838.59 1936.69
## 243 243 14/09/2026 1788.23 1826.99 1928.57
## 244 244 15/09/2026 1811.15 1846.98 1951.14
## 245 245 16/09/2026 1810.11 1847.52 1954.29
## 246 246 17/09/2026 1822.77 1866.59 1975.01
## 247 247 18/09/2026 1815.66 1860.89 1964.17
## 248 248 21/09/2026 1799.67 1845.03 1953.26
## 249 249 22/09/2026 1816.93 1856.95 1965.36
## 250 250 23/09/2026 1801.65 1850.79 1954.91
## 251 251 24/09/2026 1775.09 1827.38 1932.15
## 252 252 25/09/2026 1785.11 1836.01 1938.50
## 253 253 28/09/2026 1780.68 1821.76 1922.42
## 254 254 29/09/2026 1777.73 1813.71 1914.35
## 255 255 30/09/2026 1768.62 1806.88 1906.81
## 256 256 1/10/2026 1749.30 1790.08 1890.57
## 257 257 2/10/2026 1737.71 1772.01 1875.99
## 258 258 5/10/2026 1753.20 1785.54 1894.46
## 259 259 6/10/2026 1759.08 1790.89 1898.82
## 260 260 7/10/2026 1753.39 1784.17 1894.20
## 261 261 8/10/2026 1738.97 1764.95 1877.00
## 262 262 9/10/2026 1735.09 1766.75 1873.43
dim(tnhi)
## [1] 262 5
names(tnhi)
## [1] "STT" "Ngày" "VN.Index" "VN100" "VN30"
str(tnhi)
## 'data.frame': 262 obs. of 5 variables:
## $ STT : int 1 2 3 4 5 6 7 8 9 10 ...
## $ Ngày : chr "22/09/2025" "23/09/2025" "24/09/2025" "25/09/2025" ...
## $ VN.Index: num 1634 1635 1657 1666 1661 ...
## $ VN100 : num 1781 1781 1815 1822 1813 ...
## $ VN30 : num 1820 1821 1853 1859 1853 ...
#b) Mã hoá các chỉ số từ Ratio Scale -> Ordinal Scale
# Mã hóa VN-Index
MHVNI=VN.Index
MHVNI=replace(MHVNI,VN.Index<=median(VN.Index),1)
MHVNI=replace(MHVNI,VN.Index>median(VN.Index),2)
VNIcode=data.frame(VN.Index,MHVNI)
VNIcode
## VN.Index MHVNI
## 1 1634.45 1
## 2 1635.26 1
## 3 1657.46 1
## 4 1666.09 1
## 5 1660.70 1
## 6 1666.48 1
## 7 1661.70 1
## 8 1665.05 1
## 9 1652.71 1
## 10 1645.82 1
## 11 1695.50 1
## 12 1685.30 1
## 13 1697.83 1
## 14 1716.47 1
## 15 1747.55 1
## 16 1765.12 1
## 17 1761.06 1
## 18 1757.95 1
## 19 1766.85 1
## 20 1731.19 1
## 21 1636.43 1
## 22 1663.43 1
## 23 1678.50 1
## 24 1687.06 1
## 25 1683.18 1
## 26 1652.54 1
## 27 1680.50 1
## 28 1685.83 1
## 29 1669.57 1
## 30 1639.65 1
## 31 1617.00 1
## 32 1651.98 1
## 33 1654.89 1
## 34 1642.64 1
## 35 1599.10 1
## 36 1580.54 1
## 37 1593.61 1
## 38 1631.86 1
## 39 1631.44 1
## 40 1635.46 1
## 41 1654.42 1
## 42 1659.92 1
## 43 1649.00 1
## 44 1655.99 1
## 45 1654.93 1
## 46 1667.98 1
## 47 1660.36 1
## 48 1680.36 1
## 49 1684.32 1
## 50 1690.99 1
## 51 1701.67 1
## 52 1717.06 1
## 53 1731.77 1
## 54 1737.24 1
## 55 1741.32 1
## 56 1753.74 1
## 57 1747.17 1
## 58 1718.98 1
## 59 1698.90 1
## 60 1646.89 1
## 61 1646.01 1
## 62 1679.18 1
## 63 1673.66 1
## 64 1676.98 1
## 65 1704.31 1
## 66 1751.03 1
## 67 1772.15 1
## 68 1782.82 2
## 69 1742.85 1
## 70 1729.80 1
## 71 1754.84 1
## 72 1766.90 1
## 73 1784.49 2
## 74 1788.40 2
## 75 1816.27 2
## 76 1861.58 2
## 77 1855.56 2
## 78 1867.90 2
## 79 1877.33 2
## 80 1902.93 2
## 81 1894.44 2
## 82 1864.80 2
## 83 1879.13 2
## 84 1896.59 2
## 85 1893.78 2
## 86 1885.44 2
## 87 1882.73 2
## 88 1870.79 2
## 89 1843.72 2
## 90 1830.50 2
## 91 1802.91 2
## 92 1814.98 2
## 93 1829.04 2
## 94 1806.50 2
## 95 1813.40 2
## 96 1791.43 2
## 97 1782.56 2
## 98 1755.49 1
## 99 1754.82 1
## 100 1754.03 1
## 101 1796.85 2
## 102 1814.09 2
## 103 1824.09 2
## 104 1860.14 2
## 105 1867.62 2
## 106 1860.91 2
## 107 1879.64 2
## 108 1880.33 2
## 109 1846.10 2
## 110 1813.14 2
## 111 1818.27 2
## 112 1808.51 2
## 113 1767.84 1
## 114 1652.79 1
## 115 1676.73 1
## 116 1728.34 1
## 117 1709.61 1
## 118 1696.24 1
## 119 1693.21 1
## 120 1710.29 1
## 121 1713.83 1
## 122 1699.13 1
## 123 1647.81 1
## 124 1591.17 1
## 125 1614.77 1
## 126 1658.19 1
## 127 1644.63 1
## 128 1672.80 1
## 129 1662.54 1
## 130 1674.49 1
## 131 1702.93 1
## 132 1694.82 1
## 133 1684.04 1
## 134 1674.99 1
## 135 1677.54 1
## 136 1756.55 1
## 137 1736.68 1
## 138 1750.00 1
## 139 1758.96 1
## 140 1775.65 2
## 141 1800.65 2
## 142 1819.83 2
## 143 1817.17 2
## 144 1837.11 2
## 145 1833.48 2
## 146 1857.30 2
## 147 1870.36 2
## 148 1853.29 2
## 149 1875.84 2
## 150 1854.10 2
## 151 1854.06 2
## 152 1874.85 2
## 153 1891.20 2
## 154 1909.01 2
## 155 1915.37 2
## 156 1895.50 2
## 157 1901.10 2
## 158 1898.37 2
## 159 1925.46 2
## 160 1921.60 2
## 161 1927.94 2
## 162 1912.93 2
## 163 1913.23 2
## 164 1896.89 2
## 165 1877.13 2
## 166 1886.03 2
## 167 1884.18 2
## 168 1874.43 2
## 169 1863.67 2
## 170 1863.49 2
## 171 1844.54 2
## 172 1826.47 2
## 173 1819.01 2
## 174 1831.55 2
## 175 1838.90 2
## 176 1790.53 2
## 177 1793.05 2
## 178 1803.71 2
## 179 1798.61 2
## 180 1791.65 2
## 181 1799.31 2
## 182 1807.94 2
## 183 1806.20 2
## 184 1830.47 2
## 185 1824.53 2
## 186 1857.91 2
## 187 1869.04 2
## 188 1878.02 2
## 189 1863.07 2
## 190 1871.91 2
## 191 1854.97 2
## 192 1860.01 2
## 193 1867.21 2
## 194 1866.35 2
## 195 1862.08 2
## 196 1843.50 2
## 197 1848.25 2
## 198 1853.70 2
## 199 1840.70 2
## 200 1828.34 2
## 201 1800.54 2
## 202 1806.63 2
## 203 1782.12 2
## 204 1804.24 2
## 205 1787.45 2
## 206 1743.51 1
## 207 1730.56 1
## 208 1668.53 1
## 209 1699.38 1
## 210 1686.11 1
## 211 1669.01 1
## 212 1680.62 1
## 213 1704.68 1
## 214 1744.66 1
## 215 1735.78 1
## 216 1762.84 1
## 217 1777.23 2
## 218 1776.46 2
## 219 1764.78 1
## 220 1768.06 1
## 221 1776.77 2
## 222 1773.41 1
## 223 1793.18 2
## 224 1765.63 1
## 225 1729.08 1
## 226 1727.46 1
## 227 1732.02 1
## 228 1726.69 1
## 229 1734.24 1
## 230 1768.12 1
## 231 1788.78 2
## 232 1791.41 2
## 233 1821.32 2
## 234 1831.56 2
## 235 1832.12 2
## 236 1827.72 2
## 237 1853.08 2
## 238 1821.64 2
## 239 1830.44 2
## 240 1827.12 2
## 241 1829.23 2
## 242 1795.21 2
## 243 1788.23 2
## 244 1811.15 2
## 245 1810.11 2
## 246 1822.77 2
## 247 1815.66 2
## 248 1799.67 2
## 249 1816.93 2
## 250 1801.65 2
## 251 1775.09 1
## 252 1785.11 2
## 253 1780.68 2
## 254 1777.73 2
## 255 1768.62 1
## 256 1749.30 1
## 257 1737.71 1
## 258 1753.20 1
## 259 1759.08 1
## 260 1753.39 1
## 261 1738.97 1
## 262 1735.09 1
# Mã hóa VN100
MHVN1=VN100
MHVN1=replace(MHVN1,VN100<=median(VN100),1)
MHVN1=replace(MHVN1,VN100>median(VN100),2)
VN1code=data.frame(VN100,MHVN1)
VN1code
## VN100 MHVN1
## 1 1780.83 1
## 2 1781.08 1
## 3 1814.83 1
## 4 1821.95 1
## 5 1812.71 1
## 6 1820.38 1
## 7 1817.70 1
## 8 1826.30 1
## 9 1810.69 1
## 10 1803.15 1
## 11 1864.34 2
## 12 1850.55 1
## 13 1861.87 2
## 14 1879.62 2
## 15 1911.55 2
## 16 1935.65 2
## 17 1935.25 2
## 18 1937.17 2
## 19 1950.03 2
## 20 1908.40 2
## 21 1800.51 1
## 22 1832.98 1
## 23 1847.12 1
## 24 1854.20 1
## 25 1848.01 1
## 26 1810.65 1
## 27 1853.69 1
## 28 1859.11 1
## 29 1835.82 1
## 30 1798.62 1
## 31 1765.36 1
## 32 1812.10 1
## 33 1802.59 1
## 34 1785.91 1
## 35 1735.99 1
## 36 1715.56 1
## 37 1732.96 1
## 38 1779.78 1
## 39 1773.67 1
## 40 1780.35 1
## 41 1804.18 1
## 42 1810.72 1
## 43 1796.70 1
## 44 1804.39 1
## 45 1802.47 1
## 46 1813.34 1
## 47 1803.41 1
## 48 1826.35 1
## 49 1827.32 1
## 50 1829.88 1
## 51 1834.56 1
## 52 1847.40 1
## 53 1868.56 2
## 54 1878.72 2
## 55 1876.20 2
## 56 1880.12 2
## 57 1869.86 2
## 58 1844.94 1
## 59 1827.05 1
## 60 1767.08 1
## 61 1765.08 1
## 62 1807.62 1
## 63 1798.31 1
## 64 1804.06 1
## 65 1831.78 1
## 66 1885.64 2
## 67 1907.33 2
## 68 1922.80 2
## 69 1873.79 2
## 70 1861.57 2
## 71 1884.24 2
## 72 1900.27 2
## 73 1918.00 2
## 74 1912.69 2
## 75 1935.93 2
## 76 1973.69 2
## 77 1956.73 2
## 78 1944.22 2
## 79 1961.77 2
## 80 1976.38 2
## 81 1950.48 2
## 82 1935.22 2
## 83 1956.26 2
## 84 1971.81 2
## 85 1964.23 2
## 86 1954.03 2
## 87 1961.84 2
## 88 1956.64 2
## 89 1909.38 2
## 90 1893.38 2
## 91 1869.14 2
## 92 1881.82 2
## 93 1895.40 2
## 94 1869.17 2
## 95 1872.59 2
## 96 1861.24 2
## 97 1847.81 1
## 98 1817.49 1
## 99 1819.63 1
## 100 1826.92 1
## 101 1879.24 2
## 102 1895.81 2
## 103 1904.16 2
## 104 1927.92 2
## 105 1936.82 2
## 106 1932.71 2
## 107 1952.69 2
## 108 1948.72 2
## 109 1903.50 2
## 110 1864.78 2
## 111 1856.66 1
## 112 1850.66 1
## 113 1817.73 1
## 114 1699.60 1
## 115 1740.71 1
## 116 1796.23 1
## 117 1776.44 1
## 118 1772.56 1
## 119 1772.60 1
## 120 1788.73 1
## 121 1783.56 1
## 122 1771.74 1
## 123 1726.81 1
## 124 1664.96 1
## 125 1696.93 1
## 126 1744.40 1
## 127 1728.57 1
## 128 1759.10 1
## 129 1749.23 1
## 130 1767.05 1
## 131 1792.16 1
## 132 1781.33 1
## 133 1766.97 1
## 134 1760.97 1
## 135 1769.19 1
## 136 1856.15 1
## 137 1843.49 1
## 138 1856.76 1
## 139 1860.33 2
## 140 1878.19 2
## 141 1892.92 2
## 142 1908.08 2
## 143 1907.34 2
## 144 1927.96 2
## 145 1925.79 2
## 146 1943.91 2
## 147 1944.05 2
## 148 1930.67 2
## 149 1952.67 2
## 150 1936.31 2
## 151 1928.89 2
## 152 1944.47 2
## 153 1967.64 2
## 154 1992.48 2
## 155 1997.02 2
## 156 1978.37 2
## 157 1988.39 2
## 158 1976.98 2
## 159 1996.22 2
## 160 1982.89 2
## 161 1978.45 2
## 162 1965.64 2
## 163 1958.35 2
## 164 1952.00 2
## 165 1937.14 2
## 166 1950.39 2
## 167 1956.87 2
## 168 1951.75 2
## 169 1933.02 2
## 170 1925.55 2
## 171 1915.27 2
## 172 1894.18 2
## 173 1893.09 2
## 174 1902.15 2
## 175 1904.58 2
## 176 1859.67 2
## 177 1871.55 2
## 178 1881.42 2
## 179 1871.39 2
## 180 1865.93 2
## 181 1882.63 2
## 182 1887.98 2
## 183 1887.29 2
## 184 1899.67 2
## 185 1894.86 2
## 186 1915.25 2
## 187 1928.14 2
## 188 1940.44 2
## 189 1932.68 2
## 190 1937.79 2
## 191 1933.17 2
## 192 1931.52 2
## 193 1942.19 2
## 194 1938.99 2
## 195 1932.28 2
## 196 1914.29 2
## 197 1922.44 2
## 198 1926.12 2
## 199 1908.47 2
## 200 1891.94 2
## 201 1861.08 2
## 202 1868.47 2
## 203 1839.44 1
## 204 1864.31 2
## 205 1852.65 1
## 206 1805.16 1
## 207 1797.80 1
## 208 1743.88 1
## 209 1767.70 1
## 210 1751.49 1
## 211 1729.22 1
## 212 1745.18 1
## 213 1762.60 1
## 214 1800.95 1
## 215 1786.90 1
## 216 1825.30 1
## 217 1841.63 1
## 218 1836.22 1
## 219 1822.07 1
## 220 1823.70 1
## 221 1838.02 1
## 222 1837.69 1
## 223 1850.94 1
## 224 1824.51 1
## 225 1795.69 1
## 226 1793.86 1
## 227 1791.35 1
## 228 1788.39 1
## 229 1795.57 1
## 230 1837.77 1
## 231 1853.99 1
## 232 1851.76 1
## 233 1879.52 2
## 234 1888.06 2
## 235 1891.44 2
## 236 1874.20 2
## 237 1894.93 2
## 238 1870.98 2
## 239 1873.51 2
## 240 1869.84 2
## 241 1876.33 2
## 242 1838.59 1
## 243 1826.99 1
## 244 1846.98 1
## 245 1847.52 1
## 246 1866.59 2
## 247 1860.89 2
## 248 1845.03 1
## 249 1856.95 1
## 250 1850.79 1
## 251 1827.38 1
## 252 1836.01 1
## 253 1821.76 1
## 254 1813.71 1
## 255 1806.88 1
## 256 1790.08 1
## 257 1772.01 1
## 258 1785.54 1
## 259 1790.89 1
## 260 1784.17 1
## 261 1764.95 1
## 262 1766.75 1
# Mã hóa VN30
MHVN3=VN30
MHVN3=replace(MHVN3,VN30<=median(VN30),1)
MHVN3=replace(MHVN3,VN30>median(VN30),2)
VN3code=data.frame(VN30,MHVN3)
VN3code
## VN30 MHVN3
## 1 1819.60 1
## 2 1821.45 1
## 3 1853.48 1
## 4 1858.67 1
## 5 1852.65 1
## 6 1862.23 1
## 7 1863.13 1
## 8 1870.63 1
## 9 1859.80 1
## 10 1859.51 1
## 11 1918.97 1
## 12 1909.65 1
## 13 1922.95 1
## 14 1940.89 1
## 15 1980.57 2
## 16 2012.28 2
## 17 2013.69 2
## 18 2009.64 2
## 19 2022.27 2
## 20 1977.14 2
## 21 1870.86 1
## 22 1915.90 1
## 23 1930.88 1
## 24 1945.78 1
## 25 1944.60 1
## 26 1900.76 1
## 27 1949.28 2
## 28 1949.76 2
## 29 1925.18 1
## 30 1885.36 1
## 31 1857.64 1
## 32 1897.71 1
## 33 1886.47 1
## 34 1869.60 1
## 35 1824.71 1
## 36 1804.18 1
## 37 1821.60 1
## 38 1872.27 1
## 39 1864.23 1
## 40 1871.54 1
## 41 1893.54 1
## 42 1898.07 1
## 43 1886.20 1
## 44 1897.46 1
## 45 1899.89 1
## 46 1916.36 1
## 47 1909.60 1
## 48 1923.55 1
## 49 1921.18 1
## 50 1923.92 1
## 51 1933.56 1
## 52 1950.12 2
## 53 1971.99 2
## 54 1979.53 2
## 55 1975.50 2
## 56 1983.82 2
## 57 1973.02 2
## 58 1946.98 1
## 59 1924.29 1
## 60 1867.03 1
## 61 1869.84 1
## 62 1909.87 1
## 63 1897.95 1
## 64 1903.47 1
## 65 1933.28 1
## 66 1985.28 2
## 67 2012.87 2
## 68 2023.13 2
## 69 1976.21 2
## 70 1965.97 2
## 71 1990.66 2
## 72 2009.70 2
## 73 2030.63 2
## 74 2028.68 2
## 75 2055.96 2
## 76 2096.76 2
## 77 2074.03 2
## 78 2066.21 2
## 79 2080.24 2
## 80 2089.21 2
## 81 2067.10 2
## 82 2047.48 2
## 83 2080.35 2
## 84 2094.24 2
## 85 2085.61 2
## 86 2080.38 2
## 87 2082.35 2
## 88 2077.76 2
## 89 2032.28 2
## 90 2019.88 2
## 91 1997.13 2
## 92 2018.98 2
## 93 2029.81 2
## 94 2004.31 2
## 95 1997.69 2
## 96 1988.19 2
## 97 1969.27 2
## 98 1943.60 1
## 99 1947.75 1
## 100 1951.58 2
## 101 2000.90 2
## 102 2016.47 2
## 103 2018.64 2
## 104 2039.80 2
## 105 2050.85 2
## 106 2049.64 2
## 107 2069.82 2
## 108 2061.75 2
## 109 2010.75 2
## 110 1959.35 2
## 111 1956.53 2
## 112 1942.76 1
## 113 1904.19 1
## 114 1780.71 1
## 115 1836.96 1
## 116 1889.94 1
## 117 1859.80 1
## 118 1853.60 1
## 119 1852.99 1
## 120 1873.65 1
## 121 1868.84 1
## 122 1854.19 1
## 123 1797.99 1
## 124 1741.05 1
## 125 1770.16 1
## 126 1814.73 1
## 127 1791.25 1
## 128 1821.53 1
## 129 1811.92 1
## 130 1829.59 1
## 131 1861.84 1
## 132 1852.99 1
## 133 1837.43 1
## 134 1836.25 1
## 135 1840.96 1
## 136 1931.01 1
## 137 1915.01 1
## 138 1928.23 1
## 139 1925.66 1
## 140 1946.55 1
## 141 1961.60 2
## 142 1979.19 2
## 143 1988.11 2
## 144 2009.04 2
## 145 2007.08 2
## 146 2025.41 2
## 147 2024.74 2
## 148 2011.42 2
## 149 2041.40 2
## 150 2022.75 2
## 151 2014.62 2
## 152 2032.30 2
## 153 2053.41 2
## 154 2079.10 2
## 155 2074.06 2
## 156 2040.51 2
## 157 2053.97 2
## 158 2043.51 2
## 159 2068.62 2
## 160 2050.58 2
## 161 2046.37 2
## 162 2027.45 2
## 163 2028.94 2
## 164 2027.51 2
## 165 2010.93 2
## 166 2021.72 2
## 167 2027.90 2
## 168 2022.46 2
## 169 1999.82 2
## 170 1997.06 2
## 171 1989.71 2
## 172 1972.99 2
## 173 1974.60 2
## 174 1982.29 2
## 175 1986.28 2
## 176 1936.81 1
## 177 1951.92 2
## 178 1960.97 2
## 179 1947.28 1
## 180 1944.36 1
## 181 1962.48 2
## 182 1960.19 2
## 183 1957.17 2
## 184 1967.22 2
## 185 1963.57 2
## 186 1980.13 2
## 187 1995.12 2
## 188 2010.15 2
## 189 2004.62 2
## 190 2008.57 2
## 191 2004.29 2
## 192 1995.71 2
## 193 2013.35 2
## 194 2009.04 2
## 195 2002.56 2
## 196 1991.11 2
## 197 1995.91 2
## 198 1998.44 2
## 199 1987.11 2
## 200 1970.82 2
## 201 1939.84 1
## 202 1946.06 1
## 203 1916.05 1
## 204 1941.65 1
## 205 1931.65 1
## 206 1887.32 1
## 207 1881.32 1
## 208 1826.90 1
## 209 1845.32 1
## 210 1829.36 1
## 211 1806.50 1
## 212 1824.35 1
## 213 1849.12 1
## 214 1886.27 1
## 215 1872.07 1
## 216 1917.69 1
## 217 1927.35 1
## 218 1916.88 1
## 219 1902.79 1
## 220 1911.09 1
## 221 1925.18 1
## 222 1922.53 1
## 223 1936.46 1
## 224 1909.23 1
## 225 1876.81 1
## 226 1877.68 1
## 227 1876.14 1
## 228 1875.73 1
## 229 1887.06 1
## 230 1927.79 1
## 231 1942.00 1
## 232 1936.15 1
## 233 1970.01 2
## 234 1979.23 2
## 235 1982.96 2
## 236 1961.57 2
## 237 1984.89 2
## 238 1963.01 2
## 239 1968.52 2
## 240 1967.64 2
## 241 1976.82 2
## 242 1936.69 1
## 243 1928.57 1
## 244 1951.14 2
## 245 1954.29 2
## 246 1975.01 2
## 247 1964.17 2
## 248 1953.26 2
## 249 1965.36 2
## 250 1954.91 2
## 251 1932.15 1
## 252 1938.50 1
## 253 1922.42 1
## 254 1914.35 1
## 255 1906.81 1
## 256 1890.57 1
## 257 1875.99 1
## 258 1894.46 1
## 259 1898.82 1
## 260 1894.20 1
## 261 1877.00 1
## 262 1873.43 1
#c) Thống kê các chỉ số
# Thống kê VN-Index
summary(VN.Index)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1581 1696 1775 1769 1832 1928
var(VN.Index)
## [1] 6837.034
sd(VN.Index)
## [1] 82.68636
# Thống kê VN100
summary(VN100)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1665 1805 1859 1859 1914 1997
var(VN100)
## [1] 4677.053
sd(VN100)
## [1] 68.38899
# Thống kê VN30
summary(VN30)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1741 1894 1949 1947 2004 2097
var(VN30)
## [1] 5367.961
sd(VN30)
## [1] 73.26637
library(psych)
# Thống kê Ratio Scale
TK=data.frame(VN.Index,VN100,VN30)
describe(TK)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## VN.Index 1 262 1768.87 82.69 1775.37 1769.78 104.83 1580.54 1927.94 347.40 -0.13 -1.02 5.11
## VN100 2 262 1858.85 68.39 1859.39 1859.61 81.47 1664.96 1997.02 332.06 -0.09 -0.70 4.23
## VN30 3 262 1947.10 73.27 1948.51 1947.96 82.71 1741.05 2096.76 355.71 -0.14 -0.59 4.53
# Thống kê Ordinal
table(MHVNI)
## MHVNI
## 1 2
## 131 131
table(MHVN1)
## MHVN1
## 1 2
## 131 131
table(MHVN3)
## MHVN3
## 1 2
## 131 131
###Nhận xét:
##VN Index có giá trị trung bình là 1768.87, có mức độ phân tán lớn nhất 82.69
#Nhóm 1 có 131 quan sát(50%), nhóm 2 có 131 quan sát(50%)
##VN100 có giá trị trung bình là 1858.85, có mức độ phân tán thấp nhất 68.93
#Nhóm 1 có 131 quan sát(50%), nhóm 2 có 131 quan sát(50%)
##VN30 có giá trị trung bình là 1947.10, có mức độ phân tán 73,27
#Nhóm 1 có 131 quan sát(50%), nhóm 2 có 131 quan sát(50%)
#d) Vẽ biểu đồ
# Biểu đồ cột VN-Index
VNI.freq=table(MHVNI)
VNI.freq
## MHVNI
## 1 2
## 131 131
barplot(VNI.freq,xlab="Mức VN-Index",ylab="Số lượng",main="Biểu đồ cột VN-Index", col=c("red","yellow"))
# Biểu đồ cột VN100
VN1.freq=table(MHVN1)
VN1.freq
## MHVN1
## 1 2
## 131 131
barplot(VN1.freq,xlab="Mức VN100",ylab="Số lượng",main="Biểu đồ cột VN100", col=c("blue","pink"))
# Biểu đồ cột VN30
VN3.freq=table(MHVN3)
VN3.freq
## MHVN3
## 1 2
## 131 131
barplot(VN3.freq,xlab="Mức VN30",ylab="Số lượng",main="Biểu đồ cột VN30", col=c("orange","green"))
### Nhận xét:
##Hai nhóm của 3 biến VN Index,VN100,Vn30 đều có số lượng quan sát bằng nhau, mỗi nhóm gồm 131 quan sát.
##Mỗi nhóm chiếm 50% tổng số quan sát.
##Dữ liệu phân bố đồng đều giữa hai nhóm của các biến.
# Biểu đồ tròn VN-Index
pie(table(MHVNI), col=c("red","yellow"),main="Biểu đồ tròn VN-Index")
# Biểu đồ tròn VN100
pie(table(MHVN1), col=c("blue","pink"),main="Biểu đồ tròn VN100")
# Biểu đồ tròn VN30
pie(table(MHVN3), col=c("orange","green"),main="Biểu đồ tròn VN30")
### Nhận xét:
##Hai nhóm của 3 biến VN Index,VN100,Vn30 đều có tỷ trọng bằng nhau, mỗi nhóm chiếm 50%.
##Không có nhóm nào chiếm ưu thế về tỷ trọng.
# Biểu đồ đường
plot(VN.Index,type="o",xlab="Ngày thứ",ylab="Giá trị VN-Index",pch=19,main="Giá trị VN-Index theo ngày",col="pink")
plot(VN100,type="o",xlab="Ngày thứ",ylab="Giá trị VN100",pch=19,main="Giá trị VN100 theo ngày",col="violet")
plot(VN30,type="o",xlab="Ngày thứ",ylab="Giá trị VN30",pch=19,main="Giá trị VN30 theo ngày",col="yellow")
###Nhận xét:
##Ba chỉ số VN Index, VN100 và VN30 có xu hướng biến động cùng chiều.
##Các chỉ số có sự thay đổi qua các quan sát theo thời gian.
##Mức giá trị của VN30 nhìn chung cao hơn VN100 và VN Index.
#e) Phân tích tương quan + đồ thị tương quan các chỉ số tài chính dạng Ratio Scale
library(psych)
vars=cbind(VN.Index,VN100,VN30)
pairs.panels(vars)
library(Hmisc)
##
## Attaching package: 'Hmisc'
## The following object is masked from 'package:psych':
##
## describe
## The following objects are masked from 'package:base':
##
## format.pval, units
result=rcorr(as.matrix(vars),type="pearson")
result$r
## VN.Index VN100 VN30
## VN.Index 1.0000000 0.8767609 0.8731357
## VN100 0.8767609 1.0000000 0.9634635
## VN30 0.8731357 0.9634635 1.0000000
result$P
## VN.Index VN100 VN30
## VN.Index NA 0 0
## VN100 0 NA 0
## VN30 0 0 NA
###Nhận xét:
## VN Index có mối tương quan dương với VN100,VN30 vì hệ số tương quan Pearson của các mối quan hệ này >0
## mối tương quan giữ VN Index và VN100,VN30 có ý nghĩa thống kê vì giá trị Sig.(P-Value)của các mối quan hệ này<5%
#f) Mô hình hồi quy tính tuyến đa biến
mreg=lm(VN.Index~VN100+VN30)
summary(mreg)
##
## Call:
## lm(formula = VN.Index ~ VN100 + VN30)
##
## Residuals:
## Min 1Q Median 3Q Max
## -90.21 -41.06 13.98 30.03 57.65
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -214.3274 65.6450 -3.265 0.00124 **
## VN100 0.5988 0.1316 4.551 0.00000823 ***
## VN30 0.4469 0.1228 3.639 0.00033 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 38.94 on 259 degrees of freedom
## Multiple R-squared: 0.78, Adjusted R-squared: 0.7783
## F-statistic: 459 on 2 and 259 DF, p-value: < 2.2e-16
anova(mreg)
## Analysis of Variance Table
##
## Response: VN.Index
## Df Sum Sq Mean Sq F value Pr(>F)
## VN100 1 1371736 1371736 904.814 < 2.2e-16 ***
## VN30 1 20075 20075 13.242 0.0003304 ***
## Residuals 259 392655 1516
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
###Nhận xét:
##Phương trình hồi quy ước lượng:VN Index = -214,3274 + 0,5988 × VN100 + 0,4469 × VN30 + eᵢ
##Kết quả hồi quy cho thấy mô hình có ý nghĩa thống kê tổng thể
##với F-statistic = 459 và p-value < 0,001.
##Hệ số R² = 0,78 cho thấy hai biến VN100 và VN30 giải thích khoảng 78% sự biến thiên của VN Index
##trong khi R² hiệu chỉnh đạt 0,7783.
##Xét riêng từng biến, VN100 có tác động cùng chiều đến VN Index với hệ số hồi quy 0,5988 (p < 0,001).
##Điều này có nghĩa khi VN100 tăng 1 điểm thì VN Index tăng trung bình khoảng 0,5988 điểm, trong điều kiện VN30 không đổi.
##Tương tự, VN30 có tác động cùng chiều với hệ số hồi quy 0,4469 (p < 0,001), nghĩa là khi VN30 tăng 1 điểm thì VN Index tăng trung bình khoảng 0,4469 điểm, trong điều kiện VN100 không đổi.
##Nhìn chung, cả VN100 và VN30 đều có mối quan hệ cùng chiều và có ý nghĩa thống kê với VN Index.
##Mô hình giải thích được phần lớn sự biến thiên của VN Index
##tuy nhiên vẫn còn khoảng 22% sự biến thiên chưa được giải thích bởi các biến độc lập trong mô hình.
#g) Các phép kiểm định
# Kiểm định đa cộng tuyến tính
library(car)
## Loading required package: carData
##
## Attaching package: 'car'
## The following object is masked from 'package:psych':
##
## logit
vif(mreg)
## VN100 VN30
## 13.93961 13.93961
# Kiểm định tự tương quan
require(lmtest)
## Loading required package: lmtest
## Loading required package: zoo
##
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
##
## as.Date, as.Date.numeric
dwtest(mreg)
##
## Durbin-Watson test
##
## data: mreg
## DW = 0.033457, p-value < 2.2e-16
## alternative hypothesis: true autocorrelation is greater than 0
# Kiểm định các giả thuyết mô hình bằng đồ thị
op<-par(mfrow=c(2,2))
plot(mreg)
### Nhận xét:
#Kiểm định đa cộng tuyến: Kết quả cho thấy hệ số VIF của VN100 và VN30 đều bằng 13,93961, lớn hơn 10.
#=> Mô hình có hiện tượng đa cộng tuyến.
#Kiểm định tự tương quan: Kiểm định Durbin–Watson cho kết quả DW ≈ 0,033 và p-value < 0,05.
#=> Mô hình có hiện tượng tự tương quan dương.
# Nhận xét đồ thị:
#Đồ thị bên trái dòng 1 (Residuals vs Fitted): Các giá trị phần dư phân bố tương đối quanh đường 0 nhưng có dấu hiệu chưa phân bố hoàn toàn ngẫu nhiên.
#=> Giả định về giá trị trung bình của sai số bằng 0 có thể được xem xét, tuy nhiên cần kiểm tra thêm các giả định khác của mô hình.
#Đồ thị bên phải dòng 1 (Q-Q Residuals): Các điểm phần dư có dấu hiệu lệch khỏi đường thẳng tham chiếu, đặc biệt ở hai đầu.
#=> Phần dư có dấu hiệu chưa tuân theo phân phối chuẩn hoàn toàn.
#Đồ thị bên trái dòng 2 (Scale-Location): Các điểm phân bố chưa hoàn toàn đồng đều theo giá trị dự báo và có dấu hiệu thay đổi độ phân tán.
#=> Có khả năng xảy ra hiện tượng phương sai sai số thay đổi.
#Đồ thị bên phải dòng 2 (Residuals vs Leverage): Phần lớn các quan sát có leverage thấp, tuy nhiên cần kiểm tra thêm các điểm có ảnh hưởng lớn đến mô hình.
#=> Chưa có dấu hiệu rõ ràng về các quan sát có ảnh hưởng mạnh đến mô hình.