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