formula で，~ の右に来る変数は factor でないといけない。よくわからないエラーメッセージでしばらく悩んだ。 > library(coin)> d <- data.frame(x=c(3,1,2,4,3,2,5,4,3,4,5),+ g=c(1,1,1,1,1,2,2,2,2,2,2))> wilcox_test(x~g, data=d) 以下にエラー check…

変数 x のとる値が，1，3，6，9 のとき，対応する色 red, green, blue, brown を割り当てるx <- c(1,6,9,3,3,6,9,1,3)y <- as.integer(names(table(x)))color <- c("red", "green", "blue", "brown")sapply(x, function(z) color[which(z==y)]) # 解1color[m…

元のプログラムsystem.time({set.seed(12345)hits <- res <- 1:1000for (j in 1:1000) {for (i in 1:1000) { res[i] <- ifelse( (20 %in% sample(0:99,50,replace=TRUE)),1,0) }hits[j] <- sum(res) }hist(hits)print(summary(hits))})# Min. 1st Qu. Median…

いつの間にか（R 2.13.0 ですでに），factanal( rotation="promax") で，因子間相関係数行列が出力されるようになっているが，因子との対応がおかしいのではないか? Loadings は SS loadings の大きい順に並べ替えられているのに，因子間相関係数行列の行と…

「R による計算機統計学」が発売されていた。 「順列検定」は「並べ替え検定」に， 「確率変数を生み出す方法」は「確率変数の発生方法」に直っていた。

裏で何がどのようにやられていようと，一つの関数を呼ぶだけで描画したいとでも?以下のような関数を定義しておいて，plot(glmオブジェクト) とやれば，グラフは描かれます。 -- 河童の屁は，河童にあらず，屁である。? 2011-06-03 (金) 14:18:45plot.glm <- …

ロジスティック回帰の曲線のグラフを描きたい † うー? (2011-06-02 (木) 17:03:52) ロジスティック回帰をしています。 量xと反応の割合pについて、glm関数でロジスティック回帰式を出しました。 xとpについてプロットを行い、この図に回帰曲線を加えたいと思…