version 1.0.7

DESCRIPTION

@@ -1,6 +1,6 @@
 Package: likelihoodR
 Title: Likelihood Analyses for Common Statistical Tests
-Version: 1.0.6
+Version: 1.0.7
 Authors@R: 
     person(given = "Peter",
            family = "Cahusac",
@@ -20,11 +20,10 @@ Description: A collection of functions that calculate the log likelihood
     P. Cahusac (2020, ISBN-13 : 978-1119549802).
 License: GPL-2
 Encoding: UTF-8
-LazyData: false
-RoxygenNote: 7.1.1
+RoxygenNote: 7.1.2
 NeedsCompilation: no
-Packaged: 2021-10-29 13:38:40 UTC; cahusac
+Packaged: 2021-12-01 14:05:03 UTC; peteq
 Author: Peter Cahusac [aut, cre] (<https://orcid.org/0000-0003-4976-2834>)
 Maintainer: Peter Cahusac <peteqsac@gmail.com>
 Repository: CRAN
-Date/Publication: 2021-10-29 14:10:02 UTC
+Date/Publication: 2021-12-01 14:20:02 UTC

MD5

@@ -1,9 +1,9 @@
-75952dc6795835f76b397e2ac38217b6 *DESCRIPTION
+c127441c8f6c05f4edb1a6ee94d6b815 *DESCRIPTION
 a17a73059d9ce8de56237528a0b3c965 *NAMESPACE
 5c55bd2e51514e28c91a13acc2faf1fb *R/L_1way_ANOVA.R
 5c03a9613ee4fb2e7ef1e321b8c2a599 *R/L_1way_RM_ANOVA.R
 0262779a5a0822be8d772eed6251e38c *R/L_1way_cat.R
-22aaf99cf48f11110412a72c88e7f456 *R/L_2S_ttest.R
+2ccdbbf5772076d3162db05027806d16 *R/L_2S_ttest.R
 5bb61c75c922e04646309a490ae43aa5 *R/L_2way_Factorial_ANOVA.R
 161048edea31e3721d5f8acd5b9074ef *R/L_2way_cat.R
 c251ad544c501246bb5a3fb86f6edb12 *R/L_OR.R

R/L_2S_ttest.R

@@ -90,10 +90,11 @@ L_2S_ttest <- function(data, group, veq=0, null=0, d=0.5, alt.2=NULL, L.int=2, v
   ad <- dat[complete.cases(dat), ] # remove missing, NA or NaN, case-wise
   dat <- ad$data
   gp <- ad$group
-  ns <- length(sort(dat[gp==1])) #sort to remove NA
-  nc <- length(sort(dat[gp==0]))
-  sd1 <- sd(dat[gp==1],na.rm=TRUE)
-  sd2 <- sd(dat[gp==0],na.rm=TRUE)
+  lev <- levels(factor(gp))
+  ns <- length(sort(dat[gp==lev[1]])) #sort to remove NA
+  nc <- length(sort(dat[gp==lev[2]]))
+  sd1 <- sd(dat[gp==lev[1]],na.rm=TRUE)
+  sd2 <- sd(dat[gp==lev[2]],na.rm=TRUE)
   SD <- sqrt(sd1^2+sd2^2)
   tres0 <- t.test(dat~gp, mu = null, paired=FALSE,var.equal=veq) #
   m.obs <- unname(tres0$estimate[1]-tres0$estimate[2]) # t test does it this way