2017-03-14 6 views
1

これはよく知られているか望ましい動作かどうか疑問に思っています...?"merModLmerTest"オブジェクトの更新は "lmerMod"クラスのオブジェクトを返します

再現例:

## import lmerTest package 
library(lmerTest) 

## an object of class merModLmerTest 
m <- lmer(Informed.liking ~ Gender+Information+Product +(1|Consumer), data=ham) 

class(m) 
class(update(m,REML=FALSE)) 
+0

質問はここにオフトピックですが。これはStack Overflowのトピックに載っているはずですので、待っていればそこに移行しようとします。 – gung

答えて

2

それは、少なくとも私には明らかではありません。あなたは

library(lmerTest) 
m <- lmer(Reaction~Days+(Days|Subject),sleepstudy) 
m2 <- update(m,REML=FALSE) 
m3 <- as(m2,"merModLmerTest") 
class(m3) ## "merModLmerTest" 

... as(.,"merModLmerTest")でそれを回避することができますか、このヘルパー関数を使用することができます:Rにコーディングする方法について

update2 <- function(x,...) { 
    z <- update(x,...) 
    return(as(z,"merModLmerTest")) 
} 
update2(m,REML=FALSE) 
+0

ありがとう!それは非常に便利です! –

関連する問題