Sorry, I forgot the ltraj-summary:

 > summary(ha11ltraj)
         id    burst nb.reloc NAs          date.begin date.end
1 HA1_2011 HA1_2011       60   0 2011-05-13 21:11:33 2011-05-14 00:19:11
2 HA2_2011 HA2_2011     3619   0 2011-05-13 04:45:10 2011-05-30 10:52:24
3 HA3_2011 HA3_2011      408   0 2011-05-12 03:59:37 2011-05-12 15:11:04
4 HA4_2011 HA4_2011       33   0 2011-05-13 18:40:58 2011-05-13 20:22:38
5 HA5_2011 HA5_2011     3474   0 2011-05-13 21:17:33 2011-05-31 14:15:02
6 HA6_2011 HA6_2011      465   0 2011-05-13 16:00:32 2011-05-15 22:39:29


-------- Original-Nachricht --------
Betreff: 	[R-sig-Geo] Problems in handling several animals in BRB 
calculation
Datum: 	Sat, 12 Apr 2014 14:20:16 +0200
Von: 	Anna-Marie Corman <annajess@gmx.de>
An: 	r-sig-geo@r-project.org



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

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo




	[[alternative HTML version deleted]]

