中人网

标题: 请教dyadic data数据处理 [打印本页]
作者: chaoswang    时间: 2014-8-22 11:14
标题: 请教dyadic data数据处理
Kenny好，各位前辈好！4 M' d6 |" B8 F8 j
$ m  C# f1 a6 k
小弟最近在处理谈判数据，即dyadic data，通行的做法似乎是David Kenny的actor-parter interdependence model，这个模型主要目的是要考虑谈判双方相互的interdependence，这个模型一开始用SAS或者HLM运行的。但是我们想检测两条中介路径（并且比较两者大小），感觉用mplus比较方便一些，不知道用mplus的twolevel random是否是可以得到跟HLM一样的结果？# x: o( W: p' ~
4 u" i, n" F0 V9 z* Q+ n2 ?
另外我还有个问题：dyadic data用multilevel modeling会不会有问题？每个group只有两个人，这样within group的结果有意义么？; x, A" n8 v/ p2 {7 P

5 }: A4 g6 `4 D  ?* J谢谢大家！( M* |# O2 C+ o0 U. ~
作者: chaoswang    时间: 2014-8-25 10:22
还望大家多多关注。
作者: xinting.J    时间: 2014-8-25 14:18
chaoswang 发表于 2014-8-25 10:22
- c8 J2 l; B  C还望大家多多关注。

; ~" @; X% g6 R# Ochaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着找到reference 再回你。
: h, |& q# F1 E: G) N% C$ D& i6 T5 X" j* Q9 b3 o$ K5 P/ D
不记得在哪里看过了，用HLM时可以做dyadic data分析的，但因为每个组里只有两个数据点，存在模型识别的问题，所以需要对数据做一点处理。两种处理方法：一种是如果你用的是量表（而且题目足够多），可以把量表拆成两半，在每组中人为的构造出两个数据点，拆出的一半作为它们的测量；另一种是自己计算Y的error variance, 公式EV =(1- Reliability of Y)x Var (Y)  (根据信度计算的原理)，本来如果每组数足够多，系统会自己计算的，但是只有两个数据点，计算会有问题，只能自己先算好。. P$ `9 H9 [; l& u: |; a

) R+ Y, V8 G+ H, ]4 j9 X  X另外有一篇介绍用HLM做dyadic data的文章你可以参考，但是这篇文章里也没有提到上面的方法:
* _! Y4 }( p  z; z1 X0 Q# HCampbell, 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.
* \/ ^2 [( ]+ s6 I* \! [3 t% D4 N2 m2 b2 v
希望不会让你更糊涂了。
作者: chaoswang    时间: 2014-8-25 21:21
xinting.J 发表于 2014-8-25 14:18
% s9 D0 ~# J  [5 p; n5 _chaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着 ...

+ ?' x* K' d: s# K% |" MJane 你好，谢谢你的回答。
& Q7 ^6 k+ g5 C; m* P1 B( `1 q, N
我知道我的问题比较小众，所以大家都比较谦虚……5 g2 m( m( u5 s+ ^( {
- d8 h5 R) ]( ^- g% `( B2 h
我现在确实是用 Campbell（2002）的方法处理数据，但我现在有两点和他们不一样：1）我们加了中介变量，2）因为有中介变量，为了方便，我们用Mplus的TWOLEVEL RANDOM代替HLM来运行模型。但是前人似乎没有这种做法，所以我们不太有把握是不是能这么做。% d( R/ {3 p) ~7 ?
: I4 C5 L) e. M2 d
您说的两个数据点存在模型识别的问题我没有太明白，但确实我做模型的软件没有报错，这是怎么回事呢？T_T
: Z$ Q! k  m7 c5 j8 n; y2 O/ K
& N* {5 ]$ J) R2 K# @
8 D8 m: s6 h4 z  j4 H- {: U4 _* y: @5 M
作者: chaoswang    时间: 2014-8-25 21:36
xinting.J 发表于 2014-8-25 14:18
- L4 d7 i' N7 P  D4 {, Q2 dchaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着 ...

6 d4 m$ g& ?6 G5 _6 l8 C* FCampbell这篇文章里给出的方法（以及其他用Multilevel modeling处理dyadic data的方法)都有一个共同点，两个自变量（通常情况下是actor IV 和 partner IV）在整体上的相关为-1。因为对每个dyad的成员，成员A的partner是B，成员B的partner是A。但是他们通过将actor IV和partner IV作为within dyad variable，就将这两个自变量同时放入回归方程。
3 x: P) I  q* Y* E2 K! H; l/ r" e7 e7 @, S* {
我们加的中介变量之间的相关性也是-1，因此也想用这个方法，将他们局限在within dyad中分析（当然他们确实在dyad层面没有variance），但是不知道这样行不行……
作者: chaoswang    时间: 2014-8-27 11:03
本帖最后由 chaoswang 于 2014-8-27 11:24 编辑 . W2 o4 Q# V- M  |

: M8 p+ l% L2 y' a今天用mplus做上述模型时，得到下面这个warning:1 p/ V9 Y1 C  m% c2 E- G

/ I3 |6 o' ]' s9 a# E3 {" ETHE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES MAY NOT BE' y( X0 A) ^# h: `4 X4 |5 Q# n* u
     TRUSTWORTHY FOR SOME PARAMETERS DUE TO A NON-POSITIVE DEFINITE
, L6 B$ v8 }9 `$ y2 m     FIRST-ORDER DERIVATIVE PRODUCT MATRIX.  THIS MAY BE DUE TO THE STARTING
0 Q$ d8 H! }4 f9 Z5 ^1 ^7 z, ~: [     VALUES BUT MAY ALSO BE AN INDICATION OF MODEL NONIDENTIFICATION.  THE- I0 q2 v  Q/ C3 z4 B' m
     CONDITION NUMBER IS      -0.178D-15.  PROBLEM INVOLVING PARAMETER 2.: P) o( S+ i7 j2 S
我去网上查了，似乎需要检查模型是否可以识别（identified），但是怎么去检查呢？不是只要自由度不小于0模型都可以识别么？
0 _5 ]/ C' p9 b" m另外，我不太明白警告里的“A NON-POSITIVE DEFINITE FIRST-ORDER DERIVATIVE PRODUCT MATRIX”是什么意思……
# G9 Z! @' H' C% p- ^3 m1 |- E+ Y
( m9 X% ?) |3 {5 W. o我去查看出问题的参数 parameter 2，发现是模型其中一个中介变量的intercept，发现这个参数只是不显著而已，没什么特别的地方啊。
作者: xinting.J    时间: 2014-8-30 03:10
chaoswang 发表于 2014-8-27 11:03 5 m5 W$ w0 G! R; r5 C  a
今天用mplus做上述模型时，得到下面这个warning:
, Y! A" ~2 t4 y3 k& D: A7 X# H+ F: d/ P& N
THE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES  ...

: v7 ?+ X  c# f8 t, H! W, V应该就是模型不可识别的问题，但从这个信息看不出是哪里的问题。
( ^2 O& ]# }5 W" b- b9 V我觉得包含中介的和两个变量的直接关系原理相同，只要都是同一层面的，应该问题不大。
. ?% Z" O( o. V# e2 L1 K你有试过我前面说的方法增加每组的data point或者自己计算standard error吗？
作者: Kenneth    时间: 2014-8-31 08:02
chaoswang,
+ C& Z, p5 ]( W! w, I  ^) R) n3 l! y: h: n5 V
如果是因为起始值的问题，可以试试放弃模型中其他变量，只是把有问题的线性函数做回归，得出截距后，用该截距作为起始值再看看。
作者: chaoswang    时间: 2014-9-1 11:38
xinting.J 发表于 2014-8-30 03:10
" c5 M( o! U: D7 w" o  _, F应该就是模型不可识别的问题，但从这个信息看不出是哪里的问题。6 Y5 S4 u! G. B( b" d, c7 A
我觉得包含中介的和两个变量的直接关系 ...

2 t5 Q+ k* v+ W" ^9 W/ K我们模型就是个路径分析，没有潜变量，每个变量只有一个测量值（都是客观指标），是不是没办法用您说的两个方法呢？- ~" W0 j) I- ^2 i' y% \

, T8 s- s" _! e( y+ H我先试试Kenny说的设置初始值，谢谢您！
$ a" l% Y# u: {: Y/ d4 r+ A
作者: chaoswang    时间: 2014-9-1 11:39
Kenneth 发表于 2014-8-31 08:02 # r- B1 n, x* u4 z% u' I% b3 `! p/ n
chaoswang,
0 M) ?" Q1 l) e6 E: z
; W0 D5 Q# U+ P1 o: p如果是因为起始值的问题，可以试试放弃模型中其他变量，只是把有问题的线性函数做回归，得出截 ...

6 H7 q+ S& d' _/ R! d好的，谢谢您！我试试看！
作者: chaoswang    时间: 2014-9-11 14:02
哈哈，原来出问题的原因是我这两个中介变量的variance和intercept应该设为相等，设好之后问题就没有了，太开心了！
( Z$ h( J4 a% E$ }* r谢谢各位！



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