hello everybody, this is my first post to the list. i am facing two foolish but tedious problems: 1) i have a data frame like this: var1 var2 var3 var4 var5 var6 var7 var8 var9 1 4.12 1314 867 552 276 302 105 13 420 2 8.99 2712 2197 1492 597 271 217 27 1031 3 9.64 2263 1720 1064 588 136 385 23 1064 4 2.53 1328 876 558 226 212 212 13 465 5 3.60 2260 1808 1130 565 181 203 23 972 6 8.88 2274 1819 864 682 182 227 0 1137 7 4.92 1110 744 477 233 189 144 11 244 8 12.11 2895 2085 1419 608 550 232 29 1216 9 2.37 1438 963 604 316 273 144 14 489 10 17.30 1332 932 546 333 253 146 13 493 11 2.65 1285 977 604 308 154 103 13 501 12 4.35 2069 1221 828 290 621 186 21 600 13 2.37 2015 1471 967 403 302 181 0 806 14 2.09 1707 1280 802 410 171 188 17 563 15 0.70 2027 1318 811 365 405 223 20 527 16 7.77 2264 1471 702 634 589 181 0 860 17 15.18 1578 1073 600 347 395 95 16 631 18 24.46 1991 1573 1035 478 139 239 20 776 19 14.55 2222 1489 933 422 467 244 89 489 20 4.00 1818 1455 909 418 109 218 18 745 i need to compute the linear model lm() pairwise between the columns. i need also to compute the anova(). i would like to get into one object all the results, instead of running many lm() and anova() separately for each couple of columns. 2) i would like to export into a text file the results from aov() and TukeyHSD(). is there any suggestion? thankyou very much for your kind help. sincerely d -- Dario Greco Institute of Biotechnology - University of Helsinki Building Cultivator II P.O.Box 56 Viikinkaari 4 FIN-00014 Finland Office: +358 9 191 58951 Fax: +358 9 191 58952 Mobile: +358 44 023 5780