Fatal编程技术网
  • r/
  • R 复制交叉设计中的线性混合模型

R 复制交叉设计中的线性混合模型

r

R 复制交叉设计中的线性混合模型,r,lme4,crossover,R,Lme4,Crossover,我正在努力研究如何使用REML方法来适应复制交叉设计的模型。FDA建议的模型如上所述，有人能帮助将其编码为R编码吗？这是我的代码，我想知道它是对的还是错的 samplePK.lmer <- lmer(ykir2~1+Treatment:Sequence:Replication+ (1|Subject:Sequence:Treatment), data=samplePK, REML=TRUE) samplePK.lmer我认为公式应该是 response~trt+trt:seq:

我正在努力研究如何使用REML方法来适应复制交叉设计的模型。FDA建议的模型如上所述，有人能帮助将其编码为R编码吗？这是我的代码，我想知道它是对的还是错的

samplePK.lmer <- lmer(ykir2~1+Treatment:Sequence:Replication+
    (1|Subject:Sequence:Treatment), data=samplePK, REML=TRUE)


samplePK.lmer我认为公式应该是


response~trt+trt:seq:rep+（trt | subc:seq）


与您的规范的关键区别在于
（trt | sub:seq）
适用于类似于随机区组或随机斜率模型的模型，其中我们考虑了不同受试者/序列组合的治疗效果变化

在尝试将此模型与模拟数据进行拟合时，我遇到/注意到了一些问题：


    

  • 如果我们将其与“现代”混合模型方法相匹配，则在治疗效果（
    trt
    ）和
    trt:seq:rep
    术语之间会有一些参数混叠。在矩量法方法中，这并不重要（因为您从未明确估计参数），但它会导致对秩缺陷模型的抱怨（如果您知道自己在做什么，这是可以忽略的…）
    • 

  • 随机效应
    delta{ij}
    的平均值为
    {mu\u R，mu\u T}
    ，这似乎是错误的；这与固定效果是多余的
    mu_k
    • 


很明显，我可能对原始模型规范有错误或误解


我可能建议您尝试在
r-sig-mixed上提出后续问题-models@r-org
邮件列表，其中有广泛的读者群和广泛的专业知识（即，在混合模型主题上比在StackOverflow上更专业）
我认为公式应该是


response~trt+trt:seq:rep+（trt | subc:seq）


与您的规范的关键区别在于
（trt | sub:seq）
适用于类似于随机区组或随机斜率模型的模型，其中我们考虑了不同受试者/序列组合的治疗效果变化

在尝试将此模型与模拟数据进行拟合时，我遇到/注意到了一些问题：


    

  • 如果我们将其与“现代”混合模型方法相匹配，则在治疗效果（
    trt
    ）和
    trt:seq:rep
    术语之间会有一些参数混叠。在矩量法方法中，这并不重要（因为您从未明确估计参数），但它会导致对秩缺陷模型的抱怨（如果您知道自己在做什么，这是可以忽略的…）
    • 

  • 随机效应
    delta{ij}
    的平均值为
    {mu\u R，mu\u T}
    ，这似乎是错误的；这与固定效果是多余的
    mu_k
    • 


很明显，我可能对原始模型规范有错误或误解


我可能建议您尝试在
r-sig-mixed上提出后续问题-models@r-org
邮件列表，其中有广泛的读者群和广泛的专业知识（即，在混合模型主题上比在StackOverflow上更专业）
PS这是从这里来的吗？将此软件包与
MixedModels.jl
（Julia）和
lme4
（R）进行比较会很有趣。@BenBolker我已在R中安装了Julia软件包，但不幸的是出现了一些错误。您能帮我找到一种修复方法吗？库（“JuliaCall”）警告消息：包“JuliaCall”是在R version 3.6.3>julia\u setup（）julia version 1.5.3的位置C:\Users\User\AppData\Local\Programs\JULIA1~1.3\bin下构建的。inDL（x，as.logical（local），as.logical（now），…）中出错：无法加载共享对象“C:/Users/User/AppData/local/Programs/JULIA1~1.3/bin/libjulia.dll”：加载库故障：%1不是有效的Win32应用程序。juliacall_初始化（.julia$dll_文件）时出错：C:\Users\User\AppData\Local\Programs\JULIA1~1.3\bin\libjulia.dll-%1不是有效的Win32应用程序。对不起；Julia/R接口经常出现故障，我自己也很难跟上。PS，这是从这里开始的吗？将此软件包与
MixedModels.jl
（Julia）和
lme4
（R）进行比较会很有趣。@BenBolker我已在R中安装了Julia软件包，但不幸的是出现了一些错误。您能帮我找到一种修复方法吗？库（“JuliaCall”）警告消息：包“JuliaCall”是在R version 3.6.3>julia\u setup（）julia version 1.5.3的位置C:\Users\User\AppData\Local\Programs\JULIA1~1.3\bin下构建的。inDL（x，as.logical（local），as.logical（now），…）中出错：无法加载共享对象“C:/Users/User/AppData/local/Programs/JULIA1~1.3/bin/libjulia.dll”：加载库故障：%1不是有效的Win32应用程序。juliacall_初始化（.julia$dll_文件）时出错：C:\Users\User\AppData\Local\Programs\JULIA1~1.3\bin\libjulia.dll-%1不是有效的Win32应用程序。对不起；Julia/R界面经常被破坏，我自己也很难跟上。所以我改变了拟合模型，是的，它显示了对秩缺陷模型的抱怨。如果我收到了信息，可以继续使用模型吗？我想可以。如果能够跨软件包/平台（SAS、ReplicateBE.jl等）进行交叉检查，那就太好了……因此我已经更改了已安装的型号，是的，它显示了对排名不足型号的投诉。如果我收到了信息，可以继续使用模型吗？我想可以。跨软件包/平台（SAS、ReplicateBE.jl等）进行交叉检查以确保。。。