Acredito que o jeito mais simples seja assim

dados <- expand.grid(trat=gl(3,4), tempo=1:5)
dados$y <- rpois(nrow(dados), lambda=10)
dados$z <- 1

dados2 <- dados[rep(1:nrow(dados), dados$y),]

xtabs(y~trat+tempo, data=dados)
xtabs(z~trat+tempo, data=dados2)

À disposição.
Walmes.