Các phép toán lượng giác trên độ lệch chuẩn


14

Phép cộng, phép trừ, phép nhân và phép chia các biến ngẫu nhiên bình thường được xác định rõ, nhưng còn các phép toán lượng giác thì sao?

Ví dụ, chúng ta giả sử rằng tôi đang cố gắng tìm góc của một hình nêm hình tam giác (được mô phỏng như một tam giác góc vuông) với hai catheti có kích thước d 2 , cả hai đều được mô tả là phân phối bình thường.d1d2

Cả hai trực giác và mô phỏng cho tôi biết rằng việc phân phối kết quả là bình thường, với trung bình . Nhưng có cách nào để tính toán phân phối của góc kết quả không? Tài liệu tham khảo về nơi tôi sẽ tìm thấy câu trả lời?arctan(mean(d1)mean(d2))

(Đối với một chút bối cảnh, tôi đang làm việc với khả năng chịu thống kê của các bộ phận cơ học. Sự thôi thúc đầu tiên của tôi sẽ chỉ là mô phỏng toàn bộ quá trình, kiểm tra xem kết quả cuối cùng có bình thường hay không và tính toán độ lệch chuẩn. Nhưng tôi đang tự hỏi nếu có thể có một cách tiếp cận phân tích gọn gàng hơn.)


2
Bạn có thể xác nhận rằng (a) d1 và d2 là độ dài cạnh (và không phải là góc); (b) rằng bạn đang giả sử góc giữa chúng là một góc vuông (nếu không thì công thức atan là nghi ngờ); và (c) bạn có quan tâm đến việc phân phối một trong các góc khác của tam giác vuông này không? Ngoài ra, có lẽ, SD của mỗi phân phối độ dài nhỏ hơn nhiều so với dự kiến ​​của nó vì tam giác không nên có bất kỳ xác suất đáng kể nào về độ dài cạnh âm :-).
whuber

Chính xác. Tôi đã nói lại vấn đề để làm cho nó rõ ràng hơn một chút. Và có, SD sẽ nhỏ so với kích thước.
Bossykena

Sử dụng các công thức để nhân và cộng, bạn có thể thử mở rộng Taylor.

Cảm ơn cả hai câu trả lời xuất sắc của bạn, điều mà (theo như tôi có thể nói với chuyên môn thống kê hạn chế của mình) cả về trực quan và âm thanh.
Bossykena

Câu trả lời:


15

Trong cách giải thích này, tam giác là một hình tam giác bên phải của mặt dài Y phân phối binormally với kỳ vọng μ xμ y , độ lệch chuẩn σ xσ y , và tương quan ρ . Chúng tôi tìm kiếm sự phân phối của arctan ( Y / X ) . Cuối cùng, chuẩn hóa XY sao choXYμxμyσxσyρarctan(Y/X)XY

Y = σ y η + μ y

X=σxξ+μx
Y=σyη+μy

với η variates bình thường tiêu chuẩn với tương quan ρ . Hãy θ là một góc và để thuận tiện ghi q = tan ( θ ) . Sau đóξηρθq=tan(θ)

P[arctan(Y/X)θ]=P[YqX]

=P[σyη+μyq(σxξ+μx)

=P[σyηqσxξqμxμy]

Phía bên tay trái, là một sự kết hợp tuyến tính của normals, là bình thường, với trung bình và phương sai σ 2 y + q 2 σ 2 x - 2 q ρ σ x σ y . μyσyqμxσxσy2+q2σx22qρσxσy

Phân biệt cdf bình thường của các tham số này đối với tạo ra pdf của góc. Biểu thức khá ghê gớm, nhưng một phần quan trọng của nó là hàm mũθ

exp((μy(σy+1)μx(σx+1)tan(θ))22(2ρσxσytan(θ)+σx2+σy2+tan2(θ))),

μxμyσxσyρθ) của logarit của pdf (như thể hiện trong các phương trình (2.6) và (3.1) của tài liệu tham khảo). Tôi đề nghị một hệ thống đại số máy tính (như MatLab hoặc Mathicala) để thực hiện điều này!


1
[π,π)

1
P (Y / X q) = P (Y qX) không đúng nếu X là rv bình thường - X cũng có thể âm.
ronaf

@ronaf: thực ra, vì Y là độ dài cạnh của một tam giác vật lý, chúng ta không nên có X âm ! XYX
shabbychef

2
@ronaf: Đó là ý kiến ​​đúng. Nếu người ta sử dụng độ dài cạnh đã ký và cũng coi góc là giá trị thực (chứ không phải giá trị modulo ), thì không có sự không nhất quán với tính quy tắc trong cả hai trường hợp. Quan điểm của bạn về sự bất bình đẳng có thể sai là tuyệt vời. Tất cả những gì tôi có thể làm để đáp lại là tuyên bố rằng phương trình là một xấp xỉ xuất sắc theo các giả định được đưa ra bởi vì khả năng X hoặc Y bị âm là không đáng kể. 2π
whuber

1
@YBE Tôi đồng ý rằng "+" cuối cùng trong biểu hiện của tôi có vẻ như không thuộc về nó - nó có thể bị trượt khi tôi dọn dẹp đánh dấu TeX. Tôi không có tài liệu tham khảo vì tôi đã tự tính đạo hàm.
whuber

12

Bạn đang xem thống kê vòng tròn và đặc biệt là phân phối vòng tròn được gọi là phân phối chuẩn dự kiến .

Vì một số lý do, chủ đề này có thể hơi khó với google, nhưng hai văn bản chính về thống kê vòng tròn là Phân tích thống kê dữ liệu thông tư của Fisher và Thống kê định hướng của Mardia và Jupp.

Để biết phân tích kỹ lưỡng về phân phối bình thường dự kiến, xem trang 46 của Mardia và Jupp. Có các biểu thức dạng đóng (tích phân hàm lỗi) cho phân phối và như whuber đã đề xuất, nó trông tương tự như bình thường khi 'phương sai' của nó (cẩn thận ở đây, phương sai có nghĩa gì cho một biến ngẫu nhiên trên một vòng tròn? !) là nhỏ, tức là khi phân phối khá tập trung tại một điểm (hoặc hướng hoặc góc).

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.