Boa tarde Pessoal,

Tenho experimentos em blocos(Locais) e dentro deles gostaria de por algumas de repetições de alguns tratamentos para estimar a variância dentro do bloco, que desconfio ser grande. Por exemplo, tenho 5 Cultivares gostaria de repetir 3 delas dentro do bloco, tendo ao todo 8 parcelas, o ideal seria ter uma repetição de tudo, porem na pratica isso não é possível. 

Minhas duvidas são:

1) Estou definindo de maneira correta essas repetições dentro do modelo?

2) Como o método  de calculo dos dois modelos abaixo é igual para os efeitos fixos, por que o intercepto esta com valores diferentes?

3) Neste dois casos(m0 e m1, abaixo) o BLUP seria a média mais o efeito da cultivar? 

Obrigado pela ajuda

require(lme4)

### Simulando dados

dados <- expand.grid(Variedade=factor(c(letters[1:3],letters[1:5])),Local=factor(1:5))

dados$Rep <- factor(rep(c(2,2,2,1,1,1,1,1),t=5))

dados$Resp <- sort(rnorm(nrow(dados)))

dados

### Modelos

# Sem efeito dentro do local

m0 <- lmer(Resp~1+(1|Local)+(1|Variedade),dados)

summary(m0)

ranef(m0)

# Com efeito dentro do local

m1 <- lmer(Resp~1+(0+Rep|Local)+(1|Variedade),dados)

summary(m1)

ranef(m1)