実験計画法

Rの実験計画ライブラリの直交表の行列化

メンバー一覧の取得 attributes*1 desnumメンバーを取得 attr*2,"desnum") attrはクラスの引数で与えた要素の取得こっちのが簡単 data.matrix(oa.design(nfactors=factor_len,nlevels=2,seed=1)) data.matrix関数は「行列化+数値化」 *1:oa.design(nfactors=…

主効果でなく交絡がみえる実験

f<-function(v){ 5*v[1]+4*v[2]+1000*v[1]*v[2]*v[2]*v[2]+rnorm(1,mean=0,sd=1); } L8<-function(x=1.0,fun=f,v1,v2){ mat<-cbind(v1,v2) e1<-c(1,1,1,1,1,1,1) e2<-c(1,1,1,2,2,2,2) e3<-c(1,2,2,1,1,2,2) e4<-c(1,2,2,2,2,1,1) e5<-c(2,1,2,1,2,1,2) e6<-…