my.datafile <- tempfile() cat(file=my.datafile, " 65 a 138 66 a 331 67 a 177 68 a 446 69 a 279 70 a 317 71 a 279 72 a 443 73 a 1391 74 a 567 75 a 477 76 a 294 77 a 292 78 a 201 79 a 267 80 a 357 81 a 599 82 a 563 83 a 481 84 a 1576 85 a 1170 86 a 1217 87 a 377 88 a 431 89 a 459 65 b 815 66 b 1143 67 b 607 68 b 571 69 b 631 70 b 495 71 b 1210 72 b 987 73 b 956 74 b 859 75 b 1179 76 b 772 77 b 1224 78 b 1146 79 b 661 80 b 729 81 b 845 82 b 1166 83 b 1854 84 b 835 85 b 968 86 b 907 87 b 604 88 b 1304 89 b 559 65 c 259 66 c 202 67 c 102 68 c 214 69 c 211 70 c 330 71 c 516 72 c 178 73 c 833 74 c 265 75 c 348 76 c 236 77 c 570 78 c 674 79 c 494 80 c 454 81 c 270 82 c 238 83 c 98 84 c 268 85 c 449 86 c 562 87 c 380 88 c 392 89 c 425 ", sep=" ") options(scipen=999) # suppressing scientific notation birddata <- read.table(my.datafile, header=FALSE, col.names=c("year", "route", "bird")) attach(birddata) year <- factor(year) route <- factor(route)