Fit HMMs with fixed emissions matrix

[ahwillia]

Is there a way to do this? Essentially I'd like to pass a flag to fit(..., method="em") that would tell the code to optimize the transitions matrix and keep the emissions fixed (or vice versa!)