# NHẬP LIỆU
likert=file.choose()
datalikert=read.csv(likert,header=TRUE)
save(datalikert,file="datalikert.rda")
attach(datalikert)
is.data.frame(datalikert)
## [1] TRUE
### CÂU B
MHVIXtrungbinh=datalikert$VIX
MHVIXtrungbinh=replace(MHVIXtrungbinh,datalikert$VIX<=14,1)
MHVIXtrungbinh=replace(MHVIXtrungbinh,datalikert$VIX>14&datalikert$VIX<=17,2)
MHVIXtrungbinh=replace(MHVIXtrungbinh,datalikert$VIX>17,3)
VIX=data.frame(datalikert$VIX,MHVIXtrungbinh)
VIX
## datalikert.VIX MHVIXtrungbinh
## 1 12,650.00 1
## 2 12,750.00 1
## 3 12,950.00 1
## 4 13,050.00 1
## 5 12,800.00 1
## 6 13,150.00 1
## 7 12,850.00 1
## 8 12,650.00 1
## 9 12,900.00 1
## 10 12,600.00 1
## 11 13,250.00 1
## 12 13,800.00 1
## 13 13,600.00 1
## 14 13,750.00 1
## 15 13,600.00 1
## 16 14,000.00 2
## 17 14,050.00 2
## 18 14,100.00 2
## 19 14,250.00 2
## 20 14,350.00 2
## 21 14,100.00 2
## 22 13,800.00 1
## 23 13,500.00 1
## 24 12,650.00 1
## 25 12,476.20 1
## 26 12,523.80 1
## 27 12,714.30 1
## 28 12,809.50 1
## 29 13,238.10 1
## 30 13,619.00 1
## 31 13,523.80 1
## 32 13,285.70 1
## 33 12,952.40 1
## 34 13,047.60 1
## 35 13,333.30 1
## 36 13,238.10 1
## 37 12,761.90 1
## 38 12,381.00 1
## 39 12,428.60 1
## 40 11,619.00 1
## 41 11,714.30 1
## 42 11,190.50 1
## 43 11,857.10 1
## 44 12,285.70 1
## 45 11,523.80 1
## 46 11,904.80 1
## 47 12,190.50 1
## 48 13,095.20 1
## 49 13,381.00 1
## 50 13,285.70 1
## 51 14,000.00 2
## 52 13,952.40 1
## 53 14,714.30 2
## 54 15,285.70 2
## 55 15,714.30 2
## 56 15,714.30 2
## 57 15,523.80 2
## 58 16,238.10 2
## 59 16,095.20 2
## 60 16,285.70 2
## 61 16,095.20 2
## 62 16,000.00 2
## 63 16,047.60 2
## 64 15,904.80 2
## 65 16,142.90 2
## 66 16,285.70 2
## 67 16,571.40 2
## 68 16,619.00 2
## 69 16,857.10 2
## 70 17,095.20 3
## 71 17,142.90 3
## 72 16,761.90 2
## 73 16,238.10 2
## 74 16,047.60 2
## 75 16,238.10 2
## 76 16,285.70 2
## 77 16,047.60 2
## 78 16,904.80 2
## 79 17,047.60 3
## 80 17,047.60 3
## 81 16,571.40 2
## 82 17,047.60 3
## 83 16,857.10 2
## 84 16,761.90 2
## 85 17,476.20 3
## 86 17,761.90 3
## 87 17,666.70 3
## 88 17,666.70 3
## 89 18,190.50 3
## 90 17,904.80 3
## 91 18,333.30 3
## 92 18,095.20 3
## 93 17,857.10 3
## 94 17,952.40 3
## 95 17,904.80 3
## 96 18,142.90 3
## 97 18,428.60 3
## 98 17,714.30 3
## 99 16,857.10 2
## 100 16,761.90 2
## 101 15,666.70 2
## 102 16,047.60 2
## 103 15,904.80 2
## 104 15,571.40 2
## 105 15,809.50 2
## 106 16,047.60 2
## 107 16,238.10 2
## 108 16,238.10 2
## 109 16,571.40 2
## 110 16,666.70 2
## 111 16,952.40 2
## 112 17,333.30 3
## 113 17,285.70 3
## 114 17,190.50 3
## 115 17,190.50 3
## 116 17,142.90 3
## 117 17,238.10 3
## 118 16,142.90 2
## 119 15,095.20 2
## 120 15,428.60 2
## 121 15,809.50 2
## 122 16,190.50 2
## 123 15,857.10 2
## 124 15,571.40 2
## 125 15,952.40 2
## 126 15,238.10 2
## 127 15,428.60 2
## 128 14,904.80 2
## 129 14,619.00 2
## 130 15,476.20 2
## 131 15,428.60 2
## 132 15,809.50 2
## 133 16,095.20 2
## 134 15,428.60 2
## 135 15,666.70 2
## 136 16,000.00 2
## 137 16,285.70 2
## 138 15,238.10 2
## 139 15,666.70 2
## 140 16,809.50 2
## 141 16,488.10 2
## 142 16,970.70 2
## 143 17,533.70 3
## 144 17,212.00 3
## 145 18,096.70 3
## 146 18,217.30 3
## 147 18,378.20 3
## 148 18,659.70 3
## 149 18,297.80 3
## 150 18,136.90 3
## 151 17,614.10 3
## 152 17,855.40 3
## 153 17,051.10 3
## 154 16,769.60 2
## 155 17,051.10 3
## 156 17,855.40 3
## 157 18,177.10 3
## 158 18,257.50 3
## 159 17,614.10 3
## 160 18,177.10 3
## 161 17,694.50 3
## 162 17,734.80 3
## 163 18,177.10 3
## 164 18,297.80 3
## 165 19,544.40 3
## 166 19,785.70 3
## 167 19,584.60 3
## 168 20,549.80 3
## 169 20,911.70 3
## 170 19,866.10 3
## 171 20,187.90 3
## 172 19,866.10 3
## 173 19,624.90 3
## 174 18,378.20 3
## 175 17,212.00 3
## 176 17,694.50 3
## 177 17,573.90 3
## 178 16,890.20 2
## 179 16,890.20 2
## 180 18,096.70 3
## 181 18,418.40 3
## 182 18,338.00 3
## 183 18,338.00 3
## 184 18,498.80 3
## 185 18,981.40 3
## 186 19,222.70 3
## 187 18,901.00 3
## 188 17,935.80 3
## 189 17,694.50 3
## 190 17,895.60 3
## 191 18,418.40 3
## 192 17,533.70 3
## 193 17,533.70 3
## 194 18,820.60 3
## 195 18,498.80 3
## 196 18,901.00 3
## 197 18,941.20 3
## 198 19,624.90 3
## 199 20,147.60 3
## 200 19,544.40 3
## 201 19,665.10 3
## 202 19,705.30 3
## 203 19,705.30 3
## 204 20,228.10 3
## 205 19,705.30 3
## 206 18,418.40 3
## 207 18,941.20 3
## 208 18,740.10 3
## 209 19,745.50 3
## 210 20,107.40 3
## 211 20,911.70 3
## 212 20,911.70 3
## 213 20,750.90 3
## 214 20,107.40 3
## 215 1
VIXMH=MHVIXtrungbinh
VIXMH=replace(VIXMH,MHVIXtrungbinh==1,"<14 điểm")
VIXMH=replace(VIXMH,MHVIXtrungbinh==2,"Từ 14-17 điểm")
VIXMH=replace(VIXMH,MHVIXtrungbinh==3,"trên 17 điểm")
MHVIXtrungbinhcode=data.frame(MHVIXtrungbinh,VIXMH)
MHVIXtrungbinhcode
## MHVIXtrungbinh VIXMH
## 1 1 <14 điểm
## 2 1 <14 điểm
## 3 1 <14 điểm
## 4 1 <14 điểm
## 5 1 <14 điểm
## 6 1 <14 điểm
## 7 1 <14 điểm
## 8 1 <14 điểm
## 9 1 <14 điểm
## 10 1 <14 điểm
## 11 1 <14 điểm
## 12 1 <14 điểm
## 13 1 <14 điểm
## 14 1 <14 điểm
## 15 1 <14 điểm
## 16 2 Từ 14-17 điểm
## 17 2 Từ 14-17 điểm
## 18 2 Từ 14-17 điểm
## 19 2 Từ 14-17 điểm
## 20 2 Từ 14-17 điểm
## 21 2 Từ 14-17 điểm
## 22 1 <14 điểm
## 23 1 <14 điểm
## 24 1 <14 điểm
## 25 1 <14 điểm
## 26 1 <14 điểm
## 27 1 <14 điểm
## 28 1 <14 điểm
## 29 1 <14 điểm
## 30 1 <14 điểm
## 31 1 <14 điểm
## 32 1 <14 điểm
## 33 1 <14 điểm
## 34 1 <14 điểm
## 35 1 <14 điểm
## 36 1 <14 điểm
## 37 1 <14 điểm
## 38 1 <14 điểm
## 39 1 <14 điểm
## 40 1 <14 điểm
## 41 1 <14 điểm
## 42 1 <14 điểm
## 43 1 <14 điểm
## 44 1 <14 điểm
## 45 1 <14 điểm
## 46 1 <14 điểm
## 47 1 <14 điểm
## 48 1 <14 điểm
## 49 1 <14 điểm
## 50 1 <14 điểm
## 51 2 Từ 14-17 điểm
## 52 1 <14 điểm
## 53 2 Từ 14-17 điểm
## 54 2 Từ 14-17 điểm
## 55 2 Từ 14-17 điểm
## 56 2 Từ 14-17 điểm
## 57 2 Từ 14-17 điểm
## 58 2 Từ 14-17 điểm
## 59 2 Từ 14-17 điểm
## 60 2 Từ 14-17 điểm
## 61 2 Từ 14-17 điểm
## 62 2 Từ 14-17 điểm
## 63 2 Từ 14-17 điểm
## 64 2 Từ 14-17 điểm
## 65 2 Từ 14-17 điểm
## 66 2 Từ 14-17 điểm
## 67 2 Từ 14-17 điểm
## 68 2 Từ 14-17 điểm
## 69 2 Từ 14-17 điểm
## 70 3 trên 17 điểm
## 71 3 trên 17 điểm
## 72 2 Từ 14-17 điểm
## 73 2 Từ 14-17 điểm
## 74 2 Từ 14-17 điểm
## 75 2 Từ 14-17 điểm
## 76 2 Từ 14-17 điểm
## 77 2 Từ 14-17 điểm
## 78 2 Từ 14-17 điểm
## 79 3 trên 17 điểm
## 80 3 trên 17 điểm
## 81 2 Từ 14-17 điểm
## 82 3 trên 17 điểm
## 83 2 Từ 14-17 điểm
## 84 2 Từ 14-17 điểm
## 85 3 trên 17 điểm
## 86 3 trên 17 điểm
## 87 3 trên 17 điểm
## 88 3 trên 17 điểm
## 89 3 trên 17 điểm
## 90 3 trên 17 điểm
## 91 3 trên 17 điểm
## 92 3 trên 17 điểm
## 93 3 trên 17 điểm
## 94 3 trên 17 điểm
## 95 3 trên 17 điểm
## 96 3 trên 17 điểm
## 97 3 trên 17 điểm
## 98 3 trên 17 điểm
## 99 2 Từ 14-17 điểm
## 100 2 Từ 14-17 điểm
## 101 2 Từ 14-17 điểm
## 102 2 Từ 14-17 điểm
## 103 2 Từ 14-17 điểm
## 104 2 Từ 14-17 điểm
## 105 2 Từ 14-17 điểm
## 106 2 Từ 14-17 điểm
## 107 2 Từ 14-17 điểm
## 108 2 Từ 14-17 điểm
## 109 2 Từ 14-17 điểm
## 110 2 Từ 14-17 điểm
## 111 2 Từ 14-17 điểm
## 112 3 trên 17 điểm
## 113 3 trên 17 điểm
## 114 3 trên 17 điểm
## 115 3 trên 17 điểm
## 116 3 trên 17 điểm
## 117 3 trên 17 điểm
## 118 2 Từ 14-17 điểm
## 119 2 Từ 14-17 điểm
## 120 2 Từ 14-17 điểm
## 121 2 Từ 14-17 điểm
## 122 2 Từ 14-17 điểm
## 123 2 Từ 14-17 điểm
## 124 2 Từ 14-17 điểm
## 125 2 Từ 14-17 điểm
## 126 2 Từ 14-17 điểm
## 127 2 Từ 14-17 điểm
## 128 2 Từ 14-17 điểm
## 129 2 Từ 14-17 điểm
## 130 2 Từ 14-17 điểm
## 131 2 Từ 14-17 điểm
## 132 2 Từ 14-17 điểm
## 133 2 Từ 14-17 điểm
## 134 2 Từ 14-17 điểm
## 135 2 Từ 14-17 điểm
## 136 2 Từ 14-17 điểm
## 137 2 Từ 14-17 điểm
## 138 2 Từ 14-17 điểm
## 139 2 Từ 14-17 điểm
## 140 2 Từ 14-17 điểm
## 141 2 Từ 14-17 điểm
## 142 2 Từ 14-17 điểm
## 143 3 trên 17 điểm
## 144 3 trên 17 điểm
## 145 3 trên 17 điểm
## 146 3 trên 17 điểm
## 147 3 trên 17 điểm
## 148 3 trên 17 điểm
## 149 3 trên 17 điểm
## 150 3 trên 17 điểm
## 151 3 trên 17 điểm
## 152 3 trên 17 điểm
## 153 3 trên 17 điểm
## 154 2 Từ 14-17 điểm
## 155 3 trên 17 điểm
## 156 3 trên 17 điểm
## 157 3 trên 17 điểm
## 158 3 trên 17 điểm
## 159 3 trên 17 điểm
## 160 3 trên 17 điểm
## 161 3 trên 17 điểm
## 162 3 trên 17 điểm
## 163 3 trên 17 điểm
## 164 3 trên 17 điểm
## 165 3 trên 17 điểm
## 166 3 trên 17 điểm
## 167 3 trên 17 điểm
## 168 3 trên 17 điểm
## 169 3 trên 17 điểm
## 170 3 trên 17 điểm
## 171 3 trên 17 điểm
## 172 3 trên 17 điểm
## 173 3 trên 17 điểm
## 174 3 trên 17 điểm
## 175 3 trên 17 điểm
## 176 3 trên 17 điểm
## 177 3 trên 17 điểm
## 178 2 Từ 14-17 điểm
## 179 2 Từ 14-17 điểm
## 180 3 trên 17 điểm
## 181 3 trên 17 điểm
## 182 3 trên 17 điểm
## 183 3 trên 17 điểm
## 184 3 trên 17 điểm
## 185 3 trên 17 điểm
## 186 3 trên 17 điểm
## 187 3 trên 17 điểm
## 188 3 trên 17 điểm
## 189 3 trên 17 điểm
## 190 3 trên 17 điểm
## 191 3 trên 17 điểm
## 192 3 trên 17 điểm
## 193 3 trên 17 điểm
## 194 3 trên 17 điểm
## 195 3 trên 17 điểm
## 196 3 trên 17 điểm
## 197 3 trên 17 điểm
## 198 3 trên 17 điểm
## 199 3 trên 17 điểm
## 200 3 trên 17 điểm
## 201 3 trên 17 điểm
## 202 3 trên 17 điểm
## 203 3 trên 17 điểm
## 204 3 trên 17 điểm
## 205 3 trên 17 điểm
## 206 3 trên 17 điểm
## 207 3 trên 17 điểm
## 208 3 trên 17 điểm
## 209 3 trên 17 điểm
## 210 3 trên 17 điểm
## 211 3 trên 17 điểm
## 212 3 trên 17 điểm
## 213 3 trên 17 điểm
## 214 3 trên 17 điểm
## 215 1 <14 điểm
MHVCItrungbinh=datalikert$VCI
MHVCItrungbinh=replace(MHVCItrungbinh,datalikert$VCI<=22,1)
MHVCItrungbinh=replace(MHVCItrungbinh,datalikert$VCI>22&datalikert$VCI<=27,2)
MHVCItrungbinh=replace(MHVCItrungbinh,datalikert$VCI>27,3)
VCI=data.frame(datalikert$VCI,MHVCItrungbinh)
VCI
## datalikert.VCI MHVCItrungbinh
## 1 19,800.00 1
## 2 20,100.00 1
## 3 20,200.00 1
## 4 20,400.00 1
## 5 20,150.00 1
## 6 20,650.00 1
## 7 20,000.00 1
## 8 19,250.00 1
## 9 19,800.00 1
## 10 19,600.00 1
## 11 20,050.00 1
## 12 21,300.00 1
## 13 21,200.00 1
## 14 21,700.00 1
## 15 21,700.00 1
## 16 22,300.00 2
## 17 22,250.00 2
## 18 22,850.00 2
## 19 22,600.00 2
## 20 22,500.00 2
## 21 22,400.00 2
## 22 22,100.00 2
## 23 22,100.00 2
## 24 20,700.00 1
## 25 20,950.00 1
## 26 21,150.00 1
## 27 21,500.00 1
## 28 21,800.00 1
## 29 21,600.00 1
## 30 22,050.00 2
## 31 22,100.00 2
## 32 22,100.00 2
## 33 21,750.00 1
## 34 21,750.00 1
## 35 21,500.00 1
## 36 21,850.00 1
## 37 21,750.00 1
## 38 21,300.00 1
## 39 21,500.00 1
## 40 20,800.00 1
## 41 19,700.00 1
## 42 18,450.00 1
## 43 19,800.00 1
## 44 20,500.00 1
## 45 19,900.00 1
## 46 20,150.00 1
## 47 20,000.00 1
## 48 21,250.00 1
## 49 21,850.00 1
## 50 22,400.00 2
## 51 23,400.00 2
## 52 23,700.00 2
## 53 24,450.00 2
## 54 25,100.00 2
## 55 25,400.00 2
## 56 24,900.00 2
## 57 24,450.00 2
## 58 25,100.00 2
## 59 24,850.00 2
## 60 24,900.00 2
## 61 24,350.00 2
## 62 24,350.00 2
## 63 24,050.00 2
## 64 24,050.00 2
## 65 24,150.00 2
## 66 24,500.00 2
## 67 25,000.00 2
## 68 24,600.00 2
## 69 24,500.00 2
## 70 24,650.00 2
## 71 24,850.00 2
## 72 24,450.00 2
## 73 23,600.00 2
## 74 23,400.00 2
## 75 23,500.00 2
## 76 23,400.00 2
## 77 23,400.00 2
## 78 24,150.00 2
## 79 24,350.00 2
## 80 24,500.00 2
## 81 24,350.00 2
## 82 24,700.00 2
## 83 24,650.00 2
## 84 25,000.00 2
## 85 25,100.00 2
## 86 25,100.00 2
## 87 24,750.00 2
## 88 25,000.00 2
## 89 24,600.00 2
## 90 24,950.00 2
## 91 25,600.00 2
## 92 25,750.00 2
## 93 25,600.00 2
## 94 25,950.00 2
## 95 25,700.00 2
## 96 25,500.00 2
## 97 25,700.00 2
## 98 25,950.00 2
## 99 25,900.00 2
## 100 26,500.00 2
## 101 25,600.00 2
## 102 26,150.00 2
## 103 26,200.00 2
## 104 26,250.00 2
## 105 27,000.00 3
## 106 27,150.00 3
## 107 26,550.00 2
## 108 26,800.00 2
## 109 27,400.00 3
## 110 27,300.00 3
## 111 27,600.00 3
## 112 27,700.00 3
## 113 27,950.00 3
## 114 28,100.00 3
## 115 27,900.00 3
## 116 28,000.00 3
## 117 29,000.00 3
## 118 27,200.00 3
## 119 26,300.00 2
## 120 26,500.00 2
## 121 26,950.00 2
## 122 27,350.00 3
## 123 26,700.00 2
## 124 26,150.00 2
## 125 26,700.00 2
## 126 26,400.00 2
## 127 25,788.80 2
## 128 25,051.90 2
## 129 24,499.30 2
## 130 26,304.50 2
## 131 25,825.60 2
## 132 26,525.60 2
## 133 27,446.60 3
## 134 27,225.60 3
## 135 27,225.60 3
## 136 26,967.70 2
## 137 27,409.80 3
## 138 26,304.50 2
## 139 26,046.60 2
## 140 27,999.20 3
## 141 29,841.30 3
## 142 30,099.20 3
## 143 28,146.60 3
## 144 28,367.60 3
## 145 28,367.60 3
## 146 27,372.90 3
## 147 26,894.00 2
## 148 27,262.40 3
## 149 27,225.60 3
## 150 26,894.00 2
## 151 26,451.90 2
## 152 26,194.00 2
## 153 25,604.60 2
## 154 25,383.50 2
## 155 25,236.10 2
## 156 26,599.30 2
## 157 27,078.20 3
## 158 26,746.60 2
## 159 26,525.60 2
## 160 27,004.50 3
## 161 26,636.10 2
## 162 26,636.10 2
## 163 27,262.40 3
## 164 27,115.00 3
## 165 26,451.90 2
## 166 25,899.30 2
## 167 25,788.80 2
## 168 26,046.60 2
## 169 25,641.40 2
## 170 25,567.70 2
## 171 25,678.20 2
## 172 26,157.20 2
## 173 26,267.70 2
## 174 26,415.10 2
## 175 24,720.40 2
## 176 25,051.90 2
## 177 25,162.50 2
## 178 24,690.40 2
## 179 24,654.10 2
## 180 25,634.50 2
## 181 25,707.10 2
## 182 25,634.50 2
## 183 25,852.30 2
## 184 24,980.90 2
## 185 25,888.70 2
## 186 24,908.30 2
## 187 25,561.90 2
## 188 25,235.10 2
## 189 24,835.70 2
## 190 24,835.70 2
## 191 24,835.70 2
## 192 24,218.40 2
## 193 23,746.40 2
## 194 24,835.70 2
## 195 24,908.30 2
## 196 24,763.10 2
## 197 24,581.50 2
## 198 25,017.20 2
## 199 25,452.90 2
## 200 25,162.50 2
## 201 24,690.40 2
## 202 24,254.70 2
## 203 24,508.90 2
## 204 24,581.50 2
## 205 24,980.90 2
## 206 24,182.10 2
## 207 24,835.70 2
## 208 24,944.60 2
## 209 25,452.90 2
## 210 25,997.60 2
## 211 26,433.30 2
## 212 26,288.10 2
## 213 25,961.30 2
## 214 26,251.70 2
## 215 1
VCIMH=MHVCItrungbinh
VCIMH=replace(VCIMH,MHVCItrungbinh==1,"<22 điểm")
VCIMH=replace(VCIMH,MHVCItrungbinh==2,"Từ 22-27 điểm")
VCIMH=replace(VCIMH,MHVCItrungbinh==3,"trên 27 điểm")
MHVCItrungbinhcode=data.frame(MHVCItrungbinh,VCIMH)
MHVCItrungbinhcode
## MHVCItrungbinh VCIMH
## 1 1 <22 điểm
## 2 1 <22 điểm
## 3 1 <22 điểm
## 4 1 <22 điểm
## 5 1 <22 điểm
## 6 1 <22 điểm
## 7 1 <22 điểm
## 8 1 <22 điểm
## 9 1 <22 điểm
## 10 1 <22 điểm
## 11 1 <22 điểm
## 12 1 <22 điểm
## 13 1 <22 điểm
## 14 1 <22 điểm
## 15 1 <22 điểm
## 16 2 Từ 22-27 điểm
## 17 2 Từ 22-27 điểm
## 18 2 Từ 22-27 điểm
## 19 2 Từ 22-27 điểm
## 20 2 Từ 22-27 điểm
## 21 2 Từ 22-27 điểm
## 22 2 Từ 22-27 điểm
## 23 2 Từ 22-27 điểm
## 24 1 <22 điểm
## 25 1 <22 điểm
## 26 1 <22 điểm
## 27 1 <22 điểm
## 28 1 <22 điểm
## 29 1 <22 điểm
## 30 2 Từ 22-27 điểm
## 31 2 Từ 22-27 điểm
## 32 2 Từ 22-27 điểm
## 33 1 <22 điểm
## 34 1 <22 điểm
## 35 1 <22 điểm
## 36 1 <22 điểm
## 37 1 <22 điểm
## 38 1 <22 điểm
## 39 1 <22 điểm
## 40 1 <22 điểm
## 41 1 <22 điểm
## 42 1 <22 điểm
## 43 1 <22 điểm
## 44 1 <22 điểm
## 45 1 <22 điểm
## 46 1 <22 điểm
## 47 1 <22 điểm
## 48 1 <22 điểm
## 49 1 <22 điểm
## 50 2 Từ 22-27 điểm
## 51 2 Từ 22-27 điểm
## 52 2 Từ 22-27 điểm
## 53 2 Từ 22-27 điểm
## 54 2 Từ 22-27 điểm
## 55 2 Từ 22-27 điểm
## 56 2 Từ 22-27 điểm
## 57 2 Từ 22-27 điểm
## 58 2 Từ 22-27 điểm
## 59 2 Từ 22-27 điểm
## 60 2 Từ 22-27 điểm
## 61 2 Từ 22-27 điểm
## 62 2 Từ 22-27 điểm
## 63 2 Từ 22-27 điểm
## 64 2 Từ 22-27 điểm
## 65 2 Từ 22-27 điểm
## 66 2 Từ 22-27 điểm
## 67 2 Từ 22-27 điểm
## 68 2 Từ 22-27 điểm
## 69 2 Từ 22-27 điểm
## 70 2 Từ 22-27 điểm
## 71 2 Từ 22-27 điểm
## 72 2 Từ 22-27 điểm
## 73 2 Từ 22-27 điểm
## 74 2 Từ 22-27 điểm
## 75 2 Từ 22-27 điểm
## 76 2 Từ 22-27 điểm
## 77 2 Từ 22-27 điểm
## 78 2 Từ 22-27 điểm
## 79 2 Từ 22-27 điểm
## 80 2 Từ 22-27 điểm
## 81 2 Từ 22-27 điểm
## 82 2 Từ 22-27 điểm
## 83 2 Từ 22-27 điểm
## 84 2 Từ 22-27 điểm
## 85 2 Từ 22-27 điểm
## 86 2 Từ 22-27 điểm
## 87 2 Từ 22-27 điểm
## 88 2 Từ 22-27 điểm
## 89 2 Từ 22-27 điểm
## 90 2 Từ 22-27 điểm
## 91 2 Từ 22-27 điểm
## 92 2 Từ 22-27 điểm
## 93 2 Từ 22-27 điểm
## 94 2 Từ 22-27 điểm
## 95 2 Từ 22-27 điểm
## 96 2 Từ 22-27 điểm
## 97 2 Từ 22-27 điểm
## 98 2 Từ 22-27 điểm
## 99 2 Từ 22-27 điểm
## 100 2 Từ 22-27 điểm
## 101 2 Từ 22-27 điểm
## 102 2 Từ 22-27 điểm
## 103 2 Từ 22-27 điểm
## 104 2 Từ 22-27 điểm
## 105 3 trên 27 điểm
## 106 3 trên 27 điểm
## 107 2 Từ 22-27 điểm
## 108 2 Từ 22-27 điểm
## 109 3 trên 27 điểm
## 110 3 trên 27 điểm
## 111 3 trên 27 điểm
## 112 3 trên 27 điểm
## 113 3 trên 27 điểm
## 114 3 trên 27 điểm
## 115 3 trên 27 điểm
## 116 3 trên 27 điểm
## 117 3 trên 27 điểm
## 118 3 trên 27 điểm
## 119 2 Từ 22-27 điểm
## 120 2 Từ 22-27 điểm
## 121 2 Từ 22-27 điểm
## 122 3 trên 27 điểm
## 123 2 Từ 22-27 điểm
## 124 2 Từ 22-27 điểm
## 125 2 Từ 22-27 điểm
## 126 2 Từ 22-27 điểm
## 127 2 Từ 22-27 điểm
## 128 2 Từ 22-27 điểm
## 129 2 Từ 22-27 điểm
## 130 2 Từ 22-27 điểm
## 131 2 Từ 22-27 điểm
## 132 2 Từ 22-27 điểm
## 133 3 trên 27 điểm
## 134 3 trên 27 điểm
## 135 3 trên 27 điểm
## 136 2 Từ 22-27 điểm
## 137 3 trên 27 điểm
## 138 2 Từ 22-27 điểm
## 139 2 Từ 22-27 điểm
## 140 3 trên 27 điểm
## 141 3 trên 27 điểm
## 142 3 trên 27 điểm
## 143 3 trên 27 điểm
## 144 3 trên 27 điểm
## 145 3 trên 27 điểm
## 146 3 trên 27 điểm
## 147 2 Từ 22-27 điểm
## 148 3 trên 27 điểm
## 149 3 trên 27 điểm
## 150 2 Từ 22-27 điểm
## 151 2 Từ 22-27 điểm
## 152 2 Từ 22-27 điểm
## 153 2 Từ 22-27 điểm
## 154 2 Từ 22-27 điểm
## 155 2 Từ 22-27 điểm
## 156 2 Từ 22-27 điểm
## 157 3 trên 27 điểm
## 158 2 Từ 22-27 điểm
## 159 2 Từ 22-27 điểm
## 160 3 trên 27 điểm
## 161 2 Từ 22-27 điểm
## 162 2 Từ 22-27 điểm
## 163 3 trên 27 điểm
## 164 3 trên 27 điểm
## 165 2 Từ 22-27 điểm
## 166 2 Từ 22-27 điểm
## 167 2 Từ 22-27 điểm
## 168 2 Từ 22-27 điểm
## 169 2 Từ 22-27 điểm
## 170 2 Từ 22-27 điểm
## 171 2 Từ 22-27 điểm
## 172 2 Từ 22-27 điểm
## 173 2 Từ 22-27 điểm
## 174 2 Từ 22-27 điểm
## 175 2 Từ 22-27 điểm
## 176 2 Từ 22-27 điểm
## 177 2 Từ 22-27 điểm
## 178 2 Từ 22-27 điểm
## 179 2 Từ 22-27 điểm
## 180 2 Từ 22-27 điểm
## 181 2 Từ 22-27 điểm
## 182 2 Từ 22-27 điểm
## 183 2 Từ 22-27 điểm
## 184 2 Từ 22-27 điểm
## 185 2 Từ 22-27 điểm
## 186 2 Từ 22-27 điểm
## 187 2 Từ 22-27 điểm
## 188 2 Từ 22-27 điểm
## 189 2 Từ 22-27 điểm
## 190 2 Từ 22-27 điểm
## 191 2 Từ 22-27 điểm
## 192 2 Từ 22-27 điểm
## 193 2 Từ 22-27 điểm
## 194 2 Từ 22-27 điểm
## 195 2 Từ 22-27 điểm
## 196 2 Từ 22-27 điểm
## 197 2 Từ 22-27 điểm
## 198 2 Từ 22-27 điểm
## 199 2 Từ 22-27 điểm
## 200 2 Từ 22-27 điểm
## 201 2 Từ 22-27 điểm
## 202 2 Từ 22-27 điểm
## 203 2 Từ 22-27 điểm
## 204 2 Từ 22-27 điểm
## 205 2 Từ 22-27 điểm
## 206 2 Từ 22-27 điểm
## 207 2 Từ 22-27 điểm
## 208 2 Từ 22-27 điểm
## 209 2 Từ 22-27 điểm
## 210 2 Từ 22-27 điểm
## 211 2 Từ 22-27 điểm
## 212 2 Từ 22-27 điểm
## 213 2 Từ 22-27 điểm
## 214 2 Từ 22-27 điểm
## 215 1 <22 điểm
MHSSItrungbinh=datalikert$SSI
MHSSItrungbinh=replace(MHSSItrungbinh,datalikert$SSI<=20,1)
MHSSItrungbinh=replace(MHSSItrungbinh,datalikert$SSI>20&datalikert$SSI<=24,2)
MHSSItrungbinh=replace(MHSSItrungbinh,datalikert$SSI>24,3)
SSI=data.frame(datalikert$SSI,MHSSItrungbinh)
SSI
## datalikert.SSI MHSSItrungbinh
## 1 20,700.00 2
## 2 20,850.00 2
## 3 20,850.00 2
## 4 20,950.00 2
## 5 20,900.00 2
## 6 21,400.00 2
## 7 21,150.00 2
## 8 20,550.00 2
## 9 20,650.00 2
## 10 20,300.00 2
## 11 20,300.00 2
## 12 21,000.00 2
## 13 20,900.00 2
## 14 21,000.00 2
## 15 20,850.00 2
## 16 21,050.00 2
## 17 20,900.00 2
## 18 21,350.00 2
## 19 21,400.00 2
## 20 21,650.00 2
## 21 21,250.00 2
## 22 21,250.00 2
## 23 20,750.00 2
## 24 19,400.00 1
## 25 19,400.00 1
## 26 19,500.00 1
## 27 19,800.00 1
## 28 19,583.30 1
## 29 19,983.00 1
## 30 20,222.80 2
## 31 20,102.90 2
## 32 20,062.90 2
## 33 19,543.40 1
## 34 19,423.50 1
## 35 19,503.40 1
## 36 19,623.30 1
## 37 19,623.30 1
## 38 18,824.00 1
## 39 19,183.70 1
## 40 18,584.20 1
## 41 18,384.40 1
## 42 17,585.00 1
## 43 18,224.50 1
## 44 18,824.00 1
## 45 18,104.60 1
## 46 18,544.20 1
## 47 18,144.60 1
## 48 19,383.50 1
## 49 19,983.00 1
## 50 19,743.20 1
## 51 20,422.60 2
## 52 20,502.50 2
## 53 21,221.90 2
## 54 21,541.70 2
## 55 21,741.50 2
## 56 21,701.50 2
## 57 21,381.80 2
## 58 21,861.40 2
## 59 21,661.60 2
## 60 21,741.50 2
## 61 21,381.80 2
## 62 21,261.90 2
## 63 21,102.00 2
## 64 21,182.00 2
## 65 21,341.80 2
## 66 21,461.70 2
## 67 21,741.50 2
## 68 21,701.50 2
## 69 21,981.30 2
## 70 22,141.20 2
## 71 22,061.20 2
## 72 21,621.60 2
## 73 20,902.20 2
## 74 21,022.10 2
## 75 21,221.90 2
## 76 21,182.00 2
## 77 21,062.10 2
## 78 21,581.60 2
## 79 21,661.60 2
## 80 21,821.40 2
## 81 21,581.60 2
## 82 21,941.30 2
## 83 21,981.30 2
## 84 22,021.30 2
## 85 22,301.00 2
## 86 22,381.00 2
## 87 21,941.30 2
## 88 22,021.30 2
## 89 21,661.60 2
## 90 21,781.50 2
## 91 22,221.10 2
## 92 22,500.80 2
## 93 22,301.00 2
## 94 22,460.90 2
## 95 22,221.10 2
## 96 22,420.90 2
## 97 22,620.70 2
## 98 22,620.70 2
## 99 22,660.70 2
## 100 22,780.60 2
## 101 21,821.40 2
## 102 22,141.20 2
## 103 22,101.20 2
## 104 22,101.20 2
## 105 22,141.20 2
## 106 22,341.00 2
## 107 22,500.80 2
## 108 22,620.70 2
## 109 23,140.30 2
## 110 22,820.60 2
## 111 23,060.40 2
## 112 22,940.50 2
## 113 22,940.50 2
## 114 22,700.70 2
## 115 22,860.50 2
## 116 22,980.40 2
## 117 23,500.00 2
## 118 21,981.30 2
## 119 21,541.70 2
## 120 21,581.60 2
## 121 22,021.30 2
## 122 22,301.00 2
## 123 21,541.70 2
## 124 21,221.90 2
## 125 21,621.60 2
## 126 21,062.10 2
## 127 21,461.70 2
## 128 20,702.40 2
## 129 20,302.70 2
## 130 21,541.70 2
## 131 21,581.60 2
## 132 22,381.00 2
## 133 22,900.50 2
## 134 22,700.70 2
## 135 22,700.70 2
## 136 23,300.20 2
## 137 24,179.40 3
## 138 23,699.80 2
## 139 23,420.10 2
## 140 25,178.60 3
## 141 26,777.20 3
## 142 26,937.10 3
## 143 26,057.80 3
## 144 25,458.30 3
## 145 25,898.00 3
## 146 25,898.00 3
## 147 25,698.10 3
## 148 25,658.20 3
## 149 25,258.50 3
## 150 25,098.60 3
## 151 24,978.70 3
## 152 24,938.80 3
## 153 23,939.60 2
## 154 23,659.90 2
## 155 23,899.70 2
## 156 24,778.90 3
## 157 25,218.50 3
## 158 24,699.00 3
## 159 24,619.00 3
## 160 24,898.80 3
## 161 24,898.80 3
## 162 24,738.90 3
## 163 24,659.00 3
## 164 24,778.90 3
## 165 25,498.30 3
## 166 25,498.30 3
## 167 25,338.40 3
## 168 26,057.80 3
## 169 26,097.80 3
## 170 26,097.80 3
## 171 26,097.80 3
## 172 26,657.30 3
## 173 26,177.70 3
## 174 25,937.90 3
## 175 24,259.40 3
## 176 24,419.20 3
## 177 24,059.50 3
## 178 23,699.80 2
## 179 23,300.20 2
## 180 24,179.40 3
## 181 24,459.20 3
## 182 24,379.30 3
## 183 24,579.10 3
## 184 24,379.30 3
## 185 25,058.70 3
## 186 24,619.00 3
## 187 25,138.60 3
## 188 24,619.00 3
## 189 23,859.70 2
## 190 23,699.80 2
## 191 24,179.40 3
## 192 22,700.70 2
## 193 22,460.90 2
## 194 23,220.20 2
## 195 23,380.10 2
## 196 22,980.40 2
## 197 23,500.00 2
## 198 23,480.00 2
## 199 23,880.50 2
## 200 23,552.80 2
## 201 23,516.40 2
## 202 23,298.00 2
## 203 23,880.50 2
## 204 24,171.70 3
## 205 24,754.10 3
## 206 23,916.90 2
## 207 25,045.40 3
## 208 24,972.60 3
## 209 25,336.60 3
## 210 25,227.40 3
## 211 25,882.60 3
## 212 25,773.40 3
## 213 25,445.80 3
## 214 25,445.80 3
## 215 1
SSIMH=MHSSItrungbinh
SSIMH=replace(VCIMH,MHSSItrungbinh==1,"<20 điểm")
SSIMH=replace(VCIMH,MHSSItrungbinh==2,"Từ 20-24 điểm")
SSIMH=replace(VCIMH,MHSSItrungbinh==3,"trên 24 điểm")
MHSSItrungbinhcode=data.frame(MHSSItrungbinh,SSIMH)
MHSSItrungbinhcode
## MHSSItrungbinh SSIMH
## 1 2 <22 điểm
## 2 2 <22 điểm
## 3 2 <22 điểm
## 4 2 <22 điểm
## 5 2 <22 điểm
## 6 2 <22 điểm
## 7 2 <22 điểm
## 8 2 <22 điểm
## 9 2 <22 điểm
## 10 2 <22 điểm
## 11 2 <22 điểm
## 12 2 <22 điểm
## 13 2 <22 điểm
## 14 2 <22 điểm
## 15 2 <22 điểm
## 16 2 Từ 22-27 điểm
## 17 2 Từ 22-27 điểm
## 18 2 Từ 22-27 điểm
## 19 2 Từ 22-27 điểm
## 20 2 Từ 22-27 điểm
## 21 2 Từ 22-27 điểm
## 22 2 Từ 22-27 điểm
## 23 2 Từ 22-27 điểm
## 24 1 <22 điểm
## 25 1 <22 điểm
## 26 1 <22 điểm
## 27 1 <22 điểm
## 28 1 <22 điểm
## 29 1 <22 điểm
## 30 2 Từ 22-27 điểm
## 31 2 Từ 22-27 điểm
## 32 2 Từ 22-27 điểm
## 33 1 <22 điểm
## 34 1 <22 điểm
## 35 1 <22 điểm
## 36 1 <22 điểm
## 37 1 <22 điểm
## 38 1 <22 điểm
## 39 1 <22 điểm
## 40 1 <22 điểm
## 41 1 <22 điểm
## 42 1 <22 điểm
## 43 1 <22 điểm
## 44 1 <22 điểm
## 45 1 <22 điểm
## 46 1 <22 điểm
## 47 1 <22 điểm
## 48 1 <22 điểm
## 49 1 <22 điểm
## 50 1 Từ 22-27 điểm
## 51 2 Từ 22-27 điểm
## 52 2 Từ 22-27 điểm
## 53 2 Từ 22-27 điểm
## 54 2 Từ 22-27 điểm
## 55 2 Từ 22-27 điểm
## 56 2 Từ 22-27 điểm
## 57 2 Từ 22-27 điểm
## 58 2 Từ 22-27 điểm
## 59 2 Từ 22-27 điểm
## 60 2 Từ 22-27 điểm
## 61 2 Từ 22-27 điểm
## 62 2 Từ 22-27 điểm
## 63 2 Từ 22-27 điểm
## 64 2 Từ 22-27 điểm
## 65 2 Từ 22-27 điểm
## 66 2 Từ 22-27 điểm
## 67 2 Từ 22-27 điểm
## 68 2 Từ 22-27 điểm
## 69 2 Từ 22-27 điểm
## 70 2 Từ 22-27 điểm
## 71 2 Từ 22-27 điểm
## 72 2 Từ 22-27 điểm
## 73 2 Từ 22-27 điểm
## 74 2 Từ 22-27 điểm
## 75 2 Từ 22-27 điểm
## 76 2 Từ 22-27 điểm
## 77 2 Từ 22-27 điểm
## 78 2 Từ 22-27 điểm
## 79 2 Từ 22-27 điểm
## 80 2 Từ 22-27 điểm
## 81 2 Từ 22-27 điểm
## 82 2 Từ 22-27 điểm
## 83 2 Từ 22-27 điểm
## 84 2 Từ 22-27 điểm
## 85 2 Từ 22-27 điểm
## 86 2 Từ 22-27 điểm
## 87 2 Từ 22-27 điểm
## 88 2 Từ 22-27 điểm
## 89 2 Từ 22-27 điểm
## 90 2 Từ 22-27 điểm
## 91 2 Từ 22-27 điểm
## 92 2 Từ 22-27 điểm
## 93 2 Từ 22-27 điểm
## 94 2 Từ 22-27 điểm
## 95 2 Từ 22-27 điểm
## 96 2 Từ 22-27 điểm
## 97 2 Từ 22-27 điểm
## 98 2 Từ 22-27 điểm
## 99 2 Từ 22-27 điểm
## 100 2 Từ 22-27 điểm
## 101 2 Từ 22-27 điểm
## 102 2 Từ 22-27 điểm
## 103 2 Từ 22-27 điểm
## 104 2 Từ 22-27 điểm
## 105 2 trên 27 điểm
## 106 2 trên 27 điểm
## 107 2 Từ 22-27 điểm
## 108 2 Từ 22-27 điểm
## 109 2 trên 27 điểm
## 110 2 trên 27 điểm
## 111 2 trên 27 điểm
## 112 2 trên 27 điểm
## 113 2 trên 27 điểm
## 114 2 trên 27 điểm
## 115 2 trên 27 điểm
## 116 2 trên 27 điểm
## 117 2 trên 27 điểm
## 118 2 trên 27 điểm
## 119 2 Từ 22-27 điểm
## 120 2 Từ 22-27 điểm
## 121 2 Từ 22-27 điểm
## 122 2 trên 27 điểm
## 123 2 Từ 22-27 điểm
## 124 2 Từ 22-27 điểm
## 125 2 Từ 22-27 điểm
## 126 2 Từ 22-27 điểm
## 127 2 Từ 22-27 điểm
## 128 2 Từ 22-27 điểm
## 129 2 Từ 22-27 điểm
## 130 2 Từ 22-27 điểm
## 131 2 Từ 22-27 điểm
## 132 2 Từ 22-27 điểm
## 133 2 trên 27 điểm
## 134 2 trên 27 điểm
## 135 2 trên 27 điểm
## 136 2 Từ 22-27 điểm
## 137 3 trên 24 điểm
## 138 2 Từ 22-27 điểm
## 139 2 Từ 22-27 điểm
## 140 3 trên 24 điểm
## 141 3 trên 24 điểm
## 142 3 trên 24 điểm
## 143 3 trên 24 điểm
## 144 3 trên 24 điểm
## 145 3 trên 24 điểm
## 146 3 trên 24 điểm
## 147 3 trên 24 điểm
## 148 3 trên 24 điểm
## 149 3 trên 24 điểm
## 150 3 trên 24 điểm
## 151 3 trên 24 điểm
## 152 3 trên 24 điểm
## 153 2 Từ 22-27 điểm
## 154 2 Từ 22-27 điểm
## 155 2 Từ 22-27 điểm
## 156 3 trên 24 điểm
## 157 3 trên 24 điểm
## 158 3 trên 24 điểm
## 159 3 trên 24 điểm
## 160 3 trên 24 điểm
## 161 3 trên 24 điểm
## 162 3 trên 24 điểm
## 163 3 trên 24 điểm
## 164 3 trên 24 điểm
## 165 3 trên 24 điểm
## 166 3 trên 24 điểm
## 167 3 trên 24 điểm
## 168 3 trên 24 điểm
## 169 3 trên 24 điểm
## 170 3 trên 24 điểm
## 171 3 trên 24 điểm
## 172 3 trên 24 điểm
## 173 3 trên 24 điểm
## 174 3 trên 24 điểm
## 175 3 trên 24 điểm
## 176 3 trên 24 điểm
## 177 3 trên 24 điểm
## 178 2 Từ 22-27 điểm
## 179 2 Từ 22-27 điểm
## 180 3 trên 24 điểm
## 181 3 trên 24 điểm
## 182 3 trên 24 điểm
## 183 3 trên 24 điểm
## 184 3 trên 24 điểm
## 185 3 trên 24 điểm
## 186 3 trên 24 điểm
## 187 3 trên 24 điểm
## 188 3 trên 24 điểm
## 189 2 Từ 22-27 điểm
## 190 2 Từ 22-27 điểm
## 191 3 trên 24 điểm
## 192 2 Từ 22-27 điểm
## 193 2 Từ 22-27 điểm
## 194 2 Từ 22-27 điểm
## 195 2 Từ 22-27 điểm
## 196 2 Từ 22-27 điểm
## 197 2 Từ 22-27 điểm
## 198 2 Từ 22-27 điểm
## 199 2 Từ 22-27 điểm
## 200 2 Từ 22-27 điểm
## 201 2 Từ 22-27 điểm
## 202 2 Từ 22-27 điểm
## 203 2 Từ 22-27 điểm
## 204 3 trên 24 điểm
## 205 3 trên 24 điểm
## 206 2 Từ 22-27 điểm
## 207 3 trên 24 điểm
## 208 3 trên 24 điểm
## 209 3 trên 24 điểm
## 210 3 trên 24 điểm
## 211 3 trên 24 điểm
## 212 3 trên 24 điểm
## 213 3 trên 24 điểm
## 214 3 trên 24 điểm
## 215 1 <22 điểm
###CÂU C
library(psych)
vxi=data.frame(MHVIXtrungbinh,MHVCItrungbinh,MHSSItrungbinh)
describe(vxi)
## vars n mean sd median trimmed mad min max range skew
## MHVIXtrungbinh* 1 215 2.22 0.78 2 2.28 1.48 1 3 2 -0.41
## MHVCItrungbinh* 2 215 1.96 0.57 2 1.95 0.00 1 3 2 -0.01
## MHSSItrungbinh* 3 215 2.14 0.60 2 2.18 0.00 1 3 2 -0.06
## kurtosis se
## MHVIXtrungbinh* -1.25 0.05
## MHVCItrungbinh* 0.09 0.04
## MHSSItrungbinh* -0.35 0.04
# NHẬN XÉT:
# Cả 3 chỉ số đều có 240 quan sát
# Chỉ số SSIMH có giá trị trung bình cao nhất (22.505)và trung vị (22.221)
# Chỉ số VCIMH có giá trị trung bình (24.674)
# Chỉ số VIXMH có giá trị trung bình (16.376)
# Tính mode và phương sai VIX
library(DescTools)
##
## Attaching package: 'DescTools'
## The following objects are masked from 'package:psych':
##
## AUC, ICC, SD
Mode(MHVIXtrungbinh)
## [1] "3"
## attr(,"freq")
## [1] 94
table(MHVIXtrungbinh)
## MHVIXtrungbinh
## 1 2 3
## 46 75 94
var(MHVIXtrungbinh)
## [1] 0.6041295
# NHẬN XÉT:Mode của biến VIXMH là 3 với 94 quan sát
# Phương sai của biến VIXMH=0.6041295
# Tính mode và phương sai VCI
library(DescTools)
Mode(MHVCItrungbinh)
## [1] "2"
## attr(,"freq")
## [1] 146
table(MHVCItrungbinh)
## MHVCItrungbinh
## 1 2 3
## 39 146 30
var(MHVCItrungbinh)
## [1] 0.3206694
# NHẬN XÉT:Mode của biến VCIMH là 2 với 146 quan sát
# Phương sai của biến VCIMH=0.3206694
# Tính mode và phương sai SSI
library(DescTools)
Mode(MHSSItrungbinh)
## [1] "2"
## attr(,"freq")
## [1] 134
table(MHSSItrungbinh)
## MHSSItrungbinh
## 1 2 3
## 25 134 56
var(MHSSItrungbinh)
## [1] 0.3576179
# NHẬN XÉT:Mode của biến SSIMH là 2 với 134 quan sát
# Phương sai của biến SSIMH=0.3576179
###CÂU D
#Vẽ biểu đồ cột và tròn của cổ phiếu VIX
VIX.freq<-table(VIXMH)
barplot(VIX.freq,col=c("lightpink","lightgreen","lightyellow"),xlab="Mức giá",ylab="Số lượng",main="Biểu đồ cột phân phối mức giá cổ phiếu VIX")
pie(table(VIXMH),col=c("lightpink","lightgreen","lightyellow"),sub="Biểu đồ tròn phân phối mức giá cổ phiếu VIX")
# NHẬN XÉT
# Biểu đồ cột phân phối mức giá cổ phiếu VIX:cho thấy số lượng phiên giao dịch tập trung nhiều nhất ở nhóm trên 17 điểm.
# Biểu đồ tròn phân phối mức giá cổ phiếu VIX:nhóm giá trên 17 điểm chiếm tỷ trọng ưu thế trong tổng thể.
#Vẽ biểu đồ cột và tròn của cổ phiếu VCI
VCI.freq<-table(VCIMH)
barplot(VCI.freq,col=c("lightpink","darkgreen","lightyellow"),xlab= "Mức giá",ylab="Số lượng",main="Biểu đồ cột phân phối mức giá cổ phiếu VCI")
pie(table(VCIMH),col=c("lightpink","darkgreen","lightyellow"),sub="Biểu đồ tròn phân phối mức giá cổ phiếu VCI")
# NHẬN XÉT
# Biểu đồ cột phân phối mức giá cổ phiếu VCI: cho thấy đa số các phiên giao dịch tập trung nhiều nhất ở nhóm từ 22-27 điểm.
# Biểu đồ tròn phân phối mức giá cổ phiếu VCI:nhóm giá từ 22-27 điểm chiếm tỷ trọng ưu thế lớn nhất trong kỳ.
#Vẽ biểu đồ cột và tròn của cổ phiếu SSI
SSI.freq<-table(SSIMH)
barplot(SSI.freq,col=c("lightpink","brown","lightyellow","lightblue"),xlab="Mức giá",ylab="Số lượng",main="Biểu đồ cột phân phối mức giá cổ phiếu SSI")
pie(table(SSIMH),col=c("lightpink","brown","lightyellow","lightblue"),sub="Biểu đồ tròn phân phối mức giá cổ phiếu SSI")
# NHẬN XÉT
# Biểu đồ cột phân phối mức giá cổ phiếu SSI:cho thấy số lượng phiên giao dịch phân bổ qua các khoảng giá trị khác nhau trong mẫu nghiên cứu.
# Biểu đồ tròn phân phối mức giá cổ phiếu SSI:cơ cấu tỷ trọng đóng góp của các nhóm giá đối với cổ phiếu SSI.
# Khai báo số thứ tự ngày dựa trên số dòng của bảng dữ liệu datalikert
STT<-1:nrow(datalikert)
VIXtb<-as.numeric(gsub(",","",datalikert$VIX))
VCItb<-as.numeric(gsub(",","",datalikert$VCI))
SSItb<-as.numeric(gsub(",","",datalikert$SSI))
CleanData<-data.frame(STT,VIXtb,VCItb,SSItb)
head(CleanData)
## STT VIXtb VCItb SSItb
## 1 1 12650 19800 20700
## 2 2 12750 20100 20850
## 3 3 12950 20200 20850
## 4 4 13050 20400 20950
## 5 5 12800 20150 20900
## 6 6 13150 20650 21400
summary(CleanData)
## STT VIXtb VCItb SSItb
## Min. : 1.0 Min. :11190 Min. :18450 Min. :17585
## 1st Qu.: 54.5 1st Qu.:14762 1st Qu.:23425 1st Qu.:21182
## Median :108.0 Median :16714 Median :25000 Median :22221
## Mean :108.0 Mean :16376 Mean :24674 Mean :22505
## 3rd Qu.:161.5 3rd Qu.:18097 3rd Qu.:26297 3rd Qu.:24177
## Max. :215.0 Max. :20912 Max. :30099 Max. :26937
## NAs :1 NAs :1 NAs :1
# Vẽ biểu đồ đường kết hợp cho 3 cổ phiếu
plot(STT,VIXtb,type="o",xlab="Ngày thứ",ylab="Giá trị giao dịch VIX,VCI,SSI",pch=17,col="lightpink",lty=2,main="Biểu đồ đường biến động giá cổ phiếu VIX VCI SSI từ ngày 1 đến phiên cuối",ylim=c(min(min(VIXtb,na.rm=TRUE),min(VCItb,na.rm=TRUE),min(SSItb,na.rm=TRUE)),max(max(VIXtb,na.rm=TRUE),max(VCItb,na.rm=TRUE),max(SSItb,na.rm=TRUE))))
lines(STT,VCItb,type="o",pch=18,col="blue",lty=3)
lines(STT,SSItb,type="o",pch=19,col="red",lty=4)
legend("bottomright",legend=c("VIX","VCI","SSI"),col=c("blue","brown","red"),lty=c(2,3,4),pch=c(17,18,19))
# NHẬN XÉT
# Biểu đồ đường diễn biến xu hướng giá giao dịch của 3 mã cổ phiếu VIX, VCI và SSI qua các phiên giao dịch.
# Các đường biểu diễn cho thấy xu hướng vận động đồng pha tăng giảm đều theo thời gian.
###CÂU E
cor(VIXtb,SSItb)
## [1] NA
cor.test(VIXtb,SSItb)
##
## Pearson's product-moment correlation
##
## data: VIXtb and SSItb
## t = 25.255, df = 212, p-value < 2.2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.8284799 0.8963098
## sample estimates:
## cor
## 0.8663354
plot(SSItb~VIXtb,xlab="VIXtb",ylab="SSItb")
abline(lm(SSItb~VIXtb))
library(psych)
vars_slide<-cbind(VIXtb,VCItb,SSItb)
pairs.panels(vars_slide)
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.2.1 ✔ readr 2.2.0
## ✔ forcats 1.0.1 ✔ stringr 1.6.0
## ✔ ggplot2 4.0.3 ✔ tibble 3.3.1
## ✔ lubridate 1.9.5 ✔ tidyr 1.3.2
## ✔ purrr 1.2.2
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ ggplot2::%+%() masks psych::%+%()
## ✖ ggplot2::alpha() masks psych::alpha()
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
qplot(x=VIXtb,y=SSItb)+stat_smooth(method="lm")
## Warning: `qplot()` was deprecated in ggplot2 3.4.0.
## This warning is displayed once per session.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## `geom_smooth()` using formula = 'y ~ x'
## Warning: Removed 1 row containing non-finite outside the scale range
## (`stat_smooth()`).
## Warning: Removed 1 row containing missing values or values outside the scale range
## (`geom_point()`).
# NHẬN XÉT:
# Biểu đồ phân tán (scatter plot) kết hợp đường xu hướng tuyến tính thể hiện mối quan hệ đồng biến mạnh giữa giá cổ phiếu VIX và SSI.
# Các điểm quan sát tập trung bám sát đường hồi quy, cho thấy mức độ tương quan tuyến tính dương cao giữa hai cổ phiếu này trong suốt kỳ khảo sát.
###CÂU F
mreg=lm(SSItb~VIXtb+VCItb,data=CleanData)
summary(mreg)
##
## Call:
## lm(formula = SSItb ~ VIXtb + VCItb, data = CleanData)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1899.3 -791.8 -233.2 730.6 2991.2
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 7.025e+03 6.948e+02 10.110 < 2e-16 ***
## VIXtb 5.913e-01 4.023e-02 14.698 < 2e-16 ***
## VCItb 2.350e-01 4.017e-02 5.849 1.86e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 948.9 on 211 degrees of freedom
## (1 observation deleted due to missingness)
## Multiple R-squared: 0.7853, Adjusted R-squared: 0.7833
## F-statistic: 386 on 2 and 211 DF, p-value: < 2.2e-16
anova(mreg)
## Analysis of Variance Table
##
## Response: SSItb
## Df Sum Sq Mean Sq F value Pr(>F)
## VIXtb 1 664335447 664335447 737.753 < 2.2e-16 ***
## VCItb 1 30808929 30808929 34.214 1.862e-08 ***
## Residuals 211 190002373 900485
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#NHẬN XÉT:
# Hệ số xác định (Multiple R-squared): Đạt 0.7853 (tương đương 78.53%), cho thấy mô hình hồi quy tuyến tính với hai biến độc lập VIXtb và VCItb giải thích được 78.53% sự biến thiên của biến phụ thuộc SSItb.
# Kiểm định F (F-statistic): Giá trị F-statistic = 386 với p-value < 2.2e-16 (rất nhỏ), chứng tỏ mô hình hồi quy tuyến tính bội tổng thể có ý nghĩa thống kê và các biến độc lập thực sự tác động đến biến phụ thuộc.
# Bảng ANOVA: Cả hai biến VIXtb và VCItb đều có mức ý nghĩa Pr(>F) rất nhỏ (< 0.001, kí hiệu ***), xác nhận rằng cả hai biến này đều có đóng góp và ý nghĩa thống kê quan trọng trong việc giải thích giá cổ phiếu SSI.
###CÂU G
library(car)
## Loading required package: carData
##
## Attaching package: 'car'
##
## The following object is masked from 'package:dplyr':
##
## recode
##
## The following object is masked from 'package:purrr':
##
## some
##
## The following object is masked from 'package:DescTools':
##
## Recode
##
## The following object is masked from 'package:psych':
##
## logit
vif(mreg)
## VIXtb VCItb
## 2.065452 2.065452
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.093383, p-value < 2.2e-16
## alternative hypothesis: true autocorrelation is greater than 0
op<-par(mfrow=c(2,2))
plot(mreg)
# NHẬN XÉT:
# Đồ thị Residuals vs Fitted: Phần dư dao động tương đối đều quanh đường ngang y=0,cho thấy giả định về tính tuyến tính của mô hình được thỏa mãn.
# Đồ thị Q-Q Residuals: Các điểm phần dư chuẩn hóa bám sát khá sát đường chéo chuẩn (diagonal line), qua đó cho thấy phần dư tuân theo quy luật phân phối chuẩn.
# Đồ thị Scale-Location: Đường xu hướng tương đối phẳng và các điểm phân bố rải rác đều, chứng tỏ phương sai của phần dư thay đổi không đáng kể (giả định phương sai sai số đồng nhất được đảm bảo).
# Đồ thị Residuals vs Leverage: Hầu hết các quan sát đều nằm bên trong giới hạn khoảng cách Cook's distance, cho thấy không có điểm dữ liệu ngoại lai cực đoan nào gây ảnh hưởng quá lớn đến độ chệch của toàn bộ mô hình hồi quy.