fitted values from mmer - fixed effects only about sommer HOT 2 CLOSED

data(DT_cornhybrids)
DT_cornhybrids <- DT

hybrid2 <- DT # extract cross data
A <- GT
K1 <- A[levels(hybrid2$GCA1), levels(hybrid2$GCA1)]; dim(K1)
K2 <- A[levels(hybrid2$GCA2), levels(hybrid2$GCA2)]; dim(K2)
S <- kronecker(K1, K2) ; dim(S)
rownames(S) <- colnames(S) <- levels(hybrid2$SCA)

ans <- mmer(Yield ~ Location,
             random = ~ vs(GCA1,Gu=K1) + vs(GCA2,Gu=K2) + vs(SCA,Gu=S),
             rcov=~units,
             data=hybrid2)
pp <- predict(ans)
head(pp$fitted) # has all fixed and random effects added