Bài tập "Mô hình thống kê tuyến tính"
Bắt đầu bởi Alligator, 20-10-2005 - 20:37
#1
Đã gửi 20-10-2005 - 20:37
Do dòng đời đưa đẩy nên mình đang phải học một môn gọi là "Linear Statistical Models". Trước đây mình có học qua môn Xác suất-Thống kê nhưng cũng không mấy khi dùng nên hơi bị... quên... hết. Background của mình lại là Engineering nên mấy phần này không rành lắm. Nhờ các bạn thành thạo giúp thêm.
==== Bài tập ====
1. [1, problem 1.5, p. 20]
Suppose that http://dientuvietnam...n/mimetex.cgi?Y if the http://dientuvietnam.../mimetex.cgi?Xs have the Weibull density
,
Reference
1. McCullagh, P., and Nelder, J.A., Generalized Linear Models, 2nd Edition, Chapman & Hall/CRC, 1999.
2. Searle, S.R., Linear Models, John Wiley & Sons, Inc., 1971.
==== Bài tập ====
1. [1, problem 1.5, p. 20]
Suppose that http://dientuvietnam...n/mimetex.cgi?Y if the http://dientuvietnam.../mimetex.cgi?Xs have the Weibull density
,
Reference
1. McCullagh, P., and Nelder, J.A., Generalized Linear Models, 2nd Edition, Chapman & Hall/CRC, 1999.
2. Searle, S.R., Linear Models, John Wiley & Sons, Inc., 1971.
<span style='color:blue'>Roses are red,
violets are blue,
Fermat is dead,
but his theorem is true. </span>
violets are blue,
Fermat is dead,
but his theorem is true. </span>
#2
Đã gửi 20-10-2005 - 20:55
2. [1, problem 1.6, p. 20]
The probable error, http://dientuvietnam...n/mimetex.cgi?M is the median of http://dientuvietnam.../mimetex.cgi?Y. Find the probable errors of
1. the exponential distribution;
2. the double exponential distribution
3. the logistic distribution
4. the Cauchy distribution
5. the Normal distribution
Discuss briefly the differences between the probable error and the inter-quartile range.
The historical definition of probable error appears to be vague. Some authors take http://dientuvietnam...n/mimetex.cgi?M to be the mean; others take http://dientuvietnam.net/cgi-bin/mimetex.cgi?\tau to be a multiple (0.67) of the standard deviation.
The probable error, http://dientuvietnam...n/mimetex.cgi?M is the median of http://dientuvietnam.../mimetex.cgi?Y. Find the probable errors of
1. the exponential distribution;
2. the double exponential distribution
3. the logistic distribution
4. the Cauchy distribution
5. the Normal distribution
Discuss briefly the differences between the probable error and the inter-quartile range.
The historical definition of probable error appears to be vague. Some authors take http://dientuvietnam...n/mimetex.cgi?M to be the mean; others take http://dientuvietnam.net/cgi-bin/mimetex.cgi?\tau to be a multiple (0.67) of the standard deviation.
<span style='color:blue'>Roses are red,
violets are blue,
Fermat is dead,
but his theorem is true. </span>
violets are blue,
Fermat is dead,
but his theorem is true. </span>
#3
Đã gửi 21-10-2005 - 22:39
Bài 1 có thể làm như sau. Ký hiệu http://dientuvietnam.net/cgi-bin/mimetex.cgi?f_{\!\!X} là density function, http://dientuvietnam.net/cgi-bin/mimetex.cgi?F_{\!\!X} là distribution function của http://dientuvietnam.../mimetex.cgi?X. Ta có http://dientuvietnam.net/cgi-bin/mimetex.cgi?f_{\!\!X}(x)=F_{\!\!X}^{'}(x). Dựa trên định nghĩa http://dientuvietnam.net/cgi-bin/mimetex.cgi?F_{\!\!X}(x)=P(X<x), rồi lấy đạo hàm, ta chứng minh được rằng
Bổ đề 1a : Nếu http://dientuvietnam.net/cgi-bin/mimetex.cgi?f_{\!X_1}(z)=f_{\!X_2}(z)=\exp(-z) (normal exponential) vào , rồi tính tích phân (dùng change of variable http://dientuvietnam.net/cgi-bin/mimetex.cgi?t=\exp(-x), tích phân mới đi từ 0 đến http://dientuvietnam.net/cgi-bin/mimetex.cgi?+\infty) thì được
, đpcm.
Câu kia tương tự (Weibull)
Bổ đề 1a : Nếu http://dientuvietnam.net/cgi-bin/mimetex.cgi?f_{\!X_1}(z)=f_{\!X_2}(z)=\exp(-z) (normal exponential) vào , rồi tính tích phân (dùng change of variable http://dientuvietnam.net/cgi-bin/mimetex.cgi?t=\exp(-x), tích phân mới đi từ 0 đến http://dientuvietnam.net/cgi-bin/mimetex.cgi?+\infty) thì được
, đpcm.
Câu kia tương tự (Weibull)
Bài viết đã được chỉnh sửa nội dung bởi TieuSonTrangSi: 22-10-2005 - 01:26
Chí lớn trong thiên hạ không đựng đầy đôi mắt của giai nhân
#4
Đã gửi 21-10-2005 - 23:25
#5
Đã gửi 21-10-2005 - 23:30
#6
Đã gửi 21-10-2005 - 23:41
#7
Đã gửi 21-10-2005 - 23:43
#8
Đã gửi 21-10-2005 - 23:46
#9
Đã gửi 21-10-2005 - 23:49
#10
Đã gửi 21-10-2005 - 23:50
9. [2, problem 5, p. 224]
If T has full row rank prove that
If T has full row rank prove that
<span style='color:blue'>Roses are red,
violets are blue,
Fermat is dead,
but his theorem is true. </span>
violets are blue,
Fermat is dead,
but his theorem is true. </span>
#11
Đã gửi 22-10-2005 - 00:43
Không biết em có nhầm không, nhưng mà bài này em thấy giống vớiBài 9
#12
Đã gửi 22-10-2005 - 00:49
T chưa chắc là ma trận vuông đâu Lim.Không biết em có nhầm không, nhưng mà bài này em thấy giống vớiBài 9
<span style='color:blue'>Roses are red,
violets are blue,
Fermat is dead,
but his theorem is true. </span>
violets are blue,
Fermat is dead,
but his theorem is true. </span>
#13
Đã gửi 22-10-2005 - 00:56
T' là cái gì vậy anh ?T chưa chắc là ma trận vuông đâu Lim.
Không biết em có nhầm không, nhưng mà bài này em thấy giống vớiBài 9
#14
Đã gửi 22-10-2005 - 01:14
T' là ma trận chuyển vị (transpose) của T (hàng thành cột). Ờ ha... mà để nhân được với nhau... để anh coi lại ...
À à coi lại rồi, tổng quát T là ma trận chữ nhật, thí dụ kích thước (m x n), T' sẽ là ma trận (n x m) --> T' T = (n x m)x(m x n) = (n x n) --> vuông, lúc đó mới có thể nghịch đảo.
À à coi lại rồi, tổng quát T là ma trận chữ nhật, thí dụ kích thước (m x n), T' sẽ là ma trận (n x m) --> T' T = (n x m)x(m x n) = (n x n) --> vuông, lúc đó mới có thể nghịch đảo.
<span style='color:blue'>Roses are red,
violets are blue,
Fermat is dead,
but his theorem is true. </span>
violets are blue,
Fermat is dead,
but his theorem is true. </span>
#15
Đã gửi 22-10-2005 - 01:40
Chỗ em dùng ký hiệu cho transpose matrix , tiếng Việt mình gọi là ma trận chuyển vị àh, hay nhỉ.T' là ma trận chuyển vị (transpose) của T (hàng thành cột). Ờ ha... mà để nhân được với nhau... để anh coi lại ...
À à coi lại rồi, tổng quát T là ma trận chữ nhật, thí dụ kích thước (m x n), T' sẽ là ma trận (n x m) --> T' T = (n x m)x(m x n) = (n x n) --> vuông, lúc đó mới có thể nghịch đảo.
Bài toán 8 có thể làm tương tự không anh, lấy inverse của vế trái, nó sẽ ra vế phải.
#16
Đã gửi 22-10-2005 - 01:43
Ui da, sao mà nhiều wá xá wà xa :cry2. The probable error, http://dientuvietnam...n/mimetex.cgi?M is the median of http://dientuvietnam.../mimetex.cgi?Y. Find the probable errors of
1. the exponential distribution;
Bài 2 tôi làm như sau. Tạm xem http://dientuvietnam...n/mimetex.cgi?M là một hằng số. Sự kiện (event) http://dientuvietnam.net/cgi-bin/mimetex.cgi?\{|Y-M|\,>\,\tau\} tương đương với
http://dientuvietnam.net/cgi-bin/mimetex.cgi?F_Y cho distribution function của http://dientuvietnam...n/mimetex.cgi?Y (như trong bài 1), ta suy ra
http://dientuvietnam.net/cgi-bin/mimetex.cgi?\tau. Trong phương trình này, vì http://dientuvietnam...n/mimetex.cgi?M là median, nên http://dientuvietnam...n/mimetex.cgi?M được xác định bởi
http://dientuvietnam.net/cgi-bin/mimetex.cgi?F_Y, ta thế vào (***) để tìm http://dientuvietnam.net/cgi-bin/mimetex.cgi?M, rồi dùng (**) để tính http://dientuvietnam.net/cgi-bin/mimetex.cgi?\tau. Ví dụ, trong câu 1, ta có
http://dientuvietnam.net/cgi-bin/mimetex.cgi?F_Y(y)=1-\exp(-y).
Dùng (***), ta có http://dientuvietnam.net/cgi-bin/mimetex.cgi?\sinh\tau=\dfrac{1}{2}, tức là http://dientuvietnam.net/cgi-bin/mimetex.cgi?\tau=\mbox{argsh}\dfrac{1}{2}=\dfrac{1+\sqrt{5}}{2}.
Chí lớn trong thiên hạ không đựng đầy đôi mắt của giai nhân
#17
Đã gửi 22-10-2005 - 02:11
Nếu thuận tiện anh 2TS viết bài giải tiếng Anh được thì hay, có gì mai mốt copy qua bên trang tiếng Anh là ngon (à mà không biết bên kia có latex không).
@Lim: để chỉ transpose matrix, bên VN cũng thường dùng chữ T chỉ số trên (superscript) , trong khi sách Mỹ thì hay dùng dấu phẩy trên (prime).
<span style='color:blue'>Roses are red,
violets are blue,
Fermat is dead,
but his theorem is true. </span>
violets are blue,
Fermat is dead,
but his theorem is true. </span>
#18
Đã gửi 01-11-2005 - 08:17
[quote name='Alligator' date='Oct 21 2005, 12:50 PM']9. [2, problem 5, p. 224]
If T has full row rank prove that
http://dientuvietnam...metex.cgi?G_{1} is the unique Penrose generalized inverse of http://dientuvietnam...metex.cgi?G_{2} is the unique Penrose generalized inverse of http://dientuvietnam...mimetex.cgi?G_1 and http://dientuvietnam...mimetex.cgi?G_2 are the unique Penrose generalized inverse of , they should be the same. Therefore
If T has full row rank prove that
http://dientuvietnam...metex.cgi?G_{1} is the unique Penrose generalized inverse of http://dientuvietnam...metex.cgi?G_{2} is the unique Penrose generalized inverse of http://dientuvietnam...mimetex.cgi?G_1 and http://dientuvietnam...mimetex.cgi?G_2 are the unique Penrose generalized inverse of , they should be the same. Therefore
<span style='color:blue'>Roses are red,
violets are blue,
Fermat is dead,
but his theorem is true. </span>
violets are blue,
Fermat is dead,
but his theorem is true. </span>
1 người đang xem chủ đề
0 thành viên, 1 khách, 0 thành viên ẩn danh