デミング博士のニューエコノミクスって

デミング博士の”ニューエコノミクス”に書かれた内容と，それに関連する内容を，「マターリ」と理解するページ

分散分析（直交表）のRのプログラム【作成途中 2/5更新】

2006-02-05 18:17:25 | R　統計
　いま，分散分析を勉強するために，Rでの直交表を使った分散分析用のツールを作成中です．
　途中経過をさらしていきます．
　このページは，予告無く変更していきます．

１．プログラム（長いですし，バグありだと思います．）
　　【２月５日更新】

#
# Latin Table 
# References
#
# Washio, Yasutoshi,実験計画法入門　改訂版,
# "Introduction of  Eexperimentation Design", 2nd. Ed.,1997,
# ISBN:4-542-50330-5, P281 - 293.
#
# Funao, Nobuo, The R Tips, 2005, ISBN:4-86167-039-X
#
# RjpWiki,P Johnson tips1,http://www.okada.jp.org/RWiki/index.php?P_Johnson_tips_1
# RjpWiki,P Johnson tips2,http://www.okada.jp.org/RWiki/index.php?P_Johnson_tips_2
# 
# RjpWiki,http://www.okada.jp.org/RWiki/index.php?RjpWiki
#
# R Development Core Team,  R: A language and environment for statistical computing (Ver. 2.21), 2005,
# ISBN 3-900051-07-0, URL http://www.R-project.org.
#
# Usage:
# latinTable.L"no"(y)
# y= Data [vector]
# no=4	L4(2^3)
# no=8	L8(2^7)
# no=9	L9(3^4)
# no=12	L12(2^11)
# no=16	L16(2^15)
# no=18	L18(2^1*3^7)
# no=27	L27(3^13)
# no=32 L32(2^31)
# no=36	L36(2^11*3^12)
#
# Return [data.frame]
#

latin.table.L4<-function(y){

#		Row no	 1 2 3 4
	f01<-as.factor(c(1,1,2,2))

	f02<-as.factor(c(1,2,1,2))

	f03<-as.factor(c(1,2,2,1))

	
	return(data.frame(f01=f01,f02=f02,f03=f03,y=y))
}

latin.table.L8<-function(y){

#		Row no	 1 2  3 4  5 6  7 8
	f01<-as.factor(c(1,1, 1,1, 2,2, 2,2))

	f02<-as.factor(c(1,1, 2,2, 1,1, 2,2))

	f03<-as.factor(c(1,1, 2,2, 2,2, 1,1))

	f04<-as.factor(c(1,2, 1,2, 1,2, 1,2))

	f05<-as.factor(c(1,2, 1,2, 2,1, 2,1))

	f06<-as.factor(c(1,2, 2,1, 1,2, 2,1))

	f07<-as.factor(c(1,2, 2,1, 2,1, 1,2))

	return(data.frame(f01=f01,f02=f02,f03=f03,f04=f04,f05=f05,
		f06=f06,f07=f07,y=y))
}

latin.table.L9<-function(y){

#		Row no	 1 2 3  4 5 6  7 8 9
	f01<-as.factor(c(1,1,1, 2,2,2, 3,3,3))

	f02<-as.factor(c(1,2,3, 1,2,3, 1,2,3))

	f03<-as.factor(c(1,2,3, 2,3,1, 3,1,2))

	f04<-as.factor(c(1,2,3, 3,1,2, 2,3,1))

	return(data.frame(f01=f01,f02=f02,f03=f03,f04=f04,y=y))
}

latin.table.L12<-function(y){

#		Row no	 1 2 3 4  5 6 7 8  9 0 1 2
	f01<-as.factor(c(1,1,1,1, 1,1,2,2, 2,2,2,2))

	f02<-as.factor(c(1,1,1,2, 2,2,1,1, 1,2,2,2))

	f03<-as.factor(c(1,1,2,1, 2,2,2,2, 1,2,1,1))

	f04<-as.factor(c(1,1,2,2, 1,2,2,1, 2,1,2,1))

	f05<-as.factor(c(1,1,2,2, 2,1,1,2, 2,1,1,2))

	f06<-as.factor(c(1,2,1,1, 2,2,1,2, 2,1,2,1))

	f07<-as.factor(c(1,2,1,2, 1,2,2,2, 1,1,1,2))

	f08<-as.factor(c(1,2,1,2, 2,1,2,1, 2,2,1,1))

	f09<-as.factor(c(1,2,2,1, 1,2,1,1, 2,2,1,2))

	f10<-as.factor(c(1,2,2,1, 2,1,2,1, 1,1,2,2))

	f11<-as.factor(c(1,2,2,2, 1,1,1,2, 1,2,2,1))

	return(data.frame(f01=f01,f02=f02,f03=f03,f04=f04,f05=f05,
		f06=f06,f07=f07,f08=f08,f09=f09,f10=f10,f11=f11,y=y))
}

latin.table.L16<-function(y){

#		Row no	 1 2  3 4  5 6  7 8  9 0  1 2  3 4  5 6
	f01<-as.factor(c(1,1, 1,1, 1,1, 1,1, 2,2, 2,2, 2,2, 2,2))

	f02<-as.factor(c(1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2))

	f03<-as.factor(c(1,1, 1,1, 2,2, 2,2, 2,2, 2,2, 1,1, 1,1))

	f04<-as.factor(c(1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2))

	f05<-as.factor(c(1,1, 2,2, 1,1, 2,2, 2,2, 1,1, 2,2, 1,1))

	f06<-as.factor(c(1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1))

	f07<-as.factor(c(1,1, 2,2, 2,2, 1,1, 2,2, 1,1, 1,1, 2,2))

	f08<-as.factor(c(1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2))

	f09<-as.factor(c(1,2, 1,2, 1,2, 1,2, 2,1, 2,1, 2,1, 2,1))

	f10<-as.factor(c(1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1))

	f11<-as.factor(c(1,2, 1,2, 2,1, 2,1, 2,1, 2,1, 1,2, 1,2))

	f12<-as.factor(c(1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1))

	f13<-as.factor(c(1,2, 2,1, 1,2, 2,1, 2,1, 1,2, 2,1, 1,2))

	f14<-as.factor(c(1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2))

	f15<-as.factor(c(1,2, 2,1, 2,1, 1,2, 2,1, 1,2, 1,2, 2,1))

	return(data.frame(f01=f01,f02=f02,f03=f03,f04=f04,f05=f05,f06=f06,
		f07=f07,f08=f08,f09=f09,f10=f10,f11=f11,f12=f12,f13=f13,
		f14=f14,f15=f15,y=y))
}

latin.table.L18<-function(y){

#		Row no	 1 2 3  4 5 6  7 8 9  0 1 2  3 4 5  6 7 8
	f01<-as.factor(c(1,1,1, 1,1,1, 1,1,1, 2,2,2, 2,2,2, 2,2,2))

	f02<-as.factor(c(1,1,1, 2,2,2, 3,3,3, 1,1,1, 2,2,2, 3,3,3))

	f03<-as.factor(c(1,2,3, 1,2,3, 1,2,3, 1,2,3, 1,2,3, 1,2,3))

	f04<-as.factor(c(1,2,3, 1,2,3, 2,3,1, 3,1,2, 2,3,1, 3,1,2))

	f05<-as.factor(c(1,2,3, 2,3,1, 1,2,3, 3,1,2, 3,1,2, 2,3,1))

	f06<-as.factor(c(1,2,3, 2,3,1, 3,1,2, 2,3,1, 1,2,3, 3,1,2))

	f07<-as.factor(c(1,2,3, 3,1,2, 2,3,1, 2,3,1, 3,1,2, 1,2,3))

	f08<-as.factor(c(1,2,3, 3,1,2, 3,1,2, 1,2,3, 2,3,1, 2,3,1))

	return(data.frame(f01=f01,f02=f02,f03=f03,f04=f04,f05=f05,f06=f06,
		f07=f07,f08=f08,y=y))
}

latin.table.L27<-function(y){

#			                      1                      2
#		Row no	 1 2 3  4 5 6  7 8 9  0 1 2  3 4 5  6 7 8  9 0 1  2 3 4  5 6 7
	f01<-as.factor(c(1,1,1, 1,1,1, 1,1,1, 2,2,2, 2,2,2, 2,2,2, 3,3,3, 3,3,3, 3,3,3))

	f02<-as.factor(c(1,1,1, 2,2,2, 3,3,3, 1,1,1, 2,2,2, 3,3,3, 1,1,1, 2,2,2, 3,3,3))

	f03<-as.factor(c(1,1,1, 2,2,2, 3,3,3, 2,2,2, 3,3,3, 1,1,1, 3,3,3, 1,1,1, 2,2,2))

	f04<-as.factor(c(1,1,1, 2,2,2, 3,3,3, 3,3,3, 1,1,1, 2,2,2, 2,2,2, 3,3,3, 1,1,1))

	f05<-as.factor(c(1,2,3, 1,2,3, 1,2,3, 1,2,3, 1,2,3, 1,2,3, 1,2,3, 1,2,3, 1,2,3))

	f06<-as.factor(c(1,2,3, 1,2,3, 1,2,3, 2,3,1, 2,3,1, 2,3,1, 3,1,2, 3,1,2, 3,1,2))

	f07<-as.factor(c(1,2,3, 1,2,3, 1,2,3, 3,1,2, 3,1,2, 3,1,2, 2,3,1, 2,3,1, 2,3,1))

	f08<-as.factor(c(1,2,3, 2,3,1, 3,1,2, 1,2,3, 2,3,1, 3,1,2, 1,2,3, 2,3,1, 3,1,2))

	f09<-as.factor(c(1,2,3, 2,3,1, 3,1,2, 2,3,1, 3,1,2, 1,2,3, 3,1,2, 1,2,3, 2,3,1))

	f10<-as.factor(c(1,2,3, 2,3,1, 3,1,2, 3,1,2, 1,2,3, 2,3,1, 2,3,1, 3,1,2, 1,2,3))

	f11<-as.factor(c(1,2,3, 3,1,2, 2,3,1, 1,2,3, 3,1,2, 2,3,1, 1,2,3, 3,1,2, 2,3,1))

	f12<-as.factor(c(1,2,3, 3,1,2, 2,3,1, 2,3,1, 1,2,3, 3,1,2, 3,1,2, 2,3,1, 1,2,3))

	f13<-as.factor(c(1,2,3, 3,1,2, 2,3,1, 3,1,2, 2,3,1, 1,2,3, 2,3,1, 1,2,3, 3,1,2))

	return(data.frame(f01=f01,f02=f02,f03=f03,f04=f04,f05=f05,f06=f06,f07=f07,f08=f08,
		f09=f09,f10=f10,f11=f11,f12=f12,f13=f13,y=y))
}

latin.table.L32<-function(y){

#			                       1                        2                        3
#		Row no	 1 2  3 4  5 6  7 8  9 0  1 2  3 4  5 6  7 8  9 0  1 2  3 4  5 6  7 8  9 0  1 2
	f01<-as.factor(c(1,1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,1, 2,2, 2,2, 2,2, 2,2, 2,2, 2,2, 2,2, 2,2))

	f02<-as.factor(c(1,1, 1,1, 1,1, 1,1, 2,2, 2,2, 2,2, 2,2, 1,1, 1,1, 1,1, 1,1, 2,2, 2,2, 2,2, 2,2))

	f03<-as.factor(c(1,1, 1,1, 1,1, 1,1, 2,2, 2,2, 2,2, 2,2, 2,2, 2,2, 2,2, 2,2, 1,1, 1,1, 1,1, 1,1))

	f04<-as.factor(c(1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2))

	f05<-as.factor(c(1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 1,1))

	f06<-as.factor(c(1,1, 1,1, 2,2, 2,2, 2,2, 2,2, 1,1, 1,1, 1,1, 1,1, 2,2, 2,2, 2,2, 2,2, 1,1, 1,1))

	f07<-as.factor(c(1,1, 1,1, 2,2, 2,2, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 1,1, 1,1, 2,2, 2,2))

	f08<-as.factor(c(1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2))

	f09<-as.factor(c(1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1))

	f10<-as.factor(c(1,1, 2,2, 1,1, 2,2, 2,2, 1,1, 2,2, 1,1, 1,1, 2,2, 1,1, 2,2, 2,2, 1,1, 2,2, 1,1))

	f11<-as.factor(c(1,1, 2,2, 1,1, 2,2, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 2,2, 1,1, 1,1, 2,2, 1,1, 2,2))

	f12<-as.factor(c(1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1))

	f13<-as.factor(c(1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2))

	f14<-as.factor(c(1,1, 2,2, 2,2, 1,1, 2,2, 1,1, 1,1, 2,2, 1,1, 2,2, 2,2, 1,1, 2,2, 1,1, 1,1, 2,2))

	f15<-as.factor(c(1,1, 2,2, 2,2, 1,1, 2,2, 1,1, 1,1, 2,2, 2,2, 1,1, 1,1, 2,2, 1,1, 2,2, 2,2, 1,1))

	f16<-as.factor(c(1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2))

	f17<-as.factor(c(1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 1,2, 2,1, 2,1, 2,1, 2,1, 2,1, 2,1, 2,1, 2,1))

	f18<-as.factor(c(1,2, 1,2, 1,2, 1,2, 2,1, 2,1, 2,1, 2,1, 1,2, 1,2, 1,2, 1,2, 2,1, 2,1, 2,1, 2,1))

	f19<-as.factor(c(1,2, 1,2, 1,2, 1,2, 2,1, 2,1, 2,1, 2,1, 2,1, 2,1, 2,1, 2,1, 1,2, 1,2, 1,2, 1,2))

	f20<-as.factor(c(1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1))

	f21<-as.factor(c(1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 1,2))

	f22<-as.factor(c(1,2, 1,2, 2,1, 2,1, 2,1, 2,1, 1,2, 1,2, 1,2, 1,2, 2,1, 2,1, 2,1, 2,1, 1,2, 1,2))

	f23<-as.factor(c(1,2, 1,2, 2,1, 2,1, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 1,2, 1,2, 2,1, 2,1))

	f24<-as.factor(c(1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1))

	f25<-as.factor(c(1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2))

	f26<-as.factor(c(1,2, 2,1, 1,2, 2,1, 2,1, 1,2, 2,1, 1,2, 1,2, 2,1, 1,2, 2,1, 2,1, 1,2, 2,1, 1,2))

	f27<-as.factor(c(1,2, 2,1, 1,2, 2,1, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 2,1, 1,2, 1,2, 2,1, 1,2, 2,1))

	f28<-as.factor(c(1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2))

	f29<-as.factor(c(1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1))

	f30<-as.factor(c(1,2, 2,1, 2,1, 1,2, 2,1, 1,2, 1,2, 2,1, 1,2, 2,1, 2,1, 1,2, 2,1, 1,2, 1,2, 2,1))

	f31<-as.factor(c(1,2, 2,1, 2,1, 1,2, 2,1, 1,2, 1,2, 2,1, 2,1, 1,2, 1,2, 2,1, 1,2, 2,1, 2,1, 1,2))

	return(data.frame(f01=f01,f02=f02,f03=f03,f04=f04,f05=f05,f06=f06,f07=f07,f08=f08,f09=f09,f10=f10,
		f11=f11,f12=f12,f13=f13,f14=f14,f15=f15,f16=f16,f17=f17,f18=f18,f19=f19,f20=f20,f21=f21,
		f22=f22,f23=f23,f24=f24,f25=f25,f26=f26,f27=f27,f28=f28,f29=f29,f30=f30,f31=f31,y=y))
}


最新の画像［もっと見る］