diff --git a/tests/testthat/test-create_LDprofile.R b/tests/testthat/test-create_LDprofile.R
index 5bb46a23bfd52513b82b15cc087e961f3b643e1a..1d6cb9174776f15aff8b568bb3fc6962e1d13a10 100644
--- a/tests/testthat/test-create_LDprofile.R
+++ b/tests/testthat/test-create_LDprofile.R
@@ -14,14 +14,14 @@ df<-data.frame(
test_that("create_LDprofile calculates the LD profile correctly", {
expect_equal(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- data.frame(
- bin=c(0,0.001,0.002,0.003,0.004),
- rsq=c(0.285622427983539,0.280913978494624,0.263888888888889,0.319444444444444,NA),
- sd=c(0.270862044573862,0.201905775929377,0.321786617161322,0.142318760638328,NA),
- Beta_a=c(0.619957744381906,1.125028692019340,0.635410044952769,3.941019442363900,NA),
- Beta_b=c(1.062459890834270,2.446706389704430,1.149319432462400,8.454825333760550,NA),
- n=c(54,31,15,5,0)
- ),tolerance=0.0001)
+ data.frame(
+ bin=c(0,0.001,0.002,0.003,0.004),
+ rsq=c(0.285622427983539,0.280913978494624,0.263888888888889,0.319444444444444,NA),
+ sd=c(0.270862044573862,0.201905775929377,0.321786617161322,0.142318760638328,NA),
+ Beta_a=c(0.619957744381906,1.125028692019340,0.635410044952769,3.941019442363900,NA),
+ Beta_b=c(1.062459890834270,2.446706389704430,1.149319432462400,8.454825333760550,NA),
+ n=c(54,31,15,5,0)
+ ),tolerance=0.0001)
})
## Test the function with a different max_dist
@@ -29,28 +29,28 @@ test_that("create_LDprofile calculates the LD profile correctly", {
test_that("create_LDprofile calculates the LD profile correctly with a different max_dist", {
expect_equal(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = 0.003, beta_params = TRUE),
- data.frame(
- bin=c(0,0.001,0.002),
- rsq=c(0.285622427983539,0.280913978494624,0.263888888888889),
- sd=c(0.270862044573862,0.201905775929377,0.321786617161322),
- Beta_a=c(0.619957744381906,1.125028692019340,0.635410044952769),
- Beta_b=c(1.062459890834270,2.446706389704430,1.149319432462400),
- n=c(54,31,15)
- ),tolerance=0.0001)
+ data.frame(
+ bin=c(0,0.001,0.002),
+ rsq=c(0.285622427983539,0.280913978494624,0.263888888888889),
+ sd=c(0.270862044573862,0.201905775929377,0.321786617161322),
+ Beta_a=c(0.619957744381906,1.125028692019340,0.635410044952769),
+ Beta_b=c(1.062459890834270,2.446706389704430,1.149319432462400),
+ n=c(54,31,15)
+ ),tolerance=0.0001)
})
## Test the function with no max_dist given
test_that("create_LDprofile calculates the LD profile correctly with no max_dist supplied", {
expect_equal(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = 0.001, beta_params = TRUE),
- data.frame(
- bin=c(0,0.001,0.002,0.003),
- rsq=c(0.285622427983539,0.280913978494624,0.263888888888889,0.319444444444444),
- sd=c(0.270862044573862,0.201905775929377,0.321786617161322,0.142318760638328),
- Beta_a=c(0.619957744381906,1.125028692019340,0.635410044952769,3.941019442363900),
- Beta_b=c(1.062459890834270,2.446706389704430,1.149319432462400,8.454825333760550),
- n=c(54,31,15,5)
- ),tolerance=0.0001)
+ data.frame(
+ bin=c(0,0.001,0.002,0.003),
+ rsq=c(0.285622427983539,0.280913978494624,0.263888888888889,0.319444444444444),
+ sd=c(0.270862044573862,0.201905775929377,0.321786617161322,0.142318760638328),
+ Beta_a=c(0.619957744381906,1.125028692019340,0.635410044952769,3.941019442363900),
+ Beta_b=c(1.062459890834270,2.446706389704430,1.149319432462400,8.454825333760550),
+ n=c(54,31,15,5)
+ ),tolerance=0.0001)
})
@@ -59,14 +59,14 @@ test_that("create_LDprofile calculates the LD profile correctly with no max_dist
test_that("create_LDprofile calculates the LD profile correctly with a different bin size", {
expect_equal(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = 0.0005, beta_params = TRUE),
- data.frame(
- bin=c(0,0.0005,0.001,0.0015,0.002,0.0025,0.003,0.0035),
- rsq=c(0.238505747126437,0.340277777777778,0.283459595959596,0.274691358024691,0.215277777777778,0.361111111111111,0.288194444444444,0.444444444444445),
- sd=c(0.211600341827602,0.322468326753589,0.220527561550113,0.158590157477369,0.293808275018141,0.387895557,0.14316339,NA),
- Beta_a=c(0.916070145958307,0.637072700079744,1.046576044485340,1.909812912830260,0.775059123115346,1.088198634018290,3.789877096116000,NA),
- Beta_b=c(2.326350552394540,0.872215477086822,2.166981335251990,5.166454170350760,1.748740564135290,1.488374161884570,9.367197007381050,NA),
- n=c(29,25,22,9,10,5,4,1)
- ),tolerance=0.0001)
+ data.frame(
+ bin=c(0,0.0005,0.001,0.0015,0.002,0.0025,0.003,0.0035),
+ rsq=c(0.238505747126437,0.340277777777778,0.283459595959596,0.274691358024691,0.215277777777778,0.361111111111111,0.288194444444444,0.444444444444445),
+ sd=c(0.211600341827602,0.322468326753589,0.220527561550113,0.158590157477369,0.293808275018141,0.387895557,0.14316339,NA),
+ Beta_a=c(0.916070145958307,0.637072700079744,1.046576044485340,1.909812912830260,0.775059123115346,1.088198634018290,3.789877096116000,NA),
+ Beta_b=c(2.326350552394540,0.872215477086822,2.166981335251990,5.166454170350760,1.748740564135290,1.488374161884570,9.367197007381050,NA),
+ n=c(29,25,22,9,10,5,4,1)
+ ),tolerance=0.0001)
})
## Test the function with beta_params not specified
@@ -74,14 +74,14 @@ test_that("create_LDprofile calculates the LD profile correctly with a different
test_that("create_LDprofile calculates the LD profile correctly with beta_params not specified", {
expect_equal(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = 0.005),
- data.frame(
- bin=c(0,0.001,0.002,0.003,0.004),
- rsq=c(0.285622427983539,0.280913978494624,0.263888888888889,0.319444444444444,NA),
- sd=c(0.270862044573862,0.201905775929377,0.321786617161322,0.142318760638328,NA),
- Beta_a=c(NA,NA,NA,NA,NA),
- Beta_b=c(NA,NA,NA,NA,NA),
- n=c(54,31,15,5,0)
- ),tolerance=0.0001)
+ data.frame(
+ bin=c(0,0.001,0.002,0.003,0.004),
+ rsq=c(0.285622427983539,0.280913978494624,0.263888888888889,0.319444444444444,NA),
+ sd=c(0.270862044573862,0.201905775929377,0.321786617161322,0.142318760638328,NA),
+ Beta_a=c(NA,NA,NA,NA,NA),
+ Beta_b=c(NA,NA,NA,NA,NA),
+ n=c(54,31,15,5,0)
+ ),tolerance=0.0001)
})
## Test the function with a character matrix as x
@@ -92,14 +92,14 @@ df1[df1==2]<-"B"
test_that("create_LDprofile calculates the LD profile correctly with character matrix", {
expect_equal(create_LDprofile(dist = df1$dist, x = as.matrix(df1[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- data.frame(
- bin=c(0,0.001,0.002,0.003,0.004),
- rsq=c(0.285622427983539,0.280913978494624,0.263888888888889,0.319444444444444,NA),
- sd=c(0.270862044573862,0.201905775929377,0.321786617161322,0.142318760638328,NA),
- Beta_a=c(0.619957744381906,1.125028692019340,0.635410044952769,3.941019442363900,NA),
- Beta_b=c(1.062459890834270,2.446706389704430,1.149319432462400,8.454825333760550,NA),
- n=c(54,31,15,5,0)
- ),tolerance=0.0001)
+ data.frame(
+ bin=c(0,0.001,0.002,0.003,0.004),
+ rsq=c(0.285622427983539,0.280913978494624,0.263888888888889,0.319444444444444,NA),
+ sd=c(0.270862044573862,0.201905775929377,0.321786617161322,0.142318760638328,NA),
+ Beta_a=c(0.619957744381906,1.125028692019340,0.635410044952769,3.941019442363900,NA),
+ Beta_b=c(1.062459890834270,2.446706389704430,1.149319432462400,8.454825333760550,NA),
+ n=c(54,31,15,5,0)
+ ),tolerance=0.0001)
})
## Test all the checks
@@ -109,7 +109,7 @@ test_that("create_LDprofile calculates the LD profile correctly with character m
test_that("create_LDprofile fails when dist is non-numeric", {
expect_error(create_LDprofile(dist = paste0(df$dist,"dist"), x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- "dist must be a numeric vector or list of numeric vectors")
+ "dist must be a numeric vector or list of numeric vectors")
})
## Test the function with dists not a vector
@@ -117,7 +117,7 @@ test_that("create_LDprofile fails when dist is non-numeric", {
test_that("create_LDprofile fails when dist is not a vector", {
expect_error(create_LDprofile(dist = as.matrix(df$dist), x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- "dist must be a numeric vector or list of numeric vectors")
+ "dist must be a numeric vector or list of numeric vectors")
})
@@ -126,7 +126,7 @@ test_that("create_LDprofile fails when dist is not a vector", {
test_that("create_LDprofile fails when x is not a matrix", {
expect_error(create_LDprofile(dist = df$dist, x = df$C1, bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- "x must be a matrix or list of matrices")
+ "x must be a matrix or list of matrices")
})
## Test the function with x not having the correct amount of rows
@@ -134,7 +134,7 @@ test_that("create_LDprofile fails when x is not a matrix", {
test_that("create_LDprofile fails when the number of rows in x is not equal to the length of pos", {
expect_error(create_LDprofile(dist = df$dist, x = t(as.matrix(df[,3:7])), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- "The number of rows in x must equal the number of SNP genetic distances given in the corresponding dist")
+ "The number of rows in x must equal the number of SNP genetic distances given in the corresponding dist")
})
## Test the function with a SNP having only one allele
@@ -144,7 +144,7 @@ test_that("create_LDprofile fails when a SNP has only one allele", {
df1<-df
df1[1,3:7]<-1
expect_error(create_LDprofile(dist = df1$dist, x = as.matrix(df1[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- "SNPs must all be biallelic")
+ "SNPs must all be biallelic")
})
## Test the function with a SNP having more than two alleles
@@ -154,7 +154,7 @@ test_that("create_LDprofile fails when a SNP has more than two alleles", {
df1<-df
df1[1,7]<-3
expect_error(create_LDprofile(dist = df1$dist, x = as.matrix(df1[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- "SNPs must all be biallelic")
+ "SNPs must all be biallelic")
})
## Test the function with dists and x have a different number of elements
@@ -162,7 +162,7 @@ test_that("create_LDprofile fails when a SNP has more than two alleles", {
test_that("create_LDprofile fails when dist and x are different lengths", {
expect_error(create_LDprofile(dist = list(df$dist,df$dist), x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- "dist and x should contain the same number of elements")
+ "dist and x should contain the same number of elements")
})
@@ -171,7 +171,7 @@ test_that("create_LDprofile fails when dist and x are different lengths", {
test_that("create_LDprofile fails when bin_size is non-numeric", {
expect_error(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = "0.001cM", max_dist = 0.005, beta_params = TRUE),
- "bin_size must be a number greater than 0")
+ "bin_size must be a number greater than 0")
})
## Test the function with bin_size as negative
@@ -179,7 +179,7 @@ test_that("create_LDprofile fails when bin_size is non-numeric", {
test_that("create_LDprofile fails when bin_size is negative", {
expect_error(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = -1, max_dist = 0.005, beta_params = TRUE),
- "bin_size must be a number greater than 0")
+ "bin_size must be a number greater than 0")
})
## Test the function with max_dist as non-numeric
@@ -187,7 +187,7 @@ test_that("create_LDprofile fails when bin_size is negative", {
test_that("create_LDprofile fails when max_dist is non-numeric", {
expect_error(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = "0.005cM", beta_params = TRUE),
- "max_dist must be a number greater than 0")
+ "max_dist must be a number greater than 0")
})
## Test the function with max_dist as negative
@@ -195,7 +195,7 @@ test_that("create_LDprofile fails when max_dist is non-numeric", {
test_that("create_LDprofile fails when max_dist is negative", {
expect_error(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = -1, beta_params = TRUE),
- "max_dist must be a number greater than 0")
+ "max_dist must be a number greater than 0")
})
## Test the function with beta_params not logical
@@ -203,7 +203,7 @@ test_that("create_LDprofile fails when max_dist is negative", {
test_that("create_LDprofile fails when beta_params is not logical", {
expect_error(create_LDprofile(dist = df$dist, x = as.matrix(df[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = 1),
- "beta_params must be TRUE or FALSE")
+ "beta_params must be TRUE or FALSE")
})
## Test the function with missing value in x
@@ -212,14 +212,14 @@ df1$C1[15]<-NA
test_that("create_LDprofile calculates the LD profile correctly with missing value in x", {
expect_equal(create_LDprofile(dist = df1$dist, x = as.matrix(df1[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- data.frame(
- bin=c(0,0.001,0.002,0.003,0.004),
- rsq=c(0.302340534979424,0.298387096774194,0.256481481481481,0.297222222222222,NA),
- sd=c(0.285102946872176,0.231243979115867,0.324764476229430,0.125615767273278,NA),
- Beta_a=c(0.606846782972070,0.932888189465068,0.616307318328496,4.643089841249520,NA),
- Beta_b=c(0.942600409939103,1.692309859072240,1.139111486418360,11.012925076592600,NA),
- n=c(54,31,15,5,0)
- ),tolerance=0.0001)
+ data.frame(
+ bin=c(0,0.001,0.002,0.003,0.004),
+ rsq=c(0.302340534979424,0.298387096774194,0.256481481481481,0.297222222222222,NA),
+ sd=c(0.285102946872176,0.231243979115867,0.324764476229430,0.125615767273278,NA),
+ Beta_a=c(0.606846782972070,0.932888189465068,0.616307318328496,4.643089841249520,NA),
+ Beta_b=c(0.942600409939103,1.692309859072240,1.139111486418360,11.012925076592600,NA),
+ n=c(54,31,15,5,0)
+ ),tolerance=0.0001)
})
## Test the function with missing value in dist
@@ -228,17 +228,46 @@ df1$dist[5]<-NA
test_that("create_LDprofile calculates the LD profile correctly with missing value in dist", {
expect_equal(create_LDprofile(dist = df1$dist, x = as.matrix(df1[,3:7]), bin_size = 0.001, max_dist = 0.005, beta_params = TRUE),
- data.frame(
- bin=c(0,0.001,0.002,0.003,0.004),
- rsq=c(0.279320987654321,0.270833333333333,0.255952380952381,0.319444444444444,NA),
- sd=c(0.268871132983145,0.155580832042849,0.332406755525893,0.142318760638328,NA),
- Beta_a=c(0.634602938184746,1.570771368863860,0.611627386753874,3.941019442363900,NA),
- Beta_b=c(1.133751642734910,4.401600723631340,1.116161049339830,8.454825333760550,NA),
- n=c(45,27,14,5,0)
- ),tolerance=0.0001)
+ data.frame(
+ bin=c(0,0.001,0.002,0.003,0.004),
+ rsq=c(0.279320987654321,0.270833333333333,0.255952380952381,0.319444444444444,NA),
+ sd=c(0.268871132983145,0.155580832042849,0.332406755525893,0.142318760638328,NA),
+ Beta_a=c(0.634602938184746,1.570771368863860,0.611627386753874,3.941019442363900,NA),
+ Beta_b=c(1.133751642734910,4.401600723631340,1.116161049339830,8.454825333760550,NA),
+ n=c(45,27,14,5,0)
+ ),tolerance=0.0001)
})
## Test the function with fitdistrplus package not loaded and beta_params = TRUE
+#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)
+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),
+ data.frame(
+ bin=c(0),
+ rsq=c(0.0754751721676325),
+ sd=c(0.0928525038706074),
+ Beta_a=c(0.21155096488681),
+ Beta_b=c(2.8704915799235),
+ n=c(105)
+ ),tolerance=0.0001)
+})
+
+## Test the function when beta estimation doesn't work on the second try