楼主: tssong
683 0

[问答] 【悬赏R语言运行混合线性回归出现的结果问题】 [推广有奖]

  • 0关注
  • 0粉丝

初中生

38%

还不是VIP/贵宾

-

威望
0
论坛币
856 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
83 点
帖子
4
精华
0
在线时间
19 小时
注册时间
2020-3-12
最后登录
2021-10-28

100论坛币

请大神帮忙看下我的linear mixed-effects modeling 的结果 我有点看不明白. random effect 和fixed effect都出现了两种截距 。而且其中个变量的截距=0


我的实验涉及是在两个条件下阅读一组英文单词,看condition对于单词阅读速度的影响。是between-subjectdesign

Fixed effects: condition (2个条件) ; word(一组单词)

Random effects: word, condition

Covariate: length(单词长度)


  1. Formula: Reading Speed ~ condition+ Length + (1 | Word) + (1 | Participant) +  
  2. (1 | Length) + (condition |Word) + (1 + Participant |Length),data=reading_speed_data
复制代码
-----------------------------------------------------------------------------------------------------------------------------
R给出的结果:
REML criterion at convergence: 1246

Scaled residuals:
    Min      1Q  Median      3Q     Max --
-1.3318 -0.7162 -0.1443  0.4815  5.1412

Random effects:
Groups      Name            Variance  Std.Dev. Corr      
Word   (Intercept)     4.818e+02  21.9509            
             condition         7.214e+01   8.4934 -1.00      
Word.1  (Intercept)     0.000e+00   0.0000            
Length      (Intercept)     1.262e+04 112.3456            
             Participant66_2 1.879e-02   0.1371 -0.97      
             Participant67_3 1.638e-01   0.4047 -0.98  1.00
Length.1    (Intercept)     6.373e+03  79.8315            
Participant (Intercept)     0.000e+00   0.0000            
Residual                    8.056e+03  89.7526            
Number of obs: 108, groups:  word, 4; Length, 3; Participant, 3

Fixed effects:
            Estimate Std. Error       df t value Pr(>|t|)  
(Intercept) 281.7483   140.4396  72.6324   2.006   0.0486 *
condition     -6.1640    11.1522   3.6607  -0.553   0.6125  
Length9      -9.2571   195.8803   0.7980  -0.047   0.9713  
Length10      9.2929   195.8953   0.8434   0.047   0.9709  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
         (Intr) condition Lngth9
condition -0.176              
Length9  -0.696  0.005      
Length10 -0.696  0.007  0.498
optimizer (nloptwrap) convergence code: 0 (OK)
boundary (singular) fit: see ?isSingular

>

Q. 1 请问为什么结果random effect 和fixed effect都出现了两种截距 word 和 word.1; length和length1。而且其中word 和length 的截距都是=0. 请问是否和我的样本量很小有关呢(n=3个被试,初步运行一下计算所以测了3个人的数据)?还是我的model不对呢?


Q2. 请问这种情况下random error怎么计算呢?例如计算word 的random error,是word 和word1分开计算么?

Q3 .请问根据我的实验,我的 方程式是否是正确的呢? 非常感谢!!


a. 组间设计,两组人各自在两个不同条件下阅读同一批单词

b. 所有的单词均在两个条件内,所以有 (condition |Word) ,假设每个condition对于每个单词的影响不同。

c. 解释变量是阅读速度 reading speed 。

d. 单词的长度不同 因此将长度设为协变量。

e.单词的长度设置为了factor;word设置为了factor;condition为factor

Formula: Reading Speed ~ condition+ Length + (1 | Word) + (1 | Participant) +
(1 | Length) + (condition |Word) + (1 + Participant |Length),data=reading_speed_data


关键词:线性回归 R语言 fixed effect condition covariate 混合线性模型 R语言程序请教 混合线性模型截距
您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-24 17:37