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