Week3_EX7
CKC
2018年3月21日
# 資料讀取與整理
fL <- "http://www.stat.columbia.edu/~gelman/book/data/schiz.asc"
dta <- read.table(fL, skip = 4, col.names = paste0("T", 101:130))
head(dta)## T101 T102 T103 T104 T105 T106 T107 T108 T109 T110 T111 T112 T113 T114
## 1 312 272 350 286 268 328 298 356 292 308 296 372 396 402
## 2 354 346 384 342 302 312 322 376 306 402 320 298 308 414
## 3 256 284 320 274 324 268 370 430 314 312 362 256 342 388
## 4 260 294 306 292 264 290 272 268 344 362 330 280 354 320
## 5 204 272 250 260 314 308 246 236 208 268 272 264 308 236
## 6 590 312 286 310 778 364 318 316 316 298 344 262 274 330
## T115 T116 T117 T118 T119 T120 T121 T122 T123 T124 T125 T126 T127 T128
## 1 280 330 254 282 350 328 332 308 292 258 340 242 306 328
## 2 304 422 388 422 426 338 332 426 478 372 392 374 430 388
## 3 302 366 298 396 274 226 328 274 258 220 236 272 322 284
## 4 334 276 418 288 338 350 350 324 286 322 280 256 218 256
## 5 238 350 272 252 252 236 306 238 350 206 260 280 274 318
## 6 312 310 376 326 346 334 282 292 282 300 290 302 300 306
## T129 T130
## 1 294 272
## 2 354 368
## 3 274 356
## 4 220 356
## 5 268 210
## 6 294 444#
dtaNew <- dta %>%
mutate(schizoid = c(rep("Y", 11), rep("N", 6)), id = 1:17) %>%
gather(key = "measure", value = "ms", 1:30) %>%
mutate(measure = as.factor(measure))# 敘述統計分析
describe.by(dtaNew,dtaNew$measure)##
## Descriptive statistics by group
## group: T101
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 1.00 0.00 1 1.0 0.00 1 1 0 NaN
## ms 4 17 355.18 134.38 312 347.6 100.82 204 620 416 0.86
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.72 32.59
## --------------------------------------------------------
## group: T102
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 2.00 0.00 2 2.0 0.00 2 2 0 NaN
## ms 4 17 387.18 198.36 312 355.2 59.30 240 1014 774 2.03
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 3.36 48.11
## --------------------------------------------------------
## group: T103
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 3.00 0.00 3 3.00 0.00 3 3 0 NaN
## ms 4 17 396.82 311.70 320 328.67 62.27 230 1586 1356 3.23
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 9.48 75.60
## --------------------------------------------------------
## group: T104
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 4.00 0.00 4 4.0 0.00 4 4 0 NaN
## ms 4 17 356.12 259.05 286 299.2 47.44 222 1344 1122 3.23
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 9.45 62.83
## --------------------------------------------------------
## group: T105
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 5.00 0.00 5 5.0 0.00 5 5 0 NaN
## ms 4 17 395.53 165.51 324 382.4 85.99 210 778 568 0.99
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.38 40.14
## --------------------------------------------------------
## group: T106
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 6.00 0.00 6 6.0 0.00 6 6 0 NaN
## ms 4 17 419.65 187.90 328 402.8 65.23 254 838 584 1.09
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.40 45.57
## --------------------------------------------------------
## group: T107
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 7.00 0.00 7 7.00 0.00 7 7 0 NaN
## ms 4 17 359.53 136.57 318 340.53 77.10 232 772 540 1.62
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 2.20 33.12
## --------------------------------------------------------
## group: T108
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 8.00 0.00 8 8.0 0.00 8 8 0 NaN
## ms 4 17 342.12 80.38 340 339.2 106.75 228 500 272 0.21
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -1.15 19.49
## --------------------------------------------------------
## group: T109
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 9.00 0.00 9 9.00 0.00 9 9 0 NaN
## ms 4 17 370.82 150.24 314 356.53 50.41 208 748 540 1.18
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 0.14 36.44
## --------------------------------------------------------
## group: T110
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 10.00 0.00 10 10.00 0.00 10 10 0 NaN
## ms 4 17 400.82 198.46 334 366.13 97.85 246 1076 830 2.31
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 5.04 48.13
## --------------------------------------------------------
## group: T111
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.0
## measure* 3 17 11.00 0.00 11 11.0 0.00 11 11 0 NaN
## ms 4 17 381.53 177.94 320 348.8 62.27 256 998 742 2.4
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 5.56 43.16
## --------------------------------------------------------
## group: T112
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 12.00 0.00 12 12.00 0.00 12 12 0 NaN
## ms 4 17 356.94 150.06 314 334.93 85.99 248 796 548 1.75
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 2.12 36.39
## --------------------------------------------------------
## group: T113
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 13 0.00 13 13.00 0.00 13 13 0 NaN
## ms 4 17 406 179.97 332 378.27 83.03 260 968 708 1.77
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 2.79 43.65
## --------------------------------------------------------
## group: T114
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 14.00 0.00 14 14.0 0.00 14 14 0 NaN
## ms 4 17 461.41 340.48 364 395.2 65.23 236 1680 1444 2.71
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 6.76 82.58
## --------------------------------------------------------
## group: T115
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 15.00 0.00 15 15.00 0.00 15 15 0 NaN
## ms 4 17 432.94 343.13 318 361.47 56.34 238 1700 1462 2.94
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 7.98 83.22
## --------------------------------------------------------
## group: T116
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 16.00 0.00 16 16.0 0.00 16 16 0 NaN
## ms 4 17 387.41 145.88 344 362.4 53.37 268 882 614 2.28
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 4.96 35.38
## --------------------------------------------------------
## group: T117
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 17.00 0.00 17 17.0 0.00 17 17 0 NaN
## ms 4 17 377.41 94.26 376 373.2 97.85 254 564 310 0.55
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.64 22.86
## --------------------------------------------------------
## group: T118
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 18.00 0.00 18 18.00 0.00 18 18 0 NaN
## ms 4 17 377.76 141.89 348 365.73 109.71 216 720 504 1.07
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 0.18 34.41
## --------------------------------------------------------
## group: T119
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 19 0.00 19 19.0 0.00 19 19 0 NaN
## ms 4 17 382 124.70 338 368.4 85.99 252 716 464 1.37
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 0.96 30.24
## --------------------------------------------------------
## group: T120
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9 5.93 1 17 16 0.00
## measure* 3 17 20.00 0.00 20 20 0.00 20 20 0 NaN
## ms 4 17 322.47 57.68 324 320 41.51 226 456 230 0.29
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.26 13.99
## --------------------------------------------------------
## group: T121
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 21.00 0.00 21 21.0 0.00 21 21 0 NaN
## ms 4 17 352.82 61.00 350 348.8 62.27 264 502 238 0.73
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.02 14.80
## --------------------------------------------------------
## group: T122
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 22.00 0.00 22 22.00 0.00 22 22 0 NaN
## ms 4 17 357.65 173.75 316 324.93 68.20 232 974 742 2.58
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 6.40 42.14
## --------------------------------------------------------
## group: T123
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 23.00 0.00 23 23.00 0.00 23 23 0 NaN
## ms 4 17 482.59 415.33 292 417.87 62.27 222 1714 1492 1.83
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 2.17 100.73
## --------------------------------------------------------
## group: T124
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.0
## measure* 3 17 24.00 0.00 24 24.00 0.00 24 24 0 NaN
## ms 4 17 344.35 139.15 300 326.67 94.89 206 748 542 1.5
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 1.69 33.75
## --------------------------------------------------------
## group: T125
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 25.00 0.00 25 25.00 0.00 25 25 0 NaN
## ms 4 17 344.12 90.06 320 340.93 88.96 226 510 284 0.54
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -1.11 21.84
## --------------------------------------------------------
## group: T126
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 26.00 0.00 26 26.0 0.00 26 26 0 NaN
## ms 4 17 369.18 171.15 306 348.8 74.13 234 810 576 1.83
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 1.98 41.51
## --------------------------------------------------------
## group: T127
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 27 0.00 27 27.00 0.00 27 27 0 NaN
## ms 4 17 396 201.90 306 369.33 47.44 208 984 776 1.58
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms 1.76 48.97
## --------------------------------------------------------
## group: T128
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.0 5.93 1 17 16 0.00
## measure* 3 17 28.00 0.00 28 28.0 0.00 28 28 0 NaN
## ms 4 17 325.53 70.73 306 322.8 71.16 232 460 228 0.54
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.90 17.15
## --------------------------------------------------------
## group: T129
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9 5.93 1 17 16 0.00
## measure* 3 17 29.00 0.00 29 29 0.00 29 29 0 NaN
## ms 4 17 366.24 128.85 352 358 115.64 206 650 444 0.66
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.58 31.25
## --------------------------------------------------------
## group: T130
## vars n mean sd median trimmed mad min max range skew
## schizoid* 1 17 NaN NA NA NaN NA Inf -Inf -Inf NA
## id 2 17 9.00 5.05 9 9.00 5.93 1 17 16 0.00
## measure* 3 17 30.00 0.00 30 30.00 0.00 30 30 0 NaN
## ms 4 17 379.65 121.71 368 373.73 115.64 206 642 436 0.44
## kurtosis se
## schizoid* NA NA
## id -1.41 1.22
## measure* NaN 0.00
## ms -0.64 29.52# 畫圖
ggplot(dtaNew, aes(measure, ms, color = schizoid))+
coord_flip()+
geom_boxplot()
# 用Ancova檢定是否存在差異
m0 <- lm(ms ~ measure + schizoid, data = dtaNew)
anova(m0)## Analysis of Variance Table
##
## Response: ms
## Df Sum Sq Mean Sq F value Pr(>F)
## measure 29 638735 22025 0.8046 0.7571
## schizoid 1 4506212 4506212 164.6050 <2e-16 ***
## Residuals 479 13113064 27376
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1從結果可發現，在控制受試者接受的方法時，精神疾病的有無在反應時間上存在差異。