#Họ và tên: Lê Nguyễn Hà Nam
#Lớp: 24NH MSSV:2454020035
#Bài thực hành cuối kỳ
#a) Nhập liệu
Angnem=file.choose()
dataAngnem=read.csv(Angnem,header=TRUE)
save(dataAngnem,file="dataAngnem.rda")
attach(dataAngnem)
is.data.frame(dataAngnem)
## [1] TRUE
#b) Mã hoá các chỉ số từ Ratio Scale -> Ordinal Scale
# Mã hóa VN30
MHVN30=VN30
MHVN30=replace(MHVN30,VN30<=median(VN30),1)
MHVN30=replace(MHVN30,VN30>median(VN30),2)
VN30code=data.frame(VN30,MHVN30)
VN30code
## VN30 MHVN30
## 1 1862.23 1
## 2 1863.13 1
## 3 1870.63 1
## 4 1859.80 1
## 5 1859.51 1
## 6 1918.97 1
## 7 1909.65 1
## 8 1922.95 1
## 9 1940.89 1
## 10 1980.57 2
## 11 2012.28 2
## 12 2013.69 2
## 13 2009.64 2
## 14 2022.27 2
## 15 1977.14 2
## 16 1870.86 1
## 17 1915.90 1
## 18 1930.88 1
## 19 1945.78 1
## 20 1944.60 1
## 21 1900.76 1
## 22 1949.28 1
## 23 1949.76 1
## 24 1925.18 1
## 25 1885.36 1
## 26 1857.64 1
## 27 1897.71 1
## 28 1886.47 1
## 29 1869.60 1
## 30 1824.71 1
## 31 1804.18 1
## 32 1821.60 1
## 33 1872.27 1
## 34 1864.23 1
## 35 1871.54 1
## 36 1893.54 1
## 37 1898.07 1
## 38 1886.20 1
## 39 1897.46 1
## 40 1899.89 1
## 41 1916.36 1
## 42 1909.60 1
## 43 1923.55 1
## 44 1923.92 1
## 45 1933.56 1
## 46 1950.12 1
## 47 1971.99 2
## 48 1979.53 2
## 49 1975.50 2
## 50 1983.82 2
## 51 1973.02 2
## 52 1946.98 1
## 53 1924.29 1
## 54 1867.03 1
## 55 1869.84 1
## 56 1909.87 1
## 57 1897.95 1
## 58 1903.47 1
## 59 1933.28 1
## 60 1985.28 2
## 61 2012.87 2
## 62 2023.13 2
## 63 1976.21 2
## 64 1965.97 2
## 65 1990.66 2
## 66 2009.70 2
## 67 2030.63 2
## 68 2028.68 2
## 69 2055.96 2
## 70 2096.76 2
## 71 2074.03 2
## 72 2066.21 2
## 73 2080.24 2
## 74 2089.21 2
## 75 2067.10 2
## 76 2047.48 2
## 77 2080.35 2
## 78 2085.61 2
## 79 2080.38 2
## 80 2082.35 2
## 81 2077.76 2
## 82 2032.28 2
## 83 2019.88 2
## 84 1997.13 2
## 85 2018.98 2
## 86 2029.81 2
## 87 2004.31 2
## 88 1997.69 2
## 89 1988.19 2
## 90 1969.27 2
## 91 1943.60 1
## 92 1947.75 1
## 93 1951.58 2
## 94 2000.90 2
## 95 2016.47 2
## 96 2018.64 2
## 97 2039.80 2
## 98 2050.85 2
## 99 2049.64 2
## 100 2069.82 2
## 101 2061.75 2
## 102 2010.75 2
## 103 1959.35 2
## 104 1956.53 2
## 105 1942.76 1
## 106 1904.19 1
## 107 1780.71 1
## 108 1836.96 1
## 109 1889.94 1
## 110 1859.80 1
## 111 1853.60 1
## 112 1852.99 1
## 113 1873.65 1
## 114 1868.84 1
## 115 1854.19 1
## 116 1797.99 1
## 117 1741.05 1
## 118 1770.16 1
## 119 1814.73 1
## 120 1791.25 1
## 121 1821.53 1
## 122 1811.92 1
## 123 1829.59 1
## 124 1861.84 1
## 125 1852.99 1
## 126 1836.25 1
## 127 1840.96 1
## 128 1931.01 1
## 129 1915.01 1
## 130 1928.23 1
## 131 1925.66 1
## 132 1946.55 1
## 133 1961.60 2
## 134 1979.19 2
## 135 1988.11 2
## 136 2009.04 2
## 137 2007.08 2
## 138 2025.41 2
## 139 2024.74 2
## 140 2011.42 2
## 141 2041.40 2
## 142 2022.75 2
## 143 2014.62 2
## 144 2032.30 2
## 145 2053.41 2
## 146 2079.10 2
## 147 2074.06 2
## 148 2040.51 2
## 149 2053.97 2
## 150 2043.51 2
## 151 2068.62 2
## 152 2050.58 2
## 153 2046.37 2
## 154 2027.45 2
## 155 2028.94 2
## 156 2027.51 2
## 157 2010.93 2
## 158 2027.90 2
## 159 2022.46 2
## 160 1999.82 2
## 161 1997.06 2
## 162 1989.71 2
## 163 1972.99 2
## 164 1974.60 2
## 165 1982.29 2
## 166 1986.28 2
## 167 1936.81 1
## 168 1951.92 2
## 169 1960.97 2
## 170 1947.28 1
## 171 1944.36 1
## 172 1962.48 2
## 173 1960.19 2
## 174 1957.17 2
## 175 1967.22 2
## 176 1980.13 2
## 177 1995.12 2
## 178 2010.15 2
## 179 2004.62 2
## 180 2008.57 2
## 181 2004.29 2
## 182 1995.71 2
## 183 2013.35 2
## 184 2009.04 2
## 185 1991.11 2
## 186 1995.91 2
## 187 1998.44 2
## 188 1987.11 2
## 189 1970.82 2
## 190 1939.84 1
## 191 1946.06 1
## 192 1916.05 1
## 193 1941.65 1
## 194 1931.65 1
## 195 1887.32 1
## 196 1881.32 1
## 197 1826.90 1
## 198 1845.32 1
## 199 1829.36 1
## 200 1806.50 1
## 201 1824.35 1
## 202 1849.12 1
## 203 1886.27 1
## 204 1872.07 1
## 205 1917.69 1
## 206 1927.35 1
## 207 1916.88 1
## 208 1902.79 1
## 209 1911.09 1
## 210 1925.18 1
## 211 1922.53 1
## 212 1936.46 1
## 213 1909.23 1
## 214 1876.81 1
## 215 1877.68 1
## 216 1876.14 1
## 217 1875.73 1
## 218 1887.06 1
## 219 1927.79 1
## 220 1942.00 1
## 221 1936.15 1
## 222 1970.01 2
## 223 1979.23 2
## 224 1982.96 2
## 225 1961.57 2
## 226 1984.89 2
## 227 1968.52 2
## 228 1967.64 2
## 229 1976.82 2
## 230 1936.69 1
## 231 1928.57 1
## 232 1951.14 2
## 233 1954.29 2
## 234 1975.01 2
## 235 1964.17 2
## 236 1953.26 2
## 237 1965.36 2
## 238 1954.91 2
## 239 1932.15 1
## 240 1938.50 1
## 241 1922.42 1
## 242 1914.35 1
## 243 1906.81 1
## 244 1890.57 1
## 245 1875.99 1
## 246 1894.46 1
# MÃ HÓA VN30 SANG THANG DO ORDINAL SCALE
VN30a=MHVN30
VN30a=replace(VN30a,MHVN30==1,"Nhỏ hơn hoặc bằng 1951")
VN30a=replace(VN30a,MHVN30==2,"Lớn hơn 1951")
VN30aMH=data.frame(MHVN30,VN30a)
VN30aMH
## MHVN30 VN30a
## 1 1 Nhỏ hơn hoặc bằng 1951
## 2 1 Nhỏ hơn hoặc bằng 1951
## 3 1 Nhỏ hơn hoặc bằng 1951
## 4 1 Nhỏ hơn hoặc bằng 1951
## 5 1 Nhỏ hơn hoặc bằng 1951
## 6 1 Nhỏ hơn hoặc bằng 1951
## 7 1 Nhỏ hơn hoặc bằng 1951
## 8 1 Nhỏ hơn hoặc bằng 1951
## 9 1 Nhỏ hơn hoặc bằng 1951
## 10 2 Lớn hơn 1951
## 11 2 Lớn hơn 1951
## 12 2 Lớn hơn 1951
## 13 2 Lớn hơn 1951
## 14 2 Lớn hơn 1951
## 15 2 Lớn hơn 1951
## 16 1 Nhỏ hơn hoặc bằng 1951
## 17 1 Nhỏ hơn hoặc bằng 1951
## 18 1 Nhỏ hơn hoặc bằng 1951
## 19 1 Nhỏ hơn hoặc bằng 1951
## 20 1 Nhỏ hơn hoặc bằng 1951
## 21 1 Nhỏ hơn hoặc bằng 1951
## 22 1 Nhỏ hơn hoặc bằng 1951
## 23 1 Nhỏ hơn hoặc bằng 1951
## 24 1 Nhỏ hơn hoặc bằng 1951
## 25 1 Nhỏ hơn hoặc bằng 1951
## 26 1 Nhỏ hơn hoặc bằng 1951
## 27 1 Nhỏ hơn hoặc bằng 1951
## 28 1 Nhỏ hơn hoặc bằng 1951
## 29 1 Nhỏ hơn hoặc bằng 1951
## 30 1 Nhỏ hơn hoặc bằng 1951
## 31 1 Nhỏ hơn hoặc bằng 1951
## 32 1 Nhỏ hơn hoặc bằng 1951
## 33 1 Nhỏ hơn hoặc bằng 1951
## 34 1 Nhỏ hơn hoặc bằng 1951
## 35 1 Nhỏ hơn hoặc bằng 1951
## 36 1 Nhỏ hơn hoặc bằng 1951
## 37 1 Nhỏ hơn hoặc bằng 1951
## 38 1 Nhỏ hơn hoặc bằng 1951
## 39 1 Nhỏ hơn hoặc bằng 1951
## 40 1 Nhỏ hơn hoặc bằng 1951
## 41 1 Nhỏ hơn hoặc bằng 1951
## 42 1 Nhỏ hơn hoặc bằng 1951
## 43 1 Nhỏ hơn hoặc bằng 1951
## 44 1 Nhỏ hơn hoặc bằng 1951
## 45 1 Nhỏ hơn hoặc bằng 1951
## 46 1 Nhỏ hơn hoặc bằng 1951
## 47 2 Lớn hơn 1951
## 48 2 Lớn hơn 1951
## 49 2 Lớn hơn 1951
## 50 2 Lớn hơn 1951
## 51 2 Lớn hơn 1951
## 52 1 Nhỏ hơn hoặc bằng 1951
## 53 1 Nhỏ hơn hoặc bằng 1951
## 54 1 Nhỏ hơn hoặc bằng 1951
## 55 1 Nhỏ hơn hoặc bằng 1951
## 56 1 Nhỏ hơn hoặc bằng 1951
## 57 1 Nhỏ hơn hoặc bằng 1951
## 58 1 Nhỏ hơn hoặc bằng 1951
## 59 1 Nhỏ hơn hoặc bằng 1951
## 60 2 Lớn hơn 1951
## 61 2 Lớn hơn 1951
## 62 2 Lớn hơn 1951
## 63 2 Lớn hơn 1951
## 64 2 Lớn hơn 1951
## 65 2 Lớn hơn 1951
## 66 2 Lớn hơn 1951
## 67 2 Lớn hơn 1951
## 68 2 Lớn hơn 1951
## 69 2 Lớn hơn 1951
## 70 2 Lớn hơn 1951
## 71 2 Lớn hơn 1951
## 72 2 Lớn hơn 1951
## 73 2 Lớn hơn 1951
## 74 2 Lớn hơn 1951
## 75 2 Lớn hơn 1951
## 76 2 Lớn hơn 1951
## 77 2 Lớn hơn 1951
## 78 2 Lớn hơn 1951
## 79 2 Lớn hơn 1951
## 80 2 Lớn hơn 1951
## 81 2 Lớn hơn 1951
## 82 2 Lớn hơn 1951
## 83 2 Lớn hơn 1951
## 84 2 Lớn hơn 1951
## 85 2 Lớn hơn 1951
## 86 2 Lớn hơn 1951
## 87 2 Lớn hơn 1951
## 88 2 Lớn hơn 1951
## 89 2 Lớn hơn 1951
## 90 2 Lớn hơn 1951
## 91 1 Nhỏ hơn hoặc bằng 1951
## 92 1 Nhỏ hơn hoặc bằng 1951
## 93 2 Lớn hơn 1951
## 94 2 Lớn hơn 1951
## 95 2 Lớn hơn 1951
## 96 2 Lớn hơn 1951
## 97 2 Lớn hơn 1951
## 98 2 Lớn hơn 1951
## 99 2 Lớn hơn 1951
## 100 2 Lớn hơn 1951
## 101 2 Lớn hơn 1951
## 102 2 Lớn hơn 1951
## 103 2 Lớn hơn 1951
## 104 2 Lớn hơn 1951
## 105 1 Nhỏ hơn hoặc bằng 1951
## 106 1 Nhỏ hơn hoặc bằng 1951
## 107 1 Nhỏ hơn hoặc bằng 1951
## 108 1 Nhỏ hơn hoặc bằng 1951
## 109 1 Nhỏ hơn hoặc bằng 1951
## 110 1 Nhỏ hơn hoặc bằng 1951
## 111 1 Nhỏ hơn hoặc bằng 1951
## 112 1 Nhỏ hơn hoặc bằng 1951
## 113 1 Nhỏ hơn hoặc bằng 1951
## 114 1 Nhỏ hơn hoặc bằng 1951
## 115 1 Nhỏ hơn hoặc bằng 1951
## 116 1 Nhỏ hơn hoặc bằng 1951
## 117 1 Nhỏ hơn hoặc bằng 1951
## 118 1 Nhỏ hơn hoặc bằng 1951
## 119 1 Nhỏ hơn hoặc bằng 1951
## 120 1 Nhỏ hơn hoặc bằng 1951
## 121 1 Nhỏ hơn hoặc bằng 1951
## 122 1 Nhỏ hơn hoặc bằng 1951
## 123 1 Nhỏ hơn hoặc bằng 1951
## 124 1 Nhỏ hơn hoặc bằng 1951
## 125 1 Nhỏ hơn hoặc bằng 1951
## 126 1 Nhỏ hơn hoặc bằng 1951
## 127 1 Nhỏ hơn hoặc bằng 1951
## 128 1 Nhỏ hơn hoặc bằng 1951
## 129 1 Nhỏ hơn hoặc bằng 1951
## 130 1 Nhỏ hơn hoặc bằng 1951
## 131 1 Nhỏ hơn hoặc bằng 1951
## 132 1 Nhỏ hơn hoặc bằng 1951
## 133 2 Lớn hơn 1951
## 134 2 Lớn hơn 1951
## 135 2 Lớn hơn 1951
## 136 2 Lớn hơn 1951
## 137 2 Lớn hơn 1951
## 138 2 Lớn hơn 1951
## 139 2 Lớn hơn 1951
## 140 2 Lớn hơn 1951
## 141 2 Lớn hơn 1951
## 142 2 Lớn hơn 1951
## 143 2 Lớn hơn 1951
## 144 2 Lớn hơn 1951
## 145 2 Lớn hơn 1951
## 146 2 Lớn hơn 1951
## 147 2 Lớn hơn 1951
## 148 2 Lớn hơn 1951
## 149 2 Lớn hơn 1951
## 150 2 Lớn hơn 1951
## 151 2 Lớn hơn 1951
## 152 2 Lớn hơn 1951
## 153 2 Lớn hơn 1951
## 154 2 Lớn hơn 1951
## 155 2 Lớn hơn 1951
## 156 2 Lớn hơn 1951
## 157 2 Lớn hơn 1951
## 158 2 Lớn hơn 1951
## 159 2 Lớn hơn 1951
## 160 2 Lớn hơn 1951
## 161 2 Lớn hơn 1951
## 162 2 Lớn hơn 1951
## 163 2 Lớn hơn 1951
## 164 2 Lớn hơn 1951
## 165 2 Lớn hơn 1951
## 166 2 Lớn hơn 1951
## 167 1 Nhỏ hơn hoặc bằng 1951
## 168 2 Lớn hơn 1951
## 169 2 Lớn hơn 1951
## 170 1 Nhỏ hơn hoặc bằng 1951
## 171 1 Nhỏ hơn hoặc bằng 1951
## 172 2 Lớn hơn 1951
## 173 2 Lớn hơn 1951
## 174 2 Lớn hơn 1951
## 175 2 Lớn hơn 1951
## 176 2 Lớn hơn 1951
## 177 2 Lớn hơn 1951
## 178 2 Lớn hơn 1951
## 179 2 Lớn hơn 1951
## 180 2 Lớn hơn 1951
## 181 2 Lớn hơn 1951
## 182 2 Lớn hơn 1951
## 183 2 Lớn hơn 1951
## 184 2 Lớn hơn 1951
## 185 2 Lớn hơn 1951
## 186 2 Lớn hơn 1951
## 187 2 Lớn hơn 1951
## 188 2 Lớn hơn 1951
## 189 2 Lớn hơn 1951
## 190 1 Nhỏ hơn hoặc bằng 1951
## 191 1 Nhỏ hơn hoặc bằng 1951
## 192 1 Nhỏ hơn hoặc bằng 1951
## 193 1 Nhỏ hơn hoặc bằng 1951
## 194 1 Nhỏ hơn hoặc bằng 1951
## 195 1 Nhỏ hơn hoặc bằng 1951
## 196 1 Nhỏ hơn hoặc bằng 1951
## 197 1 Nhỏ hơn hoặc bằng 1951
## 198 1 Nhỏ hơn hoặc bằng 1951
## 199 1 Nhỏ hơn hoặc bằng 1951
## 200 1 Nhỏ hơn hoặc bằng 1951
## 201 1 Nhỏ hơn hoặc bằng 1951
## 202 1 Nhỏ hơn hoặc bằng 1951
## 203 1 Nhỏ hơn hoặc bằng 1951
## 204 1 Nhỏ hơn hoặc bằng 1951
## 205 1 Nhỏ hơn hoặc bằng 1951
## 206 1 Nhỏ hơn hoặc bằng 1951
## 207 1 Nhỏ hơn hoặc bằng 1951
## 208 1 Nhỏ hơn hoặc bằng 1951
## 209 1 Nhỏ hơn hoặc bằng 1951
## 210 1 Nhỏ hơn hoặc bằng 1951
## 211 1 Nhỏ hơn hoặc bằng 1951
## 212 1 Nhỏ hơn hoặc bằng 1951
## 213 1 Nhỏ hơn hoặc bằng 1951
## 214 1 Nhỏ hơn hoặc bằng 1951
## 215 1 Nhỏ hơn hoặc bằng 1951
## 216 1 Nhỏ hơn hoặc bằng 1951
## 217 1 Nhỏ hơn hoặc bằng 1951
## 218 1 Nhỏ hơn hoặc bằng 1951
## 219 1 Nhỏ hơn hoặc bằng 1951
## 220 1 Nhỏ hơn hoặc bằng 1951
## 221 1 Nhỏ hơn hoặc bằng 1951
## 222 2 Lớn hơn 1951
## 223 2 Lớn hơn 1951
## 224 2 Lớn hơn 1951
## 225 2 Lớn hơn 1951
## 226 2 Lớn hơn 1951
## 227 2 Lớn hơn 1951
## 228 2 Lớn hơn 1951
## 229 2 Lớn hơn 1951
## 230 1 Nhỏ hơn hoặc bằng 1951
## 231 1 Nhỏ hơn hoặc bằng 1951
## 232 2 Lớn hơn 1951
## 233 2 Lớn hơn 1951
## 234 2 Lớn hơn 1951
## 235 2 Lớn hơn 1951
## 236 2 Lớn hơn 1951
## 237 2 Lớn hơn 1951
## 238 2 Lớn hơn 1951
## 239 1 Nhỏ hơn hoặc bằng 1951
## 240 1 Nhỏ hơn hoặc bằng 1951
## 241 1 Nhỏ hơn hoặc bằng 1951
## 242 1 Nhỏ hơn hoặc bằng 1951
## 243 1 Nhỏ hơn hoặc bằng 1951
## 244 1 Nhỏ hơn hoặc bằng 1951
## 245 1 Nhỏ hơn hoặc bằng 1951
## 246 1 Nhỏ hơn hoặc bằng 1951
# Mã hóa VN100
MHVN100=VN100
MHVN100=replace(MHVN100,VN100<=median(VN100),1)
MHVN100=replace(MHVN100,VN100>median(VN100),2)
VN100code=data.frame(VN100,MHVN100)
VN100code
## VN100 MHVN100
## 1 1820.38 1
## 2 1817.70 1
## 3 1826.30 1
## 4 1810.69 1
## 5 1803.15 1
## 6 1864.34 2
## 7 1850.55 1
## 8 1861.87 2
## 9 1879.62 2
## 10 1911.55 2
## 11 1935.65 2
## 12 1935.25 2
## 13 1937.17 2
## 14 1950.03 2
## 15 1908.40 2
## 16 1800.51 1
## 17 1832.98 1
## 18 1847.12 1
## 19 1854.20 1
## 20 1848.01 1
## 21 1810.65 1
## 22 1853.69 1
## 23 1859.11 1
## 24 1835.82 1
## 25 1798.62 1
## 26 1765.36 1
## 27 1812.10 1
## 28 1802.59 1
## 29 1785.91 1
## 30 1735.99 1
## 31 1715.56 1
## 32 1732.96 1
## 33 1779.78 1
## 34 1773.67 1
## 35 1780.35 1
## 36 1804.18 1
## 37 1810.72 1
## 38 1796.70 1
## 39 1804.39 1
## 40 1802.47 1
## 41 1813.34 1
## 42 1803.41 1
## 43 1826.35 1
## 44 1829.88 1
## 45 1834.56 1
## 46 1847.40 1
## 47 1868.56 2
## 48 1878.72 2
## 49 1876.20 2
## 50 1880.12 2
## 51 1869.86 2
## 52 1844.94 1
## 53 1827.05 1
## 54 1767.08 1
## 55 1765.08 1
## 56 1807.62 1
## 57 1798.31 1
## 58 1804.06 1
## 59 1831.78 1
## 60 1885.64 2
## 61 1907.33 2
## 62 1922.80 2
## 63 1873.79 2
## 64 1861.57 2
## 65 1884.24 2
## 66 1900.27 2
## 67 1918.00 2
## 68 1912.69 2
## 69 1935.93 2
## 70 1973.69 2
## 71 1956.73 2
## 72 1944.22 2
## 73 1961.77 2
## 74 1976.38 2
## 75 1950.48 2
## 76 1935.22 2
## 77 1956.26 2
## 78 1964.23 2
## 79 1954.03 2
## 80 1961.84 2
## 81 1956.64 2
## 82 1909.38 2
## 83 1893.38 2
## 84 1869.14 2
## 85 1881.82 2
## 86 1895.40 2
## 87 1869.17 2
## 88 1872.59 2
## 89 1861.24 2
## 90 1847.81 1
## 91 1817.49 1
## 92 1819.63 1
## 93 1826.92 1
## 94 1879.24 2
## 95 1895.81 2
## 96 1904.16 2
## 97 1927.92 2
## 98 1936.82 2
## 99 1932.71 2
## 100 1952.69 2
## 101 1948.72 2
## 102 1903.50 2
## 103 1864.78 2
## 104 1856.66 1
## 105 1850.66 1
## 106 1817.73 1
## 107 1699.60 1
## 108 1740.71 1
## 109 1796.23 1
## 110 1776.44 1
## 111 1772.56 1
## 112 1772.60 1
## 113 1788.73 1
## 114 1783.56 1
## 115 1771.74 1
## 116 1726.81 1
## 117 1664.96 1
## 118 1696.93 1
## 119 1744.40 1
## 120 1728.57 1
## 121 1759.10 1
## 122 1749.23 1
## 123 1767.05 1
## 124 1792.16 1
## 125 1781.33 1
## 126 1760.97 1
## 127 1769.19 1
## 128 1856.15 1
## 129 1843.49 1
## 130 1856.76 1
## 131 1860.33 1
## 132 1878.19 2
## 133 1892.92 2
## 134 1908.08 2
## 135 1907.34 2
## 136 1927.96 2
## 137 1925.79 2
## 138 1943.91 2
## 139 1944.05 2
## 140 1930.67 2
## 141 1952.67 2
## 142 1936.31 2
## 143 1928.89 2
## 144 1944.47 2
## 145 1967.64 2
## 146 1992.48 2
## 147 1997.02 2
## 148 1978.37 2
## 149 1988.39 2
## 150 1976.98 2
## 151 1996.22 2
## 152 1982.89 2
## 153 1978.45 2
## 154 1965.64 2
## 155 1958.35 2
## 156 1952.00 2
## 157 1937.14 2
## 158 1956.87 2
## 159 1951.75 2
## 160 1933.02 2
## 161 1925.55 2
## 162 1915.27 2
## 163 1894.18 2
## 164 1893.09 2
## 165 1902.15 2
## 166 1904.58 2
## 167 1859.67 1
## 168 1871.55 2
## 169 1881.42 2
## 170 1871.39 2
## 171 1865.93 2
## 172 1882.63 2
## 173 1887.98 2
## 174 1887.29 2
## 175 1899.67 2
## 176 1915.25 2
## 177 1928.14 2
## 178 1940.44 2
## 179 1932.68 2
## 180 1937.79 2
## 181 1933.17 2
## 182 1931.52 2
## 183 1942.19 2
## 184 1938.99 2
## 185 1914.29 2
## 186 1922.44 2
## 187 1926.12 2
## 188 1908.47 2
## 189 1891.94 2
## 190 1861.08 2
## 191 1868.47 2
## 192 1839.44 1
## 193 1864.31 2
## 194 1852.65 1
## 195 1805.16 1
## 196 1797.80 1
## 197 1743.88 1
## 198 1767.70 1
## 199 1751.49 1
## 200 1729.22 1
## 201 1745.18 1
## 202 1762.60 1
## 203 1800.95 1
## 204 1786.90 1
## 205 1825.30 1
## 206 1841.63 1
## 207 1836.22 1
## 208 1822.07 1
## 209 1823.70 1
## 210 1838.02 1
## 211 1837.69 1
## 212 1850.94 1
## 213 1824.51 1
## 214 1795.69 1
## 215 1793.86 1
## 216 1791.35 1
## 217 1788.39 1
## 218 1795.57 1
## 219 1837.77 1
## 220 1853.99 1
## 221 1851.76 1
## 222 1879.52 2
## 223 1888.06 2
## 224 1891.44 2
## 225 1874.20 2
## 226 1894.93 2
## 227 1873.51 2
## 228 1869.84 2
## 229 1876.33 2
## 230 1838.59 1
## 231 1826.99 1
## 232 1846.98 1
## 233 1847.52 1
## 234 1866.59 2
## 235 1860.89 1
## 236 1845.03 1
## 237 1856.95 1
## 238 1850.79 1
## 239 1827.38 1
## 240 1836.01 1
## 241 1821.76 1
## 242 1813.71 1
## 243 1806.88 1
## 244 1790.08 1
## 245 1772.01 1
## 246 1785.54 1
# MÃ HÓA VN100 SANG THANG DO ORDINAL SCALE
VN100a=MHVN100
VN100a=replace(VN100a,MHVN100==1,"Nhỏ hơn hoặc bằng 1861")
VN100a=replace(VN100a,MHVN100==2,"Lớn hơn 1861")
VN100aMH=data.frame(MHVN100,VN100a)
VN100aMH
## MHVN100 VN100a
## 1 1 Nhỏ hơn hoặc bằng 1861
## 2 1 Nhỏ hơn hoặc bằng 1861
## 3 1 Nhỏ hơn hoặc bằng 1861
## 4 1 Nhỏ hơn hoặc bằng 1861
## 5 1 Nhỏ hơn hoặc bằng 1861
## 6 2 Lớn hơn 1861
## 7 1 Nhỏ hơn hoặc bằng 1861
## 8 2 Lớn hơn 1861
## 9 2 Lớn hơn 1861
## 10 2 Lớn hơn 1861
## 11 2 Lớn hơn 1861
## 12 2 Lớn hơn 1861
## 13 2 Lớn hơn 1861
## 14 2 Lớn hơn 1861
## 15 2 Lớn hơn 1861
## 16 1 Nhỏ hơn hoặc bằng 1861
## 17 1 Nhỏ hơn hoặc bằng 1861
## 18 1 Nhỏ hơn hoặc bằng 1861
## 19 1 Nhỏ hơn hoặc bằng 1861
## 20 1 Nhỏ hơn hoặc bằng 1861
## 21 1 Nhỏ hơn hoặc bằng 1861
## 22 1 Nhỏ hơn hoặc bằng 1861
## 23 1 Nhỏ hơn hoặc bằng 1861
## 24 1 Nhỏ hơn hoặc bằng 1861
## 25 1 Nhỏ hơn hoặc bằng 1861
## 26 1 Nhỏ hơn hoặc bằng 1861
## 27 1 Nhỏ hơn hoặc bằng 1861
## 28 1 Nhỏ hơn hoặc bằng 1861
## 29 1 Nhỏ hơn hoặc bằng 1861
## 30 1 Nhỏ hơn hoặc bằng 1861
## 31 1 Nhỏ hơn hoặc bằng 1861
## 32 1 Nhỏ hơn hoặc bằng 1861
## 33 1 Nhỏ hơn hoặc bằng 1861
## 34 1 Nhỏ hơn hoặc bằng 1861
## 35 1 Nhỏ hơn hoặc bằng 1861
## 36 1 Nhỏ hơn hoặc bằng 1861
## 37 1 Nhỏ hơn hoặc bằng 1861
## 38 1 Nhỏ hơn hoặc bằng 1861
## 39 1 Nhỏ hơn hoặc bằng 1861
## 40 1 Nhỏ hơn hoặc bằng 1861
## 41 1 Nhỏ hơn hoặc bằng 1861
## 42 1 Nhỏ hơn hoặc bằng 1861
## 43 1 Nhỏ hơn hoặc bằng 1861
## 44 1 Nhỏ hơn hoặc bằng 1861
## 45 1 Nhỏ hơn hoặc bằng 1861
## 46 1 Nhỏ hơn hoặc bằng 1861
## 47 2 Lớn hơn 1861
## 48 2 Lớn hơn 1861
## 49 2 Lớn hơn 1861
## 50 2 Lớn hơn 1861
## 51 2 Lớn hơn 1861
## 52 1 Nhỏ hơn hoặc bằng 1861
## 53 1 Nhỏ hơn hoặc bằng 1861
## 54 1 Nhỏ hơn hoặc bằng 1861
## 55 1 Nhỏ hơn hoặc bằng 1861
## 56 1 Nhỏ hơn hoặc bằng 1861
## 57 1 Nhỏ hơn hoặc bằng 1861
## 58 1 Nhỏ hơn hoặc bằng 1861
## 59 1 Nhỏ hơn hoặc bằng 1861
## 60 2 Lớn hơn 1861
## 61 2 Lớn hơn 1861
## 62 2 Lớn hơn 1861
## 63 2 Lớn hơn 1861
## 64 2 Lớn hơn 1861
## 65 2 Lớn hơn 1861
## 66 2 Lớn hơn 1861
## 67 2 Lớn hơn 1861
## 68 2 Lớn hơn 1861
## 69 2 Lớn hơn 1861
## 70 2 Lớn hơn 1861
## 71 2 Lớn hơn 1861
## 72 2 Lớn hơn 1861
## 73 2 Lớn hơn 1861
## 74 2 Lớn hơn 1861
## 75 2 Lớn hơn 1861
## 76 2 Lớn hơn 1861
## 77 2 Lớn hơn 1861
## 78 2 Lớn hơn 1861
## 79 2 Lớn hơn 1861
## 80 2 Lớn hơn 1861
## 81 2 Lớn hơn 1861
## 82 2 Lớn hơn 1861
## 83 2 Lớn hơn 1861
## 84 2 Lớn hơn 1861
## 85 2 Lớn hơn 1861
## 86 2 Lớn hơn 1861
## 87 2 Lớn hơn 1861
## 88 2 Lớn hơn 1861
## 89 2 Lớn hơn 1861
## 90 1 Nhỏ hơn hoặc bằng 1861
## 91 1 Nhỏ hơn hoặc bằng 1861
## 92 1 Nhỏ hơn hoặc bằng 1861
## 93 1 Nhỏ hơn hoặc bằng 1861
## 94 2 Lớn hơn 1861
## 95 2 Lớn hơn 1861
## 96 2 Lớn hơn 1861
## 97 2 Lớn hơn 1861
## 98 2 Lớn hơn 1861
## 99 2 Lớn hơn 1861
## 100 2 Lớn hơn 1861
## 101 2 Lớn hơn 1861
## 102 2 Lớn hơn 1861
## 103 2 Lớn hơn 1861
## 104 1 Nhỏ hơn hoặc bằng 1861
## 105 1 Nhỏ hơn hoặc bằng 1861
## 106 1 Nhỏ hơn hoặc bằng 1861
## 107 1 Nhỏ hơn hoặc bằng 1861
## 108 1 Nhỏ hơn hoặc bằng 1861
## 109 1 Nhỏ hơn hoặc bằng 1861
## 110 1 Nhỏ hơn hoặc bằng 1861
## 111 1 Nhỏ hơn hoặc bằng 1861
## 112 1 Nhỏ hơn hoặc bằng 1861
## 113 1 Nhỏ hơn hoặc bằng 1861
## 114 1 Nhỏ hơn hoặc bằng 1861
## 115 1 Nhỏ hơn hoặc bằng 1861
## 116 1 Nhỏ hơn hoặc bằng 1861
## 117 1 Nhỏ hơn hoặc bằng 1861
## 118 1 Nhỏ hơn hoặc bằng 1861
## 119 1 Nhỏ hơn hoặc bằng 1861
## 120 1 Nhỏ hơn hoặc bằng 1861
## 121 1 Nhỏ hơn hoặc bằng 1861
## 122 1 Nhỏ hơn hoặc bằng 1861
## 123 1 Nhỏ hơn hoặc bằng 1861
## 124 1 Nhỏ hơn hoặc bằng 1861
## 125 1 Nhỏ hơn hoặc bằng 1861
## 126 1 Nhỏ hơn hoặc bằng 1861
## 127 1 Nhỏ hơn hoặc bằng 1861
## 128 1 Nhỏ hơn hoặc bằng 1861
## 129 1 Nhỏ hơn hoặc bằng 1861
## 130 1 Nhỏ hơn hoặc bằng 1861
## 131 1 Nhỏ hơn hoặc bằng 1861
## 132 2 Lớn hơn 1861
## 133 2 Lớn hơn 1861
## 134 2 Lớn hơn 1861
## 135 2 Lớn hơn 1861
## 136 2 Lớn hơn 1861
## 137 2 Lớn hơn 1861
## 138 2 Lớn hơn 1861
## 139 2 Lớn hơn 1861
## 140 2 Lớn hơn 1861
## 141 2 Lớn hơn 1861
## 142 2 Lớn hơn 1861
## 143 2 Lớn hơn 1861
## 144 2 Lớn hơn 1861
## 145 2 Lớn hơn 1861
## 146 2 Lớn hơn 1861
## 147 2 Lớn hơn 1861
## 148 2 Lớn hơn 1861
## 149 2 Lớn hơn 1861
## 150 2 Lớn hơn 1861
## 151 2 Lớn hơn 1861
## 152 2 Lớn hơn 1861
## 153 2 Lớn hơn 1861
## 154 2 Lớn hơn 1861
## 155 2 Lớn hơn 1861
## 156 2 Lớn hơn 1861
## 157 2 Lớn hơn 1861
## 158 2 Lớn hơn 1861
## 159 2 Lớn hơn 1861
## 160 2 Lớn hơn 1861
## 161 2 Lớn hơn 1861
## 162 2 Lớn hơn 1861
## 163 2 Lớn hơn 1861
## 164 2 Lớn hơn 1861
## 165 2 Lớn hơn 1861
## 166 2 Lớn hơn 1861
## 167 1 Nhỏ hơn hoặc bằng 1861
## 168 2 Lớn hơn 1861
## 169 2 Lớn hơn 1861
## 170 2 Lớn hơn 1861
## 171 2 Lớn hơn 1861
## 172 2 Lớn hơn 1861
## 173 2 Lớn hơn 1861
## 174 2 Lớn hơn 1861
## 175 2 Lớn hơn 1861
## 176 2 Lớn hơn 1861
## 177 2 Lớn hơn 1861
## 178 2 Lớn hơn 1861
## 179 2 Lớn hơn 1861
## 180 2 Lớn hơn 1861
## 181 2 Lớn hơn 1861
## 182 2 Lớn hơn 1861
## 183 2 Lớn hơn 1861
## 184 2 Lớn hơn 1861
## 185 2 Lớn hơn 1861
## 186 2 Lớn hơn 1861
## 187 2 Lớn hơn 1861
## 188 2 Lớn hơn 1861
## 189 2 Lớn hơn 1861
## 190 2 Lớn hơn 1861
## 191 2 Lớn hơn 1861
## 192 1 Nhỏ hơn hoặc bằng 1861
## 193 2 Lớn hơn 1861
## 194 1 Nhỏ hơn hoặc bằng 1861
## 195 1 Nhỏ hơn hoặc bằng 1861
## 196 1 Nhỏ hơn hoặc bằng 1861
## 197 1 Nhỏ hơn hoặc bằng 1861
## 198 1 Nhỏ hơn hoặc bằng 1861
## 199 1 Nhỏ hơn hoặc bằng 1861
## 200 1 Nhỏ hơn hoặc bằng 1861
## 201 1 Nhỏ hơn hoặc bằng 1861
## 202 1 Nhỏ hơn hoặc bằng 1861
## 203 1 Nhỏ hơn hoặc bằng 1861
## 204 1 Nhỏ hơn hoặc bằng 1861
## 205 1 Nhỏ hơn hoặc bằng 1861
## 206 1 Nhỏ hơn hoặc bằng 1861
## 207 1 Nhỏ hơn hoặc bằng 1861
## 208 1 Nhỏ hơn hoặc bằng 1861
## 209 1 Nhỏ hơn hoặc bằng 1861
## 210 1 Nhỏ hơn hoặc bằng 1861
## 211 1 Nhỏ hơn hoặc bằng 1861
## 212 1 Nhỏ hơn hoặc bằng 1861
## 213 1 Nhỏ hơn hoặc bằng 1861
## 214 1 Nhỏ hơn hoặc bằng 1861
## 215 1 Nhỏ hơn hoặc bằng 1861
## 216 1 Nhỏ hơn hoặc bằng 1861
## 217 1 Nhỏ hơn hoặc bằng 1861
## 218 1 Nhỏ hơn hoặc bằng 1861
## 219 1 Nhỏ hơn hoặc bằng 1861
## 220 1 Nhỏ hơn hoặc bằng 1861
## 221 1 Nhỏ hơn hoặc bằng 1861
## 222 2 Lớn hơn 1861
## 223 2 Lớn hơn 1861
## 224 2 Lớn hơn 1861
## 225 2 Lớn hơn 1861
## 226 2 Lớn hơn 1861
## 227 2 Lớn hơn 1861
## 228 2 Lớn hơn 1861
## 229 2 Lớn hơn 1861
## 230 1 Nhỏ hơn hoặc bằng 1861
## 231 1 Nhỏ hơn hoặc bằng 1861
## 232 1 Nhỏ hơn hoặc bằng 1861
## 233 1 Nhỏ hơn hoặc bằng 1861
## 234 2 Lớn hơn 1861
## 235 1 Nhỏ hơn hoặc bằng 1861
## 236 1 Nhỏ hơn hoặc bằng 1861
## 237 1 Nhỏ hơn hoặc bằng 1861
## 238 1 Nhỏ hơn hoặc bằng 1861
## 239 1 Nhỏ hơn hoặc bằng 1861
## 240 1 Nhỏ hơn hoặc bằng 1861
## 241 1 Nhỏ hơn hoặc bằng 1861
## 242 1 Nhỏ hơn hoặc bằng 1861
## 243 1 Nhỏ hơn hoặc bằng 1861
## 244 1 Nhỏ hơn hoặc bằng 1861
## 245 1 Nhỏ hơn hoặc bằng 1861
## 246 1 Nhỏ hơn hoặc bằng 1861
# Mã hóa Dow Jones
MHDJ=Dow.Jones
MHDJ=replace(MHDJ,Dow.Jones<=median(Dow.Jones),1)
MHDJ=replace(MHDJ,Dow.Jones>median(Dow.Jones),2)
DJcode=data.frame(Dow.Jones,MHDJ)
DJcode
## Dow.Jones MHDJ
## 1 46316.07 1
## 2 46397.89 1
## 3 46441.10 1
## 4 46519.72 1
## 5 46758.28 1
## 6 46694.97 1
## 7 46602.98 1
## 8 46601.78 1
## 9 46358.42 1
## 10 45479.60 1
## 11 46067.58 1
## 12 46270.46 1
## 13 46253.31 1
## 14 45952.24 1
## 15 46190.61 1
## 16 46706.58 1
## 17 46924.74 1
## 18 46590.41 1
## 19 46734.61 1
## 20 47207.12 1
## 21 47544.59 1
## 22 47706.37 1
## 23 47632.00 1
## 24 47522.12 1
## 25 47562.87 1
## 26 47336.68 1
## 27 47085.24 1
## 28 47311.00 1
## 29 46912.30 1
## 30 46987.10 1
## 31 47368.63 1
## 32 47927.96 1
## 33 48254.82 1
## 34 47457.22 1
## 35 47147.48 1
## 36 46590.24 1
## 37 46091.74 1
## 38 46138.77 1
## 39 45752.26 1
## 40 46245.41 1
## 41 46448.27 1
## 42 47112.45 1
## 43 47427.12 1
## 44 47716.42 1
## 45 47289.33 1
## 46 47474.46 1
## 47 47882.90 1
## 48 47850.94 1
## 49 47954.99 1
## 50 47739.32 1
## 51 47560.29 1
## 52 48057.75 1
## 53 48704.01 1
## 54 48458.05 1
## 55 48416.56 1
## 56 48114.26 1
## 57 47885.97 1
## 58 47951.85 1
## 59 48134.89 1
## 60 48362.68 1
## 61 48442.41 1
## 62 48731.16 1
## 63 48731.16 1
## 64 48710.97 1
## 65 48461.93 1
## 66 48367.06 1
## 67 48063.29 1
## 68 48977.18 1
## 69 49462.08 2
## 70 48996.08 1
## 71 49266.11 1
## 72 49504.07 2
## 73 49590.20 2
## 74 49191.99 1
## 75 49149.63 1
## 76 49442.44 2
## 77 49359.33 2
## 78 48488.59 1
## 79 49077.23 1
## 80 49384.01 2
## 81 49098.71 1
## 82 49412.40 2
## 83 49003.41 1
## 84 49015.60 1
## 85 49071.56 1
## 86 48892.47 1
## 87 49407.66 2
## 88 49240.99 1
## 89 49501.30 2
## 90 48908.72 1
## 91 50115.67 2
## 92 50135.87 2
## 93 50188.14 2
## 94 50121.40 2
## 95 49451.98 2
## 96 49500.93 2
## 97 48804.06 1
## 98 49174.50 1
## 99 49482.15 2
## 100 49499.20 2
## 101 48977.92 1
## 102 48904.78 1
## 103 48501.27 1
## 104 48739.41 1
## 105 47954.74 1
## 106 47501.55 1
## 107 47740.80 1
## 108 47706.51 1
## 109 47417.27 1
## 110 46677.85 1
## 111 46558.47 1
## 112 46946.41 1
## 113 46993.26 1
## 114 46225.15 1
## 115 46021.43 1
## 116 45577.47 1
## 117 46208.47 1
## 118 46124.06 1
## 119 46429.49 1
## 120 45960.11 1
## 121 45166.64 1
## 122 45216.14 1
## 123 46341.51 1
## 124 46565.74 1
## 125 46504.67 1
## 126 46669.88 1
## 127 46584.46 1
## 128 47909.92 1
## 129 48185.80 1
## 130 47916.57 1
## 131 48218.25 1
## 132 48535.99 1
## 133 48463.72 1
## 134 48578.72 1
## 135 49447.43 2
## 136 49442.56 2
## 137 49149.38 1
## 138 49490.03 2
## 139 49310.32 2
## 140 49230.71 1
## 141 49141.93 1
## 142 48861.81 1
## 143 48941.90 1
## 144 49298.25 1
## 145 49910.59 2
## 146 49596.97 2
## 147 49609.16 2
## 148 49704.47 2
## 149 49760.56 2
## 150 49693.20 2
## 151 50063.46 2
## 152 49526.17 2
## 153 49686.12 2
## 154 49363.88 2
## 155 50009.35 2
## 156 50285.66 2
## 157 50579.70 2
## 158 50461.68 2
## 159 50644.28 2
## 160 50668.97 2
## 161 51032.46 2
## 162 51078.88 2
## 163 51307.79 2
## 164 50687.07 2
## 165 51561.93 2
## 166 50866.78 2
## 167 50786.01 2
## 168 50872.11 2
## 169 49918.78 2
## 170 50848.75 2
## 171 51202.26 2
## 172 51671.03 2
## 173 51999.67 2
## 174 51492.55 2
## 175 51564.70 2
## 176 51712.71 2
## 177 51666.84 2
## 178 51848.90 2
## 179 51920.62 2
## 180 51876.11 2
## 181 52182.74 2
## 182 52319.20 2
## 183 52305.24 2
## 184 52900.07 2
## 185 53055.91 2
## 186 52925.15 2
## 187 52348.39 2
## 188 52487.41 2
## 189 52637.01 2
## 190 52498.64 2
## 191 52508.27 2
## 192 52658.64 2
## 193 52552.97 2
## 194 52146.42 2
## 195 51839.26 2
## 196 52224.64 2
## 197 52218.58 2
## 198 51711.65 2
## 199 51947.25 2
## 200 52210.08 2
## 201 52747.32 2
## 202 51594.14 2
## 203 52208.06 2
## 204 52485.03 2
## 205 53178.41 2
## 206 54085.88 2
## 207 54349.12 2
## 208 53885.10 2
## 209 54036.93 2
## 210 53975.98 2
## 211 53791.85 2
## 212 53770.27 2
## 213 53839.99 2
## 214 53732.41 2
## 215 53459.78 2
## 216 53343.40 2
## 217 53463.05 2
## 218 52759.21 2
## 219 53277.01 2
## 220 53417.16 2
## 221 53577.40 2
## 222 53463.88 2
## 223 53569.44 2
## 224 53559.99 2
## 225 53686.11 2
## 226 53414.25 2
## 227 52786.07 2
## 228 52380.66 2
## 229 52064.10 2
## 230 52573.29 2
## 231 52421.20 2
## 232 52093.11 2
## 233 51461.90 2
## 234 51778.04 2
## 235 51682.64 2
## 236 52048.83 2
## 237 51863.69 2
## 238 51511.59 2
## 239 51349.98 2
## 240 51828.62 2
## 241 51481.51 2
## 242 51349.92 2
## 243 50906.05 2
## 244 50926.56 2
## 245 51176.96 2
## 246 51025.44 2
# MÃ HÓA Dow Jones SANG THANG DO ORDINAL SCALE
Dow.Jones1=MHDJ
Dow.Jones1=replace(Dow.Jones1,MHDJ==1,"Nhỏ hơn hoặc bằng 49304")
Dow.Jones1=replace(Dow.Jones1,MHDJ==2,"Lớn hơn 49304")
Dow.Jones1MH=data.frame(MHDJ,Dow.Jones1)
Dow.Jones1MH
## MHDJ Dow.Jones1
## 1 1 Nhỏ hơn hoặc bằng 49304
## 2 1 Nhỏ hơn hoặc bằng 49304
## 3 1 Nhỏ hơn hoặc bằng 49304
## 4 1 Nhỏ hơn hoặc bằng 49304
## 5 1 Nhỏ hơn hoặc bằng 49304
## 6 1 Nhỏ hơn hoặc bằng 49304
## 7 1 Nhỏ hơn hoặc bằng 49304
## 8 1 Nhỏ hơn hoặc bằng 49304
## 9 1 Nhỏ hơn hoặc bằng 49304
## 10 1 Nhỏ hơn hoặc bằng 49304
## 11 1 Nhỏ hơn hoặc bằng 49304
## 12 1 Nhỏ hơn hoặc bằng 49304
## 13 1 Nhỏ hơn hoặc bằng 49304
## 14 1 Nhỏ hơn hoặc bằng 49304
## 15 1 Nhỏ hơn hoặc bằng 49304
## 16 1 Nhỏ hơn hoặc bằng 49304
## 17 1 Nhỏ hơn hoặc bằng 49304
## 18 1 Nhỏ hơn hoặc bằng 49304
## 19 1 Nhỏ hơn hoặc bằng 49304
## 20 1 Nhỏ hơn hoặc bằng 49304
## 21 1 Nhỏ hơn hoặc bằng 49304
## 22 1 Nhỏ hơn hoặc bằng 49304
## 23 1 Nhỏ hơn hoặc bằng 49304
## 24 1 Nhỏ hơn hoặc bằng 49304
## 25 1 Nhỏ hơn hoặc bằng 49304
## 26 1 Nhỏ hơn hoặc bằng 49304
## 27 1 Nhỏ hơn hoặc bằng 49304
## 28 1 Nhỏ hơn hoặc bằng 49304
## 29 1 Nhỏ hơn hoặc bằng 49304
## 30 1 Nhỏ hơn hoặc bằng 49304
## 31 1 Nhỏ hơn hoặc bằng 49304
## 32 1 Nhỏ hơn hoặc bằng 49304
## 33 1 Nhỏ hơn hoặc bằng 49304
## 34 1 Nhỏ hơn hoặc bằng 49304
## 35 1 Nhỏ hơn hoặc bằng 49304
## 36 1 Nhỏ hơn hoặc bằng 49304
## 37 1 Nhỏ hơn hoặc bằng 49304
## 38 1 Nhỏ hơn hoặc bằng 49304
## 39 1 Nhỏ hơn hoặc bằng 49304
## 40 1 Nhỏ hơn hoặc bằng 49304
## 41 1 Nhỏ hơn hoặc bằng 49304
## 42 1 Nhỏ hơn hoặc bằng 49304
## 43 1 Nhỏ hơn hoặc bằng 49304
## 44 1 Nhỏ hơn hoặc bằng 49304
## 45 1 Nhỏ hơn hoặc bằng 49304
## 46 1 Nhỏ hơn hoặc bằng 49304
## 47 1 Nhỏ hơn hoặc bằng 49304
## 48 1 Nhỏ hơn hoặc bằng 49304
## 49 1 Nhỏ hơn hoặc bằng 49304
## 50 1 Nhỏ hơn hoặc bằng 49304
## 51 1 Nhỏ hơn hoặc bằng 49304
## 52 1 Nhỏ hơn hoặc bằng 49304
## 53 1 Nhỏ hơn hoặc bằng 49304
## 54 1 Nhỏ hơn hoặc bằng 49304
## 55 1 Nhỏ hơn hoặc bằng 49304
## 56 1 Nhỏ hơn hoặc bằng 49304
## 57 1 Nhỏ hơn hoặc bằng 49304
## 58 1 Nhỏ hơn hoặc bằng 49304
## 59 1 Nhỏ hơn hoặc bằng 49304
## 60 1 Nhỏ hơn hoặc bằng 49304
## 61 1 Nhỏ hơn hoặc bằng 49304
## 62 1 Nhỏ hơn hoặc bằng 49304
## 63 1 Nhỏ hơn hoặc bằng 49304
## 64 1 Nhỏ hơn hoặc bằng 49304
## 65 1 Nhỏ hơn hoặc bằng 49304
## 66 1 Nhỏ hơn hoặc bằng 49304
## 67 1 Nhỏ hơn hoặc bằng 49304
## 68 1 Nhỏ hơn hoặc bằng 49304
## 69 2 Lớn hơn 49304
## 70 1 Nhỏ hơn hoặc bằng 49304
## 71 1 Nhỏ hơn hoặc bằng 49304
## 72 2 Lớn hơn 49304
## 73 2 Lớn hơn 49304
## 74 1 Nhỏ hơn hoặc bằng 49304
## 75 1 Nhỏ hơn hoặc bằng 49304
## 76 2 Lớn hơn 49304
## 77 2 Lớn hơn 49304
## 78 1 Nhỏ hơn hoặc bằng 49304
## 79 1 Nhỏ hơn hoặc bằng 49304
## 80 2 Lớn hơn 49304
## 81 1 Nhỏ hơn hoặc bằng 49304
## 82 2 Lớn hơn 49304
## 83 1 Nhỏ hơn hoặc bằng 49304
## 84 1 Nhỏ hơn hoặc bằng 49304
## 85 1 Nhỏ hơn hoặc bằng 49304
## 86 1 Nhỏ hơn hoặc bằng 49304
## 87 2 Lớn hơn 49304
## 88 1 Nhỏ hơn hoặc bằng 49304
## 89 2 Lớn hơn 49304
## 90 1 Nhỏ hơn hoặc bằng 49304
## 91 2 Lớn hơn 49304
## 92 2 Lớn hơn 49304
## 93 2 Lớn hơn 49304
## 94 2 Lớn hơn 49304
## 95 2 Lớn hơn 49304
## 96 2 Lớn hơn 49304
## 97 1 Nhỏ hơn hoặc bằng 49304
## 98 1 Nhỏ hơn hoặc bằng 49304
## 99 2 Lớn hơn 49304
## 100 2 Lớn hơn 49304
## 101 1 Nhỏ hơn hoặc bằng 49304
## 102 1 Nhỏ hơn hoặc bằng 49304
## 103 1 Nhỏ hơn hoặc bằng 49304
## 104 1 Nhỏ hơn hoặc bằng 49304
## 105 1 Nhỏ hơn hoặc bằng 49304
## 106 1 Nhỏ hơn hoặc bằng 49304
## 107 1 Nhỏ hơn hoặc bằng 49304
## 108 1 Nhỏ hơn hoặc bằng 49304
## 109 1 Nhỏ hơn hoặc bằng 49304
## 110 1 Nhỏ hơn hoặc bằng 49304
## 111 1 Nhỏ hơn hoặc bằng 49304
## 112 1 Nhỏ hơn hoặc bằng 49304
## 113 1 Nhỏ hơn hoặc bằng 49304
## 114 1 Nhỏ hơn hoặc bằng 49304
## 115 1 Nhỏ hơn hoặc bằng 49304
## 116 1 Nhỏ hơn hoặc bằng 49304
## 117 1 Nhỏ hơn hoặc bằng 49304
## 118 1 Nhỏ hơn hoặc bằng 49304
## 119 1 Nhỏ hơn hoặc bằng 49304
## 120 1 Nhỏ hơn hoặc bằng 49304
## 121 1 Nhỏ hơn hoặc bằng 49304
## 122 1 Nhỏ hơn hoặc bằng 49304
## 123 1 Nhỏ hơn hoặc bằng 49304
## 124 1 Nhỏ hơn hoặc bằng 49304
## 125 1 Nhỏ hơn hoặc bằng 49304
## 126 1 Nhỏ hơn hoặc bằng 49304
## 127 1 Nhỏ hơn hoặc bằng 49304
## 128 1 Nhỏ hơn hoặc bằng 49304
## 129 1 Nhỏ hơn hoặc bằng 49304
## 130 1 Nhỏ hơn hoặc bằng 49304
## 131 1 Nhỏ hơn hoặc bằng 49304
## 132 1 Nhỏ hơn hoặc bằng 49304
## 133 1 Nhỏ hơn hoặc bằng 49304
## 134 1 Nhỏ hơn hoặc bằng 49304
## 135 2 Lớn hơn 49304
## 136 2 Lớn hơn 49304
## 137 1 Nhỏ hơn hoặc bằng 49304
## 138 2 Lớn hơn 49304
## 139 2 Lớn hơn 49304
## 140 1 Nhỏ hơn hoặc bằng 49304
## 141 1 Nhỏ hơn hoặc bằng 49304
## 142 1 Nhỏ hơn hoặc bằng 49304
## 143 1 Nhỏ hơn hoặc bằng 49304
## 144 1 Nhỏ hơn hoặc bằng 49304
## 145 2 Lớn hơn 49304
## 146 2 Lớn hơn 49304
## 147 2 Lớn hơn 49304
## 148 2 Lớn hơn 49304
## 149 2 Lớn hơn 49304
## 150 2 Lớn hơn 49304
## 151 2 Lớn hơn 49304
## 152 2 Lớn hơn 49304
## 153 2 Lớn hơn 49304
## 154 2 Lớn hơn 49304
## 155 2 Lớn hơn 49304
## 156 2 Lớn hơn 49304
## 157 2 Lớn hơn 49304
## 158 2 Lớn hơn 49304
## 159 2 Lớn hơn 49304
## 160 2 Lớn hơn 49304
## 161 2 Lớn hơn 49304
## 162 2 Lớn hơn 49304
## 163 2 Lớn hơn 49304
## 164 2 Lớn hơn 49304
## 165 2 Lớn hơn 49304
## 166 2 Lớn hơn 49304
## 167 2 Lớn hơn 49304
## 168 2 Lớn hơn 49304
## 169 2 Lớn hơn 49304
## 170 2 Lớn hơn 49304
## 171 2 Lớn hơn 49304
## 172 2 Lớn hơn 49304
## 173 2 Lớn hơn 49304
## 174 2 Lớn hơn 49304
## 175 2 Lớn hơn 49304
## 176 2 Lớn hơn 49304
## 177 2 Lớn hơn 49304
## 178 2 Lớn hơn 49304
## 179 2 Lớn hơn 49304
## 180 2 Lớn hơn 49304
## 181 2 Lớn hơn 49304
## 182 2 Lớn hơn 49304
## 183 2 Lớn hơn 49304
## 184 2 Lớn hơn 49304
## 185 2 Lớn hơn 49304
## 186 2 Lớn hơn 49304
## 187 2 Lớn hơn 49304
## 188 2 Lớn hơn 49304
## 189 2 Lớn hơn 49304
## 190 2 Lớn hơn 49304
## 191 2 Lớn hơn 49304
## 192 2 Lớn hơn 49304
## 193 2 Lớn hơn 49304
## 194 2 Lớn hơn 49304
## 195 2 Lớn hơn 49304
## 196 2 Lớn hơn 49304
## 197 2 Lớn hơn 49304
## 198 2 Lớn hơn 49304
## 199 2 Lớn hơn 49304
## 200 2 Lớn hơn 49304
## 201 2 Lớn hơn 49304
## 202 2 Lớn hơn 49304
## 203 2 Lớn hơn 49304
## 204 2 Lớn hơn 49304
## 205 2 Lớn hơn 49304
## 206 2 Lớn hơn 49304
## 207 2 Lớn hơn 49304
## 208 2 Lớn hơn 49304
## 209 2 Lớn hơn 49304
## 210 2 Lớn hơn 49304
## 211 2 Lớn hơn 49304
## 212 2 Lớn hơn 49304
## 213 2 Lớn hơn 49304
## 214 2 Lớn hơn 49304
## 215 2 Lớn hơn 49304
## 216 2 Lớn hơn 49304
## 217 2 Lớn hơn 49304
## 218 2 Lớn hơn 49304
## 219 2 Lớn hơn 49304
## 220 2 Lớn hơn 49304
## 221 2 Lớn hơn 49304
## 222 2 Lớn hơn 49304
## 223 2 Lớn hơn 49304
## 224 2 Lớn hơn 49304
## 225 2 Lớn hơn 49304
## 226 2 Lớn hơn 49304
## 227 2 Lớn hơn 49304
## 228 2 Lớn hơn 49304
## 229 2 Lớn hơn 49304
## 230 2 Lớn hơn 49304
## 231 2 Lớn hơn 49304
## 232 2 Lớn hơn 49304
## 233 2 Lớn hơn 49304
## 234 2 Lớn hơn 49304
## 235 2 Lớn hơn 49304
## 236 2 Lớn hơn 49304
## 237 2 Lớn hơn 49304
## 238 2 Lớn hơn 49304
## 239 2 Lớn hơn 49304
## 240 2 Lớn hơn 49304
## 241 2 Lớn hơn 49304
## 242 2 Lớn hơn 49304
## 243 2 Lớn hơn 49304
## 244 2 Lớn hơn 49304
## 245 2 Lớn hơn 49304
## 246 2 Lớn hơn 49304
#c) Thống kê các chỉ số
# Thống kê VN30
summary(VN30)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1741 1898 1951 1950 2006 2097
var(VN30)
## [1] 5236.864
sd(VN30)
## [1] 72.36618
# Thống kê VN100
summary(VN100)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1665 1808 1861 1861 1915 1997
var(VN100)
## [1] 4643.104
sd(VN100)
## [1] 68.14032
# Thống kê Dow Jones
summary(Dow.Jones)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 45167 47580 49304 49547 51670 54349
var(Dow.Jones)
## [1] 5572499
sd(Dow.Jones)
## [1] 2360.614
library(psych)
# Thống kê Ratio Scale
TKRS=data.frame(VN30,VN100,Dow.Jones)
describe(TKRS)
## vars n mean sd median trimmed mad min max
## VN30 1 246 1949.54 72.37 1950.63 1950.92 79.57 1741.05 2096.76
## VN100 2 246 1860.51 68.14 1860.99 1861.81 79.67 1664.96 1997.02
## Dow.Jones 3 246 49546.70 2360.61 49304.29 49482.34 2962.82 45166.64 54349.12
## range skew kurtosis se
## VN30 355.71 -0.20 -0.50 4.61
## VN100 332.06 -0.15 -0.63 4.34
## Dow.Jones 9182.48 0.19 -1.08 150.51
#Nhận xét: cả 3 chỉ số đều có n = 246
#Đối với chỉ số VN30: Giá trị trung bình (mean) là 1949.54 gần tương đương với trung vị (median) là 1950.63 cho thấy phân phối tương đối đối xứng; độ lệch chuẩn (sd) là 72.37 với khoảng biến thiên (range) là 355.71 cho thấy độ phân tán ở mức thấp
#Đối với chỉ số VN100: Giá trị trung bình (mean) là 1860.51 gần tương đương với trung vị (median) là 1860.99 cho thấy phân phối tương đối đối xứng; độ lệch chuẩn (sd) là 68.14 với khoảng biến thiên (range) là 332.06 cho thấy độ phân tán ở mức thấp
#Đối với chỉ số Dow.Jones: Giá trị trung bình (mean) là 49546.70 gần tương đương với trung vị (median) là 49304.29 cho thấy phân phối tương đối đối xứng; độ lệch chuẩn (sd) là 2360.61 với khoảng biến thiên (range) là 9182.48 cho thấy độ phân tán ở mức cao
# Thống kê Ordinal Scale
TKOS=data.frame(MHVN30,MHVN100,MHDJ)
describe(TKOS)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## MHVN30 1 246 1.5 0.5 1.5 1.5 0.74 1 2 1 0 -2.01 0.03
## MHVN100 2 246 1.5 0.5 1.5 1.5 0.74 1 2 1 0 -2.01 0.03
## MHDJ 3 246 1.5 0.5 1.5 1.5 0.74 1 2 1 0 -2.01 0.03
#Nhận xét: Cả ba chỉ số đều có n = 246 mẫu
# Đối với chỉ số MHVN30: Giá trị trung bình (mean) là 1.5 gần tương đương với trung vị (median) là 1.5 cho thấy phân phối tương đối đối xứng; độ lệch chuẩn (sd) là 0.5 với khoảng biến thiên (range) là 1
# Đối với chỉ số MHVN100: Giá trị trung bình (mean) là 1.5 gần tương đương với trung vị (median) là 1.5 cho thấy phân phối tương đối đối xứng; độ lệch chuẩn (sd) là 0.5 với khoảng biến thiên (range) là 1
# Đối với chỉ số MHDJ: Giá trị trung bình (mean) là 1.5 gần tương đương với trung vị (median) là 1.5 cho thấy phân phối tương đối đối xứng; độ lệch chuẩn (sd) là 0.5 với khoảng biến thiên (range) là 1
#d) Vẽ đồ thị dưới dạng Ordinal Scale
# Đồ thị cột VN30
VN30.freq<-table(VN30aMH)
VN30.freq
## VN30a
## MHVN30 Lớn hơn 1951 Nhỏ hơn hoặc bằng 1951
## 1 0 123
## 2 123 0
barplot(VN30.freq,xlab="Mức VN30",ylab="Số lượng",main="Đồ thị cột VN30",col=c("lightgreen","lightblue"))
#Nhận xét: Nhóm "Lớn hơn 1951" có 123 quan sát, nhóm "Nhỏ hơn hoặc bằng 1951" có 123 quan sát.
# Đồ thị cột VN100
VN100.freq<-table(VN100aMH)
VN100.freq
## VN100a
## MHVN100 Lớn hơn 1861 Nhỏ hơn hoặc bằng 1861
## 1 0 123
## 2 123 0
barplot(VN100.freq,xlab="Mức VN100",ylab="Số lượng",main="Đồ thị cột VN100",col=c("lightyellow","lightpink"))
#Nhận xét: Nhóm "Lớn hơn 1861" có 123 quan sát, nhóm "Nhỏ hơn hoặc bằng 1861" có 123 quan sát.
# Đồ thị cột Dow Jones
DJ.freq<-table(Dow.Jones1MH)
DJ.freq
## Dow.Jones1
## MHDJ Lớn hơn 49304 Nhỏ hơn hoặc bằng 49304
## 1 0 123
## 2 123 0
barplot(DJ.freq,xlab="Mức Dow Jones",ylab="Số lượng",main="Đồ thị cột Dow Jones",col=c("orange","brown"))
#Nhận xét: Nhóm "Lớn hơn 49304" có 123 quan sát, nhóm "Nhỏ hơn hoặc bằng 49304" có 123 quan sát.
# Đồ thị tròn VN30
pie(table(VN30a),col=c("yellow","blue"),main="Đồ thị tròn VN30")
#Nhận xét: Nhóm "Lớn hơn 1951" và nhóm "Nhỏ hơn hoặc bằng 1951" chiếm tỷ trọng bằng nhau (mỗi nhóm chiếm 50% tổng số quan sát)
# Đồ thị tròn VN100
pie(table(VN100a),col=c("lightgreen","pink"),main="Đồ thị tròn VN100")
#Nhận xét: Nhóm "Lớn hơn 1861" và nhóm "Nhỏ hơn hoặc bằng 1861" chiếm tỷ trọng bằng nhau (mỗi nhóm chiếm 50% tổng số quan sát)
# Đồ thị tròn Dow Jones
pie(table(Dow.Jones1),col=c("red","yellow"),main="Đồ thị tròn Dow Jones")
#Nhận xét: Nhóm "Lớn hơn 49304" và nhóm "Nhỏ hơn hoặc bằng 49304" chiếm tỷ trọng bằng nhau (mỗi nhóm chiếm 50% tổng số quan sát)
# Đồ thị đường dưới dạng Ratio Scale
plot(STT,VN30,type="o",xlab="Ngày thứ",ylab="Giá trị VN30",pch=19,main="Giá trị VN30 theo ngày")
#Nhận xét: Biến VN30 có sự biến động liên tục qua 246 ngày, dao động trong khoảng từ ~1.740 đến đỉnh cao nhất ~2.100, với các nhịp tăng trưởng xen kẽ những đợt điều chỉnh giảm mạnh như đợt sụt giảm sâu xuống ~1.740 ở khoảng ngày 115-120, và đạt các mức đỉnh quan trọng ở khoảng ngày 70-80 và ngày 150
plot
## function (x, y, ...)
## UseMethod("plot")
## <bytecode: 0x0000017e8fe9d8c0>
## <environment: namespace:base>
#Nhận xét: Biến VN100 có sự biến động liên tục qua 246 ngày, dao động trong khoảng từ ~1.660 đến đỉnh cao nhất gần 2.000 (đạt được ở khoảng ngày 145-150), với các nhịp tăng trưởng xen kẽ những đợt điều chỉnh giảm mạnh như đợt sụt giảm sâu xuống mức thấp nhất ~1.660 ở khoảng ngày 115-120, thể hiện rõ xu hướng trồi sụt theo từng chu kỳ giao dịch
plot(STT,Dow.Jones,type="o",xlab="Ngày thứ",ylab="Giá trị Dow Jones",pch=19,main="Giá trị Dow Jones theo ngày")
#Nhận xét: Biến Dow Jones có sự biến động qua 246 ngày, dao động trong khoảng từ ~45.000 đến đỉnh cao nhất vượt ngưỡng 54.000 (tại khoảng ngày 200-210), trải qua nhịp giảm mạnh xuống vùng đáy ~45.000 ở khoảng ngày 120 trước khi bật tăng mạnh mẽ theo xu hướng đi lên và có nhịp điều chỉnh nhẹ về cuối giai đoạn
#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(VN30,VN100,Dow.Jones)
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
## VN30 VN100 Dow.Jones
## VN30 1.0000000 0.9665728 0.2052141
## VN100 0.9665728 1.0000000 0.2024598
## Dow.Jones 0.2052141 0.2024598 1.0000000
result$P
## VN30 VN100 Dow.Jones
## VN30 NA 0.000000000 0.001208861
## VN100 0.000000000 NA 0.001410908
## Dow.Jones 0.001208861 0.001410908 NA
#Nhận xét
#VN30 với VN100: Biến VN30 có biến tương quan dương với VN100 vì hệ số Pearson > 0 (r = 0.97); Mối tương quan giữa VN30 với VN100 có ý nghĩa thống kê vì Sig bé hơn 5%
#VN30 với Dow.Jones: Biến VN30 có biến tương quan dương với Dow.Jones vì hệ số Pearson > 0 (r = 0.21); Mối tương quan giữa VN30 với Dow.Jones có ý nghĩa thống kê vì Sig bé hơn 5%
#VN100 với Dow.Jones: Biến VN100 có biến tương quan dương với Dow.Jones vì hệ số Pearson > 0 (r = 0.20); Mối tương quan giữa VN100 với Dow.Jones có ý nghĩa thống kê vì Sig bé hơn 5%
#f) Mô hình hồi quy tính tuyến đa biến
mreg=lm(VN30~VN100+Dow.Jones)
summary(mreg)
##
## Call:
## lm(formula = VN30 ~ VN100 + Dow.Jones)
##
## Residuals:
## Min 1Q Median 3Q Max
## -45.218 -13.329 -2.694 13.351 47.755
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.859e+01 3.753e+01 0.762 0.447
## VN100 1.024e+00 1.782e-02 57.472 <2e-16 ***
## Dow.Jones 3.044e-04 5.145e-04 0.592 0.555
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 18.62 on 243 degrees of freedom
## Multiple R-squared: 0.9344, Adjusted R-squared: 0.9338
## F-statistic: 1729 on 2 and 243 DF, p-value: < 2.2e-16
anova(mreg)
## Analysis of Variance Table
##
## Response: VN30
## Df Sum Sq Mean Sq F value Pr(>F)
## VN100 1 1198689 1198689 3458.52 <2e-16 ***
## Dow.Jones 1 121 121 0.35 0.5547
## Residuals 243 84221 347
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Nhận xét:
#Phương trình hồi quy: VN30 = 28.59 + 1.024*VN100 + 0.0003*Dow.Jones + ei
#Kết quả hồi quy cho thấy mô hình có ý nghĩa thống kê tổng thể với giá trị F-statistic bằng 1729 và mức ý nghĩa (p < 0,001).
#Hệ số R2 = 0,9344; cho thấy các biến VN100 và Dow.Jones giải thích khoảng 93,44% sự biến thiên của VN30, trong khi R2 hiệu chỉnh đạt 0,9338.
#Xét riêng từng biến, VN100 có tác động cùng chiều đến VN30 với hệ số hồi quy 1.024, điều này có nghĩa khi VN100 tăng 1 đơn vị thì VN30 tăng 1.024 đơn vị và mối quan hệ này có ý nghĩa thống kê rất chặt chẽ ở mức 1% (p < 0,001).
#Riêng biến Dow.Jones có p-value = 0.555 (> 0,05) nên không đạt ý nghĩa thống kê trong mô hình hồi quy này.
#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 Dow.Jones
## 1.042742 1.042742
# 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.054688, 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 tự tương quan cho thấy giá trị DW = 0.054688 không nằm trong khoảng chấp nhận (1,5 < DW < 2,5) => Mô hình có hiện tượng tự tương quan
#Đồ thị bên trái dòng 1, vẽ phần dư ei và giá trị dự báo: Đồ thị này cho thấy các giá trị phần dư phân tán tập trung quanh đường y = 0 cho nên giả định về giá trị trung bình của ei = 0 là có thể chấp nhận được
#Đồ thị bên phải dòng 1, vẽ giá trị phần dư và giá trị lý vọng dựa vào phân phối chuẩn: Kết quả cho thấy các điểm phần dư tập trung rất gần các giá trị trên đường chuẩn tuy nhiên có sự lệch nhẹ giữa 2 đầu, do đó phần dư ei tuân thủ theo quy luật phân phối chuẩn
#Đồ thị bên trái dòng 2, vẽ căn bậc hai của sai số chuẩn phần dư và các giá trị: Đồ thị này cho thấy độ phân tán thay đổi theo các giá trị dự báo, do đó phần phương sai tuân thủ theo quy luật phân phối chuẩn
#Đồ thị bên phải dòng 2, đa số các quan sát có độ đòn bẩy nhỏ, kết quả cho thấy các giá trị đều nằm trong giới hạn đường ranh giới khoảng cách Cook và có rất ít giá trị ngoại lai vượt qua đường ranh giới đó, do đó mô hình ổn định và chi phối