Description Usage Arguments Details See Also Examples
View source: R/namespace.R View source: R/plot3D_dif.R
plot3D_dif
differentiates boundary and interior spatial locations in a spatial domain (a collection of
spatial locations with their coordinates). Since TSCS method is only capable of interpolation but not
extrapolation, it is necessary to highlight the difference between interior spatial locations and system boundary.
1 2 |
coords |
data frame; should only contain the three variables: X coordinate, Y coordinate and Z coordinate. Each row uniquely denotes a spatial location. (coordinates must be numeric) |
h1 |
numeric; side length of the unit cubic grid in X coordinate direction (horizontal). |
h2 |
numeric; side length of the unit cubic grid in Y coordinate direction (horizontal). |
v |
numeric; side length of the unit cubic grid in Z coordinate direction (vertical). |
xlab |
a label for the x axis, defaults to the name of X coordinate. |
ylab |
a label for the y axis, defaults to the name of Y coordinate. |
zlab |
a label for the z axis, defaults to the name of Z coordinate. |
title |
a main title for the plot. |
cex |
numeric; size of point to be plotted for each spatial location. (default: 3) |
The resulting plot is interactive, where the red points are interior spatial locations while the black points denote system boundary.
plot3D_dif
is exclusive to 3D rectangular grid system. Similarly, if you want to fathom how
this package handles 2D rectangular grid system, please refer to plot_dif
.
plot_dif
, plot3D_NA
, plot3D_map
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | ## Not run:
## TSCS spatial interpolation procedure:
basis <- tscsRegression3D(data = data, h1 = 3.75, h2 = 2.5, v = 5, alpha = 0.01);
basis$percentage
est <- tscsEstimate3D(matrix = basis$coef_matrix, newdata = newdata, h1 = 3.75, h2 = 2.5, v = 5);
str(est)
## comparison of estimates and true values:
plot_compare(est = est$estimate[,4], true = true)
index <- appraisal_index(est = est$estimate[,4], true = true);
index
## data visualization:
plot3D_dif(data = data[,1:3], h1 = 3.75, h2 = 2.5, v = 5)
plot3D_NA(newdata = newdata)
plot3D_map(newdata = newdata)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.