中人网

标题: 非對稱效果的估計 [打印本页]

作者: chienhsin    时间: 2011-6-17 06:04
标题: 非對稱效果的估計
Kenny你好:
% @8 f, m. Z1 K7 M; B) X- q2 A6 f0 S0 h- d$ D7 j
我最近在一篇文章內看到作者如此敘述他的方法:% M  ~. C; z6 {' J$ x. o% |8 r% m
先說明資料結構,每一筆資料有Y  A  B
. C  r4 ?; \! c/ n0 A; D作者要估計 (A-B) 差是正值或負值時對Y的效果。
' d/ T3 U; {3 e( A9 K/ d. U" g- C5 T" b1 d5 L, ]' D9 @- Y5 l, i
we fitted a random-effects piecewise logit model on the data set using the difference between A and B as an independent variable and accounting for potentially different slopes when the A-B difference was positive versus negative.+ b$ ?7 A2 z$ i8 k8 r( e
- U) p9 E- h5 J  j' N5 J+ u
分析後,結果同時顯示了當A-B為正、A-B為負,之結果。# e& @) X; E" P4 X! u) L

" V' m' P$ j7 c5 L& I* ?& A. D我不懂的是作者如何處理A-B,A-B不是正值就是負值,或是0,當A-B是正值有資料的時候,# v, ]8 W1 U9 M5 p/ L6 M6 \4 D
負值就是missing value,vice versa, 因此在同一筆資料上,怎麼同時會有(A-B)為正值、為負值的資料?
) Z4 |8 l/ f/ N+ [/ q7 F5 r/ a
( ^/ k5 N. f- p, J9 K我本來想作者是將資料分兩半分別估計,但他的結果卻只顯示一個intercept,其他係數也是一個、R-square也一個,所以也不是分開估計。: h+ A( \0 W7 L* w* q

- }  i$ t: V8 ]* Z, ~$ t. h謝謝Kenny的幫忙。; {7 b  r" u9 x  @# R+ H' f
# {; u6 w, i; i, X1 i, Q; M8 z
' c1 v: J: z  H+ Z2 ]
Chien-Hsin: W" M8 Z3 @/ t" x& ~$ e0 g9 k

  j, |! q: v8 i* `$ z5 \# o 本帖最后由 chienhsin 于 2011-6-17 06:46 编辑
# X# u+ W; C( W' Y* c* D6 O4 I2 C- A7 r4 n5 Y

作者: Kenneth    时间: 2011-6-17 13:04
回复 1楼 chienhsin 的帖子" D1 n1 q( U: J( i
Chienhsin, 你给我的资料不够深入,所以我依你的描述上网找了一下。我猜我找到你所说的文章了。
. \! k" {6 k# K- {# G+ n我猜他们用的是一种特别的回归,叫做 piecewise regression (random effect and logit 在这里不重要。logit 是因为因变量是一个虚拟变量;random effect 是他们选择了的估计方法)。我在网上看了一下,piecewise regression 好像是特别用来解决一些问题,问题的性质是当 x<a 时做一个回归分析估计(a 是一个常数);当 x>a 时做一个回归分析估计。但是两组数据是“同时估计的”,所以只有一个截距。我想最简单的理解是用同一组数据作出两个回归分析(就是两条回归直线)当x<a 时一条;当x<a 时另外一条。但是它们有相同的截距,而且两条回归是同时估计的。不知道我讲的对不对,因为我从来没有在管理的研究中见过。
! C  x3 ~" O8 B" C, h/ j. [6 }0 _8 t" t
   
作者: chienhsin    时间: 2011-6-17 20:54
回复 2楼 Kenneth 的帖子2 S$ u9 c9 c6 _5 {. u

0 H  c" o# d' T! H( PKenny您好:( ]& K+ X+ l- ]( S; m# y3 m# [5 r

: \* a. i6 x6 S: G6 g我根據您的指示,也找了相關的資料,原來在我的迴歸分析的課本也有簡短提到,我的課本把這種方法稱為spline regression.
* _! p1 }( |* c8 L) Q) z. m# P- Q! a0 M  t  s1 ]
現在我依照我原來的問題把模式設定一下:- Y* U; w& H/ b) R5 c6 `7 i
* B+ d. `/ U: A' N- x
Let  X = (A-B)3 M, B! K* i! R9 L  ?3 v% `
     A1 = max (X-0, 0)         <<因為該文章要研究的是(A-B) > 0 和 < 0 的效果差異>>) i  Q" r: t" P/ m0 w! Q( S2 Z+ z
& C% D  _& A- _, s( k7 T# x
     Y = b0 + b1X + b2A1 + e
7 H9 @. ]2 E- d
6 h: l1 h; [- ?( U" M& `& {, ]   for (A-B) > 0
. l# q* _* q. b! A" v, _       Y = b0 + b1(A-B) + b2(A-B)
. m* ?- k2 X5 N) e- d8 N) ^          = b0 + (b1+b2) (A-B)
6 W; ^) w' A" x1 Y( c6 L: y' n7 I
, C& ~. I* c  @  i   for (A-B) <=0
- n6 t4 v& J0 v) C7 n! Z+ X       Y = b0 + b1(A-B)
7 C4 v( j" C5 p& u        . N8 a% b- |6 z% G, ^; V8 d; ~
比較效果的差異:
# F5 o  p" Q/ k. M, O+ y( ~9 v(A-B)負值時: 效果為b1
9 Z# v$ e. P0 `1 D! i( X: Z(A-B)正值時: 效果為b1+b2! W% q5 W1 M$ f3 ?  W" L

4 }0 y! h" X. @- {4 o我想應該是這個樣子。
9 [. t0 ~  }. L$ @% O但還有個不明瞭的地方,原文的結果報告分別對) i7 m* }# H, n8 P8 `
(A-B)正值的效果、(A-B)負值的效果 分別報告顯著性,
6 c0 `+ d3 j4 Y) t* P# [但從上面的model來看,我僅能檢定b1, b2的顯著性,並無法檢定(b1+b2)的顯著性?
5 \1 x- h; w2 a
6 f2 k+ }' K6 T- ~' vKenny有任何看法嗎?; u$ Y7 H3 a4 R* a* ?9 B

5 }* @2 B8 M* q
9 A1 C& K  m  `0 l, H: i/ \8 ]$ s. T. e9 x

5 G8 A$ e$ {+ S0 h) }- q! l4 _* L   
作者: Kenneth    时间: 2011-6-18 11:29
回复 3楼 chienhsin 的帖子  f; ?' {8 T) N6 z8 c" [0 g% @% ^
chienhsin,我猜你说的是“可以检验b1的显著性,和(b1+b2)的显著性”吧。
0 {- V; L) E5 E4 N& j- p. w一般当我们知道b1和b2的标准差时,只要作一点假设(比如 multivariate normality 等),(b1+b2)的标准差是可以用公式估计出来的。无论如何,这已经是一个统计学的问题了。我不是念统计的,不知道答案。我也试过在网上替你找,可是找到的也不多。反正我觉得这不是一个常见的问题,要用的时候才花时间去找吧。
5 F( B9 w1 e# W& ^9 |4 ?( v! k# o- D6 i! z, n
   
作者: chienhsin    时间: 2011-6-19 07:09
回复 4楼 Kenneth 的帖子
* r$ _5 w8 \6 k7 B" t1 W. E% M
- d0 u9 X& Z7 e1 m* y( X+ m1 ^. D% J, F# p3 b, K$ B, i" p3 B
    Kenny謝謝你。




欢迎光临 中人网 (http://bbs.chinahrd.net/) Powered by Discuz! X2.5