search

dbspitz / migrateR

Shared tools for analyzing animal movement behavior in R
6 stars 0 forks source link

Error when constraining migration #12

Open malenacandino opened 3 years ago

malenacandino commented 3 years ago

Hi, I'm trying to use top.bhs.elev2<-topmvmt(bhs.elev2,mrho = 30 , mdelta = 100) to constrain migration to be at least a month long, and at least 100mt in altitudinal difference, and for some reason, when I get the output, one of my individuals with a rho=8.57 is still being classified as migratory. I've tried removing the mdelta and leaving just the mrho=30, and it still doesn't work. DO you know what the problem could be?

Thank you,

Malena

dbspitz commented 3 years ago

Hi Malena,

Thanks for bringing this to my attention! I wasn't aware of this issue, but it will soon be moot. I recently overhauled this function to increase its flexibility. The new version will be included in v1.1.0 (currently in review). In the meantime, hopefully you can handle this manually. The threshold arguments for 'topmvmt' just exclude models from consideration. Inelegant, but perhaps something like:

lowRho <- which(sapply(topmvmt(yourElevMod), function(z) coef(z)["rho"] < 30)) yourElevMod2 <- yourElevMod yourElevMod2[lowRho] <- lapply(yourElevMod[lowRho], function(z){ @. <- @*.**@*.***)%in%"migrant")]. # remove migrant model return(z) }) topmvmt(yourElevMod2)

Best,

Derek Spitz PhD Diaper Changes & Quantitative Ecology CityOfKitesAndCrows.com http://cityofkitesandcrows.com Folsom, CA (831) 737-3120

On Tue, Mar 9, 2021 at 2:04 PM malenacandino @.***> wrote:

Hi, I'm trying to use top.bhs.elev2<-topmvmt(bhs.elev2,mrho = 30 , mdelta = 100) to constrain migration to be at least a month long, and at least 100mt in altitudinal difference, and for some reason, when I get the output, one of my individuals with a rho=8.57 is still being classified as migratory. I've tried removing the mdelta and leaving just the mrho=30, and it still doesn't work. DO you know what the problem could be?

Thank you,

Malena

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/dbspitz/migrateR/issues/12, or unsubscribe https://github.com/notifications/unsubscribe-auth/ACTZKUSY4BVY44GHO563AMTTC2EMZANCNFSM4Y4SJNMQ .