Dear list,

I am using the package adehabitatLT to create null models using the function bellow. I set my treatment function to NULL and I included a constraint function based on an MCP of all my animal locations.

my.CRW <- NMs.randomCRW(na.omit(inter3), rangles = TRUE, rdist = TRUE, fixedStart = TRUE,
            x0 = NULL, rx = NULL, ry = NULL, treatment.func = NULL,
            treatment.par = NULL, constraint.func = consfun,
            constraint.par = GHA_const, nrep = 15)

my.sim<- testNM(my.CRW)

When I was ordering the data by dates to extract my environmental variables, I realised that no matter the number of nrep I perform, my results (locations of the null model) for a specific date are all the same. I don’t think this is right, correct? The idea of generating the crw is to have a possible distribution where the animal could have gone. Although the results are different than the original coordinates, instead of generating 15 different trajectories it is just generating one. Here an example of my output:

        x       y       date    Id      CRW_Id
1       4775351 1059107 13/05/2003 20:27         O66A   1
50      4775351 1059107 13/05/2003 20:27         O66A   2
99      4775351 1059107 13/05/2003 20:27         O66A   3
148     4775351 1059107 13/05/2003 20:27         O66A   4
197     4775351 1059107 13/05/2003 20:27         O66A   5
246     4775351 1059107 13/05/2003 20:27         O66A   6
295     4775351 1059107 13/05/2003 20:27         O66A   7
344     4775351 1059107 13/05/2003 20:27         O66A   8
393     4775351 1059107 13/05/2003 20:27         O66A   9
442     4775351 1059107 13/05/2003 20:27         O66A   10


Could you help me undertand that? Any possible solution?

Thanks a lot for your help!

Best,

Ana


--

PhD student

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

Reply via email to