library(HSAUR3)## Loading required package: toolsdta <- HSAUR3::backpain
list(dta)## [[1]]
## ID status driver suburban
## 1 1 case yes yes
## 2 1 control yes no
## 3 2 case yes yes
## 4 2 control yes yes
## 5 3 case yes no
## 6 3 control yes yes
## 7 4 case no no
## 8 4 control no no
## 9 5 case yes no
## 10 5 control yes yes
## 11 6 case yes yes
## 12 6 control yes no
## 13 7 case yes yes
## 14 7 control no no
## 15 8 case yes no
## 16 8 control no no
## 17 9 case yes yes
## 18 9 control yes no
## 19 10 case no yes
## 20 10 control no yes
## 21 11 case yes yes
## 22 11 control yes no
## 23 12 case yes no
## 24 12 control yes yes
## 25 13 case yes yes
## 26 13 control yes yes
## 27 14 case yes yes
## 28 14 control yes yes
## 29 15 case yes yes
## 30 15 control yes yes
## 31 16 case yes no
## 32 16 control yes no
## 33 17 case yes no
## 34 17 control yes yes
## 35 18 case yes yes
## 36 18 control yes yes
## 37 19 case yes yes
## 38 19 control yes yes
## 39 20 case yes yes
## 40 20 control no no
## 41 21 case yes yes
## 42 21 control no no
## 43 22 case yes no
## 44 22 control yes no
## 45 23 case no no
## 46 23 control yes yes
## 47 24 case yes no
## 48 24 control yes yes
## 49 25 case yes yes
## 50 25 control yes yes
## 51 26 case no no
## 52 26 control yes no
## 53 27 case yes yes
## 54 27 control yes yes
## 55 28 case yes yes
## 56 28 control yes yes
## 57 29 case yes no
## 58 29 control yes yes
## 59 30 case yes no
## 60 30 control yes yes
## 61 31 case yes no
## 62 31 control yes yes
## 63 32 case yes no
## 64 32 control yes yes
## 65 33 case yes yes
## 66 33 control yes yes
## 67 34 case no no
## 68 34 control yes no
## 69 35 case yes yes
## 70 35 control yes no
## 71 36 case yes yes
## 72 36 control yes yes
## 73 37 case yes no
## 74 37 control yes no
## 75 38 case yes yes
## 76 38 control yes yes
## 77 39 case yes yes
## 78 39 control yes yes
## 79 40 case yes no
## 80 40 control yes yes
## 81 41 case yes yes
## 82 41 control yes no
## 83 42 case yes yes
## 84 42 control yes yes
## 85 43 case yes yes
## 86 43 control no no
## 87 44 case yes no
## 88 44 control no no
## 89 45 case no yes
## 90 45 control no no
## 91 46 case no no
## 92 46 control yes yes
## 93 47 case yes no
## 94 47 control yes no
## 95 48 case yes no
## 96 48 control no no
## 97 49 case yes no
## 98 49 control yes no
## 99 50 case no no
## 100 50 control no no
## 101 51 case yes yes
## 102 51 control yes yes
## 103 52 case yes yes
## 104 52 control yes yes
## 105 53 case yes yes
## 106 53 control yes yes
## 107 54 case no no
## 108 54 control no no
## 109 55 case no no
## 110 55 control no no
## 111 56 case yes no
## 112 56 control yes yes
## 113 57 case no no
## 114 57 control no no
## 115 58 case yes yes
## 116 58 control yes yes
## 117 59 case yes yes
## 118 59 control yes yes
## 119 60 case yes no
## 120 60 control yes no
## 121 61 case yes yes
## 122 61 control yes yes
## 123 62 case yes no
## 124 62 control yes yes
## 125 63 case no yes
## 126 63 control no no
## 127 64 case yes yes
## 128 64 control yes yes
## 129 65 case yes yes
## 130 65 control yes no
## 131 66 case yes no
## 132 66 control yes yes
## 133 67 case yes yes
## 134 67 control yes yes
## 135 68 case yes yes
## 136 68 control yes no
## 137 69 case yes yes
## 138 69 control yes yes
## 139 70 case yes yes
## 140 70 control yes yes
## 141 71 case yes yes
## 142 71 control yes no
## 143 72 case yes yes
## 144 72 control no no
## 145 73 case yes yes
## 146 73 control yes no
## 147 74 case yes no
## 148 74 control yes yes
## 149 75 case yes yes
## 150 75 control yes yes
## 151 76 case yes yes
## 152 76 control yes no
## 153 77 case yes no
## 154 77 control yes yes
## 155 78 case yes no
## 156 78 control yes yes
## 157 79 case yes yes
## 158 79 control no no
## 159 80 case yes yes
## 160 80 control yes no
## 161 81 case no no
## 162 81 control no no
## 163 82 case yes yes
## 164 82 control yes yes
## 165 83 case yes no
## 166 83 control yes yes
## 167 84 case yes yes
## 168 84 control no yes
## 169 85 case yes yes
## 170 85 control no no
## 171 86 case yes no
## 172 86 control yes no
## 173 87 case yes no
## 174 87 control yes no
## 175 88 case yes yes
## 176 88 control yes no
## 177 89 case yes yes
## 178 89 control yes no
## 179 90 case yes yes
## 180 90 control yes yes
## 181 91 case no no
## 182 91 control yes no
## 183 92 case yes yes
## 184 92 control no no
## 185 93 case yes no
## 186 93 control yes yes
## 187 94 case yes yes
## 188 94 control no yes
## 189 95 case yes yes
## 190 95 control no no
## 191 96 case yes no
## 192 96 control no no
## 193 97 case yes yes
## 194 97 control yes no
## 195 98 case yes yes
## 196 98 control no yes
## 197 99 case yes yes
## 198 99 control yes yes
## 199 100 case yes no
## 200 100 control yes yes
## 201 101 case yes no
## 202 101 control yes no
## 203 102 case yes no
## 204 102 control yes no
## 205 103 case yes yes
## 206 103 control yes yes
## 207 104 case no no
## 208 104 control no no
## 209 105 case yes yes
## 210 105 control yes yes
## 211 106 case no no
## 212 106 control yes no
## 213 107 case yes no
## 214 107 control no no
## 215 108 case yes no
## 216 108 control no no
## 217 109 case yes yes
## 218 109 control yes yes
## 219 110 case yes yes
## 220 110 control yes no
## 221 111 case yes no
## 222 111 control no no
## 223 112 case yes yes
## 224 112 control no no
## 225 113 case no no
## 226 113 control yes no
## 227 114 case yes yes
## 228 114 control yes yes
## 229 115 case yes yes
## 230 115 control no yes
## 231 116 case yes yes
## 232 116 control yes yes
## 233 117 case yes yes
## 234 117 control yes yes
## 235 118 case yes no
## 236 118 control yes no
## 237 119 case yes yes
## 238 119 control no no
## 239 120 case yes yes
## 240 120 control yes yes
## 241 121 case no yes
## 242 121 control yes no
## 243 122 case yes yes
## 244 122 control yes no
## 245 123 case yes yes
## 246 123 control yes yes
## 247 124 case yes no
## 248 124 control no yes
## 249 125 case yes yes
## 250 125 control yes no
## 251 126 case yes yes
## 252 126 control yes no
## 253 127 case yes yes
## 254 127 control yes no
## 255 128 case yes yes
## 256 128 control no no
## 257 129 case no no
## 258 129 control yes no
## 259 130 case yes yes
## 260 130 control yes yes
## 261 131 case yes no
## 262 131 control yes no
## 263 132 case yes yes
## 264 132 control yes no
## 265 133 case yes no
## 266 133 control yes no
## 267 134 case yes yes
## 268 134 control yes yes
## 269 135 case no no
## 270 135 control yes no
## 271 136 case yes no
## 272 136 control yes no
## 273 137 case yes yes
## 274 137 control yes yes
## 275 138 case yes no
## 276 138 control yes no
## 277 139 case no no
## 278 139 control no no
## 279 140 case yes yes
## 280 140 control yes yes
## 281 141 case no no
## 282 141 control yes no
## 283 142 case yes yes
## 284 142 control yes yes
## 285 143 case no yes
## 286 143 control yes yes
## 287 144 case yes yes
## 288 144 control no no
## 289 145 case yes no
## 290 145 control no no
## 291 146 case yes no
## 292 146 control yes yes
## 293 147 case yes no
## 294 147 control yes no
## 295 148 case yes yes
## 296 148 control no no
## 297 149 case yes no
## 298 149 control yes yes
## 299 150 case yes yes
## 300 150 control yes no
## 301 151 case yes no
## 302 151 control yes no
## 303 152 case yes yes
## 304 152 control yes yes
## 305 153 case yes no
## 306 153 control yes no
## 307 154 case yes yes
## 308 154 control yes yes
## 309 155 case yes yes
## 310 155 control yes yes
## 311 156 case yes no
## 312 156 control yes no
## 313 157 case yes yes
## 314 157 control yes yes
## 315 158 case yes yes
## 316 158 control no no
## 317 159 case yes no
## 318 159 control yes yes
## 319 160 case no no
## 320 160 control yes yes
## 321 161 case yes yes
## 322 161 control yes yes
## 323 162 case yes no
## 324 162 control yes yes
## 325 163 case yes yes
## 326 163 control yes no
## 327 164 case yes yes
## 328 164 control yes yes
## 329 165 case no no
## 330 165 control yes yes
## 331 166 case no no
## 332 166 control yes no
## 333 167 case yes yes
## 334 167 control yes yes
## 335 168 case yes no
## 336 168 control yes yes
## 337 169 case yes no
## 338 169 control yes yes
## 339 170 case yes yes
## 340 170 control yes yes
## 341 171 case yes no
## 342 171 control no no
## 343 172 case yes yes
## 344 172 control yes no
## 345 173 case yes yes
## 346 173 control yes no
## 347 174 case yes no
## 348 174 control no no
## 349 175 case yes yes
## 350 175 control no no
## 351 176 case no no
## 352 176 control yes no
## 353 177 case yes no
## 354 177 control no no
## 355 178 case yes yes
## 356 178 control yes yes
## 357 179 case yes no
## 358 179 control no no
## 359 180 case yes no
## 360 180 control no no
## 361 181 case yes yes
## 362 181 control yes yes
## 363 182 case yes no
## 364 182 control yes yes
## 365 183 case yes no
## 366 183 control yes yes
## 367 184 case yes yes
## 368 184 control no no
## 369 185 case yes yes
## 370 185 control no no
## 371 186 case yes yes
## 372 186 control yes yes
## 373 187 case yes yes
## 374 187 control yes no
## 375 188 case yes yes
## 376 188 control no no
## 377 189 case yes yes
## 378 189 control yes yes
## 379 190 case yes yes
## 380 190 control yes yes
## 381 191 case yes yes
## 382 191 control yes no
## 383 192 case yes yes
## 384 192 control no no
## 385 193 case yes yes
## 386 193 control yes no
## 387 194 case yes yes
## 388 194 control yes yes
## 389 195 case yes yes
## 390 195 control yes yes
## 391 196 case yes yes
## 392 196 control no no
## 393 197 case yes yes
## 394 197 control yes no
## 395 198 case yes yes
## 396 198 control yes no
## 397 199 case yes yes
## 398 199 control no no
## 399 200 case no no
## 400 200 control yes yes
## 401 201 case no no
## 402 201 control yes yes
## 403 202 case yes yes
## 404 202 control yes yes
## 405 203 case yes yes
## 406 203 control yes yes
## 407 204 case yes no
## 408 204 control no no
## 409 205 case yes yes
## 410 205 control yes yes
## 411 206 case yes no
## 412 206 control yes no
## 413 207 case yes yes
## 414 207 control yes yes
## 415 208 case no yes
## 416 208 control no yes
## 417 209 case no no
## 418 209 control no no
## 419 210 case yes yes
## 420 210 control yes no
## 421 211 case yes no
## 422 211 control yes yes
## 423 212 case yes yes
## 424 212 control yes yes
## 425 213 case yes no
## 426 213 control yes yes
## 427 214 case yes yes
## 428 214 control yes yes
## 429 215 case yes yes
## 430 215 control yes yes
## 431 216 case no no
## 432 216 control yes yes
## 433 217 case yes yes
## 434 217 control yes no##看資料
help(backpain)## starting httpd help server ... donestr(dta)## 'data.frame': 434 obs. of 4 variables:
## $ ID : Factor w/ 217 levels "1","2","3","4",..: 1 1 2 2 3 3 4 4 5 5 ...
## $ status : Factor w/ 2 levels "case","control": 1 2 1 2 1 2 1 2 1 2 ...
## $ driver : Factor w/ 2 levels "no","yes": 2 2 2 2 2 2 1 1 2 2 ...
## $ suburban: Factor w/ 2 levels "no","yes": 2 1 2 2 1 2 1 1 1 2 ...##只將driver,suburban 算total
with(dta, table(driver,suburban))## suburban
## driver no yes
## no 73 13
## yes 127 221##加入status分類
with(dta, table(driver,suburban,status ))## , , status = case
##
## suburban
## driver no yes
## no 26 6
## yes 64 121
##
## , , status = control
##
## suburban
## driver no yes
## no 47 7
## yes 63 100##到這即可以看出各個總數
drivers <- c("no", "no", "yes", "yes")
suburbans <- c("no", "yes", "no", "yes")
cases <- c(26, 6, 64, 121)
controls <- c(47, 7, 63, 100)
totals <- c(73, 13, 127, 221)
Summarize <- data.frame(driver = drivers, suburban =suburbans , case = cases, control = controls, total = totals)
Summarize## driver suburban case control total
## 1 no no 26 47 73
## 2 no yes 6 7 13
## 3 yes no 64 63 127
## 4 yes yes 121 100 221##以下只是進行嘗試
##將整體資料以status分類driver和suburban
dss <- aggregate(cbind(driver, suburban) ~ status, data=dta, FUN= length)
dss## status driver suburban
## 1 case 217 217
## 2 control 217 217##將整體資料以driver分類status
sss <- aggregate(cbind(status) ~ driver, data=dta, FUN= length)
sss## driver status
## 1 no 86
## 2 yes 348yss <- aggregate(cbind(status) ~ suburban, data=dta, FUN= length)
yss## suburban status
## 1 no 200
## 2 yes 234