mychisq.test()の検証 2 母分散が同じがどうかの検定をしてみる。 s <- function(){ sample(1:3,20,replace=TRUE) } s.t <- function(){ table(s())}#テーブル作成 x<-c() for(i in 1:100)x<-c(x,mychisq.test(s.t())) y<-c() for(i in 1:100)y<-c(y,chisq.t…

いよいよ正規のマイコースプログラム期間開始である。気を引き締めてかかりたいところである。 mychisq.test()の自作 とりあえず求めるものは出来上がった。 mat1<-function(n,N){ if (N==0) return (matrix(0,ncol=n)) else if(n==1)return(matrix(N)) else…

やはり家だとあまり働かないので、できるだけ研究室に行くようにしよう。先生の助言と試行錯誤の繰り返しにより、行列をつくる関数は完成した。 mat1<-function(n,N){ if (N==0) return (matrix(0,ncol=n)) else if(n==1)return(matrix(N)) else{ tmp.list<-…