diff --git a/R/tf.R b/R/tf.R index 9dada5ff..69ba5249 100644 --- a/R/tf.R +++ b/R/tf.R @@ -77,10 +77,7 @@ #' # plotting the TAI using log2 transformed expression levels #' # and performing the Flat Line Test to obtain the p-value #' PlotSignature(ExpressionSet = tf(PhyloExpressionSetExample, log2), -#' permutations = 10, -#' type = "l", -#' lwd = 5, -#' TestStatistic = "FlatLineTest") +#' permutations = 1000) #' #' #' diff --git a/docs/pkgdown.yml b/docs/pkgdown.yml index 13954456..9b340123 100644 --- a/docs/pkgdown.yml +++ b/docs/pkgdown.yml @@ -8,7 +8,7 @@ articles: Intermediate: Intermediate.html Introduction: Introduction.html Phylostratigraphy: Phylostratigraphy.html -last_built: 2024-10-03T15:57Z +last_built: 2024-10-03T16:07Z urls: reference: https://drostlab.github.io/myTAI/reference article: https://drostlab.github.io/myTAI/articles diff --git a/docs/reference/EarlyConservationTest.html b/docs/reference/EarlyConservationTest.html index aa0a0f62..740d98a5 100644 --- a/docs/reference/EarlyConservationTest.html +++ b/docs/reference/EarlyConservationTest.html @@ -191,10 +191,10 @@

Examples#> [ Computing variances of permuted transcriptome signatures ... ] #> #> $p.value -#> [1] 0.9998591 +#> [1] 0.9998289 #> #> $std.dev -#> [1] 0.05393337 0.05301197 0.05284169 0.05161798 0.05042234 0.05343254 0.05679332 +#> [1] 0.05369442 0.05209210 0.05093189 0.04996095 0.04901245 0.05003862 0.05481401 #> #> $lillie.test #> [1] NA @@ -217,11 +217,11 @@

Examples custom.perm.matrix = custom_perm_matrix) #> The phylotranscriptomic pattern may not follow an early conservation pattern (low-mid-high or low-high-high). #> $p.value -#> [1] 0.9999901 +#> [1] 0.9999693 #> #> $std.dev #> Zygote Quadrant Globular Heart Torpedo Bent Mature -#> 0.05246516 0.05312139 0.04978308 0.04955870 0.04989766 0.05093588 0.05233829 +#> 0.05635929 0.05285043 0.04957050 0.05001778 0.04811156 0.05039466 0.05989830 #> #> $lillie.test #> [1] NA diff --git a/docs/reference/FlatLineTest.html b/docs/reference/FlatLineTest.html index 8c61cd2b..4f824460 100644 --- a/docs/reference/FlatLineTest.html +++ b/docs/reference/FlatLineTest.html @@ -166,17 +166,17 @@

Examples#> #> -> We recommended using at least 20000 permutations to achieve a sufficient permutation test. #> $p.value -#> [1] 3.331422e-16 +#> [1] 8.488052e-10 #> #> $std.dev -#> [1] 0.05303712 0.05209588 0.04971351 0.04849724 0.05027523 0.04967010 0.05683321 +#> [1] 0.05830159 0.05826165 0.05508265 0.05585175 0.05526261 0.05392054 0.06382353 #> #> $ks.test #> #> Exact one-sample Kolmogorov-Smirnov test #> #> data: filtered_vars -#> D = 0.04664, p-value = 0.9788 +#> D = 0.062749, p-value = 0.8118 #> alternative hypothesis: two-sided #> #> @@ -197,18 +197,18 @@

Examples#> #> -> We recommended using at least 20000 permutations to achieve a sufficient permutation test. #> $p.value -#> [1] 4.115314e-44 +#> [1] 3.048669e-14 #> #> $std.dev #> Zygote Quadrant Globular Heart Torpedo Bent Mature -#> 0.06078755 0.05841072 0.05539504 0.05449748 0.05275256 0.05482176 0.06069415 +#> 0.05154281 0.05174183 0.05101208 0.04891800 0.04824081 0.05254812 0.05079183 #> #> $ks.test #> #> Exact one-sample Kolmogorov-Smirnov test #> #> data: filtered_vars -#> D = 0.073016, p-value = 0.7466 +#> D = 0.046018, p-value = 0.9785 #> alternative hypothesis: two-sided #> #> diff --git a/docs/reference/LateConservationTest.html b/docs/reference/LateConservationTest.html index 94a46ebf..d5485ffc 100644 --- a/docs/reference/LateConservationTest.html +++ b/docs/reference/LateConservationTest.html @@ -190,10 +190,10 @@

Examples#> [ Computing variances of permuted transcriptome signatures ... ] #> #> $p.value -#> [1] 0.9999981 +#> [1] 0.9999991 #> #> $std.dev -#> [1] 0.05344308 0.05246857 0.05222011 0.05122469 0.04998389 0.05116365 0.05225663 +#> [1] 0.05343829 0.05299432 0.05104522 0.04950552 0.04919799 0.05051301 0.05572345 #> #> $lillie.test #> [1] NA @@ -216,11 +216,11 @@

Examples custom.perm.matrix = custom_perm_matrix) #> The phylotranscriptomic pattern may not follow a late conservation pattern (high-mid-low or high-high-low). #> $p.value -#> [1] 0.999999 +#> [1] 0.9999949 #> #> $std.dev #> Zygote Quadrant Globular Heart Torpedo Bent Mature -#> 0.05257291 0.05522342 0.05261010 0.05158845 0.05288796 0.05216044 0.05553098 +#> 0.05747708 0.05585285 0.05300770 0.05001848 0.04777103 0.04897164 0.05483565 #> #> $lillie.test #> [1] NA diff --git a/docs/reference/PairwiseTest.html b/docs/reference/PairwiseTest.html index 876f38dd..f1e26eb9 100644 --- a/docs/reference/PairwiseTest.html +++ b/docs/reference/PairwiseTest.html @@ -201,10 +201,10 @@

Examples#> [ Computing variances of permuted transcriptome signatures ... ] #> #> $p.value -#> [1] 0.01971706 +#> [1] 0.02584971 #> #> $std.dev -#> [1] 0.05614200 0.05508994 0.05334782 0.05143772 0.05121888 0.05221399 0.05667816 +#> [1] 0.05408029 0.05251915 0.05223826 0.05091721 0.04997079 0.05260522 0.05694218 #> #> $lillie.test #> [1] NA @@ -223,10 +223,10 @@

Examples#> [ Computing variances of permuted transcriptome signatures ... ] #> #> $p.value -#> [1] 0.973595 +#> [1] 0.9729322 #> #> $std.dev -#> [1] 0.05555567 0.05472251 0.05274843 0.05209187 0.05074195 0.05208817 0.05569443 +#> [1] 0.05744036 0.05659292 0.05449721 0.05255368 0.05186078 0.05305327 0.05646092 #> #> $lillie.test #> [1] NA @@ -245,10 +245,10 @@

Examples#> [ Computing variances of permuted transcriptome signatures ... ] #> #> $p.value -#> [1] 0.3856208 +#> [1] 0.4174285 #> #> $std.dev -#> [1] 0.05568058 0.05449264 0.05283244 0.05168838 0.05033825 0.05185507 0.05587415 +#> [1] 0.05238544 0.05178879 0.05066156 0.04913529 0.04877068 0.05127811 0.05492211 #> #> $lillie.test #> [1] NA @@ -270,11 +270,11 @@

Examples altHypothesis = "greater", custom.perm.matrix = custom_perm_matrix) #> $p.value -#> [1] 0.04031736 +#> [1] 0.02825196 #> #> $std.dev #> Zygote Quadrant Globular Heart Torpedo Bent Mature -#> 0.05655089 0.05535914 0.05549875 0.05317636 0.04994245 0.05009926 0.05629529 +#> 0.05294514 0.04831905 0.04589758 0.04496874 0.04385149 0.04804513 0.05538722 #> #> $lillie.test #> [1] NA diff --git a/docs/reference/PlotPattern-2.png b/docs/reference/PlotPattern-2.png index e98168a2..e23344be 100644 Binary files a/docs/reference/PlotPattern-2.png and b/docs/reference/PlotPattern-2.png differ diff --git a/docs/reference/PlotPattern-3.png b/docs/reference/PlotPattern-3.png index 79353b92..26bb7d7d 100644 Binary files a/docs/reference/PlotPattern-3.png and b/docs/reference/PlotPattern-3.png differ diff --git a/docs/reference/PlotPattern-4.png b/docs/reference/PlotPattern-4.png index 2d31db8d..f97bdecb 100644 Binary files a/docs/reference/PlotPattern-4.png and b/docs/reference/PlotPattern-4.png differ diff --git a/docs/reference/PlotPattern-5.png b/docs/reference/PlotPattern-5.png index 438f1ed5..e9b2cac7 100644 Binary files a/docs/reference/PlotPattern-5.png and b/docs/reference/PlotPattern-5.png differ diff --git a/docs/reference/PlotPattern.html b/docs/reference/PlotPattern.html index 100ca52d..3fd434d3 100644 --- a/docs/reference/PlotPattern.html +++ b/docs/reference/PlotPattern.html @@ -216,7 +216,7 @@

Examples#> [ Computing variances of permuted transcriptome signatures ... ] #> #> -#> Total runtime of your permutation test: 0.024 seconds. +#> Total runtime of your permutation test: 0.015 seconds. #> #> -> We recommended using at least 20000 permutations to achieve a sufficient permutation test. diff --git a/docs/reference/PlotSignature-1.png b/docs/reference/PlotSignature-1.png index e76c908c..c4057188 100644 Binary files a/docs/reference/PlotSignature-1.png and b/docs/reference/PlotSignature-1.png differ diff --git a/docs/reference/ReductiveHourglassTest.html b/docs/reference/ReductiveHourglassTest.html index 9302abf1..3c332b89 100644 --- a/docs/reference/ReductiveHourglassTest.html +++ b/docs/reference/ReductiveHourglassTest.html @@ -203,10 +203,10 @@

Examples#> [ Computing variances of permuted transcriptome signatures ... ] #> #> $p.value -#> [1] 1.53632e-08 +#> [1] 2.25838e-08 #> #> $std.dev -#> [1] 0.05620014 0.05490952 0.05323798 0.05270313 0.05260379 0.05389659 0.05483547 +#> [1] 0.05360555 0.05261389 0.05140542 0.04998371 0.04954327 0.05199497 0.05697204 #> #> $lillie.test #> [1] NA @@ -228,11 +228,11 @@

Examples modules = list(early = 1:2, mid = 3:5, late = 6:7), custom.perm.matrix = custom_perm_matrix) #> $p.value -#> [1] 6.915043e-11 +#> [1] 1.25244e-09 #> #> $std.dev #> Zygote Quadrant Globular Heart Torpedo Bent Mature -#> 0.05491497 0.05523176 0.05663849 0.05451059 0.05593353 0.05465526 0.05362093 +#> 0.05413093 0.05438733 0.05356645 0.05182419 0.04987965 0.04852364 0.05162828 #> #> $lillie.test #> [1] NA diff --git a/docs/reference/ReverseHourglassTest.html b/docs/reference/ReverseHourglassTest.html index b8f36f77..b5c69b12 100644 --- a/docs/reference/ReverseHourglassTest.html +++ b/docs/reference/ReverseHourglassTest.html @@ -204,7 +204,7 @@

Examples#> [1] 1 #> #> $std.dev -#> [1] 0.05425229 0.05269226 0.05260948 0.05143537 0.05058456 0.05246225 0.05647622 +#> [1] 0.05543708 0.05423243 0.05240352 0.05121000 0.05051120 0.05210643 0.05712842 #> #> $lillie.test #> [1] NA @@ -231,7 +231,7 @@

Examples#> #> $std.dev #> Zygote Quadrant Globular Heart Torpedo Bent Mature -#> 0.05322027 0.05442798 0.05494614 0.05256053 0.05193650 0.05330424 0.06116217 +#> 0.05179833 0.04870956 0.05064932 0.04873651 0.04951337 0.04931970 0.05426546 #> #> $lillie.test #> [1] NA diff --git a/docs/reference/Rplot004.png b/docs/reference/Rplot004.png index 701f7387..ecb93410 100644 Binary files a/docs/reference/Rplot004.png and b/docs/reference/Rplot004.png differ diff --git a/docs/reference/Rplot005.png b/docs/reference/Rplot005.png index bd2b7018..eec23e2e 100644 Binary files a/docs/reference/Rplot005.png and b/docs/reference/Rplot005.png differ diff --git a/docs/reference/ecScore.html b/docs/reference/ecScore.html index 2cd829b6..2568ffbc 100644 --- a/docs/reference/ecScore.html +++ b/docs/reference/ecScore.html @@ -144,10 +144,10 @@

Examples#> [=========================================] 100% #> [ Computing variances of permuted transcriptome signatures ... ] #> -#> 1 2 3 4 5 6 -#> 0.012846111 0.005238153 0.018602515 -0.006213682 0.018621329 -0.030893780 -#> 7 8 9 10 -#> -0.020161766 -0.002020718 -0.057479369 0.014861284 +#> 1 2 3 4 5 6 +#> -0.05681366 -0.02292928 -0.05911419 -0.04578738 -0.02225435 -0.02986441 +#> 7 8 9 10 +#> 0.01308162 -0.03649770 0.04031555 -0.04323258 # get warning if the expected pattern isn't followed ec_score <- ecScore(age_vals = TAIs,early = 1:2,mid = 3:5,late = 6:7,profile.warn=TRUE) diff --git a/docs/reference/lcScore.html b/docs/reference/lcScore.html index 931fec43..9a26264b 100644 --- a/docs/reference/lcScore.html +++ b/docs/reference/lcScore.html @@ -145,9 +145,9 @@

Examples#> [ Computing variances of permuted transcriptome signatures ... ] #> #> 1 2 3 4 5 6 -#> 0.059241551 0.004520677 -0.045917715 -0.019148566 -0.062604806 0.015434327 +#> -0.002674790 -0.051208266 -0.035363643 0.010016502 -0.052145721 0.038411729 #> 7 8 9 10 -#> -0.024608654 0.013527231 0.046585976 -0.038995459 +#> 0.009693024 0.006257383 0.014944768 0.024002544 # get warning if the expected pattern isn't followed lc_score <- lcScore(age_vals = TAIs,early = 1:2,mid = 3:5,late = 6:7,profile.warn=TRUE) diff --git a/docs/reference/pairScore.html b/docs/reference/pairScore.html index a6687ee7..3a156e4d 100644 --- a/docs/reference/pairScore.html +++ b/docs/reference/pairScore.html @@ -152,10 +152,10 @@

Examples#> [=========================================] 100% #> [ Computing variances of permuted transcriptome signatures ... ] #> -#> 1 2 3 4 5 6 -#> 0.009680977 0.008971025 0.017052149 -0.026414837 -0.022327672 -0.019311822 -#> 7 8 9 10 -#> 0.011474037 0.013211218 -0.026666449 -0.006799789 +#> 1 2 3 4 5 +#> 0.0094227374 0.0102870656 -0.0022944993 -0.0045077845 -0.0006856504 +#> 6 7 8 9 10 +#> 0.0113790720 -0.0178215366 -0.0011520399 0.0040270946 -0.0192558721 diff --git a/docs/reference/rcpp_boottei_parallel.html b/docs/reference/rcpp_boottei_parallel.html index f8d3eb23..fa053025 100644 --- a/docs/reference/rcpp_boottei_parallel.html +++ b/docs/reference/rcpp_boottei_parallel.html @@ -105,106 +105,106 @@

Examples## get permutations rcpp_boottei_parallel(spmat, ps, 100, 1) #> [,1] [,2] [,3] [,4] [,5] [,6] [,7] -#> [1,] 3.645446 3.651768 3.620124 3.608758 3.637346 3.633831 3.571788 -#> [2,] 3.499268 3.515727 3.496707 3.502995 3.507991 3.498898 3.499141 -#> [3,] 3.594584 3.607869 3.607120 3.609239 3.626909 3.647488 3.596822 -#> [4,] 3.632186 3.620927 3.629074 3.620988 3.586829 3.586578 3.631172 -#> [5,] 3.543491 3.552468 3.576062 3.585150 3.562097 3.522033 3.470307 -#> [6,] 3.533574 3.529404 3.537714 3.524943 3.531196 3.485195 3.507280 -#> [7,] 3.460863 3.472227 3.502521 3.480732 3.499382 3.503085 3.423462 -#> [8,] 3.571896 3.577924 3.549391 3.551568 3.589345 3.580445 3.558516 -#> [9,] 3.526182 3.501362 3.558576 3.571258 3.553030 3.591332 3.509143 -#> [10,] 3.548379 3.535988 3.553974 3.546582 3.565752 3.537104 3.494325 -#> [11,] 3.502389 3.502307 3.498237 3.482524 3.456566 3.455512 3.453558 -#> [12,] 3.562758 3.601485 3.594406 3.576632 3.556065 3.505993 3.472705 -#> [13,] 3.550936 3.578955 3.583660 3.573659 3.566466 3.530076 3.464357 -#> [14,] 3.523194 3.532284 3.507375 3.492898 3.481974 3.507776 3.555786 -#> [15,] 3.522301 3.514896 3.516837 3.507826 3.470673 3.428265 3.477633 -#> [16,] 3.429522 3.449266 3.439216 3.449021 3.468490 3.481611 3.484937 -#> [17,] 3.526747 3.525294 3.527188 3.533996 3.533393 3.497366 3.451280 -#> [18,] 3.645592 3.626552 3.629186 3.628807 3.617861 3.626939 3.534400 -#> [19,] 3.441509 3.454730 3.478969 3.455942 3.446786 3.485955 3.419221 -#> [20,] 3.449906 3.436778 3.442976 3.450777 3.436731 3.461867 3.474465 -#> [21,] 3.495212 3.484285 3.490290 3.501601 3.478363 3.491501 3.481568 -#> [22,] 3.505570 3.525965 3.516849 3.502507 3.470997 3.481730 3.439369 -#> [23,] 3.482492 3.468911 3.431116 3.449573 3.488423 3.509763 3.506970 -#> [24,] 3.499275 3.522743 3.567748 3.538810 3.528852 3.525148 3.463542 -#> [25,] 3.577320 3.587264 3.563869 3.554274 3.547897 3.533799 3.518891 -#> [26,] 3.535295 3.514061 3.486415 3.487491 3.498266 3.450997 3.504080 -#> [27,] 3.555566 3.566067 3.613960 3.619739 3.605961 3.588347 3.504527 -#> [28,] 3.536479 3.532588 3.512328 3.519097 3.504487 3.531890 3.598914 -#> [29,] 3.486869 3.488627 3.447366 3.451761 3.455831 3.436589 3.526933 -#> [30,] 3.486012 3.512842 3.519886 3.525156 3.517947 3.501520 3.494023 -#> [31,] 3.676238 3.661867 3.655696 3.636203 3.629330 3.615341 3.567290 -#> [32,] 3.560019 3.554434 3.553726 3.564909 3.566243 3.570990 3.519393 -#> [33,] 3.536122 3.542346 3.495193 3.504444 3.509680 3.511550 3.548710 -#> [34,] 3.516141 3.517627 3.519985 3.532826 3.530599 3.524582 3.562375 -#> [35,] 3.610742 3.611470 3.602182 3.599834 3.590742 3.568898 3.635098 -#> [36,] 3.599611 3.595344 3.545331 3.521092 3.511664 3.528591 3.645715 -#> [37,] 3.524590 3.491591 3.497356 3.492014 3.498298 3.493367 3.521059 -#> [38,] 3.589652 3.602260 3.541010 3.538722 3.561140 3.570521 3.539542 -#> [39,] 3.450576 3.457478 3.462637 3.478455 3.470134 3.484990 3.510191 -#> [40,] 3.491403 3.515123 3.504112 3.519487 3.527230 3.529208 3.434870 -#> [41,] 3.593598 3.582949 3.559509 3.547191 3.515431 3.508664 3.531339 -#> [42,] 3.454080 3.454547 3.469972 3.474528 3.475927 3.475123 3.442066 -#> [43,] 3.510217 3.519449 3.511687 3.509526 3.537192 3.502677 3.479539 -#> [44,] 3.441379 3.453080 3.483341 3.486791 3.494492 3.485325 3.445336 -#> [45,] 3.520706 3.516193 3.525002 3.535245 3.543963 3.530468 3.515073 -#> [46,] 3.418429 3.400631 3.425715 3.426269 3.449621 3.495444 3.522857 -#> [47,] 3.477020 3.470499 3.498758 3.509182 3.503547 3.514217 3.534582 -#> [48,] 3.614044 3.602531 3.550431 3.548250 3.540275 3.590380 3.572735 -#> [49,] 3.441267 3.438989 3.422310 3.424470 3.432621 3.435423 3.415842 -#> [50,] 3.567303 3.582539 3.595035 3.595780 3.580009 3.619646 3.579885 -#> [51,] 3.462449 3.485133 3.482384 3.480607 3.501607 3.513321 3.466450 -#> [52,] 3.557778 3.559668 3.535102 3.541450 3.557253 3.561062 3.489085 -#> [53,] 3.490262 3.502281 3.499272 3.503953 3.496561 3.515999 3.530291 -#> [54,] 3.617220 3.608119 3.592863 3.593035 3.562296 3.543559 3.544767 -#> [55,] 3.493181 3.466666 3.481847 3.467948 3.455848 3.493233 3.551186 -#> [56,] 3.533322 3.523451 3.521306 3.530813 3.522295 3.550346 3.478215 -#> [57,] 3.655644 3.646393 3.601653 3.635158 3.658428 3.663336 3.620346 -#> [58,] 3.481044 3.467874 3.506616 3.525207 3.546892 3.531147 3.490811 -#> [59,] 3.527161 3.535428 3.537136 3.514556 3.516811 3.541548 3.562811 -#> [60,] 3.572474 3.587010 3.620845 3.596092 3.578512 3.568401 3.561188 -#> [61,] 3.538232 3.572055 3.574998 3.567350 3.556701 3.499014 3.509944 -#> [62,] 3.552016 3.535769 3.503354 3.486769 3.459523 3.455173 3.478853 -#> [63,] 3.448474 3.433673 3.462286 3.484433 3.495386 3.511972 3.491216 -#> [64,] 3.552500 3.527062 3.529135 3.535664 3.554273 3.530658 3.550121 -#> [65,] 3.428225 3.439597 3.446539 3.468003 3.494723 3.473914 3.505546 -#> [66,] 3.468444 3.461832 3.462503 3.457843 3.466520 3.435897 3.475808 -#> [67,] 3.519352 3.513604 3.558361 3.571447 3.572705 3.530051 3.524838 -#> [68,] 3.498583 3.539275 3.506780 3.496986 3.511239 3.508571 3.411983 -#> [69,] 3.551930 3.545460 3.504974 3.511980 3.523216 3.538530 3.564309 -#> [70,] 3.441760 3.434375 3.450638 3.445432 3.462001 3.515143 3.560481 -#> [71,] 3.427147 3.437101 3.440209 3.426109 3.443521 3.452792 3.528299 -#> [72,] 3.630735 3.609708 3.593156 3.589666 3.573685 3.545902 3.601362 -#> [73,] 3.527562 3.532522 3.572743 3.548216 3.529802 3.579851 3.511185 -#> [74,] 3.495636 3.477183 3.494507 3.494594 3.490886 3.458137 3.511664 -#> [75,] 3.538741 3.561163 3.561751 3.565442 3.549192 3.546592 3.512101 -#> [76,] 3.483484 3.491073 3.465041 3.482243 3.477392 3.495157 3.558825 -#> [77,] 3.459694 3.457583 3.477406 3.474231 3.464080 3.482597 3.497245 -#> [78,] 3.573974 3.579897 3.576827 3.575533 3.550819 3.547339 3.569561 -#> [79,] 3.548139 3.562031 3.560744 3.565260 3.572391 3.532333 3.501065 -#> [80,] 3.505597 3.502591 3.468336 3.464719 3.477044 3.488752 3.419755 -#> [81,] 3.596064 3.558968 3.505616 3.503557 3.525487 3.500294 3.487716 -#> [82,] 3.501704 3.495716 3.529450 3.528173 3.509569 3.513799 3.503477 -#> [83,] 3.525786 3.526887 3.496330 3.497191 3.503050 3.492553 3.464577 -#> [84,] 3.483831 3.490026 3.487445 3.491418 3.505683 3.484523 3.452159 -#> [85,] 3.567584 3.570437 3.581038 3.591426 3.588990 3.581742 3.547324 -#> [86,] 3.589437 3.551511 3.519456 3.531001 3.536938 3.564402 3.500823 -#> [87,] 3.498695 3.524810 3.477090 3.479568 3.475929 3.486445 3.470598 -#> [88,] 3.550673 3.553670 3.546858 3.529219 3.544539 3.538689 3.497476 -#> [89,] 3.537869 3.541829 3.522230 3.522641 3.514686 3.486155 3.508998 -#> [90,] 3.638532 3.594089 3.582845 3.570046 3.570448 3.528380 3.583235 -#> [91,] 3.564175 3.550971 3.572009 3.566859 3.553373 3.559304 3.551182 -#> [92,] 3.534428 3.535832 3.523268 3.531259 3.546027 3.561872 3.539346 -#> [93,] 3.616726 3.630904 3.631727 3.623956 3.608749 3.598033 3.557954 -#> [94,] 3.549996 3.549212 3.549319 3.565265 3.562696 3.577095 3.472850 -#> [95,] 3.471457 3.498334 3.434733 3.428171 3.438463 3.414542 3.489118 -#> [96,] 3.611630 3.630401 3.633882 3.620879 3.628869 3.623517 3.578291 -#> [97,] 3.607310 3.607964 3.620935 3.614474 3.585991 3.589948 3.487840 -#> [98,] 3.545539 3.542808 3.540429 3.528892 3.513406 3.541346 3.560043 -#> [99,] 3.529872 3.548568 3.529503 3.549718 3.553388 3.573103 3.528795 -#> [100,] 3.517955 3.537782 3.549419 3.530608 3.511852 3.485675 3.475535 +#> [1,] 3.485065 3.500548 3.523695 3.513210 3.518488 3.536586 3.505940 +#> [2,] 3.450433 3.456191 3.459595 3.454438 3.437657 3.480261 3.446454 +#> [3,] 3.534142 3.534034 3.548797 3.550227 3.552603 3.554733 3.543461 +#> [4,] 3.536101 3.533254 3.513446 3.524565 3.517994 3.557340 3.536497 +#> [5,] 3.508651 3.518543 3.504199 3.540596 3.549397 3.552073 3.533044 +#> [6,] 3.558369 3.546603 3.531305 3.539373 3.552368 3.537420 3.447822 +#> [7,] 3.543463 3.539038 3.590499 3.559016 3.559545 3.572261 3.511127 +#> [8,] 3.450836 3.446856 3.440262 3.426459 3.421018 3.428527 3.506781 +#> [9,] 3.578727 3.576340 3.555575 3.558988 3.563689 3.556704 3.520515 +#> [10,] 3.569295 3.550210 3.514895 3.517598 3.539634 3.556474 3.541706 +#> [11,] 3.507978 3.525927 3.505090 3.510493 3.514511 3.457838 3.438887 +#> [12,] 3.521401 3.520601 3.483204 3.490950 3.493933 3.505095 3.493845 +#> [13,] 3.473782 3.475366 3.498906 3.519962 3.553547 3.536754 3.465897 +#> [14,] 3.495975 3.494553 3.522059 3.510915 3.491990 3.486616 3.478194 +#> [15,] 3.439812 3.418574 3.448768 3.460049 3.480807 3.480883 3.514986 +#> [16,] 3.537144 3.511894 3.513971 3.531042 3.530224 3.571280 3.514913 +#> [17,] 3.457717 3.446786 3.458322 3.483095 3.498497 3.517633 3.543894 +#> [18,] 3.455554 3.447912 3.470496 3.462529 3.452899 3.494469 3.529821 +#> [19,] 3.454914 3.469083 3.484552 3.491104 3.454225 3.452984 3.536013 +#> [20,] 3.531732 3.547022 3.584379 3.581904 3.559391 3.546595 3.488368 +#> [21,] 3.615637 3.569811 3.541041 3.537724 3.536676 3.523952 3.601040 +#> [22,] 3.597227 3.607700 3.546949 3.550177 3.546260 3.548395 3.609069 +#> [23,] 3.425372 3.423062 3.436914 3.427645 3.411498 3.442659 3.427586 +#> [24,] 3.650197 3.627605 3.645724 3.651598 3.631273 3.660346 3.586527 +#> [25,] 3.455602 3.472342 3.464089 3.464316 3.456929 3.437475 3.457195 +#> [26,] 3.473583 3.495222 3.518632 3.498197 3.507131 3.513286 3.480314 +#> [27,] 3.437233 3.422187 3.386025 3.382569 3.382434 3.352647 3.359759 +#> [28,] 3.444755 3.449036 3.452384 3.452610 3.460757 3.540128 3.569811 +#> [29,] 3.523749 3.533837 3.565080 3.535494 3.508461 3.494340 3.505571 +#> [30,] 3.454379 3.436866 3.458962 3.470221 3.469511 3.493466 3.563412 +#> [31,] 3.506344 3.548313 3.525084 3.498452 3.496660 3.475995 3.471965 +#> [32,] 3.402017 3.382286 3.395447 3.411462 3.426799 3.440598 3.465712 +#> [33,] 3.535107 3.495623 3.500223 3.503116 3.491167 3.515757 3.501457 +#> [34,] 3.545782 3.527204 3.533275 3.541225 3.546722 3.557578 3.553596 +#> [35,] 3.469558 3.466290 3.471240 3.470325 3.492974 3.496060 3.532502 +#> [36,] 3.439787 3.463764 3.465573 3.465791 3.460079 3.451410 3.435245 +#> [37,] 3.592363 3.609051 3.577068 3.579376 3.617450 3.614728 3.532373 +#> [38,] 3.442015 3.444140 3.451190 3.463300 3.487857 3.502834 3.448861 +#> [39,] 3.401813 3.404089 3.385019 3.384933 3.389660 3.382516 3.454582 +#> [40,] 3.477130 3.477000 3.519266 3.505307 3.484813 3.491189 3.430860 +#> [41,] 3.537599 3.509745 3.520441 3.525631 3.517660 3.519152 3.469915 +#> [42,] 3.502789 3.519938 3.540519 3.533407 3.516727 3.519334 3.549323 +#> [43,] 3.522362 3.512204 3.508342 3.514974 3.505605 3.498172 3.548005 +#> [44,] 3.419094 3.449131 3.452715 3.461094 3.435231 3.447928 3.450133 +#> [45,] 3.452982 3.447575 3.461756 3.476219 3.465800 3.461766 3.529605 +#> [46,] 3.530784 3.535163 3.539798 3.526743 3.563731 3.574288 3.470460 +#> [47,] 3.546035 3.549349 3.571674 3.580288 3.578451 3.585698 3.485385 +#> [48,] 3.443652 3.436089 3.509961 3.489863 3.488037 3.495170 3.615633 +#> [49,] 3.513142 3.508967 3.479432 3.494024 3.516489 3.536835 3.516466 +#> [50,] 3.460031 3.452279 3.485304 3.500410 3.497081 3.499209 3.520941 +#> [51,] 3.563983 3.536105 3.554214 3.568197 3.566877 3.592618 3.496183 +#> [52,] 3.513813 3.507023 3.525951 3.522930 3.510582 3.493686 3.520616 +#> [53,] 3.503404 3.499828 3.519098 3.508959 3.513663 3.538344 3.502143 +#> [54,] 3.487687 3.477153 3.448287 3.456206 3.447132 3.423523 3.506686 +#> [55,] 3.483835 3.472546 3.479919 3.502852 3.507656 3.497641 3.493258 +#> [56,] 3.544301 3.551873 3.548581 3.549295 3.528093 3.508926 3.607338 +#> [57,] 3.466352 3.488251 3.500020 3.497063 3.491411 3.472133 3.467444 +#> [58,] 3.437950 3.467349 3.476617 3.462564 3.470358 3.439007 3.440881 +#> [59,] 3.513794 3.513371 3.544162 3.563606 3.584393 3.616472 3.538583 +#> [60,] 3.395633 3.407416 3.447300 3.448800 3.473976 3.456841 3.413727 +#> [61,] 3.521370 3.527778 3.510665 3.516791 3.533563 3.530460 3.590304 +#> [62,] 3.435900 3.439003 3.462339 3.469575 3.476532 3.446932 3.489802 +#> [63,] 3.460701 3.478696 3.503909 3.497206 3.481968 3.477703 3.434062 +#> [64,] 3.564472 3.552119 3.529966 3.540339 3.570445 3.574167 3.557788 +#> [65,] 3.665218 3.641421 3.645443 3.643485 3.614275 3.566026 3.482296 +#> [66,] 3.549158 3.545901 3.551212 3.539276 3.535103 3.511192 3.516926 +#> [67,] 3.517841 3.520738 3.498495 3.488284 3.496990 3.462758 3.521473 +#> [68,] 3.498910 3.525966 3.510460 3.525591 3.540729 3.531492 3.526949 +#> [69,] 3.491820 3.502947 3.475760 3.492799 3.488865 3.470330 3.469068 +#> [70,] 3.509015 3.516993 3.534962 3.522264 3.537645 3.544526 3.505181 +#> [71,] 3.554149 3.534371 3.569309 3.579121 3.565091 3.533619 3.607779 +#> [72,] 3.513397 3.515119 3.497527 3.500893 3.509220 3.500741 3.483599 +#> [73,] 3.492957 3.487800 3.515289 3.500386 3.487112 3.542580 3.513061 +#> [74,] 3.535984 3.531182 3.543331 3.561751 3.556978 3.577525 3.533953 +#> [75,] 3.499591 3.488090 3.483900 3.493926 3.500501 3.531731 3.499708 +#> [76,] 3.396401 3.375157 3.406631 3.402845 3.432941 3.487822 3.540748 +#> [77,] 3.600023 3.620185 3.581089 3.578346 3.568079 3.571804 3.594126 +#> [78,] 3.551101 3.541187 3.542033 3.552960 3.553655 3.534751 3.534068 +#> [79,] 3.546061 3.538312 3.539328 3.518645 3.492269 3.497597 3.471863 +#> [80,] 3.531713 3.527833 3.513508 3.517726 3.513249 3.563071 3.475138 +#> [81,] 3.531367 3.530205 3.542648 3.525238 3.495363 3.525207 3.517780 +#> [82,] 3.474114 3.503793 3.512647 3.493299 3.485726 3.477402 3.437767 +#> [83,] 3.385961 3.407081 3.454082 3.452492 3.447172 3.448199 3.397078 +#> [84,] 3.486110 3.508970 3.548654 3.524254 3.521531 3.543952 3.545779 +#> [85,] 3.411092 3.428893 3.454040 3.465658 3.459026 3.466486 3.456636 +#> [86,] 3.604336 3.577133 3.557696 3.542845 3.538461 3.568805 3.632620 +#> [87,] 3.531275 3.516473 3.529812 3.506306 3.501487 3.521841 3.518759 +#> [88,] 3.513079 3.499215 3.551186 3.546389 3.545731 3.561184 3.566536 +#> [89,] 3.538608 3.544846 3.564917 3.525106 3.532875 3.558406 3.547311 +#> [90,] 3.599651 3.616330 3.558358 3.568562 3.552894 3.532971 3.514893 +#> [91,] 3.470882 3.482628 3.508899 3.520075 3.526829 3.571617 3.521389 +#> [92,] 3.545927 3.548086 3.535224 3.559724 3.580593 3.612495 3.586157 +#> [93,] 3.452810 3.452795 3.445646 3.451387 3.453661 3.475940 3.499907 +#> [94,] 3.425632 3.417743 3.455482 3.465863 3.461739 3.451724 3.546973 +#> [95,] 3.463940 3.484618 3.490187 3.484471 3.486134 3.513660 3.480674 +#> [96,] 3.469190 3.474024 3.526860 3.511217 3.509974 3.468255 3.456542 +#> [97,] 3.463207 3.465013 3.520461 3.526601 3.555582 3.562655 3.526191 +#> [98,] 3.500710 3.500088 3.473551 3.478605 3.479792 3.494792 3.528789 +#> [99,] 3.421468 3.421371 3.434253 3.445164 3.439634 3.424682 3.453375 +#> [100,] 3.462148 3.449993 3.431413 3.435270 3.443673 3.455121 3.479393