Re: MLE of a restricted RV
- From: Scott Seidman <namdiesttocs@xxxxxxxxxxxxxx>
- Date: 22 Aug 2007 12:49:05 GMT
Outlier <MTBrenneman@xxxxxxxxx> wrote in news:1187742184.834976.297520
@q3g2000prf.googlegroups.com:
I want to find the MLE for an angle, a. The sine of the angle is
normally distributed s.t.
X ~ N( sin(a0), s^2)
Let me try to recast this a bit, and see if that offers any insight. I'm
assuming your data is not circular.
Instead of thinking about the angle itself, try thinking about this as a
vector of length 1 and angle theta in the complex plane (you can cast the
problem this way with Euler's formula). Now, the Imaginary component of
this vector is sin(theta), and is distributed normally (as defined by
your problem). The real part, when the problem is constrained this way,
is a simple trig relationship to the imaginary part).
I'm dealing with a problem much like this in neuroscience-- except the
magnitude of the vector is also free to change-- and I think I have a
valid approach, which I'm writing up soon for Neuroscience Methods. I'm
doing my stats on the slope of a zero-intercept orthogonal least squares
fit to the data. The slope is related to the angle, of course, by an
arctan.
I have a prelim version of this "in press" in appedix form in
Experimental Brain Research,
http://www.springerlink.com/content/100473/?Content+Status=Accepted
doi number: 10.1007/s00221-007-1072-3 Author: Seidman
This approach uses an ordinary least squares approach, but the orthogonal
least squares is much more correct. Feel free to contact me if you need
more detail-- just remember to reverse the first part of my email addy.
--
Scott
Reverse name to reply
.
- Follow-Ups:
- Re: MLE of a restricted RV
- From: Outlier
- Re: MLE of a restricted RV
- References:
- MLE of a restricted RV
- From: Outlier
- MLE of a restricted RV
- Prev by Date: Re: MLE of a restricted RV
- Next by Date: Two part model question - How to combine models?
- Previous by thread: Re: MLE of a restricted RV
- Next by thread: Re: MLE of a restricted RV
- Index(es):
Relevant Pages
|
Loading