相关疑难解决方法(0)

如何制作一个很好的R可重复的例子

在与同事讨论性能,教学,发送错误报告或在邮件列表上搜索指导时,以及在Stack Overflow上,通常会询问可重复的示例并始终提供帮助.

您有什么建议创建优秀示例的提示?如何以文本格式粘贴中的数据结构?您还应该包含哪些其他信息?

在另外还有其他招数来使用dput(),dump()structure()?你什么时候应该包括library()require()声明?其中保留字应避免一个,此外c,df,data等?

怎样才能成为一位伟大的重复的例子?

r r-faq

2474
推荐指数
23
解决办法
28万
查看次数

从`lmermod`对象中提取结果

我想从一个lmermod对象中提取一些结果

require(lme4)
(fm1 <- lmer(Yield ~ 1|Batch, Dyestuff))
Run Code Online (Sandbox Code Playgroud)

这产生了

Linear mixed model fit by REML ['lmerMod']
Formula: Yield ~ 1 | Batch
   Data: Dyestuff
REML criterion at convergence: 319.6543
Random effects:
 Groups   Name        Std.Dev.
 Batch    (Intercept) 42.00   
 Residual             49.51   
Number of obs: 30, groups:  Batch, 6
Fixed Effects:
(Intercept)  
       1527 
Run Code Online (Sandbox Code Playgroud)

特别是,我想提取2个随机效应的标准偏差,即42.00和49.51.我想可能有一个内置的方法来做到这一点,但我找不到它,所以我认为我应该能够通过使用找到它str(fm1),这产生了这个:

Formal class 'lmerMod' [package "lme4"] with 13 slots
  ..@ resp   :Reference class 'lmerResp' [package "lme4"] with 9 fields
  .. ..$ Ptr    :<externalptr> 
  .. ..$ mu …
Run Code Online (Sandbox Code Playgroud)

r lme4

5
推荐指数
1
解决办法
2735
查看次数

标签 统计

r ×2

lme4 ×1

r-faq ×1