[R] Suggestion on how to improve efficiency when using MASS:::hubers on high-dimensional arrays

2007-01-19 Thread Benilton Carvalho
Hi Everyone, Given the scenario I have, I was wondering if anyone would be able to give me a hind on how to get the results from hubers() in a more efficient way. I have an outcome on an array [N x S x D]. I also have a factor (levels 1,2,3) stored on a matrix N x S. My objective is to get

Re: [R] Suggestion on how to improve efficiency when using MASS:::hubers on high-dimensional arrays

2007-01-19 Thread Prof Brian Ripley
The usual advice would seem to apply here: - profile (the real application, not a toy example) to see where the bottlenecks are. - rewrite those in C. It is quite possible that hubers would benefit in your problem by being rewritten in C. However, there is a reason why it was not. It is a