Hjælp til programmering i R - fejl ved plotting af koordinater
Hej Folkens,Jeg har et kæmpe problem ved følgende kode. Jeg får følgende fejlkode: "Error in if (antipodal(p1, p2)) { : missing value where TRUE/FALSE needed" og jeg ved ikke hvad jeg skal gøre. Har søgt nettet, men kan ikke finde noget brugbar hjælp.
Det jeg gerne vil er at plotte liner fra ét postnummer til et andet via koordinater. Jeg har afgrænset mig til europa, så jeg skulle være ude over det der antipodal problem.....men det er jeg så ikke alligevel!
Nogen der kan hjælpe???
>
> library(maps)
> library(geosphere)
>
>
> postcodes <- read.csv2("Postcodes.csv", header=TRUE)
> deliveries <- read.csv2("Deliveries.csv", header=TRUE, as.is=TRUE)
> xlim <- c(-24.3, 43.4)
> ylim <- c(31.4, 71.7)
> map("world", col="#f2f2f2", fill=TRUE, bg="white", lwd=0.05, xlim=xlim, ylim=ylim)
> fsub <- deliveries[deliveries$ShipPostcode == "DK-4600",]
> for (j in 1:length(fsub$postcodes)) {
+ ShipPostcode <- postcodes[postcodes$DW_postcode == fsub[j,]$ShipPostcode,]
+ DeliveryPostcode <- postcodes[postcodes$DW_postcode == fsub[j,]$DeliveryPostcode,]
+
+ inter <- gcIntermediate(c(ShipPostcode[1,]$Long, ShipPostcode[1,]$Lat),c(DeliveryPostcode[1,]$Long, DeliveryPostcode[1,]$Lat), n=100, addStartEnd=TRUE,)
+
+ lines(inter, col="red", lwd=0.5)
+ }