[R-br] Criação de uma função complexa com dplyr

Cesar Rabak cesar.rabak em gmail.com
Ter Jun 28 17:52:46 -03 2022 

A sua descrição está difícil de acompanhar porque é colocado entremeado nas
regras termos que necessitam melhor definição para entendê-las:


   - considerando "Class_1_3" para "Class_50"  ;
   - ClassType não mudar por nest durante o tempo   ;
   - AT entrar mas não sair do ClassType  ;
   - o último ClassType antes de AT  ;
   - último ClassType e o atual ClassType para nest  ;
   - Se atrasar qualquer ClassType em  ;
   - "Class_1_3" para "Class_50"  ;
   - valor de AT no tempo 0  ;
   - classe não presente  ;
   - último ClassType para nest  ;

Ademais, parece que sua formulação de velocidade e mesmo a saída desejável
indica, que existam velocidades diferentes para o mesmo nest e classe.

Por certo esses termos todos acima devem fazer todo o sentido para vc mas
para podermos ajudá-lo, essas definições devem ser quebradas nas operações
que elas implicam, aí naturalmente convergir-se-á para a formulação final
desejada.

Por exemplo, se ClassType implica em ordem, então necessita-se de listas
ordenadas dessas classe por nest?

Parece que a formulação para o cálculo de velocidade é mais simples quando
Se AT for >= 50 && Se o ninho for >= 1000  ->  a velocidade é c(0,diff(AT))/
c(0,diff(days)) por days para cada nest, Que se entendida corretamente sua
explicação,  significa que apenas as linhas que tenham classe "Class_0_1",
seriam usadas para esses cálculos.

Não entendo porque é catenado um zero no num e dem dos cálculos de
velocidade . . .

Supondo as condições para esse cálculo para classe e ninho, você consegue
escrever essa spec na linguagem R ou tem problemas com a sintaxe para
produzir o resultado?

HTH


On Tue, Jun 28, 2022 at 8:49 AM Alexandre Dos Santos por (R-br) <
r-br em listas.c3sl.ufpr.br> wrote:

> Prezados Membros do R-br,
>
> Meu dados de exemplo são:
>
> library(dplyr)
>
> my.ds <- read.csv("
> https://raw.githubusercontent.com/Leprechault/trash/main/test_ants.csv")
>
> my.ds
>
> # days nest       AT
>
> # 1     0    2  10.9200
>
> # 2    47    2  22.8600
>
> # 3    76    2  23.2400
>
> # 4     0    3   0.1400
>
> # 5    47    3   0.4800
>
> # 6    76    3   0.8300
>
> # 7   118    3   0.8300
>
> # 8   160    3   0.9400
>
> # 9   193    3   0.9400
>
> # 10  227    3   0.9400
>
> # 11  262    3   0.9400
>
> # 12  306    3   0.9400
>
> # 13  355    3  11.9300
>
> # 14  396    3  12.8100
>
> # 15  450    3  29.3600
>
> # 16    0    4   1.0000
>
> # 17   76    4   1.5600
>
> # 18  160    4   2.8800
>
> # 19  193    4   2.8800
>
> # 20  227    4   2.8800
>
> # 21  262    4   2.8800
>
> # 22  306    4   2.8800
>
> # 23  355    4  17.0500
>
> # 24    0    7  11.7100
>
> # 25   47    7  24.7900
>
> # 26   76    7  33.1200
>
> # 27    0   10   4.1400
>
> # 28   47   10   6.3000
>
> # 29   76   10  13.9400
>
> # 30  118   10  13.3000
>
> # 31  160   10  29.4200
>
> # 32  193   10  29.6000
>
> # 33  227   10  45.1500
>
> # 34  262   10  45.1500
>
> # 35  306   10  45.1500
>
> # 36  355   10  70.5300
>
> # 37  396   10 109.2600
>
> # 38  450   10 149.2000
>
> # 39    0 1066   0.0289
>
> # 40   29 1066   0.0870
>
> # 41   70 1066   0.1254
>
> # 42  112 1066   0.1254
>
> # 43  146 1066   0.1254
>
> # 44  180 1066   0.1254
>
> # 45  215 1066   0.1350
>
> # 46  259 1066   0.1350
>
> # 47    0 1067   0.0520
>
> # 48   29 1067   0.1254
>
> # 49   70 1067   0.1440
>
> # 50  112 1067   0.2160
>
> # 51  146 1067   0.2538
>
> # 52  180 1067   0.2538
>
> # 53  215 1067   0.9600
>
> # 54  259 1067   0.9600
>
> # 55  349 1067   0.9600
>
> # Class size
>
> my.ds$ClassType <- cut(my.ds$AT,breaks=c(-Inf,1,2.9,8.9,24.9,49.9,Inf),
> right=FALSE,
>
>                           labels=c("Class_0_1","Class_1_3","Class_3_9",
> "Class_9_25","Class_25_50","Class_50"))
>
>
> Eu gostaria de calular a variável velocity com as seguintes regras:
>
> Se nest for < 1000 - considerando "Class_1_3" para "Class_50"
>
> 1.1) Se AT for < 50
>
> 1.1.1) Se ClassType não mudar por nest durante o tempo ou AT entrar mas
> não sair do ClassType, então não calcule velocity (no meu caso nest 2 e 7
> ).,
>
> 1.1.2) Procurando o último ClassType antes de AT e calcula a `velocity`
> entre o último ClassType e o atual ClassType para nest
>
> 1.1.3) Se atrasar qualquer ClassType em "Class_1_3" para "Class_50" em uma
> classe maior que o valor de AT no tempo 0, a classe não presente é o valor
> de velocity do último ClassType para nest.
> 1.2) Se AT for >= 50
>
> 1.2.1) velocity é c(0,diff(AT))/c(0,diff(days)) por days para cada nest
> Se o ninho for >= 1000 - apenas considerando "Class_0_1"
>
> 2.1) a velocidade é c(0,diff(AT))/c(0,diff(days)) por days para cada nest
> também.
>
> Meu output desejável é:
>
> # nest  ClassType   velocity
>
> # 3 Classe_ 1_3 0,224285714
>
> # 3 Classe_ 3_9 0,224285714
>
> # 3 Classe_ 9_25 0,224285714
>
> # 3 Classe_ 9_25 0,021463415
>
> # 3 Classe_ 9_25 0,306481481
>
> # 4 Classe_ 1_3 0,0075
>
> # 4 Classe_ 1_3 0,015714286
>
> # 4 Classe_ 1_3 0
>
> # 4 Classe_ 1_3 0
>
> # 4 Classe_ 1_3 0
>
> # 4 Classe_ 1_3 0
>
> # 4 Classe_ 1_3 0,289183673
>
> # 4 Classe_ 3_9 0,289183673
>
> # 10 Classe_9_25 0,263448276
>
> # 10 Classe_9_25 -0,015238095
>
> # 10 Classe_9_25 0,383809524
>
> # 10 Class_25_50 0,383809524
>
> # 10 Classe_25_50 0,005454545
>
> # 10 Class_25_50 0,457352941
>
> # 10 Classe_25_50 0
>
> # 10 Classe_25_50 0
>
> # 10 Class_25_50 0,517959184
>
> # 10 Class_50 0,517959184
>
> # 10 Class_50 0,944634146
>
> # 10 Class_50 0,73962963
>
> # 1066 Classe_0_1 0,002003448
>
> # 1066 Classe_0_1 0,000936585
>
> # 1066 Classe_0_1 0
>
> # 1066 Classe_0_1 0
>
> # 1066 Classe_0_1 0
>
> # 1066 Classe_0_1 0,000274286
>
> # 1066 Classe_0_1 0
>
> # 1067 Classe_0_1 0,002531034
>
> # 1067 Classe_0_1 0,000453659
>
> # 1067 Classe_0_1 0,001714286
>
> # 1067 Classe_0_1 0,001111765
>
> # 1067 Classe_0_1 0
>
> # 1067 Classe_0_1 0.020177143
>
> # 1067 Classe_0_1 0
>
> # 1067 Classe_0_1 0
>
> Por favor, alguém poderia me ajudar com isso?
>
>
>
> --
>
> --
> Alexandre dos Santos
> Geotechnologies and Spatial Statistics applied to Forest Entomology
> Instituto Federal de Mato Grosso (IFMT) - Campus Caceres
> Caixa Postal 244 (PO Box)
> Avenida dos Ramires, s/n - Vila Real
> Caceres - MT - CEP 78201-380 (ZIP code)
> Phone: (+55) 65 99686-6970 / (+55) 65 3221-2674
> Lattes CV: http://lattes.cnpq.br/1360403201088680
> OrcID: orcid.org/0000-0001-8232-6722
> ResearchGate: www.researchgate.net/profile/Alexandre_Santos10
> Publons: https://publons.com/researcher/3085587/alexandre-dos-santos/
> --
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20220628/7cd50a50/attachment.htm>

Mais detalhes sobre a lista de discussão R-br