分散0の列を取り除く

piyo <- hoge[,apply(hoge,2,var) != 0]

というのがあったが，実際に何が行われるかはともかくとして，見かけ上は，
piyo <- hoge[,sd(hoge) > 0]
でよい。ここで，
piyo <- hoge[,var(hoge) > 0]
としてはいけない。

黄金四角形（黄金長方形）を R で描く

長い R プログラムも見たが，簡単に描くことができる。

phi <- (1+sqrt(5))/2
radian <- 90:180*pi/180
xy <- rbind(c(cos(radian), -1, 0, 0, -1), c(sin(radian), 1, 1, 0, 0))
par(mar=c(0, 0, 0, 0))
plot(c(-1, phi-1), c(0, 1), type="n", axes=FALSE, xlab="", ylab="", asp=1)
for (m in 1:11) {
    lines(t(xy))
    xy <- (matrix(c(0, -1, 1, 0), 2)%*%xy+c(0, 1/phi))/phi
}

実行結果は以下の通り