Compos_rebanho/script/EvolucaoVitruvio.org

#+Title:Evolução de rebanho leiteiro em PRV
#+AUTHOR: Adenor Wendling
#+DATE: 08-31-2023
#+OPTIONS: creator:nil timestamp:nil
#+PROPERTY: header-args:R :results output :session *R2* :cmdline :tangle yes
#+PROPERTY: header-args:R+ :exports both
#+SLUG: orgmode-excel
#+TAGS: tabelas, org-mode, composição ideal
#+SETUPFILE: https://fniessen.github.io/org-html-themes/org/theme-readtheorg.setup
#+HTML_HEAD: <style> #content{max-width:1800px;}</style>
#+HTML_HEAD: <style> p{max-width:800px;}</style>
#+HTML_HEAD: <style> li{max-width:800px;}</style>
#+HTML_HEAD: <style> pre.src{background: #f5f5f5;color:black;}</style>



** TODO

Continuar planilha
#+startup 2023-08-02


Vou estudar o emacs calc, tentando elaborar uma planilha para evolução do



Vou utilizar os manuais on line <https://orgmode.org/manual/Field-and-range-formulas.html>
e <https://orgmode.org/org.html#Built_002din-Table-Editor>

** Índices

IP  = intervalo entre partos (12 meses) (setq IP "12")

PP = primeiro parto (24 meses) (setq PP "24")

PL = Período de lactação (10 meses)(setq PL "10")

CL = capacidade de lotação (20UGM) (setq CL "20")
IndMor = Indice de mortalidade de bezerras
IRP = Índice de renovação do plantel - substituição de vacas velhas ou ruins


** Relação UGM e categoria animal

#+Caption: UGM de referência de acordo com a categoria animal
#+NAME: Tab1
| Cataegoria animal     |  UGM | abrev |
|-----------------------+------+-------|
| Vacas em lactação     |    1 | vaca  |
| Vacas secas           |    1 | vase  |
| Bezerras até 12 meses | 0.25 | Beze  |
| Novilhas até 18 meses | 0.75 | No18  |
| Novilhas até 24 meses | 0.75 | No24  |
| Outros                |    1 | Outr  |
|-----------------------+------+-------|

Os dados acima são padronizados, sendo que, por padrão, 1 UGM equivale a 500 kg de peso vivo. Desta forma, pode ser necessário ajustar as UGMs e acordo com o peso dos animais de cada fazenda.

** Composição do Rebanho

A composição do rebanho de cada estabelecimento varia conforme os objetivos, índices e manejo adotado. De modo geral busca-se uma composição com maior participação de categorias produtivas, ou seja, em caso de atividade leiteira, busca-se maior número de animais na categoria *vacas em lactação*.
Os dados da tabela 2 refletem a composição de um rebanho com finalidade de produção de leite, considerando um manejo reprodutivo e de alimentação excelentes.

#+Caption: Composição ideal em porcentagem , número de animais e UGM, com base em 10 vacas em lactação.
:PROPERTIES:
:TABLE_EXPORT_FILE: ~/Documentos/GITs/Compos_rebanho/data/composicao.csv
:TABLE_EXPORT_FORMAT: orgtbl-to-csv
:END:
#+NAME: C_ideal
|   | Categoria de animais  | Cabeças (N) | Cabeças (%) | UGM (N) | UGM (%) |
|---+-----------------------+-------------+-------------+---------+---------|
| # | Vacas em lactação     |          20 |        60.6 |      20 |    75.0 |
| # | Vacas secas           |           2 |         6.7 |       2 |     7.8 |
| # | Bezerras ate 12 meses |           4 |        13.3 |      1. |     3.9 |
| # | Novilhas até 18 meses |           2 |         6.7 |      1. |     3.9 |
| # | Novilhas até 24 meses |           2 |         6.7 |     1.5 |     5.9 |
| # | Outros animais        |             |         0.0 |       0 |     0.0 |
| # | Total                 |          30 |       100.0 |    25.5 |   100.0 |
| ^ |                       |        Tcab |             |    Tugm |         |
|---+-----------------------+-------------+-------------+---------+---------|
#+TBLFM: $4=$3*100/$Tcab;%.1f::@2$5=@2$3*remote(Tab1,@2$2)::@8$3=@2$3+@3$3+@4$3+@5$3+@6$3+@7$3::@3$5=@3$3*remote(Tab1, @3$2)::@4$5=@4$3*remote(Tab1, @4$2)::@5$5=@5$3*remote(Tab1, @5$2)::@6$5=@6$3*remote(Tab1, @6$2)::@7$5=@7$3*remote(Tab1, @7$2)::@8$5=vsum(@2..@7)::$6=$5*100/$Tugm;%.1f


A composição do rebanho é dinâmico, com espectativa de crescimento ao longo do tempo, até alcançar o limite de suporte das pastagens
ou de outro limitante definido pela infraestrutura ou de trabalho da fazenda. A espectativa de crescimento se deve à previsão de 1 parto anual
por vaca, que gera um incremento anual de 50% de fêmeas em relação ao número de vacas, e do descarte de 25% das vacas.
Quando for atingido o limite de vacas na fazenda, pode-se vender o excedente de fêmeas em idade de terneiras, novilhas prenhas, ou vacas,
de acordo com a conveniência da fazenda. Entretanto, o limite é medido em UGM, o que também vai limitar a opção pela venda de animais maiores.
Com base nos índices definidos em conjunto com a fazenda, estimamos que a evolução do rebanho de apresente conforme a tabela 3, ressalvando que
esses valores não são absolutos, pois podem ocorrer imprevistos que alterem a evolução. Da mesma forma, havendo alteração nos índices esperados,
haverá impacto na evolução do rebanho, o que não caracteriza necessariamente erro de previsão ou de manejo.



**  Tabela inicial, com o dimensionamento atual e previsão de entradas e saídas de animais

Os dados desta tabela deve ser alterados conforme necessário para alcançar a evolução do rebanho desejada.
As tabelas de composição de rebanho serão calculadas automáticamente, com base nos dados anotados na tabela <Tab0>.
#+Caption: Entradas e saídas de animais do rebanho durante os primeiros 6 anos (até estabilizar)
#+NAME:Tab0
|   | Constante              | ano 1 | ano 2 | ano 3 | ano 4 | ano 5 | ano .. | observação   |
|---+------------------------+-------+-------+-------+-------+-------+--------+--------------|
|   | Partos/vaca/ano        | 1     | 1     | 1     | 1     | 1     | 1      |              |
| ^ |                        | PV1   | PV2   | PV3   | PV4   | PV5   | PV6    |              |
|---+------------------------+-------+-------+-------+-------+-------+--------+--------------|
|   | Vacas Existentes       | 00    | -     | -     | -     | -     | -      | - a calcular |
| ^ |                        | VE1   |       |       |       |       |        |              |
|   | Novilhas exist.        | 00    | -     | -     | -     | -     | -      | - a calcular |
| ^ |                        | NE1   |       |       |       |       |        |              |
|   | Bezerras  existentes   | 0     | -     | -     | -     | -     | -      | - a calcular |
| ^ |                        | TE1   |       |       |       |       |        |              |
|---+------------------------+-------+-------+-------+-------+-------+--------+--------------|
|   | Novilhas >- vacas      | 0     | 0     | 0     | 0     | 0     | 0      | NÃO ALTERAR  |
| ^ |                        | TN1   | TN2   | TN3   | TN4   | TN5   | TN6    |              |
|   | % Bezerras -> novilhas | 0.4   | 0.4   | 0.75  | 0.7   | 0.7   | 0.7    | % S_Nascidos |
| ^ |                        | BA1   | BA2   | BA3   | BA4   | BA5   | BA6    |              |
|---+------------------------+-------+-------+-------+-------+-------+--------+--------------|
|   | Vacas a vender         | 0     | 0     | 0     | 2     | 5     | 4      |              |
| ^ |                        | VV1   | VV2   | VV3   | VV4   | VV5   | VV6    |              |
|   | Novilhas a vender      | 0     | 0     | 0     | 0     | 0     | 0      |              |
| ^ |                        | NV1   | NV2   | NV3   | NV4   | NV5   | NV6    |              |
|   | Bezerras a vender      | 0     | 4     | 10    | 11    | 11    | 11     |              |
| ^ |                        | BV1   | BV2   | BV3   | BV4   | BV5   | BV6    |              |
|---+------------------------+-------+-------+-------+-------+-------+--------+--------------|
|   | Vacas a comprar        | 0     | 0     | 0     | 0     | 0     | 0      |              |
| ^ |                        | VC1   | VC2   | VC3   | VC4   | VC5   | VC6    |              |
|   | Novilhas a comprar     | 20    | 0     | 0     | 0     | 0     | 0      |              |
| ^ |                        | NC1   | NC2   | NC3   | NC4   | NC5   | NC6    |              |
|   | Bezerras a compras     | 0     | 0     | 0     | 0     | 0     | 0      |              |
| ^ |                        | BC1   | BC2   | BC3   | BC4   | BC5   | BC6    |              |
|---+------------------------+-------+-------+-------+-------+-------+--------+--------------|

Com a evolução estimada através dos indicadores:
. 1 parto por vaca por ano;
. Primeiro parto aos 18 meses
. Descarte de 25% das vacas a partir do 3º ano
. Capacidade de suporte para 40 vacas
. Aproveitamento de 40% dos nascidos para repsição
Temos o seguinte:
>> Iniciamos com 20 vacas
>> Compramos 5 novilhas cobertas
>> No segundo ano serão compradas 2 novilhas


** Composição anual do rebanho da fazenda
*** Ano 1
#+Caption: Composição do rebanho na fazenda VITRUVIO (ANO 1)
#+NAME: Tab3_1
|   | Categoria       | ExttAnt | TIR | Compra | Venda | Saldo | TIC | ExttFin |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
| # | VACAS    >24    | 0       | 0   | 0      | 0     | 0     | 0   | 0       |
| ^ |                 | EAV     | TRV | COV    | VEV   | SAV   | TCV | EFV     |
| # | NOVILHAS 12<>24 | 0       | 0   | 20     | 0     | 20    | 20  | 0       |
| ^ |                 | EAN     | TRN | CON    | VEN   | SAN   | TCN | EFN     |
| #  | BEZERRAS  <12   | 0       | 0   | 0      | 0     | 0     | 0   | 0       |
| ^ |                 | EAB     | FRB | COB    | VEB   | SAB   | TCB | EFB     |
| # | Outros          |         | 0   | 0      | 0     |       | 0   | 0       |
| ^ |                 | EAO     | FRO | COO    | VEO   | SAO   | TCO | EFO     |
#+TBLFM: $EAV=remote(Tab0, $VE1)::$EFV=$SAV-$TCV::$SAV=$EAV+$TRV+$COV-$VEV::$TCN=($EAN+$TRN+$CON-$VEV)::$VEV=remote(Tab0, $VV1):: $EAN=remote(Tab0, $NE1)::$TRN=remote(Tab0, $TN1)::$CON=remote(Tab0, $NC1)::$VEN=remote(Tab0, $NV1)::$SAN=$EAN+$TRN+$CON-$VEN::$EFN=$SAN-$TCN::$EAB=remote(Tab0, $TE1)::$FRB=$EFV*remote(Tab0, $BA1);%.0f::$SAB=$EAB+$FRB+$COB-$VEB::$TCB=$SAB::$EFB=$SAB-$TCB


Na tabela 4 consta o número de animais para cada categoria, com as ocorrencias do ano. No caso das vacas, não há alteração, visto que não há entradas, nem saídas, e o número inicial é o mesmo do final. No caso das novilhas, há movimentação, com a transferência interna de 5 (que passam de bezerras para novilhas), compra de 2, e transferência cedendo 5 para serem vacas no próximo ano. Assim, o ano inicia com 1 novilha, e termina com 3.


*** Ano 2

#+Caption: Composição do rebanho na fazenda VITRUVIO (ANO 2)
#+NAME: Tab3_2
|   | CATEGORIA       | EXTTANT | TIR | COMPRA | VENDA | SALDO | TIC | EXTTFIN |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
| # | VACAS  >24      | 0       | 20  | 0      | 0     | 20    | 0   | 20      |
| ^ |                 | EAV     | TRV | COV    | VEV   | SAV   | TCV | EFV     |
| # | Novilhas 12<>24 | 0       | 0   | 0      | 0     | 0     | 0   | 0       |
| ^ |                 | EAN     | TRN | CON    | VEN   | SAN   | TCN | EFN     |
| # | Bezerras   <12  | 0       | 8   | 0      | 4     | 4     | 4   | 0       |
| ^ |                 | EAB     | FRB | COB    | VEB   | SAB   | TCB | EFB     |
| # | outros          |         |     |        |       |       |     |         |
| ^ |                 | EAO     | FRO | COO    | VEO   | SAO   | TCO | EFO     |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
#+TBLFM: $EAV=remote(Tab3_1, $EFV)::$TRV=remote(Tab3_1, $TCN)::$COV=remote(Tab0,$VC2)::$VEV=remote(Tab0, $VV2)::$SAV=$EAV+$TRV+$COV-$VEV::$EFV=$SAV::$EAN=remote(Tab3_1, $EFN)::$TRN=remote(Tab3_1, $TCB)::$CON=remote(Tab0,$NC2)::$VEN=remote(Tab0, $NV2)::$SAN=$EAN+$TRN+$CON-$VEN::$TCN=$SAN::$EFN=$SAN-$TCN::$EAB=remote(Tab3_1, $EFB)::$FRB=$EFV*remote(Tab0, $BA2);%.0f::$COB=remote(Tab0,$BC2)::$VEB=remote(Tab0,$BV2)::$SAB=$EAB+$FRB+$COB-$VEB::$TCB=$SAB::$EFB=$SAB-$TCB

 Comentários sobre tabela 5
*** Ano 3

Na tabela 3.3 podemos observar...
#+Caption: Composição do rebanho na fazenda VITRUVIO (ANO 3)
#+NAME: Tab3_3
|   | CATEGORIA       | EXTTANT | TIR | COMPRA | VENDA | SALDO | TIC | EXTTFIN |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
| # | VACAS  >24      | 20      | 0   | 0      | 0     | 20    | 0   | 20      |
| ^ |                 | EAV     | TRV | COV    | VEV   | SAV   | TCV | EFV     |
| # | Novilhas 12<>24 | 0       | 4   | 0      | 0     | 4     | 4   | 0       |
| ^ |                 | EAN     | TRN | CON    | VEN   | SAN   | TCN | EFN     |
| # | Bezerras   <12  | 0       | 15  | 0      | 10    | 5     | 5   | 0       |
| ^ |                 | EAB     | FRB | COB    | VEB   | SAB   | TCB | EFB     |
| # | outros          |         |     |        |       |       |     |         |
| ^ |                 | EAO     | FRO | COO    | VEO   | SAO   | TCO | EFO     |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
#+TBLFM: $EAV=remote(Tab3_2, $EFV)::$TRV=remote(Tab3_2, $TCN)::$COV=remote(Tab0,$VC3)::$VEV=remote(Tab0, $VV3)::$SAV=$EAV+$TRV+$COV-$VEV::$EFV=$SAV::$EAN=remote(Tab3_2, $EFN)::$TRN=remote(Tab3_2, $TCB)::$CON=remote(Tab0,$NC3)::$VEN=remote(Tab0, $NV3)::$SAN=$EAN+$TRN+$CON-$VEN::$TCN=$SAN::$EFN=$SAN-$TCN::$EAB=remote(Tab3_2, $EFB)::$FRB=$EFV*remote(Tab0, $BA3);%.0f::$COB=remote(Tab0,$BC3)::$VEB=remote(Tab0,$BV3)::$SAB=$EAB+$FRB+$COB-$VEB::$TCB=$SAB::$EFB=$SAB-$TCB
Comentários sobre a tabela 6

*** Ano 4

#+Caption: Composição do rebanho na fazenda VITRUVIO (ANO 4)
#+NAME: Tab3_4
|   | CATEGORIA       | EXTTANT | TIR | COMPRA | VENDA | SALDO | TIC | EXTTFIN |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
| # | VACAS  >24      | 20      | 4   | 0      | 2     | 22    | 0   | 22      |
| ^ |                 | EAV     | TRV | COV    | VEV   | SAV   | TCV | EFV     |
| # | Novilhas 12<>24 | 0       | 5   | 0      | 0     | 5     | 5   | 0       |
| ^ |                 | EAN     | TRN | CON    | VEN   | SAN   | TCN | EFN     |
| # | Bezerras   <12  | 0       | 15  | 0      | 11    | 4     | 4   | 0       |
| ^ |                 | EAB     | FRB | COB    | VEB   | SAB   | TCB | EFB     |
| # | outros          | 1       |     |        |       |       |     |         |
| ^ |                 | EAO     | FRO | COO    | VEO   | SAO   | TCO | EFO     |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
#+TBLFM: $EAV=remote(Tab3_3, $EFV)::$TRV=remote(Tab3_3, $TCN)::$COV=remote(Tab0,$VC4)::$VEV=remote(Tab0, $VV4)::$SAV=$EAV+$TRV+$COV-$VEV::$EFV=$SAV::$EAN=remote(Tab3_3, $EFN)::$TRN=remote(Tab3_3, $TCB)::$CON=remote(Tab0,$NC4)::$VEN=remote(Tab0, $NV4)::$SAN=$EAN+$TRN+$CON-$VEN::$TCN=$SAN::$EFN=$SAN-$TCN::$EAB=remote(Tab3_3, $EFB)::$FRB=$EFV*remote(Tab0, $BA4);%.0f::$COB=remote(Tab0,$BC4)::$VEB=remote(Tab0,$BV4)::$SAB=$EAB+$FRB+$COB-$VEB::$TCB=$SAB::$EFB=$SAB-$TCB

 Comentários da tabela 7
*** Ano 5

#+NAME: Tab3_5
#+Caption: Composição do rebanho na fazenda VITRUVIO (ANO 5)
|   | CATEGORIA       | EXTTANT | TIR | COMPRA | VENDA | SALDO | TIC | EXTTFIN |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
| # | VACAS  >24      | 22      | 5   | 0      | 5     | 22    | 0   | 22      |
| ^ |                 | EAV     | TRV | COV    | VEV   | SAV   | TCV | EFV     |
| # | Novilhas 12<>24 | 0       | 4   | 0      | 0     | 4     | 4   | 0       |
| ^ |                 | EAN     | TRN | CON    | VEN   | SAN   | TCN | EFN     |
| # | Bezerras   <12  | 0       | 15  | 0      | 11    | 4     | 4   | 0       |
| ^ |                 | EAB     | FRB | COB    | VEB   | SAB   | TCB | EFB     |
| # | outros          | 1       |     |        |       | 1     |     | 1       |
| ^ |                 | EAO     | FRO | COO    | VEO   | SAO   | TCO | EFO     |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
|   |                 |         |     |        |       |       |     |         |
#+TBLFM: $EAV=remote(Tab3_4, $EFV)::$TRV=remote(Tab3_4, $TCN)::$COV=remote(Tab0,$VC5)::$VEV=remote(Tab0, $VV5)::$SAV=$EAV+$TRV+$COV-$VEV::$EFV=$SAV::$EAN=remote(Tab3_4, $EFN)::$TRN=remote(Tab3_4, $TCB)::$CON=remote(Tab0,$NC5)::$VEN=remote(Tab0, $NV5)::$SAN=$EAN+$TRN+$CON-$VEN::$TCN=$SAN::$EFN=$SAN-$TCN::$EAB=remote(Tab3_4, $EFB)::$FRB=$EFV*remote(Tab0, $BA5);%.0f::$COB=remote(Tab0,$BC5)::$VEB=remote(Tab0,$BV4)::$SAB=$EAB+$FRB+$COB-$VEB::$TCB=$SAB::$EFB=$SAB-$TCB

Comentários da tabela 8

*** Ano 6
e, finalmente, na tabela 9, a composição do rebanho no sexto ano após a implementeação e adoção do PRV na propriedade VITRUVIO.


#+NAME:Tab3_6
#+caption: Composição do rebanho na fazenda VITRUVIO (ANO 6)
|   | CATEGORIA       | EXTTANT | TIR | COMPRA | VENDA | SALDO | TIC | EXTTFIN |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
| # | VACAS  >24      | 22      | 4   | 0      | 4     | 22    | 0   | 22      |
| ^ |                 | EAV     | TRV | COV    | VEV   | SAV   | TCV | EFV     |
| # | Novilhas 12<>24 | 0       | 4   | 0      | 0     | 4     | 4   | 0       |
| ^ |                 | EAN     | TRN | CON    | VEN   | SAN   | TCN | EFN     |
| # | Bezerras   <12  | 0       | 15  | 0      | 11    | 4     | 4   | 0       |
| ^ |                 | EAB     | FRB | COB    | VEB   | SAB   | TCB | EFB     |
| # | outros          | 1       |     |        |       |       |     |         |
| ^ |                 | EAO     | FRO | COO    | VEO   | SAO   | TCO | EFO     |
|---+-----------------+---------+-----+--------+-------+-------+-----+---------|
#+TBLFM: $EAV=remote(Tab3_5, $EFV)::$TRV=remote(Tab3_5, $TCN)::$COV=remote(Tab0,$VC6)::$VEV=remote(Tab0, $VV6)::$SAV=$EAV+$TRV+$COV-$VEV::$EFV=$SAV::$EAN=remote(Tab3_5, $EFN)::$TRN=remote(Tab3_5, $TCB)::$CON=remote(Tab0,$NC6)::$VEN=remote(Tab0, $NV6)::$SAN=$EAN+$TRN+$CON-$VEN::$TCN=$SAN::$EFN=$SAN-$TCN::$EAB=remote(Tab3_5, $EFB)::$FRB=$EFV*remote(Tab0, $BA6);%.0f::$COB=remote(Tab0,$BC6)::$VEB=remote(Tab0,$BV6)::$SAB=$EAB+$FRB+$COB-$VEB::$TCB=$SAB::$EFB=$SAB-$TCB



** Resumo da evolução

Esta tabela é a que efetivamente interessa para ser apresentada ao cliente, ou para análise técnica, pois apresenta o resumo (resultado da evolução) proposto para a fazenda.

Tabela com os dados resumidos referente aos animais existentes anualmente
Os dados da tabela abaixo resultam da estimativa de evolução de rebanho, de acordo com os dados da tabela de dados.

*** Entradas e saídas
#+NAME:TABELA5
#+CAPTION:  Resumo das entradas e saídas de animais nos primeiros anos do projeto
:PROPERTIES:
:TABLE_EXPORT_FILE: ~/Documentos/GITs/Compos_rebanho/data/ResumoEntraSaiVitruvio.csv
:TABLE_EXPORT_FORMAT: orgtbl-to-csv
:END:
|   | Categoria          | 2004 | 2005 | 2006 | 2007 | 2008 | 2009 |
|---+--------------------+------+------+------+------+------+------|
| # | Vacas entradas     |    0 |    0 |    0 |    0 |    0 |    0 |
| # | Novilhas entradas  |   20 |    0 |    0 |    0 |    0 |    0 |
| # | Terneiras entradas |    0 |    0 |    0 |    0 |    0 |    0 |
|---+--------------------+------+------+------+------+------+------|
| # | Total entradas     |   20 |    0 |    0 |    0 |    0 |    0 |
|---+--------------------+------+------+------+------+------+------|
| # | Vacas saídas       |    0 |    0 |    0 |    2 |    5 |    4 |
| # | Novilhas saídas    |    0 |    0 |    0 |    0 |    0 |    0 |
| # | Terneiras saídas   |    0 |    4 |   10 |   11 |   11 |   11 |
|---+--------------------+------+------+------+------+------+------|
| # | Total saídas       |    0 |    4 |   10 |   13 |   16 |   15 |
| # |                    |      |      |      |      |      |      |
|---+--------------------+------+------+------+------+------+------|
#+TBLFM:@2$3=remote(Tab3_1, $COV)::@2$4=remote(Tab3_2, $COV)::@2$5=remote(Tab3_3, $COV)::@2$6=remote(Tab3_4, $COV)::@2$7=remote(Tab3_5, $COV)::@2$8=remote(Tab3_6, $COV)::@3$3=remote(Tab3_1, $CON)::@3$4=remote(Tab3_2, $CON)::@3$5=remote(Tab3_3, $CON)::@3$6=remote(Tab3_4, $CON)::@3$7=remote(Tab3_5, $CON)::@3$8=remote(Tab3_6, $CON)::@4$3=remote(Tab3_1, $COB)::@4$4=remote(Tab3_2, $COB)::@4$5=remote(Tab3_3, $COB)::@4$6=remote(Tab3_4, $COB)::@4$7=remote(Tab3_5, $COB)::@4$8=remote(Tab3_6, $COB)::@5$3=@2$3+@3$3+@4$3::@5$4=@2$4+@3$4+@4$4::@5$5=@2$5+@3$5+@4$5::@5$6=@2$6+@3$6+@4$6::@5$7=@2$7+@3$7+@4$7::@5$8=@2$8+@3$8+@4$8::@6$3=remote(Tab3_1, $VEV)::@6$4=remote(Tab3_2, $VEV)::@6$5=remote(Tab3_3, $VEV)::@6$6=remote(Tab3_4, $VEV)::@6$7=remote(Tab3_5, $VEV)::@6$8=remote(Tab3_6, $VEV)::@7$3=remote(Tab3_1, $VEN)::@7$4=remote(Tab3_2, $VEN)::@7$5=remote(Tab3_3, $VEN)::@7$6=remote(Tab3_4, $VEN)::@7$7=remote(Tab3_5, $VEN)::@7$8=remote(Tab3_6, $VEN)::@8$3=remote(Tab3_1, $VEB)::@8$4=remote(Tab3_2, $VEB)::@8$5=remote(Tab3_3, $VEB)::@8$6=remote(Tab3_4, $VEB)::@8$7=remote(Tab3_5, $VEB)::@8$8=remote(Tab3_6, $VEB)::@9$3=@6$3+@7$3+@8$3::@9$4=@6$4+@7$4+@8$4::@9$5=@6$5+@7$5+@8$5::@9$6=@6$6+@7$6+@8$6::@9$7=@6$7+@7$3+@8$7::@9$8=@6$8+@7$8+@8$8

*** Animais existentes
#+NAME:TABELA4
#+CAPTION:  Resumo da evolução do rebanho (Animais existentes anualmente) nos primeiros anos do projeto
:PROPERTIES:
:TABLE_EXPORT_FILE: ~/Documentos/GITs/Compos_rebanho/data/ResumoEvolucapVitruvio.csv
:TABLE_EXPORT_FORMAT: orgtbl-to-csv
:END:
|   | Categoria | 2004 | 2005 |  2006 |  2007 | 2008 | 2009 |
|---+-----------+------+------+-------+-------+------+------|
| # | Vacas     |    0 |   20 |    20 |    22 |   22 |   22 |
| # | Novilhas  |   20 |    0 |     4 |     5 |    4 |    4 |
| # | Terneiras |    0 |    4 |     5 |     4 |    4 |    4 |
|---+-----------+------+------+-------+-------+------+------|
| # | T_cab     |   20 |   24 |    29 |    31 |   30 |   30 |
| # | T_UGMs    |  15. |  21. | 24.25 | 26.75 |  26. |  26. |
| # |           |      |      |       |       |      |      |
#+TBLFM:@2$3=remote(Tab3_1, $SAV)::@2$4=remote(Tab3_2, $SAV)::@2$5=remote(Tab3_3, $SAV)::@2$6=remote(Tab3_4, $SAV)::@2$7=remote(Tab3_5, $SAV)::@2$8=remote(Tab3_6, $SAV)::@3$3=remote(Tab3_1, $SAN)::@3$4=remote(Tab3_2, $SAN)::@3$5=remote(Tab3_3, $SAN)::@3$6=remote(Tab3_4, $SAN)::@3$7=remote(Tab3_5, $SAN)::@3$8=remote(Tab3_6, $SAN)::@4$3=remote(Tab3_1, $SAB)::@4$4=remote(Tab3_2, $SAB)::@4$5=remote(Tab3_3, $SAB)::@4$6=remote(Tab3_4, $SAB)::@4$7=remote(Tab3_5, $SAB)::@4$8=remote(Tab3_6, $SAB)::@5$3=@2$3+@3$3+@4$3::@5$4=@2$4+@3$4+@4$4::@5$5=@2$5+@3$5+@4$5::@5$6=@2$6+@3$6+@4$6:: @5$7=@2$7+@3$7+@4$7::@5$8=@2$8+@3$8+@4$8::@6$3=(@2$3*remote(Tab1,@2$2)+@3$3*remote(tab1,@6$2)+@4$3*remote(tab1,@4$2))::@6$4=(@2$4*remote(Tab1,@2$2)+@3$4*remote(tab1,@6$2)+@4$4*remote(tab1,@4$2))::@6$5=(@2$5*remote(Tab1,@2$2)+@3$5*remote(tab1,@6$2)+@4$5*remote(tab1,@4$2))::@6$6=(@2$6*remote(Tab1,@2$2)+@3$6*remote(tab1,@6$2)+@4$6*remote(tab1,@4$2))::@6$7=(@2$7*remote(Tab1,@2$2)+@3$7*remote(tab1,@6$2)+@4$7*remote(tab1,@4$2))::@6$8=(@2$8*remote(Tab1,@2$2)+@3$8*remote(tab1,@6$2)+@4$8*remote(tab1,@4$2))

*** TODO Resolver exportar ou criar data.tabela [%]

- Para continuar, é necessário exportar a Tabela acima. Faço pelo menu Table - Export to file. Nome do Arquivo a exportar = Tabela$.csv.>>


** Gráfico resumo
#+NAME: Ler data.table Resumo
#+begin_src R :exports results
  library(data.table)
  library(tidyr)
  library(dplyr)

  REV <- fread("../data/ResumoEvolucaoVitruvio.csv",  header=TRUE)

  REv <- REV %>%
      pivot_longer(cols = c("2004", "2005", "2006", "2007", "2008", "2009"),
                   names_to = "Ano",
                   values_to = "Quantidade",
                   values_drop_na = TRUE )

  REv$Ano <- as.numeric(REv$Ano)

  str(REv)
#+end_src

#+NAME: grafico resumo
#+begin_src R :grafResum.png :results graphics file

  library (ggplot2)

  ResEvolVit <- ggplot(REv, aes(Ano, Quantidade, colour=Categoria)) +
      geom_line()

  ggsave("ResEvolVit.png", ResEvolVit, path="../data/", width=7, height=3,
         units="in")

#+end_src