PRML 12章 確率的主成分分析を試す - 木曜不足PCAを試す に続いて確率的主成分分析(Probability Principal Component Analysis)。 解析的に解けてしまって、閉形式の解がわかっているので実装としてはたいしておもしろくない(いや、いいことなんですけどね)。 M <- 2; directory <- "."; argv <- commandArgs(T); if (length(argv)>0) directory <- commandArgs(T)[1]; if (length(argv)>1) M <- as.integer(commandArgs(T)[2]); oilflow <- as.matrix(read.table(sprintf("%s/DataTrn.txt", directory))); oilflow.labels <- read.table(sprintf("%s/DataT
