中人网

标题: 请教dyadic data数据处理 [打印本页]
作者: chaoswang    时间: 2014-8-22 11:14
标题: 请教dyadic data数据处理
Kenny好，各位前辈好！) u! p: q* a- j

8 O8 \' d: o; P7 e0 w* x# U0 y小弟最近在处理谈判数据，即dyadic data，通行的做法似乎是David Kenny的actor-parter interdependence model，这个模型主要目的是要考虑谈判双方相互的interdependence，这个模型一开始用SAS或者HLM运行的。但是我们想检测两条中介路径（并且比较两者大小），感觉用mplus比较方便一些，不知道用mplus的twolevel random是否是可以得到跟HLM一样的结果？: j, {5 e# [) o- i% h( s3 T+ P# R$ k

6 e. y" @- s' v( q" M另外我还有个问题：dyadic data用multilevel modeling会不会有问题？每个group只有两个人，这样within group的结果有意义么？
) S' z8 g8 z) O2 z
$ J- k) |* u. }9 e谢谢大家！
' z/ W! d2 T7 t4 \6 e6 P8 e' z
作者: chaoswang    时间: 2014-8-25 10:22
还望大家多多关注。
作者: xinting.J    时间: 2014-8-25 14:18
chaoswang 发表于 2014-8-25 10:22
: b+ h& i7 J% I1 N7 ?6 L+ Z7 }还望大家多多关注。

# w* h2 n5 `1 t3 P, I7 Jchaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着找到reference 再回你。$ z! R+ I' b7 v, w- F( t
5 E* ?! `1 o0 T$ A9 k
不记得在哪里看过了，用HLM时可以做dyadic data分析的，但因为每个组里只有两个数据点，存在模型识别的问题，所以需要对数据做一点处理。两种处理方法：一种是如果你用的是量表（而且题目足够多），可以把量表拆成两半，在每组中人为的构造出两个数据点，拆出的一半作为它们的测量；另一种是自己计算Y的error variance, 公式EV =(1- Reliability of Y)x Var (Y)  (根据信度计算的原理)，本来如果每组数足够多，系统会自己计算的，但是只有两个数据点，计算会有问题，只能自己先算好。  U. e1 e& @$ R4 K+ |' N

- H, b: j3 q/ T: }另外有一篇介绍用HLM做dyadic data的文章你可以参考，但是这篇文章里也没有提到上面的方法:1 o) q( r/ P2 u( f& u, Q
Campbell, 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.4 I4 h0 F% S% k( R- |& ^

9 g8 I& x! `% e. b希望不会让你更糊涂了。
作者: chaoswang    时间: 2014-8-25 21:21
xinting.J 发表于 2014-8-25 14:18 . [- Q1 i3 ?' N1 ]' x  Y
chaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着 ...
. ?& G! C: O" D' X  O, a# U+ F1 p
Jane 你好，谢谢你的回答。
) l, d- k( Z: V' z6 M0 y
0 V$ c, W0 n4 ^& |! |4 J我知道我的问题比较小众，所以大家都比较谦虚……
' F- C1 K7 T7 B  {. P8 x9 ~2 D) a/ o  E2 N9 ^
我现在确实是用 Campbell（2002）的方法处理数据，但我现在有两点和他们不一样：1）我们加了中介变量，2）因为有中介变量，为了方便，我们用Mplus的TWOLEVEL RANDOM代替HLM来运行模型。但是前人似乎没有这种做法，所以我们不太有把握是不是能这么做。( m3 q+ N7 y' L3 W

! T$ ~+ d; e% s4 O( D- L4 _$ Y您说的两个数据点存在模型识别的问题我没有太明白，但确实我做模型的软件没有报错，这是怎么回事呢？T_T
9 ^+ q9 ~' s6 i; m6 M+ t; U. R. E4 p4 c5 K+ A) P

4 ?* h4 a3 S' D
# O. I: ]' r3 R+ e
作者: chaoswang    时间: 2014-8-25 21:36
xinting.J 发表于 2014-8-25 14:18
0 m) p5 ]" x+ X: i# B) Hchaswang, 昨天本来想回你，但是因为只有一个印象，没有找到文献，我自己也没有做过，不敢贸然建议，想着 ...

8 k9 w5 Z+ `/ o. I/ t( z( LCampbell这篇文章里给出的方法（以及其他用Multilevel modeling处理dyadic data的方法)都有一个共同点，两个自变量（通常情况下是actor IV 和 partner IV）在整体上的相关为-1。因为对每个dyad的成员，成员A的partner是B，成员B的partner是A。但是他们通过将actor IV和partner IV作为within dyad variable，就将这两个自变量同时放入回归方程。
0 J0 [, \6 O' g: I3 M
0 Y. X% ^/ L; w& \我们加的中介变量之间的相关性也是-1，因此也想用这个方法，将他们局限在within dyad中分析（当然他们确实在dyad层面没有variance），但是不知道这样行不行……
作者: chaoswang    时间: 2014-8-27 11:03
本帖最后由 chaoswang 于 2014-8-27 11:24 编辑
: ?, R/ X6 G3 r) C2 K. G& x. n5 m- C: d) K$ W
今天用mplus做上述模型时，得到下面这个warning:; F" j8 U* {% |2 v0 z, f/ ^

. [( _3 |" F$ V& h  N0 J4 ]THE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES MAY NOT BE
6 K% Y8 q4 E" D- o7 n5 O  h, ]) d     TRUSTWORTHY FOR SOME PARAMETERS DUE TO A NON-POSITIVE DEFINITE% ]# k% P. z& G8 |+ J) T" ]/ d
     FIRST-ORDER DERIVATIVE PRODUCT MATRIX.  THIS MAY BE DUE TO THE STARTING
4 o/ h( {$ K* c- ^+ ^% Y$ H     VALUES BUT MAY ALSO BE AN INDICATION OF MODEL NONIDENTIFICATION.  THE
  ]3 `9 h' |: q, r) k     CONDITION NUMBER IS      -0.178D-15.  PROBLEM INVOLVING PARAMETER 2.
$ Y! \0 c. D& J6 T% ^2 x我去网上查了，似乎需要检查模型是否可以识别（identified），但是怎么去检查呢？不是只要自由度不小于0模型都可以识别么？
7 a0 u9 b% x' s+ z  O: D( Y另外，我不太明白警告里的“A NON-POSITIVE DEFINITE FIRST-ORDER DERIVATIVE PRODUCT MATRIX”是什么意思……& X: @+ d/ c2 |7 v: f  L5 r
! b9 x3 Z8 n& e
我去查看出问题的参数 parameter 2，发现是模型其中一个中介变量的intercept，发现这个参数只是不显著而已，没什么特别的地方啊。
作者: xinting.J    时间: 2014-8-30 03:10
chaoswang 发表于 2014-8-27 11:03 % }8 ]6 s  @1 X- Z2 Z2 c+ E! \. ]* f
今天用mplus做上述模型时，得到下面这个warning:+ {; V4 Z' g# O+ \" y
$ R! j9 S, P+ l4 Y
THE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES  ...

( l3 p9 o, g* n3 n! z应该就是模型不可识别的问题，但从这个信息看不出是哪里的问题。2 f" ~" v6 C0 l
我觉得包含中介的和两个变量的直接关系原理相同，只要都是同一层面的，应该问题不大。
" }( g6 x* [; Q0 x/ l4 w你有试过我前面说的方法增加每组的data point或者自己计算standard error吗？
作者: Kenneth    时间: 2014-8-31 08:02
chaoswang,
! q- s! P/ k0 D' w  y9 d$ j) [. T, L- p
如果是因为起始值的问题，可以试试放弃模型中其他变量，只是把有问题的线性函数做回归，得出截距后，用该截距作为起始值再看看。
作者: chaoswang    时间: 2014-9-1 11:38
xinting.J 发表于 2014-8-30 03:10 # Z" m6 q4 B# L; f" C3 V8 |  V
应该就是模型不可识别的问题，但从这个信息看不出是哪里的问题。
$ ?# G0 v0 V8 i, P: j- M$ G+ X9 I" D我觉得包含中介的和两个变量的直接关系 ...

' k) w- C% f; g) k5 A我们模型就是个路径分析，没有潜变量，每个变量只有一个测量值（都是客观指标），是不是没办法用您说的两个方法呢？
5 U9 A4 U' m; J% E9 {7 j
; H9 ~, `' a. i" ?9 u我先试试Kenny说的设置初始值，谢谢您！# g4 e7 O+ ]- K" t* M" H
作者: chaoswang    时间: 2014-9-1 11:39
Kenneth 发表于 2014-8-31 08:02 * S0 r8 E4 ?; _" u% {
chaoswang,
0 T* g+ }4 ]4 u! d, W# A% l1 o5 C% D) U& M4 P, u/ R0 y! E
如果是因为起始值的问题，可以试试放弃模型中其他变量，只是把有问题的线性函数做回归，得出截 ...

' J) y. X) s/ R# B- }" H( d好的，谢谢您！我试试看！
作者: chaoswang    时间: 2014-9-11 14:02
哈哈，原来出问题的原因是我这两个中介变量的variance和intercept应该设为相等，设好之后问题就没有了，太开心了！
6 A$ I% b+ s+ ?& u9 H2 P$ ]谢谢各位！



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