library(terra)
## terra version 1.1.2
##
## Attaching package: 'terra'
## The following object is masked from 'package:knitr':
##
##     spin


longitude <- c(-72.875691, -72.869782, -72.862467, -72.858593, -72.850833, -72.849032, -72.858058)
latitude <- c(2.212863, 2.219539, 2.22007, 2.227792, 2.230666, 2.216927, 2.2097)
lonlat <- cbind(longitude, latitude)


pts <- vect(lonlat)
class(pts)
[1] "SpatVector"
attr(,"package")
[1] "terra"
pts
 class       : SpatVector 
 geometry    : points 
 dimensions  : 7, 0  (geometries, attributes)
 extent      : -72.87569, -72.84903, 2.2097, 2.230666  (xmin, xmax, ymin, ymax)
 coord. ref. :  
geom(pts)
     geom part         x        y hole
[1,]    1    1 -72.87569 2.212863    0
[2,]    2    1 -72.86978 2.219539    0
[3,]    3    1 -72.86247 2.220070    0
[4,]    4    1 -72.85859 2.227792    0
[5,]    5    1 -72.85083 2.230666    0
[6,]    6    1 -72.84903 2.216927    0
[7,]    7    1 -72.85806 2.209700    0
crdref <- "+proj=longlat +datum=WGS84"
pts <- vect(lonlat, crs=crdref)
pts
 class       : SpatVector 
 geometry    : points 
 dimensions  : 7, 0  (geometries, attributes)
 extent      : -72.87569, -72.84903, 2.2097, 2.230666  (xmin, xmax, ymin, ymax)
 coord. ref. : +proj=longlat +datum=WGS84 +no_defs 
crs(pts)
[1] "GEOGCRS[\"unknown\",\n    DATUM[\"World Geodetic System 1984\",\n        ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n            LENGTHUNIT[\"metre\",1]],\n        ID[\"EPSG\",6326]],\n    PRIMEM[\"Greenwich\",0,\n        ANGLEUNIT[\"degree\",0.0174532925199433],\n        ID[\"EPSG\",8901]],\n    CS[ellipsoidal,2],\n        AXIS[\"longitude\",east,\n            ORDER[1],\n            ANGLEUNIT[\"degree\",0.0174532925199433,\n                ID[\"EPSG\",9122]]],\n        AXIS[\"latitude\",north,\n            ORDER[2],\n            ANGLEUNIT[\"degree\",0.0174532925199433,\n                ID[\"EPSG\",9122]]]]"

# Generate random precipitation values, same quantity as points
precipvalue <- runif(nrow(lonlat), min=100, max=1000)
df <- data.frame(ID=1:nrow(lonlat), precip=precipvalue)

#Combine Spatial vector with the Dataframe
ptv <- vect(lonlat, atts=df, crs=crdref)

ptv
 class       : SpatVector 
 geometry    : points 
 dimensions  : 7, 2  (geometries, attributes)
 extent      : -72.87569, -72.84903, 2.2097, 2.230666  (xmin, xmax, ymin, ymax)
 coord. ref. : +proj=longlat +datum=WGS84 +no_defs 
LS0tDQp0aXRsZTogIkNvb3JkcyBpbnRvIGRhdGFmcmFtZSINCm91dHB1dDogaHRtbF9ub3RlYm9vaw0KLS0tDQoNCg0KDQpgYGB7cn0NCmxpYnJhcnkodGVycmEpDQojIyB0ZXJyYSB2ZXJzaW9uIDEuMS4yDQojIw0KIyMgQXR0YWNoaW5nIHBhY2thZ2U6ICd0ZXJyYScNCiMjIFRoZSBmb2xsb3dpbmcgb2JqZWN0IGlzIG1hc2tlZCBmcm9tICdwYWNrYWdlOmtuaXRyJzoNCiMjDQojIyAgICAgc3Bpbg0KDQoNCmxvbmdpdHVkZSA8LSBjKC03Mi44NzU2OTEsIC03Mi44Njk3ODIsIC03Mi44NjI0NjcsIC03Mi44NTg1OTMsIC03Mi44NTA4MzMsIC03Mi44NDkwMzIsIC03Mi44NTgwNTgpDQpsYXRpdHVkZSA8LSBjKDIuMjEyODYzLCAyLjIxOTUzOSwgMi4yMjAwNywgMi4yMjc3OTIsIDIuMjMwNjY2LCAyLjIxNjkyNywgMi4yMDk3KQ0KbG9ubGF0IDwtIGNiaW5kKGxvbmdpdHVkZSwgbGF0aXR1ZGUpDQoNCg0KcHRzIDwtIHZlY3QobG9ubGF0KQ0KY2xhc3MocHRzKQ0KDQpwdHMNCg0KZ2VvbShwdHMpDQoNCmNyZHJlZiA8LSAiK3Byb2o9bG9uZ2xhdCArZGF0dW09V0dTODQiDQpwdHMgPC0gdmVjdChsb25sYXQsIGNycz1jcmRyZWYpDQpwdHMNCg0KY3JzKHB0cykNCg0KDQoNCg0KDQoNCg0KYGBgDQoNCg0KYGBge3J9DQoNCiMgR2VuZXJhdGUgcmFuZG9tIHByZWNpcGl0YXRpb24gdmFsdWVzLCBzYW1lIHF1YW50aXR5IGFzIHBvaW50cw0KcHJlY2lwdmFsdWUgPC0gcnVuaWYobnJvdyhsb25sYXQpLCBtaW49MTAwLCBtYXg9MTAwMCkNCmRmIDwtIGRhdGEuZnJhbWUoSUQ9MTpucm93KGxvbmxhdCksIHByZWNpcD1wcmVjaXB2YWx1ZSkNCg0KI0NvbWJpbmUgU3BhdGlhbCB2ZWN0b3Igd2l0aCB0aGUgRGF0YWZyYW1lDQpwdHYgPC0gdmVjdChsb25sYXQsIGF0dHM9ZGYsIGNycz1jcmRyZWYpDQoNCnB0dg0KDQoNCg0KYGBgDQoNCg0KYGBge3J9DQoNCg0KDQpgYGANCg0KDQpgYGB7cn0NCg0KYGBgDQoNCg0KYGBge3J9DQoNCmBgYA0KDQoNCmBgYHtyfQ0KYGBgDQoNCg==