Here's the fix

[willpearse]

These are the small changes.

I made the dimnames and colnames easier to track, and I made your loop use proper indexing