[R-sig-Geo] Problems in handling several animals in BRB calculation

Anna-Marie Corman annajess at gmx.de
Sat Apr 12 14:20:16 CEST 2014


Dear all,

I have a problem with the BRB calculation of my seabird tracking data. 
When I run the model R always has a fatal error and is aborted. Is there 
anything wrong in my ltraj-class or with the "D=" term in the BRB code? 
Does anyone see my beginner's mistake(s)?

Thank you again.

Best regards,
Anna

#Overview of my data (6 different individuals)
 > head(dcol[,c(1:11,17,25:27)])
         date day month year hour minute second time_sec time_diff latitude
1 13.05.2011  13     5 2011   21     11     33    76293       186 54.70463
2 13.05.2011  13     5 2011   21     14     37    76477       184 54.70479
3 13.05.2011  13     5 2011   21     17     40    76660       183 54.70548
4 13.05.2011  13     5 2011   21     20     44    76844       184 54.71007
5 13.05.2011  13     5 2011   21     23     48    77028       184 54.72099
6 13.05.2011  13     5 2011   21     26     52    77212       184 54.73233
   longitude  bird_id                 cet   xKoord  yKoord
1  8.349938 HA1_2011 2011-05-13 21:11:33 458111.9 6062117
2  8.349722 HA1_2011 2011-05-13 21:14:37 458098.1 6062135
3  8.344133 HA1_2011 2011-05-13 21:17:40 457738.7 6062216
4  8.338260 HA1_2011 2011-05-13 21:20:44 457365.1 6062730
5  8.320427 HA1_2011 2011-05-13 21:23:48 456227.9 6063956
6  8.312453 HA1_2011 2011-05-13 21:26:52 455726.7 6065223

ha11time <- dcol[, c("cet")]
ha11XY <- dcol[,c("xKoord", "yKoord")]
ha11ltraj <- as.ltraj(ha11XY, ha11time, id=dcol$bird_id,typeII=T)

##diffusion parameter
ha11D<-BRB.D(ha11ltraj,Tmax=420,Lmin=20)

##BRB
ha11_ud<-BRB(ha11ltraj,D=ha11_likD,type="UD",Tmax=420,tau=120,Lmin=20,radius=100, 
maxt=10*60,hmin=100,grid=500,extent=5)
##fatal error occurs



More information about the R-sig-Geo mailing list