так что у вас уже есть х и у,
затем вы помещаете их в фрейм данных и, используя библиотеку sf (простые функции), можете превратить его в пространственный объект, здесь он получится без проекции, поскольку это произвольный набор данных, а не что-то географическое, в противном случае вам следует объявление st_set_crs () с кодом ниже:
library(sf)
x <- runif(1000,min=0,max=10)
lambda=2*x+0.2*x*sin(x)
y <- rpois(1000,lam
bda)
df <- data.frame(x=x, y=y, z=runif(1000)) %>% st_as_sf(coords=c("x", "y"))
> df
Simple feature collection with 1000 features and 0 fields
geometry type: POINT
dimension: XY
bbox: xmin: 0.005045172 ymin: 0 xmax: 9.994533 ymax: 30
epsg (SRID): NA
proj4string: NA
First 10 features:
geometry
1 POINT (8.375505 20)
2 POINT (0.08116931 0)
3 POINT (3.786693 5)
4 POINT (7.68517 17)
5 POINT (9.363003 25)
6 POINT (5.114014 9)
7 POINT (5.70659 12)
8 POINT (9.936392 22)
9 POINT (9.164108 15)
10 POINT (7.524004 19)
plot(df)