2

我正在尝试将 Milliken 的论文中的分析（https://sci-hub.tw/10.1016/s0169-7161(03)22007-1，第 8 节）从 SAS 代码复制到 R。我很困惑老实说。这是一个裂图重复测量设计，其中相关结构是复合对称结构。下面是数据和 SAS 代码及其结果。

数据

library(magrittr)
library(tidyr)
library(dplyr)
dta <- data.frame(
  tmp = c(rep(900, 3), rep(1000, 3), rep(1100, 3)),
  posit = rep(c("top", "mid", "bot"), 3),
  lot_1 = c(189, 211, 178, 213, 220, 197, 194, 212, 189),
  lot_2 = c(195, 206, 162, 199, 230, 198, 215, 208, 193),
  lot_3 = c(183, 210, 173, 189, 228, 202, 194, 201, 180),
  lot_4 = c(187, 223, 181, 183, 221, 168, 232, 215, 192),
  lot_5 = c(173, 191, 149, 202, 213, 151, 190, 198, 182)
)


dta <- dta %>% 
  tidyr::pivot_longer(., cols = c(lot_1, lot_2, lot_3, lot_4, lot_5), 
               names_to = "Lot") %>% 
  dplyr::mutate(Lot = as.factor(Lot), 
                tmp = as.factor(tmp), 
                lot_tmp = as.factor(paste0(Lot, "-", tmp)))

SAS 代码

proc mixed data = dta cl covtest ic; 
class Posit temp lot; 
model thick = temp Posit Posit*temp/ddfm = kr; random lot; 
repeated posit/type = cs subject = lot*temp r rcorr

SAS 输出
在此处输入图像描述

R代码尝试


## this works but isn't doing the same thing as above 
library(nlme)

m1 <- lme(
  value ~ temp + posit + temp:posit, 
  random = ~ 1 | lot ,
  correlation = corCompSymm(form=~1|lot),
  data = dta, method = "REML"
)

在这一点上，我被困在如何向posit因子中添加重复结构。

感谢您的帮助！

4

0 回答 0