Isso pode ser feito de muitas formas, via muitos pacotes e funções. Dentre tantas, uma possibilidade seria essa. Certamente não é a melhor mas serve.

str(GPD)

xtabs(~TRAT+ANIMAL, data=GPD)
xtabs(~DATA+ANIMAL, data=GPD)

require(reshape)
require(plyr)

db <- cast(data=GPD, ANIMAL~DATA, value="Pesos")

db$mGPD <- apply(db[,-1], MARGIN=1,
                 function(x){
                     mean(diff(na.omit(x)))
                 })

À disposição.

Walmes.

​

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

-- 
Fernando Antônio de Souza
Zootecnista,Dsc. Nutrição e alimentação animal (UFMG)
Lattes: http://lattes.cnpq.br/6519538815038307
e-mail:nandodesouza@gmail.com
(82)9618-6457 (tim) 
(82)8113-8781 (vivo)
(82)3313-8781 (net)
(31)9784-5273 (vivo-MG)