请教dyadic data数据处理

chaoswang

Kenny好，各位前辈好！
5 p8 d+ e( V, N6 c2 }
小弟最近在处理谈判数据，即dyadic data，通行的做法似乎是David Kenny的actor-parter interdependence model，这个模型主要目的是要考虑谈判双方相互的interdependence，这个模型一开始用SAS或者HLM运行的。但是我们想检测两条中介路径（并且比较两者大小），感觉用mplus比较方便一些，不知道用mplus的twolevel random是否是可以得到跟HLM一样的结果？
/ W! h3 P7 G* ~- ^
  `+ T$ m+ w9 U) T/ o6 ~" f* b" R! S另外我还有个问题：dyadic data用multilevel modeling会不会有问题？每个group只有两个人，这样within group的结果有意义么？

8 u5 i. u0 z5 F0 w: s谢谢大家！
; Y: ]% S4 ~3 v+ ^8 F

chaoswang

还望大家多多关注。

xinting.J

chaoswang 发表于 2014-8-25 10:22
* k$ ^: P4 T4 ?  |还望大家多多关注。

chaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着找到reference 再回你。

5 E1 }% i$ D' o7 T& [: J; {) }不记得在哪里看过了，用HLM时可以做dyadic data分析的，但因为每个组里只有两个数据点，存在模型识别的问题，所以需要对数据做一点处理。两种处理方法：一种是如果你用的是量表（而且题目足够多），可以把量表拆成两半，在每组中人为的构造出两个数据点，拆出的一半作为它们的测量；另一种是自己计算Y的error variance, 公式EV =(1- Reliability of Y)x Var (Y)  (根据信度计算的原理)，本来如果每组数足够多，系统会自己计算的，但是只有两个数据点，计算会有问题，只能自己先算好。
% }, z4 ^8 L" [( v9 H* O
, Z5 n1 d5 ]: a6 E+ L另外有一篇介绍用HLM做dyadic data的文章你可以参考，但是这篇文章里也没有提到上面的方法:
; |% R6 M* `1 DCampbell, L., & Kashy, D. A. (2002). Estimating actor, partner, and interaction effects for dyadic data using PROC MIXED and HLM: A user–friendly guide. Personal Relationships, 9(3), 327-342.

. G' a! w) ?1 S, B* y希望不会让你更糊涂了。

chaoswang

xinting.J 发表于 2014-8-25 14:18
chaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着 ...

Jane 你好，谢谢你的回答。

我知道我的问题比较小众，所以大家都比较谦虚……

- Q8 Y" g) h2 Z& M1 q/ l我现在确实是用 Campbell（2002）的方法处理数据，但我现在有两点和他们不一样：1）我们加了中介变量，2）因为有中介变量，为了方便，我们用Mplus的TWOLEVEL RANDOM代替HLM来运行模型。但是前人似乎没有这种做法，所以我们不太有把握是不是能这么做。

您说的两个数据点存在模型识别的问题我没有太明白，但确实我做模型的软件没有报错，这是怎么回事呢？T_T



1 E  }6 i5 _/ Z( }8 M

chaoswang

xinting.J 发表于 2014-8-25 14:18
. h6 M& V/ Q$ O0 jchaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着 ...

Campbell这篇文章里给出的方法（以及其他用Multilevel modeling处理dyadic data的方法)都有一个共同点，两个自变量（通常情况下是actor IV 和 partner IV）在整体上的相关为-1。因为对每个dyad的成员，成员A的partner是B，成员B的partner是A。但是他们通过将actor IV和partner IV作为within dyad variable，就将这两个自变量同时放入回归方程。

我们加的中介变量之间的相关性也是-1，因此也想用这个方法，将他们局限在within dyad中分析（当然他们确实在dyad层面没有variance），但是不知道这样行不行……

chaoswang

- B  m! E# Y# X5 ]& _1 K& u
6 B- ?9 W1 X( Z5 n  S! j) V$ v2 O今天用mplus做上述模型时，得到下面这个warning:

& P' |# _7 o1 i" W$ HTHE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES MAY NOT BE
     TRUSTWORTHY FOR SOME PARAMETERS DUE TO A NON-POSITIVE DEFINITE
0 V( w0 z2 A* E2 Z& Y, O     FIRST-ORDER DERIVATIVE PRODUCT MATRIX.  THIS MAY BE DUE TO THE STARTING
" t) ~5 @( U! |2 Z6 a3 l     VALUES BUT MAY ALSO BE AN INDICATION OF MODEL NONIDENTIFICATION.  THE
; a9 X) o8 m" O: S) ]+ Z5 L     CONDITION NUMBER IS      -0.178D-15.  PROBLEM INVOLVING PARAMETER 2.
我去网上查了，似乎需要检查模型是否可以识别（identified），但是怎么去检查呢？不是只要自由度不小于0模型都可以识别么？
另外，我不太明白警告里的“A NON-POSITIVE DEFINITE FIRST-ORDER DERIVATIVE PRODUCT MATRIX”是什么意思……
' o. q  o, N4 y
我去查看出问题的参数 parameter 2，发现是模型其中一个中介变量的intercept，发现这个参数只是不显著而已，没什么特别的地方啊。

xinting.J

chaoswang 发表于 2014-8-27 11:03
% a$ I1 R4 }% E) Q3 x% W) `9 h今天用mplus做上述模型时，得到下面这个warning:

THE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES  ...

应该就是模型不可识别的问题，但从这个信息看不出是哪里的问题。
我觉得包含中介的和两个变量的直接关系原理相同，只要都是同一层面的，应该问题不大。
你有试过我前面说的方法增加每组的data point或者自己计算standard error吗？

Kenneth

chaoswang,
/ g  ]0 s/ W. ~% J8 L4 N
' o! A9 ]4 p9 L# ~如果是因为起始值的问题，可以试试放弃模型中其他变量，只是把有问题的线性函数做回归，得出截距后，用该截距作为起始值再看看。

chaoswang

xinting.J 发表于 2014-8-30 03:10
应该就是模型不可识别的问题，但从这个信息看不出是哪里的问题。
我觉得包含中介的和两个变量的直接关系 ...

2 a/ u% t7 M& t7 _$ F  \' }我们模型就是个路径分析，没有潜变量，每个变量只有一个测量值（都是客观指标），是不是没办法用您说的两个方法呢？
2 f$ C2 _  H) J
我先试试Kenny说的设置初始值，谢谢您！
( m6 V$ k3 k7 W, C/ C0 }5 \' a, r, o

chaoswang

Kenneth 发表于 2014-8-31 08:02
chaoswang,

7 S# d5 O+ y) J9 ~' V" {( H6 P2 w如果是因为起始值的问题，可以试试放弃模型中其他变量，只是把有问题的线性函数做回归，得出截 ...

好的，谢谢您！我试试看！