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@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.
