Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

12-05-2012 02:54 PM

Thank you, Rick. It is good enough. I got it.

My purpose is to find an information matrix at convergence. To get it, I need to know the relationship between vec(A) and vec(A`), which is only one part of the information matrix.

- Mark as New
- Bookmark
- Subscribe
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to renc

12-05-2012 02:54 PM

In SAS/IML, the SHAPECOL function is used to generate vec(A). You can use the SHAPE function to generate vec(A`) directly, without using any transpose operators.

Explicitly forming a permutation matrix is usually inefficient because you are using an nxn matrix to represent a permutation on N elements. I avoid permutation matrices whenever possible and use the permutation as a subscript index to map a vector to a permutation of iteself. However, if you insists on it, you can generate a permutation matrix from a permutation vector by using the function in this blog post: Automating the Great Christmas Gift Exchange - The DO Loop