Displaying 1 result from an estimated 1 matches for "rel_dom".
Did you mean:
get_dom
2008 Jul 25
0
Error in vector("double", length)
...newx=(xcoord-min(xcoord))/1000
newy=(ycoord-min(ycoord))/1000
plot(newx,newy,pch=20,xlab="X coordinate (km)",ylab="Y coordinate (km)")
#set of new locations to be used from here on
tr.locs=cbind(newx,newy)
#transform proportions using an angular transformation
rel.tr=asin(sqrt(rel_dom/100))
#set as geodata
rd<-as.geodata(cbind(tr.locs, rel.tr))
#calculate empirical variogram
vg.deft<-variog(rd,uvec=10)
#define distance-based neighborhood based on function written in separate code
nb.tr=dist.neighbors(tr.locs,2)
-----------------------------------------------------------...