A report on my Oja Median problem:
Gabor Grothendieck suggested replacing my loopy cofactor function with:
cofactors2 - function(x) {
x - t(rbind(1,cbind(x,1)))
p - nrow(x)
solve( x, (seq(p)==p)+0) * det(x)
}
Which is much better, especially if the dimension p is
I am shocked and dismayed (and the term hasn't even started yet ;-) )
that none of you have turned in the weekend homework problem that
I assigned last Friday. At the risk of further embarrassment I am
posting an answer in the hope that this will inspire someone to suggest
some improvements. In