lz是纯新手,最近自学R,之前一直在用spss和statistica分析数据,origin画图。
[data]
a b c d e f g h i j k group
0,111990203 0,012029078 6,752896212 0,536668872 0,014514174 6,904861351 0,0010866 1,48195239 3,813457638 42,73482249 0,54978455 A
0,106987804 0,011929605 6,772133004 0,537274784 0,013319431 6,924173778 0,003145276 1,442204758 3,905086072 42,95747412 0,552349665 A
0,123271942 0,011868791 6,769000625 0,535468297 0,017791536 6,988105614 0 1,475574003 3,807341925 43,20342948 0,547337088 A
0,136544078 0,012096461 6,80874592 0,54144094 0,013886531 6,816121156 0,022240525 1,412292059 3,761038753 42,00924996 0,575777926 A
0,17041966 0,011763584 6,779839357 0,540402173 0,012841541 6,961579033 0,009756848 1,467799269 3,892201845 42,91439005 0,561922604 A
0,145134722 0,011826287 6,861098105 0,550934642 0,012575961 7,051123952 0,015 1,534185619 3,860057932 42,87682947 0,577760929 A
0,074892467 0,012947401 6,853003794 0,514794771 0,021568814 6,927803849 0,031522079 1,136410043 3,968567789 43,20986033 0,559264251 B
0,078819173 0,01474719 6,851823037 0,509480053 0,01245698 6,915329536 0,012093416 1,119547485 3,999153126 44,21876365 0,536320659 B
0,077954012 0,016427751 6,770323109 0,491429815 0,022 6,943968377 0,016667733 1,221195898 4,050066745 44,64933894 0,5245253 B
0,078896449 0,020054913 6,843967914 0,478009218 0,012 6,993476358 0,018377406 1,219648893 4,000852513 44,4190703 0,516441537 B
0,071789516 0,013940767 6,982805628 0,535979562 0,022 7,006614305 0,011242685 1,213174402 4,049117207 44,56148595 0,561163014 B
0,078902745 0,013757741 6,893327514 0,503568608 0,012 6,856468495 0,009260867 1,164955382 3,994288232 43,30296737 0,526587216 B
0,091563544 0,013994374 6,910195353 0,4967731 0,02168352 7,105350784 0,026563871 1,292579273 4,089605478 45,03777244 0,537331344 C
0,072898687 0,013142807 7,005630093 0,491715077 0,022 7,105225934 0,022150197 1,171007854 4,064519618 44,9057476 0,527008081 C
0,07360525 0,020528199 6,966739375 0,479652301 0,022 7,05431912 0,010441761 1,214950819 4,040445422 44,83254558 0,510622261 C
0,071547346 0,014399085 7,109453608 0,495838803 0,022 7,495403106 0,023800954 1,476620461 4,146614226 45,2566688 0,534038842 C
0,074214218 0,014510562 7,007595336 0,478979815 0,022 6,97473462 0,014674392 1,084874096 4,035295468 44,42370436 0,508164769 C
0,071790559 0,018786096 6,985984366 0,504207148 0,022 6,951052039 0,016214618 1,114232835 4,054896766 45,02378895 0,539207862 C
[/data]
这是数据,第一列的a~k是相互独立的元素名称,group是分组依据,想对其进行ANOVA分析,每次anova只针对a~k中的一列在group的分组下,比较A B C三组之间的差异(当然之前要进行方差齐性检验,posthoc Tukey等),基本的对一列数据分析我已经会了,比如aov(a ~ group, data),我之前想到的笨方法是在excel里把a批量替换成bcdef等等,然后把code拷到R上进行运算。
现在,我自己尝试写了下for循环,但苦于不知道怎样表达aov()中的 eg. a ~ group ,我尝试用get(colnames(data)), 返回的是FALSE,如果去掉get(),那就不是变量名了。也搜索了很多资料,暂时不知怎么解决。
我知道这个问题对大家很简单,我想通过这个例子再深入理解一下R和编程思维。请大家多多指教。