Advertisement
Guest User

Untitled

a guest
Nov 21st, 2014
142
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.35 KB | None | 0 0
  1. WELL.ID X Y BENZENE
  2. 1 MW-02 268.8155 282.83 0.00150
  3. 2 IW-06 271.6961 377.01 0.00050
  4. 3 IW-07 251.0236 300.41 0.01040
  5. 4 IW-08 278.9238 300.37 0.03190
  6. 5 MW-10 281.4008 414.15 2.04000
  7. 6 MW-12 391.3973 449.40 0.01350
  8. 7 MW-13 309.5307 335.55 0.01940
  9. 8 MW-15 372.8967 370.04 0.01620
  10. 9 MW-17 250.0000 428.04 0.01900
  11. 10 MW-24 424.4025 295.69 0.00780
  12. 11 MW-28 419.3205 250.00 0.00100
  13. 12 MW-29 352.9197 277.27 0.00031
  14. 13 MW-31 309.3174 370.92 0.17900
  15.  
  16. setwd("C:/.....")
  17. getwd()
  18.  
  19. require(geoR)
  20. require(ggplot2)
  21.  
  22.  
  23. a <- read.table("krigbenz_loc.csv", sep = ",", header = TRUE)
  24. b <- data.matrix(a)
  25. c <- as.geodata(b)
  26.  
  27. x.range <- as.integer(range(a[,2]))
  28. y.range <- as.integer(range(a[,3]))
  29. x = seq(from=x.range[1], to=x.range[2], by=1)
  30. y = seq(from=y.range[1], to=y.range[2], by=1)
  31. length(x)
  32. length(y)
  33. xv <- rep(x,length(y))
  34. yv <- rep(y, each=length(x))
  35. in_mat <- as.matrix(cbind(xv, yv))
  36.  
  37. q <- ksline(b, cov.model="exp", cov.pars=c(10,3.33), nugget=0, locations=in_mat)
  38.  
  39. Error in geodata$coords : $ operator is invalid for atomic vectors
  40.  
  41. q <- ksline(c, cov.model="exp", cov.pars=c(10,3.33), nugget=0, locations=in_mat)
  42.  
  43. image(q, val=q$predict)
  44.  
  45. Error in eval(x$call$geodata, envir = attr(x, "parent.env"))$borders :
  46. object of type 'builtin' is not subsettable
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement