diff --git a/tests/testthat/test-create_LDprofile.R b/tests/testthat/test-create_LDprofile.R index 051ae0d665c40dbbd6d97666f135ef94a67635a0..27a587899faa5ed5c03f2c05e97f8894550a4702 100644 --- a/tests/testthat/test-create_LDprofile.R +++ b/tests/testthat/test-create_LDprofile.R @@ -243,31 +243,37 @@ test_that("create_LDprofile calculates the LD profile correctly with missing val #Tested manually ## Test the function when beta estimation doesn't work the first try -x1<-matrix(c(1,1,1,2,1,2,1,1,1,2,1,2,2,2,2, - 1,1,1,2,1,2,2,2,2,1,1,2,1,1,2, - 1,2,1,2,1,1,1,2,1,1,2,1,2,2,2, - 1,2,1,2,1,1,2,1,2,2,1,2,2,1,2, - 1,2,1,2,1,2,1,2,1,2,1,1,2,1,2, - 1,2,2,1,1,2,1,2,2,1,2,1,1,1,2, - 1,1,2,2,1,1,2,1,2,2,2,1,1,1,2, - 1,1,2,1,1,1,2,2,1,2,2,1,1,1,2, - 1,1,1,1,1,2,2,2,2,1,2,1,2,2,2, - 1,2,2,2,2,1,2,1,2,1,2,2,2,1,2, - 1,2,1,1,1,1,2,1,1,1,2,1,1,2,2, - 1,1,1,2,2,1,2,1,2,1,1,2,2,1,2, - 1,1,1,1,1,1,1,2,2,2,1,2,2,1,2, - 1,2,2,1,2,1,1,1,2,1,2,1,2,1,2, - 1,1,2,1,1,1,2,2,2,1,1,1,2,1,2),byrow=TRUE,nrow=15) +x1<-matrix(c(1,1,2,1,2,1,1,2,1,2,2,1,2,1,2, + 1,1,1,1,1,2,2,2,2,1,1,1,2,2,2, + 1,1,1,1,1,1,2,1,2,2,2,1,1,2,2, + 1,2,2,1,1,2,1,1,2,1,1,2,1,2,2, + 1,2,2,1,2,1,1,2,2,1,2,2,1,2,2, + 1,2,2,2,2,1,2,1,1,1,1,1,1,1,2, + 1,2,2,1,1,1,1,2,2,2,2,2,1,1,2, + 1,2,1,2,1,1,1,2,1,1,1,2,2,2,2, + 1,2,2,2,1,1,2,1,1,1,2,2,2,1,2, + 1,1,2,1,2,1,1,1,1,1,2,1,2,1,2, + 1,1,1,2,1,2,2,1,2,1,1,1,2,1,2, + 1,2,2,2,2,2,1,2,1,2,1,1,1,1,2, + 1,2,1,2,1,1,1,1,2,2,2,1,2,1,2, + 1,2,1,2,2,1,2,1,2,1,2,1,2,1,2, + 1,2,2,2,2,1,2,2,1,1,2,2,2,2,2, + 1,1,2,2,2,1,2,2,1,1,1,1,1,2,2, + 1,2,1,2,2,2,1,2,2,1,1,1,1,2,2, + 1,2,1,2,1,2,2,2,1,2,1,2,2,1,2, + 1,1,1,2,1,2,2,2,1,2,2,2,1,2,2, + 1,1,1,2,2,2,1,1,1,2,1,2,2,1,2),byrow=TRUE,nrow=20) test_that("create_LDprofile calculates the LD profile correctly when beta calculation fails on first attempt", { - expect_equal(create_LDprofile(dist = rep(0,15), x = x1, bin_size = 0.001, max_dist = 0.001, beta_params = TRUE), + expect_equal(create_LDprofile(dist = rep(0,20), x = x1, bin_size = 0.001, beta_params = TRUE), data.frame( bin=c(0), - rsq=c(0.0754751721676325), - sd=c(0.0928525038706074), - Beta_a=c(0.795140679587351), - Beta_b=c(9.0658860153924), - n=c(105) + rsq=c(0.0665701124504383), + sd=c(0.0932577303343809), + Beta_a=c(0.641838135131359), + Beta_b=c(8.47114586606102), + n=c(190) ),tolerance=0.0001) }) ## Test the function when beta estimation doesn't work on the second try +## Tested manually