Hi everyone,
I've searched a few R lists for an answer to this dilemma but it's got me
stumped. I have a dataset not unlike the Puechabon dataset that's included with
the Adehabitat package:
> data(puechabon)
> p = puechabon$locs
> p = p[-c(94:120,53:90,37:50,7:29),]
> p
Name Age Sex X Y Date
1 Brock 2 1 699889 3161559 930701
2 Brock 2 1 700046 3161541 930703
3 Brock 2 1 698840 3161033 930706
4 Brock 2 1 699809 3161496 930707
5 Brock 2 1 698627 3160941 930708
6 Brock 2 1 698719 3160989 930709
30 Brock 2 1 699207 3161099 930831
31 Calou 2 1 700098 3161648 930703
32 Calou 2 1 700310 3161323 930704
33 Calou 2 1 699877 3161471 930707
34 Calou 2 1 700003 3161360 930708
35 Calou 2 1 699820 3161387 930709
36 Calou 2 1 699778 3161328 930711
52 Chou 3 2 699640 3158702 920802
53 Chou 3 2 699689 3158857 920803
92 Jean 2 1 699448 3158689 930705
93 Jean 2 1 699500 3158684 930706
94 Jean 2 1 699409 3158794 930707
I'm trying to create an MCP for each animal that has 5 or more relocation
points (e.g. Brock & Calou) and export each MCP to a CSV file. I've tried the
following, but I get an error message:
> summary(p$Name)>=5
Brock Calou Chou Jean
TRUE TRUE FALSE FALSE
> for (i in 1:nrow(p)){
+ if(summary(p$Name)>=5){
+ mcp.p = mcp(pxy, pid, percent=95)
+ write.table(mcp.p, sep=",", file=paste(pid,".csv"), col.names=NA)
+ }
+ }
Error in mcp(pxy, pid, percent = 95) :
At least 5 relocations are required to fit an home range
In addition: Warning message:
In if (summary(p$Name) >= 5) { :
the condition has length > 1 and only the first element will be used
I'm not an expert in R and would appreciate any help.
Regards,
_________________________________
Hakim Abdi, M.Sc.
Research Intern
Conservation GIS/Remote Sensing Lab
Smithsonian Conservation Biology Institute
1500 Remount Road
Front Royal, VA 22630
phone: +1 540 635 6578
mobile: +1 747 224 7006
fax: +1 540 635 6506 (Attn:GIS Lab)
email: ab...@si.edu
http://nationalzoo.si.edu/SCBI/ConservationGIS/
[[alternative HTML version deleted]]
_______________________________________________
R-sig-ecology mailing list
R-sig-ecology@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
