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/ --