DoE.baseで直交表の自動生成(実験計画法)
library("DoE.base") name<-c("Taro","Jiro","Hanako","Honda") gender<-c("male","female") year<-c("1964","1852") amada<-c("1","2","3") oaTable<-oa.design(factor.names=list(name=name,genger=gender,year=year,amada=amada),seed=1); print(oaTable)
出力例
name genger year amada
1 Jiro male 1852 1
2 Honda male 1852 1
3 Hanako female 1964 2
4 Taro female 1852 2
5 Hanako male 1964 2
6 Honda female 1964 1
7 Hanako female 1852 1
8 Taro male 1852 2
9 Hanako male 1852 3
10 Taro male 1964 1
11 Jiro male 1964 3
12 Taro female 1964 3
13 Taro female 1852 3
14 Jiro female 1852 2
15 Honda male 1852 3
16 Jiro female 1964 1
17 Taro male 1964 1
18 Hanako female 1852 1
19 Honda female 1852 3
20 Honda male 1964 2
21 Jiro female 1964 3
22 Honda female 1964 2
23 Hanako male 1964 3
24 Jiro male 1852 2
class=design, type= oa
参考)Rで実験計画法 後編
http://www.slideshare.net/itoyan110/r-14261638