Como especificar efeitos aleatórios no lme?

8

Pesquisei este online por horas, mas nenhuma postagem online é o que estou procurando. Minha pergunta é muito fácil de implementar no procedimento misto SAS Proc, mas não sei como fazê-lo nos pacotes lme e / ou lmer. Suponha, eu tenho um modelo, , em que é fixo, mas e $y = \mu + \alpha + \beta +\alpha\beta + e$ $\alpha$ $\beta$ $\alpha\beta$ são aleatórios. Meu código R é

 f1 = lme(y ~ factor(a), data = mydata,
     random = list(factor(b) = ~ 1, factor(a):factor(b) = ~ 1))

Erro: inesperado =em:

 f1 = lme(y ~ factor(a), data = mydata,
          random = list(factor(a) =

Alguém poderia me dizer como especificar esses efeitos aleatórios no lme? Muito obrigado antecipadamente

r mixed-model cross-validation biostatistics

— Tu.2
fonte

É útil usar dputpara obter o código necessário para recriar seus dados. A partir do comentário que você deixou, o resultado é

structure(list(method = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,  1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("1", "2"), class = "factor"),  day = structure(c(1L, 1L, 2L, 2L, 3L, 3L, 4L, 4L, 1L, 1L,  2L, 2L, 3L, 3L, 4L, 4L), .Label = c("1", "2", "3", "4"), class = "factor"),  level = c(142.3, 144, 134.9, 146.3, 148.6, 156.5, 152, 151.4,  142.9, 147.4, 125.9, 127.6, 135.5, 138.9, 142.9, 142.3)), .Names = c("method",  "day", "level"), row.names = c(NA, -16L), class = "data.frame")

— Aaron deixou o Stack Overflow

12

Tente isso, é uma maneira padrão de fazer uma plotagem dividida. A notação /significa que o método está aninhado em dia.

lme(level~method, random=~1|day/method, data=d)

— Aaron deixou Stack Overflow
fonte

Olá Aaron, muito obrigado. Sua saída R é exatamente a mesma que a saída SAS e o manual. Mas por que precisamos usar "aninhado" na sintaxe R. Porque no livro, afirma claramente que

α β

$\alpha\beta$ é um termo de interação aleatória e eu também uso o termo de interação aleatória no SAS. Você poderia me dizer por que / como especificar um efeito de interação aleatória no lme, se for possível? muito obrigado Tu.2

— Tu.2

3

Sua pergunta não é sobre a sintaxe R, é sobre o que significa aninhamento. Aninhar B em A (com A / B) cria duas variáveis, A e a interação entre A e B, que é exatamente o que você descreve.

— Aaron deixou o Stack Overflow

Oi, esta é uma ótima explicação. Muito obrigado.

— Tu.2

2

Ajudaria muito se você fornecesse um data.frame. Agora não está claro o que é um fator de agrupamento. Eu julgo que é $\beta$ . Em seguida, na lmenotação, seu modelo deve ser escrito da seguinte maneira:

lme(y~a,random=~a|b, data=mydata)

— mpiktas
fonte

Olá mpiktas, obrigado pela sua resposta, mas o resultado R é muito diferente do procedimento misto proc SAS. A saída do SAS é a mesma com a resposta em um livro didático (Kuehk RO. 1999. Design de experimentos: princípios estatísticos de design e análise de pesquisa, 2ª edição). O conjunto de dados é o nível do dia do método 1 1 142,3 1 1 144,0 1 2 134,9 1 2 146,3 1 3 148,6 1 3 156,5 1 4 152,0 1 4 151,4 2 1 142,9 2 1 147,4 2 2 125,9 2 2 127,6 2 3 135,5 2 3 138,9 2 4 142,9 2 4 142,3

— Tu.2 17/05