From 2a7c5df03b3be58610e2032073d3e39e4a238d28 Mon Sep 17 00:00:00 2001 From: "Emerson M. Del Ponte" Date: Thu, 12 Oct 2023 19:41:06 -0300 Subject: [PATCH] update runs and doublets --- .DS_Store | Bin 6148 -> 10244 bytes R/doublets.test.R | 5 +---- R/oruns.test.R | 4 ++-- 3 files changed, 3 insertions(+), 6 deletions(-) diff --git a/.DS_Store b/.DS_Store index 115bc2ab2576a2975de6bdcebce34cfda05b7a55..69e04d94eae67ea339b9ef259b94cfcb5e49a671 100644 GIT binary patch literal 10244 zcmeHMU2GIp6h3ELU~Z?-De_kcxKOZwtL&EkQo(JvNXt(kY@wz6l-Zr3oj5zQ?96t7 zQd6U-G5*A;F^U?az8FY+Fc>915(R%AG{oQNgTDA^G@AIN=gyri%a-DUNi{k*nS1U% z=bn4#eD|BVbMF8EcICAgKo|fdy0}!TDZ4@9c79DsAzw3(B$7X*3Qjg<+9|4EnbJB? z1PBBO1PBBO1PBBO{2LIUHJc|%7Ba{K1OfyC1SS#S?}sp5T*iGlA*A=}pv>C>kgO!N zm#D1p0?Dv18ftt z=~8nCd2NsCj->o@pKY6yWi_?m&A@T(agzHW2eXOdtStZ{(Avx z%3^$>G(dAWmQ353Y)PBR;4^l2w|{scy;bwZ)us!VI!k>swsuo{c4Yd#N0T|rw5?(4 zkKZRB(}pxSC4hf;1Hh|h3UEe%*Nu%eHaAAXjcqXwcr4b~9tpR`nx`Pn5v4RO;%CV? z8KLQnCIqUX=ZX=2S)>$Wa&_lX z)uvT5ou@{6r)egQQ@Yciw#|ITat2)egv&L9j-7K8o^rq-U#=Tc=5znuS=X>FXTP3v zNHQj!o4aU5eM2PH(h=W!bzY59%h-a2{iB9sq)h#Q>ZEl`Gpyl#1w(U328o=OYaBE5 z9AmYiwSAVF(XTFEQdcW4lb090MB8h=l2R`XmUMKU@r_3wJj*XWAHdU0Z+nNcpc8e+wdNI0GH`x@MHKC zzJ~ANNB9|jfnVV__yZ+mScTPChpVt2S7Q|0aT9jnX6(WpxD$8bZtTOucm#(~#bM0i zam=BMr|}Fvgb(AR_%yzVFX7Ah3ciW&;JbJcFX1Qn8Ga+~FjRa&>=N~2j^9$Ll&<~> zX6VQb>3Yw$p2R;x*RSSsj?TYj!NNtiu3mTN#bEPEmcq5465vE88Z~m@l zALSOLIh;WfXUEd8yh>I$bwqWI+#o9swVct1XX4-y9-!Ka9Y&)~Df!q>`J*IVS9ka)h7b!E##1!QPF zwn`RYQpwP;Ky`Bs~EkGbZAV45M zAV6Ry5vUM@h57ma*6IKM&*UuyV+9BV2+S}7SlONI?xI19TjcWb`B{6IuA_AE#M?~> z>A6sL-H!mO@Z)&zq=0Px8G+{)Cq4I{5Yls}y5~-E=6?nR`~S&Ps9^v9KidEQ352^x A82|tP delta 105 zcmZn(XfcprU|?W$DortDU=RQ@Ie-{Mvv5r;6q~50$jG%ZU^g=(*Jd7p62{H*#oL({ rvvY6=G6R(WfdDs(44dP5<}d>Ql(P_h diff --git a/R/doublets.test.R b/R/doublets.test.R index ab04b31..4849b06 100644 --- a/R/doublets.test.R +++ b/R/doublets.test.R @@ -38,7 +38,7 @@ doublets.test <- function(x) { ZDb = (Db - EDb) / SDb # pvalue is the p-value of the Z-score - pvalue <- (2 * pnorm(abs(ZDb), lower.tail = FALSE)) + pvalue <- (pnorm(abs(ZDb), lower.tail = FALSE)) # result is "aggregation or clustering" if the Z-score is greater than 1.64, otherwise it is "randomness" result <- ifelse(abs(ZDb) >= 1.64, @@ -47,7 +47,4 @@ doublets.test <- function(x) { # Return the results as a list return(list(Db = Db, EDb = EDb, ZDb = ZDb, pvalue = pvalue, result = result)) - - - } diff --git a/R/oruns.test.R b/R/oruns.test.R index 87e0e18..4aefd57 100644 --- a/R/oruns.test.R +++ b/R/oruns.test.R @@ -32,11 +32,11 @@ oruns.test <- function(x) { Z = (U - EU) / sU # pvalue is the p-value of the Z-score - pvalue <- (2 * pnorm(abs(Z), lower.tail = FALSE)) + pvalue <- (pnorm(abs(Z), lower.tail = FALSE)) # result is "clustering" if the Z-score is less than 1.64, otherwise it is "randomness" result <- ifelse(Z < 1.64, - c("clustering"), + c("aggregation or clustering"), c("randomness")) # Return the results as a list