Estou fazendo o modelo linear de efeitos mistos ajustado pelo REML no pacote nlme. E estes são os códigos que funcionam para mim:

# Linear mixed-effects model fit by REML (intercept and not slope)
x <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker)
summary(x)

# Linear mixed-effects model fit by REML (slope and no intercept)
x1 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~IV3-1|speaker)
summary(x1)

# Linear mixed-effects model fit by REML (slope and intercept)
x2 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~IV3|speaker)
summary(x2)

#nested random effect
x5 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker/item)
summary(x5)

O que realmente gostaria de fazer é ter um modelo com alto-falante e item como efeitos aleatórios separadamente. Eu tentei usar esta fórmula:

x4 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker + 1|item)

No entanto, esta fórmula me fornece a seguinte mensagem de aviso:

Warning message:
In Ops.factor(speaker, 1) : + not meaningful for factors

Você tem alguma idéia do que isso significa? E como encaixar o alto-falante e o item como efeitos aleatórios separadamente?