help(formula) ## Create a formula for a model with a large number of variables: xnam <- paste("x", 1:25, sep="") (fmla <- as.formula(paste("y ~ ", paste(xnam, collapse= "+")))) k<-4 # 例えば xnam <- paste("Fs[,", 1:k, "]",sep="") (fmla <- …

１元配置 # 1元配置 Ng<-5 Ns<-sample(5:10,Ng,replace=TRUE) Fs<-c() for(i in 1:Ng){ Fs<-c(Fs,rep(i,Ns[i])) } Vs<-rnorm(sum(Ns)) D<-cbind(Fs,Vs) print(D) plot(Vs ~ Fs) myfit <- lm(Vs ~ Fs) anova(myfit) 2元配置 各水準の繰返し数によって「分散の…

実験計画法入門byNAG フィッシャーによる実験計画法 実験計画法では研究を２つに分ける 計画性を持ってデータを集まる段階 集めたデータを解析する段階 実験計画立案時のポイント 目的の確認(目的を表す数値である『特性値』を定める) 因子の選定(因子として…