Hola,
Aparte de los errores en tu código, la forma más fácil y eficiente
que se me ocurre es emplear la función loccoords del paquete geoR. El
código quedaría así:
ligustro <-read.table("arboles.txt", header=T)
bosquetotal<-read.table("bosquecompleto.txt", header=T)
library(geoR)
# ?loccoords: returns a n x N matrix with distances between data points
and prediction locations.
distalig <- loccoords(bosquetotal[,c('x','y')], ligustro[,c('x','y')])
Un saludo, Rubén.
P.D. Entiendo que quieres hacer algún tipo de análisis espacial (el
problema parece de procesos puntuales), mira los paquetes que ya
implementan herramientas en R
(http://cran.r-project.org/web/views/Spatial.html), o el libro Applied
Spatial Data Analysis with R (http://www.asdar-book.org).
El 23/04/2015 a las 17:37, Carlos J. Gil Bellosta escribió:
Hola, ¿qué tal?
En el doble bucle de tu código tienes los índices i y j
intercambiados. Además, buscas las columnas 3 y 4 de bosque, que solo
tiene 2.
distalig<-matrix (NA, nrow (lig), nrow (bosque))
for (i in 1: nrow(lig)){
for (j in 1: nrow (bosque)){
distalig[i,j]<-dist(lig[i,1],lig [i,2],bosque [j,1],bosque[j,2])
}
}
en cambio, "funciona". También lo hace algo así como
a <- outer(lig[,1], bosque[,1], "-")
b <- outer(lig[,2], bosque[,2], "-")
distalig <- sqrt(a^2 + b^2)
y da los mismos resultados.
Un saludo,
Carlos J. Gil Bellosta
http://www.datanalytics.com
El día 23 de abril de 2015, 17:22, Priscila Ana Powell
<[email protected]> escribió:
¡Hola!
Estoy tratando de hacer una matriz que contenga los valores de distancias
entre distintos elementos (con sus respectivas coordenadas x e y), pero
salta un error, y no sé como corregirlo. ¿alguna idea?
muchas gracias!
saludos!!
Priscila
--
Dra. Priscila Ana Powell
Instituto de Ecología Regional
Facultad de Ciencias Naturales e Instituto Miguel Lillo
Universidad Nacional de Tucumán
Argentina
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
