[R-sig-eco] adehabitatLT merge different bursts for same ID
g.cerritelli
giulia.cerritelli at biologia.unipi.it
Thu Feb 2 16:00:20 CET 2017
Dear list,
I'm trying to simulate the migratory movement of an animal, following
the work of Bunnefeld et al 2011.
I started from a dataset of dates and used the function simm.mou and
simm.bb in adehabitatLT to simulate five distinct movement types: the
stay in a wintering area, the migration to reproductive areas, the stay
in this sites, the return migration and finally the stay in wintering
areas.
Now I would like to merge this five simulations into one and so obtain a
single object ltraj with all the information on the whole movement. For
this purpose, I assigned to the five simulation the same ID but
different burst and than used the string c() to put all the simulations
in one object ltraj. However the five bursts remained separated.
Does anyone know how can I actually merge in adehabitatLT this tracks to
be one?
da = as.POSIXct(strptime(as.character(data$timestamp),
tz="",format="%d/%m/%Y %H:%M:%S"))
da[1:10]
[1] "2014-04-29 00:00:00 CEST" "2014-04-29 00:19:00 CEST"
[3] "2014-04-29 00:37:00 CEST" "2014-04-29 00:55:00 CEST"
[5] "2014-04-29 01:14:00 CEST" "2014-04-29 01:32:00 CEST"
[7] "2014-04-29 01:51:00 CEST" "2014-04-29 02:09:00 CEST"
[9] "2014-04-29 02:27:00 CEST" "2014-04-29 02:46:00 CEST"
hr<-simm.mou(da[1:70],b = c(0,
0),a=diag(c(0.001,0.001)),sigma=diag(c(100,100)),
x0 = c(0, 0),id = "Prova", burst = "HR1")
bb<-simm.bb(da[71:110],
begin=c(0,0),end=c(60000,60000),id="Prova",burst="migr")
hr2<-simm.mou(da[111:210],b = c(60000,
60000),a=diag(c(0.001,0.001)),sigma=diag(c(100,100)),
x0 = c(60000, 60000),id = "Prova",burst="HR2")
bb2<-simm.bb(da[211:240],
begin=c(60000,60000),end=c(0,0),id="Prova",burst="migr2")
hr3<-simm.mou(da[241:350],b = c(0,
0),a=diag(c(0.001,0.001)),sigma=diag(c(100,100)),
x0 = c(0, 0),id = "Prova", burst = "HR3")
migr<- c(hr,bb,hr2,bb2,hr3)
ld(migr)
x y date dx dy dist
dt
68 1273.749 7132.147 2014-04-29 20:31:00 -4814.308 -9767.770
10889.761 1080
69 -3540.560 -2635.623 2014-04-29 20:49:00 -2708.557 -1208.668
2966.001 1140
70 -6249.117 -3844.291 2014-04-29 21:08:00 NA NA
NA NA
111 0.000 0.000 2014-04-29 21:26:00 1444.220 1450.462
2046.854 1080
210 1444.220 1450.462 2014-04-29 21:44:00 1545.417 1533.493
2177.134 1140
310 2989.636 2983.956 2014-04-29 22:03:00 1505.393 1525.449
2143.176 1080
R2n abs.angle rel.angle id burst
pkey
68 52489961 -2.0287292 -3.1366093166 Prova HR1 Prova.2014-04-29
20:31:00
69 19482070 -2.7218694 -0.6931401962 Prova HR1 Prova.2014-04-29
20:49:00
70 53830039 NA NA Prova HR1 Prova.2014-04-29
21:08:00
111 0 0.7875548 NA Prova migr Prova.2014-04-29
21:26:00
210 4189611 0.7815256 -0.0060292387 Prova migr Prova.2014-04-29
21:44:00
310 17841916 0.7920154 0.0104897686 Prova migr Prova.2014-04-29
22:03:00
Thank you in advance for your time.
Giulia
More information about the R-sig-ecology
mailing list