Pessoal, preciso de uma ajuda...
Meus dados são esses "a21", nele tenho a observação de 219 vacas (n_vaca), que estão numeradas de maneira aleatória. As vacas foram separadas em 3 grupos, onde em cada grupo aplicaram um tratamento diferente. Esse tratamento foi aplicado em forma de injeções (17 injeções) e entre uma injeção e outra, temos 14 dias de intervalo. A variável de interesse é "producao"
> head(a21)
n_vaca n_lactacao trat pri_multiparas n_injecao dia_ciclo dia_exp producao saude
8 1731 3 1 2 0 1 -13 45.1 healthy
9 1731 3 1 2 0 2 -12 47.8 healthy
10 1731 3 1 2 0 3 -11 45.5 healthy
11 1731 3 1 2 0 4 -10 45.2 healthy
12 1731 3 1 2 0 5 -9 47.4 enrollment
13 1731 3 1 2 0 6 -8 48.3 enrollment
# Aqui estou criando um subconjunto, onde tenho todas as informações da vaca 1731:
> (vaca1731<- subset(a21, a21$n_vaca==1731)); dim(vaca1731)
n_vaca n_lactacao trat pri_multiparas n_injecao dia_ciclo dia_exp producao saude
8 1731 3 1 2 0 1 -13 45.1 healthy
9 1731 3 1 2 0 2 -12 47.8 healthy
10 1731 3 1 2 0 3 -11 45.5 healthy
11 1731 3 1 2 0 4 -10 45.2 healthy
12 1731 3 1 2 0 5 -9 47.4 enrollment
13 1731 3 1 2 0 6 -8 48.3 enrollment
14 1731 3 1 2 0 7 -7 46.4 enrollment
15 1731 3 1 2 0 8 -6 47.9 enrollment
16 1731 3 1 2 0 9 -5 50.3 enrollment
17 1731 3 1 2 0 10 -4 50.8 enrollment
18 1731 3 1 2 0 11 -3 47.9 enrollment
19 1731 3 1 2 0 12 -2 46.2 heat
20 1731 3 1 2 0 13 -1 50.8 healthy
21 1731 3 1 2 0 14 0 45.8 healthy
22 1731 3 1 2 1 1 1 50.5 healthy
23 1731 3 1 2 1 2 2 49.0 healthy
24 1731 3 1 2 1 3 3 48.8 healthy
25 1731 3 1 2 1 4 4 53.9 healthy
26 1731 3 1 2 1 5 5 45.4 healthy
27 1731 3 1 2 1 6 6 38.6 healthy
28 1731 3 1 2 1 7 7 41.3 healthy
29 1731 3 1 2 1 8 8 50.2 healthy
30 1731 3 1 2 1 9 9 50.2 healthy
31 1731 3 1 2 1 10 10 50.4 healthy
32 1731 3 1 2 1 11 11 49.9 healthy
33 1731 3 1 2 1 12 12 50.1 healthy
34 1731 3 1 2 1 13 13 52.0 healthy
35 1731 3 1 2 1 14 14 47.8 healthy
36 1731 3 1 2 2 1 15 50.0 healthy
37 1731 3 1 2 2 2 16 49.7 healthy
38 1731 3 1 2 2 3 17 50.4 healthy
39 1731 3 1 2 2 4 18 48.9 healthy
40 1731 3 1 2 2 5 19 49.5 healthy
41 1731 3 1 2 2 6 20 48.8 healthy
42 1731 3 1 2 2 7 21 56.3 healthy
43 1731 3 1 2 2 8 22 52.1 healthy
44 1731 3 1 2 2 9 23 51.6 healthy
45 1731 3 1 2 2 10 24 50.5 healthy
46 1731 3 1 2 2 11 25 46.0 healthy
47 1731 3 1 2 2 12 26 53.7 healthy
48 1731 3 1 2 2 13 27 51.7 healthy
49 1731 3 1 2 2 14 28 49.1 healthy
50 1731 3 1 2 3 1 29 51.9 healthy
51 1731 3 1 2 3 2 30 51.4 healthy
...
256 1731 3 1 2 17 11 235 21.2 lamness
257 1731 3 1 2 17 12 236 19.4 lamness
258 1731 3 1 2 17 13 237 23.5 lamness
259 1731 3 1 2 17 14 238 19.9 lamness
[1] 251 9
# Para simplificar, considere que temos apenas 10 vacas:
#--------------------------------------------
> vacas <- c(1731, 2090, 2652, 2945, 2946, 3246, 3413, 3593, 3878, 3882) # é a numeração das vacas
# O que eu quero fazer com esses dados é o seguinte:
v1<- subset(a21, a21$n_vaca==vacas[1]); dim(v1)
v2<- subset(a21, a21$n_vaca==vacas[2]); dim(v2)
v3<- subset(a21, a21$n_vaca==vacas[3]); dim(v3)
v4<- subset(a21, a21$n_vaca==vacas[4]); dim(v4)
v5<- subset(a21, a21$n_vaca==vacas[5]); dim(v5)
v6<- subset(a21, a21$n_vaca==vacas[6]); dim(v6)
v7<- subset(a21, a21$n_vaca==vacas[7]); dim(v7)
v8<- subset(a21, a21$n_vaca==vacas[8]); dim(v8)
v9<- subset(a21, a21$n_vaca==vacas[9]); dim(v9)
v10<- subset(a21, a21$n_vaca==vacas[10]); dim(v10)
b<- c(
tapply(v1$producao, v1$dia_ciclo, mean),
tapply(v2$producao, v2$dia_ciclo, mean),
tapply(v3$producao, v3$dia_ciclo, mean),
tapply(v4$producao, v4$dia_ciclo, mean),
tapply(v5$producao, v5$dia_ciclo, mean),
tapply(v6$producao, v6$dia_ciclo, mean),
tapply(v7$producao, v7$dia_ciclo, mean),
tapply(v8$producao, v8$dia_ciclo, mean),
tapply(v9$producao, v9$dia_ciclo, mean),
tapply(v10$producao, v10$dia_ciclo, mean))
#--------------------------------------------
Ainda estou engatinhando em questão de programação no R, como poderia fazer um looping com essas linhas de comando, pois preciso fazer isso para as 219 vacas!
desde já, muito obrigada!
Simone
***---***---***---***---***---***---***---***---***---***---***---***---***
* Simone Daniela Sartorio
* Professora Adjunta I no CCA/UFSCar, Araras/SP.
* Doutora e Mestre em Estatística e Experimentação Agronômica - ESALQ/USP.
* Licenciada em Matemática - UNESP/Rio Claro.
***---***---***---***---***---***---***---***---***---***---***---***---***