1
1
library(spatialrisk )
2
- context(" Concentration " )
2
+ context(" concentration " )
3
3
4
4
test_that(" value is equal to concentration for small radius" , {
5
5
conc <- concentration(Groningen [1 ,], Groningen , value = amount , radius = 0.00001 , display_progress = FALSE )
@@ -10,15 +10,18 @@ test_that("sum of value column obtained from points_in_circle fn is equal to con
10
10
11
11
# Check first observation
12
12
uit1 <- concentration(Groningen [1 ,], Groningen , value = amount , radius = 200 , display_progress = FALSE )
13
- expect_equal(uit1 $ concentration , sum(points_in_circle(Groningen , lon_center = uit1 $ lon , lat_center = uit1 $ lat , radius = 200 )$ amount ))
13
+ sum1 <- points_in_circle(Groningen , lon_center = uit1 $ lon , lat_center = uit1 $ lat , radius = 200 )$ amount
14
+ expect_equal(uit1 $ concentration , sum(sum1 ))
14
15
15
16
# Check second observation
16
17
uit2 <- concentration(Groningen [2 ,], Groningen , value = amount , radius = 200 , display_progress = FALSE )
17
- expect_equal(uit2 $ concentration , sum(points_in_circle(Groningen , lon_center = uit2 $ lon , lat_center = uit2 $ lat , radius = 200 )$ amount ))
18
+ sum2 <- points_in_circle(Groningen , lon_center = uit2 $ lon , lat_center = uit2 $ lat , radius = 200 )$ amount
19
+ expect_equal(uit2 $ concentration , sum(sum2 ))
18
20
19
21
# Check last observation
20
22
uit3 <- concentration(Groningen [nrow(Groningen ),], Groningen , value = amount , radius = 200 , display_progress = FALSE )
21
- expect_equal(uit3 $ concentration , sum(points_in_circle(Groningen , lon_center = uit3 $ lon , lat_center = uit3 $ lat , radius = 200 )$ amount ))
23
+ sum3 <- points_in_circle(Groningen , lon_center = uit3 $ lon , lat_center = uit3 $ lat , radius = 200 )$ amount
24
+ expect_equal(uit3 $ concentration , sum(sum3 ))
22
25
})
23
26
24
27
test_that(" error should be returned for nonpositive radius" , {
0 commit comments