==19863== Memcheck, a memory error detector ==19863== Copyright (C) 2002-2017, and GNU GPL'd, by Julian Seward et al. ==19863== Using Valgrind-3.13.0 and LibVEX; rerun with -h for copyright info ==19863== Command: /data/blackswan/ripley/R/R-devel-vg/bin/exec/R --vanilla ==19863== R Under development (unstable) (2018-05-10 r74708) -- "Unsuffered Consequences" Copyright (C) 2018 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > pkgname <- "pvrank" > source(file.path(R.home("share"), "R", "examples-header.R")) > options(warn = 1) > library('pvrank') Loading required package: EnvStats Attaching package: ‘EnvStats’ ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x40300CB: ??? ==19863== by 0x15BBCF67: ??? ==19863== by 0x15BBCF67: ??? ==19863== by 0x15BBCF7A: ??? ==19863== by 0x1FFEFF63BF: ??? ==19863== The following objects are masked from ‘package:stats’: predict, predict.lm ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x403FC03: ??? ==19863== by 0x107FB7E7: ??? ==19863== by 0x107FB7E7: ??? ==19863== by 0x107FB7F4: ??? ==19863== by 0x1FFEFF63BF: ??? ==19863== The following object is masked from ‘package:base’: print.default Loading required package: Rmpfr Loading required package: gmp Attaching package: ‘gmp’ ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x403FC03: ??? ==19863== by 0x12147527: ??? ==19863== by 0x12147527: ??? ==19863== by 0x12147550: ??? ==19863== by 0x1FFEFF465F: ??? ==19863== The following objects are masked from ‘package:base’: %*%, apply, crossprod, matrix, tcrossprod C code of R package 'Rmpfr': GMP using 64 bits per limb Attaching package: ‘Rmpfr’ ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x403FC03: ??? ==19863== by 0x135D8B57: ??? ==19863== by 0x135D8B57: ??? ==19863== by 0x135D8B72: ??? ==19863== by 0x1FFEFF63BF: ??? ==19863== The following objects are masked from ‘package:stats’: dbinom, dnorm, dpois, pnorm The following objects are masked from ‘package:base’: cbind, pmax, pmin, rbind > > base::assign(".oldSearch", base::search(), pos = 'CheckExEnv') > base::assign(".old_wd", base::getwd(), pos = 'CheckExEnv') > cleanEx() > nameEx("Atar") > ### * Atar > > flush(stderr()); flush(stdout()) > > ### Name: Atar > ### Title: Total backward and forward output linkage indices > ### Aliases: Atar > ### Keywords: datasets > > ### ** Examples > > data(Atar) > head(Atar) TBL TFL 1_2 0.9187 1.0374 5 0.8083 0.6228 10_14 0.9582 0.7100 15_37 1.2192 2.5258 40_41 1.1290 1.3010 45 1.1406 0.6703 > > > > cleanEx() > nameEx("Berk") > ### * Berk > > flush(stderr()); flush(stdout()) > > ### Name: Berk > ### Title: Effect of outliers > ### Aliases: Berk > ### Keywords: datasets > > ### ** Examples > > data(Berk) > head(Berk) Births Deaths 1 142 228 2 173 253 3 130 230 4 122 242 5 120 231 6 111 213 > > > > cleanEx() > nameEx("BlaAlt") > ### * BlaAlt > > flush(stderr()); flush(stdout()) > > ### Name: BlaAlt > ### Title: Pairs of measurements of forced expiratory volume > ### Aliases: BlaAlt > ### Keywords: datasets > > ### ** Examples > > data(BlaAlt) > head(BlaAlt) Fev1 Fev2 1 1.19 1.37 2 1.33 1.32 3 1.35 1.40 4 1.36 1.25 5 1.38 1.29 6 1.38 1.37 > > > > cleanEx() > nameEx("Dalyww") > ### * Dalyww > > flush(stderr()); flush(stdout()) > > ### Name: Dalyww > ### Title: Dark contrasts: the paradox of high rates of suicide in happy > ### places > ### Aliases: Dalyww > ### Keywords: datasets > > ### ** Examples > > data(Dalyww) > head(Dalyww) SRU ASHR RCLS ACLS AL 11.98 0.134 0.0000000 0.0000000 AK 23.60 1.955 0.0185221 -0.0130092 AZ 15.32 0.702 0.0493716 0.0017149 AR 13.14 0.522 0.0099526 -0.0173219 CA 9.40 1.634 -0.0104240 -0.0758851 CO 17.33 0.742 0.0594714 -0.0270579 > > > > cleanEx() > nameEx("DietFish") > ### * DietFish > > flush(stderr()); flush(stdout()) > > ### Name: DietFish > ### Title: Comparison of total diets of banded killifish and mummichog > ### Aliases: DietFish > ### Keywords: datasets > > ### ** Examples > > data(DietFish) > head(DietFish) occ_D occ_H N_virens 4.8 0.0 Cladocera 3.4 7.7 Canthocampidae 2.0 4.6 Eurytemora_affinis 15.2 1.5 Orthocyclops 11.8 13.1 Ostracoda 16.6 27.1 > > > > cleanEx() > nameEx("Dizytwin") > ### * Dizytwin > > flush(stderr()); flush(stdout()) > > ### Name: Dizytwin > ### Title: Latitude and dizygotic twinning rates > ### Aliases: Dizytwin > ### Keywords: datasets > > ### ** Examples > > data(Dizytwin) > head(Dizytwin) Latitude DZT_Rate Portugal 40 6.5 Greece 40 8.8 Spain 41 5.9 Bulgaria 42 7.0 Italy 44 8.6 France 47 7.1 > > > > cleanEx() > nameEx("Femurs") > ### * Femurs > > flush(stderr()); flush(stdout()) > > ### Name: Femurs > ### Title: Radiographic-caudolateral curvilinear osteophyte > ### Aliases: Femurs > ### Keywords: datasets > > ### ** Examples > > data(Femurs) > head(Femurs) CCO Age BW Gender 1 5.67 52 17.4 1 2 6.15 120 16.8 1 3 7.41 120 16.8 1 4 8.25 68 14.3 1 5 8.32 84 18.5 1 6 9.28 120 19.4 1 > > > > cleanEx() > nameEx("Franzen") > ### * Franzen > > flush(stderr()); flush(stdout()) > > ### Name: Franzen > ### Title: Environmental attitudes data > ### Aliases: Franzen > ### Keywords: datasets > > ### ** Examples > > data(Franzen) > head(Franzen) MECISSP PPP J 27.4 26.5 Fin 27.3 24.6 Ch 27.1 25.0 Cdn 26.5 27.3 S 26.3 23.8 Nl 26.2 26.2 > > > > cleanEx() > nameEx("FrigateB") > ### * FrigateB > > flush(stderr()); flush(stdout()) > > ### Name: FrigateB > ### Title: Magnificent frigatebird > ### Aliases: FrigateB > ### Keywords: datasets > > ### ** Examples > > data(FrigateB) > head(FrigateB) Vol Frq 1 1760 529 2 2040 566 3 2440 473 4 2550 461 5 2730 465 6 2740 532 > > > > cleanEx() > nameEx("GGFR") > ### * GGFR > > flush(stderr()); flush(stdout()) > > ### Name: VGGFR > ### Title: Vianelli density > ### Aliases: VGGFR > > ### ** Examples > > > # Density curve of a VGGFR model > VGGFR(2, 12, add=FALSE, lwd=2, lty=5, col="darkgreen", ylim=c(0,2), Main="", np=201) Parameters of VGGFR: 2 12 Std. Deviation: 0.19245 Excess Kurtosis: -0.206897 Ordinate at the mode: 2.014753 $st.dev [1] 0.19245 $kurt [1] -0.206897 $oam [1] 2.014753 > ##### > # > a<-ranktes(0.5, 28, "r4", "vg",FALSE, "two", FALSE) > b<-VGGFR(a$Lambda, add = FALSE, lwd = 2, lty = 5, col = "blue", ylim=c(0,2.5),np = 201) Parameters of VGGFR: 1.5 7.55238588 Std. Deviation: 0.193182 Excess Kurtosis: 0.150978 Ordinate at the mode: 2.286617 > ##### > # > # Lambert's semicircular distribution of errors (1760,1765). > # Given a probability distribution, the value with the higher probability density is > # deemed to be more probable than the value with the lower probability density. > # > VGGFR(2,0.5,col="red",ylim=c(0,0.75),Main="Lambert's distribution of errors") Parameters of VGGFR: 2 0.5 Std. Deviation: 0.5 Excess Kurtosis: -1 Ordinate at the mode: 0.63662 $st.dev [1] 0.5 $kurt [1] -1 $oam [1] 0.63662 > # > ##### > # Pearson type II used as an approximation to the null distribution of the Fisher-Yates > # rank correlation. Fieller, E. C. and Pearson, E. S. (1961). Tests for rank correlation > # coefficients: II. Biometrika, 48, 29-40. > n<-10 > VGGFR(2, (n-4)/2, add=FALSE, lwd=2, lty=5, col="magenta2", ylim=c(0,1.1), Main="", np=201) Parameters of VGGFR: 2 3 Std. Deviation: 0.333333 Excess Kurtosis: -0.545455 Ordinate at the mode: 1.09375 $st.dev [1] 0.333333 $kurt [1] -0.545455 $oam [1] 1.09375 > abline(h=0);abline(v=0,lty=2,lwd=2,col="pink2") > ##### > # > # Save and use the results > res<-VGGFR(1.5,5.5,add = FALSE, lwd = 2, lty = 1, col = "blue", ylim=c(0,2.5),np = 201) Parameters of VGGFR: 1.5 5.5 Std. Deviation: 0.227761 Excess Kurtosis: 0.010096 Ordinate at the mode: 1.896675 > res$kurt-res$oam/res$st.dev [1] -8.31738324 > ##### > # > # A family of symmetrical beta densities > VGGFR(2,1,col="black",ylim=c(0,1.4),Main="Symmetrical beta densities") Parameters of VGGFR: 2 1 Std. Deviation: 0.447214 Excess Kurtosis: -0.857143 Ordinate at the mode: 0.75 $st.dev [1] 0.447214 $kurt [1] -0.857143 $oam [1] 0.75 > La<-seq(1,6,0.5) > for (L1 in La){VGGFR(2,L1,add=TRUE, lwd = 1, lty = 1, col=gray(L1/6))} Parameters of VGGFR: 2 1 Std. Deviation: 0.447214 Excess Kurtosis: -0.857143 Ordinate at the mode: 0.75 Parameters of VGGFR: 2 1.5 Std. Deviation: 0.408248 Excess Kurtosis: -0.75 Ordinate at the mode: 0.848826 Parameters of VGGFR: 2 2 Std. Deviation: 0.377964 Excess Kurtosis: -0.666667 Ordinate at the mode: 0.9375 Parameters of VGGFR: 2 2.5 Std. Deviation: 0.353553 Excess Kurtosis: -0.6 Ordinate at the mode: 1.018592 Parameters of VGGFR: 2 3 Std. Deviation: 0.333333 Excess Kurtosis: -0.545455 Ordinate at the mode: 1.09375 Parameters of VGGFR: 2 3.5 Std. Deviation: 0.316228 Excess Kurtosis: -0.5 Ordinate at the mode: 1.164105 Parameters of VGGFR: 2 4 Std. Deviation: 0.301511 Excess Kurtosis: -0.461538 Ordinate at the mode: 1.230469 Parameters of VGGFR: 2 4.5 Std. Deviation: 0.288675 Excess Kurtosis: -0.428571 Ordinate at the mode: 1.29345 Parameters of VGGFR: 2 5 Std. Deviation: 0.27735 Excess Kurtosis: -0.4 Ordinate at the mode: 1.353516 Parameters of VGGFR: 2 5.5 Std. Deviation: 0.267261 Excess Kurtosis: -0.375 Ordinate at the mode: 1.411036 Parameters of VGGFR: 2 6 Std. Deviation: 0.258199 Excess Kurtosis: -0.352941 Ordinate at the mode: 1.466309 > ##### > # > # A family of GGFR curves > VGGFR(1,2, lwd = 1, lty = 1,col="black",ylim=c(0,5)) Parameters of VGGFR: 1 2 Std. Deviation: 0.316228 Excess Kurtosis: -0.142857 Ordinate at the mode: 1.5 $st.dev [1] 0.316228 $kurt [1] -0.142857 $oam [1] 1.5 > La<-seq(1,6,0.5);Lg<-seq(0,1,1/12) > for (L1 in La){ + c2<-gray(Lg, alpha= 2/6) + for (L2 in seq(1,12,1)){ + VGGFR(L1,L2,add=TRUE, lwd = 1, lty = 1, col=c2[L2]) + }} Parameters of VGGFR: 1 1 Std. Deviation: 0.408248 Excess Kurtosis: -0.6 Ordinate at the mode: 1 Parameters of VGGFR: 1 2 Std. Deviation: 0.316228 Excess Kurtosis: -0.142857 Ordinate at the mode: 1.5 Parameters of VGGFR: 1 3 Std. Deviation: 0.258199 Excess Kurtosis: 0.214286 Ordinate at the mode: 2 Parameters of VGGFR: 1 4 Std. Deviation: 0.218218 Excess Kurtosis: 0.5 Ordinate at the mode: 2.5 Parameters of VGGFR: 1 5 Std. Deviation: 0.188982 Excess Kurtosis: 0.733333 Ordinate at the mode: 3 Parameters of VGGFR: 1 6 Std. Deviation: 0.166667 Excess Kurtosis: 0.927273 Ordinate at the mode: 3.5 Parameters of VGGFR: 1 7 Std. Deviation: 0.149071 Excess Kurtosis: 1.090909 Ordinate at the mode: 4 Parameters of VGGFR: 1 8 Std. Deviation: 0.13484 Excess Kurtosis: 1.230769 Ordinate at the mode: 4.5 Parameters of VGGFR: 1 9 Std. Deviation: 0.123091 Excess Kurtosis: 1.351648 Ordinate at the mode: 5 Parameters of VGGFR: 1 10 Std. Deviation: 0.113228 Excess Kurtosis: 1.457143 Ordinate at the mode: 5.5 Parameters of VGGFR: 1 11 Std. Deviation: 0.104828 Excess Kurtosis: 1.55 Ordinate at the mode: 6 Parameters of VGGFR: 1 12 Std. Deviation: 0.09759 Excess Kurtosis: 1.632353 Ordinate at the mode: 6.5 Parameters of VGGFR: 1.5 1 Std. Deviation: 0.430331 Excess Kurtosis: -0.756923 Ordinate at the mode: 0.833333 Parameters of VGGFR: 1.5 2 Std. Deviation: 0.351364 Excess Kurtosis: -0.476538 Ordinate at the mode: 1.111111 Parameters of VGGFR: 1.5 3 Std. Deviation: 0.30089 Excess Kurtosis: -0.283355 Ordinate at the mode: 1.358025 Parameters of VGGFR: 1.5 4 Std. Deviation: 0.26536 Excess Kurtosis: -0.14223 Ordinate at the mode: 1.584362 Parameters of VGGFR: 1.5 5 Std. Deviation: 0.238754 Excess Kurtosis: -0.034644 Ordinate at the mode: 1.79561 Parameters of VGGFR: 1.5 6 Std. Deviation: 0.217952 Excess Kurtosis: 0.050081 Ordinate at the mode: 1.995123 Parameters of VGGFR: 1.5 7 Std. Deviation: 0.20116 Excess Kurtosis: 0.118526 Ordinate at the mode: 2.185134 Parameters of VGGFR: 1.5 8 Std. Deviation: 0.18727 Excess Kurtosis: 0.17497 Ordinate at the mode: 2.367229 Parameters of VGGFR: 1.5 9 Std. Deviation: 0.175554 Excess Kurtosis: 0.222313 Ordinate at the mode: 2.542579 Parameters of VGGFR: 1.5 10 Std. Deviation: 0.165514 Excess Kurtosis: 0.262592 Ordinate at the mode: 2.712084 Parameters of VGGFR: 1.5 11 Std. Deviation: 0.156796 Excess Kurtosis: 0.297277 Ordinate at the mode: 2.876453 Parameters of VGGFR: 1.5 12 Std. Deviation: 0.149143 Excess Kurtosis: 0.327458 Ordinate at the mode: 3.036256 Parameters of VGGFR: 2 1 Std. Deviation: 0.447214 Excess Kurtosis: -0.857143 Ordinate at the mode: 0.75 Parameters of VGGFR: 2 2 Std. Deviation: 0.377964 Excess Kurtosis: -0.666667 Ordinate at the mode: 0.9375 Parameters of VGGFR: 2 3 Std. Deviation: 0.333333 Excess Kurtosis: -0.545455 Ordinate at the mode: 1.09375 Parameters of VGGFR: 2 4 Std. Deviation: 0.301511 Excess Kurtosis: -0.461538 Ordinate at the mode: 1.230469 Parameters of VGGFR: 2 5 Std. Deviation: 0.27735 Excess Kurtosis: -0.4 Ordinate at the mode: 1.353516 Parameters of VGGFR: 2 6 Std. Deviation: 0.258199 Excess Kurtosis: -0.352941 Ordinate at the mode: 1.466309 Parameters of VGGFR: 2 7 Std. Deviation: 0.242536 Excess Kurtosis: -0.315789 Ordinate at the mode: 1.571045 Parameters of VGGFR: 2 8 Std. Deviation: 0.229416 Excess Kurtosis: -0.285714 Ordinate at the mode: 1.669235 Parameters of VGGFR: 2 9 Std. Deviation: 0.218218 Excess Kurtosis: -0.26087 Ordinate at the mode: 1.761971 Parameters of VGGFR: 2 10 Std. Deviation: 0.208514 Excess Kurtosis: -0.24 Ordinate at the mode: 1.850069 Parameters of VGGFR: 2 11 Std. Deviation: 0.2 Excess Kurtosis: -0.222222 Ordinate at the mode: 1.934163 Parameters of VGGFR: 2 12 Std. Deviation: 0.19245 Excess Kurtosis: -0.206897 Ordinate at the mode: 2.014753 Parameters of VGGFR: 2.5 1 Std. Deviation: 0.460566 Excess Kurtosis: -0.925714 Ordinate at the mode: 0.7 Parameters of VGGFR: 2.5 2 Std. Deviation: 0.398862 Excess Kurtosis: -0.787429 Ordinate at the mode: 0.84 Parameters of VGGFR: 2.5 3 Std. Deviation: 0.35887 Excess Kurtosis: -0.704132 Ordinate at the mode: 0.952 Parameters of VGGFR: 2.5 4 Std. Deviation: 0.330113 Excess Kurtosis: -0.648474 Ordinate at the mode: 1.0472 Parameters of VGGFR: 2.5 5 Std. Deviation: 0.30808 Excess Kurtosis: -0.60866 Ordinate at the mode: 1.130976 Parameters of VGGFR: 2.5 6 Std. Deviation: 0.29046 Excess Kurtosis: -0.578768 Ordinate at the mode: 1.206374 Parameters of VGGFR: 2.5 7 Std. Deviation: 0.275928 Excess Kurtosis: -0.555501 Ordinate at the mode: 1.27531 Parameters of VGGFR: 2.5 8 Std. Deviation: 0.263658 Excess Kurtosis: -0.536877 Ordinate at the mode: 1.339076 Parameters of VGGFR: 2.5 9 Std. Deviation: 0.253108 Excess Kurtosis: -0.521631 Ordinate at the mode: 1.39859 Parameters of VGGFR: 2.5 10 Std. Deviation: 0.243901 Excess Kurtosis: -0.508921 Ordinate at the mode: 1.454534 Parameters of VGGFR: 2.5 11 Std. Deviation: 0.235768 Excess Kurtosis: -0.498164 Ordinate at the mode: 1.507426 Parameters of VGGFR: 2.5 12 Std. Deviation: 0.228512 Excess Kurtosis: -0.48894 Ordinate at the mode: 1.557673 Parameters of VGGFR: 3 1 Std. Deviation: 0.471405 Excess Kurtosis: -0.975 Ordinate at the mode: 0.666667 Parameters of VGGFR: 3 2 Std. Deviation: 0.41574 Excess Kurtosis: -0.869805 Ordinate at the mode: 0.777778 Parameters of VGGFR: 3 3 Std. Deviation: 0.379517 Excess Kurtosis: -0.808942 Ordinate at the mode: 0.864198 Parameters of VGGFR: 3 4 Std. Deviation: 0.353311 Excess Kurtosis: -0.769285 Ordinate at the mode: 0.936214 Parameters of VGGFR: 3 5 Std. Deviation: 0.333105 Excess Kurtosis: -0.741401 Ordinate at the mode: 0.998628 Parameters of VGGFR: 3 6 Std. Deviation: 0.316846 Excess Kurtosis: -0.720728 Ordinate at the mode: 1.054108 Parameters of VGGFR: 3 7 Std. Deviation: 0.303357 Excess Kurtosis: -0.704789 Ordinate at the mode: 1.104303 Parameters of VGGFR: 3 8 Std. Deviation: 0.291905 Excess Kurtosis: -0.692126 Ordinate at the mode: 1.150316 Parameters of VGGFR: 3 9 Std. Deviation: 0.282007 Excess Kurtosis: -0.681823 Ordinate at the mode: 1.19292 Parameters of VGGFR: 3 10 Std. Deviation: 0.273328 Excess Kurtosis: -0.673276 Ordinate at the mode: 1.232684 Parameters of VGGFR: 3 11 Std. Deviation: 0.265627 Excess Kurtosis: -0.666073 Ordinate at the mode: 1.270038 Parameters of VGGFR: 3 12 Std. Deviation: 0.258726 Excess Kurtosis: -0.659919 Ordinate at the mode: 1.305317 Parameters of VGGFR: 3.5 1 Std. Deviation: 0.480384 Excess Kurtosis: -1.011765 Ordinate at the mode: 0.642857 Parameters of VGGFR: 3.5 2 Std. Deviation: 0.429669 Excess Kurtosis: -0.928922 Ordinate at the mode: 0.734694 Parameters of VGGFR: 3.5 3 Std. Deviation: 0.396566 Excess Kurtosis: -0.882446 Ordinate at the mode: 0.804665 Parameters of VGGFR: 3.5 4 Std. Deviation: 0.372509 Excess Kurtosis: -0.852726 Ordinate at the mode: 0.862141 Parameters of VGGFR: 3.5 5 Std. Deviation: 0.353872 Excess Kurtosis: -0.832091 Ordinate at the mode: 0.911406 Parameters of VGGFR: 3.5 6 Std. Deviation: 0.338807 Excess Kurtosis: -0.816931 Ordinate at the mode: 0.954806 Parameters of VGGFR: 3.5 7 Std. Deviation: 0.326254 Excess Kurtosis: -0.805323 Ordinate at the mode: 0.993778 Parameters of VGGFR: 3.5 8 Std. Deviation: 0.315554 Excess Kurtosis: -0.79615 Ordinate at the mode: 1.02927 Parameters of VGGFR: 3.5 9 Std. Deviation: 0.306271 Excess Kurtosis: -0.788718 Ordinate at the mode: 1.061945 Parameters of VGGFR: 3.5 10 Std. Deviation: 0.298102 Excess Kurtosis: -0.782576 Ordinate at the mode: 1.092287 Parameters of VGGFR: 3.5 11 Std. Deviation: 0.29083 Excess Kurtosis: -0.777414 Ordinate at the mode: 1.120658 Parameters of VGGFR: 3.5 12 Std. Deviation: 0.284294 Excess Kurtosis: -0.773015 Ordinate at the mode: 1.14734 Parameters of VGGFR: 4 1 Std. Deviation: 0.48795 Excess Kurtosis: -1.04 Ordinate at the mode: 0.625 Parameters of VGGFR: 4 2 Std. Deviation: 0.441367 Excess Kurtosis: -0.972991 Ordinate at the mode: 0.703125 Parameters of VGGFR: 4 3 Std. Deviation: 0.410891 Excess Kurtosis: -0.936304 Ordinate at the mode: 0.761719 Parameters of VGGFR: 4 4 Std. Deviation: 0.388664 Excess Kurtosis: -0.913181 Ordinate at the mode: 0.809326 Parameters of VGGFR: 4 5 Std. Deviation: 0.371381 Excess Kurtosis: -0.897281 Ordinate at the mode: 0.849792 Parameters of VGGFR: 4 6 Std. Deviation: 0.357362 Excess Kurtosis: -0.88568 Ordinate at the mode: 0.885201 Parameters of VGGFR: 4 7 Std. Deviation: 0.345642 Excess Kurtosis: -0.876843 Ordinate at the mode: 0.916815 Parameters of VGGFR: 4 8 Std. Deviation: 0.335621 Excess Kurtosis: -0.869887 Ordinate at the mode: 0.945465 Parameters of VGGFR: 4 9 Std. Deviation: 0.326902 Excess Kurtosis: -0.864271 Ordinate at the mode: 0.971728 Parameters of VGGFR: 4 10 Std. Deviation: 0.319209 Excess Kurtosis: -0.85964 Ordinate at the mode: 0.996021 Parameters of VGGFR: 4 11 Std. Deviation: 0.312344 Excess Kurtosis: -0.855758 Ordinate at the mode: 1.018658 Parameters of VGGFR: 4 12 Std. Deviation: 0.306158 Excess Kurtosis: -0.852455 Ordinate at the mode: 1.03988 Parameters of VGGFR: 4.5 1 Std. Deviation: 0.494413 Excess Kurtosis: -1.062201 Ordinate at the mode: 0.611111 Parameters of VGGFR: 4.5 2 Std. Deviation: 0.451335 Excess Kurtosis: -1.006835 Ordinate at the mode: 0.679012 Parameters of VGGFR: 4.5 3 Std. Deviation: 0.423099 Excess Kurtosis: -0.977114 Ordinate at the mode: 0.72931 Parameters of VGGFR: 4.5 4 Std. Deviation: 0.402447 Excess Kurtosis: -0.958598 Ordinate at the mode: 0.769827 Parameters of VGGFR: 4.5 5 Std. Deviation: 0.386342 Excess Kurtosis: -0.945963 Ordinate at the mode: 0.804041 Parameters of VGGFR: 4.5 6 Std. Deviation: 0.373242 Excess Kurtosis: -0.936793 Ordinate at the mode: 0.833821 Parameters of VGGFR: 4.5 7 Std. Deviation: 0.362262 Excess Kurtosis: -0.929836 Ordinate at the mode: 0.860291 Parameters of VGGFR: 4.5 8 Std. Deviation: 0.352851 Excess Kurtosis: -0.924377 Ordinate at the mode: 0.884188 Parameters of VGGFR: 4.5 9 Std. Deviation: 0.344644 Excess Kurtosis: -0.91998 Ordinate at the mode: 0.90602 Parameters of VGGFR: 4.5 10 Std. Deviation: 0.337388 Excess Kurtosis: -0.916363 Ordinate at the mode: 0.926154 Parameters of VGGFR: 4.5 11 Std. Deviation: 0.330899 Excess Kurtosis: -0.913335 Ordinate at the mode: 0.944864 Parameters of VGGFR: 4.5 12 Std. Deviation: 0.325042 Excess Kurtosis: -0.910763 Ordinate at the mode: 0.962361 Parameters of VGGFR: 5 1 Std. Deviation: 0.5 Excess Kurtosis: -1.08 Ordinate at the mode: 0.6 Parameters of VGGFR: 5 2 Std. Deviation: 0.459933 Excess Kurtosis: -1.033455 Ordinate at the mode: 0.66 Parameters of VGGFR: 5 3 Std. Deviation: 0.433629 Excess Kurtosis: -1.008873 Ordinate at the mode: 0.704 Parameters of VGGFR: 5 4 Std. Deviation: 0.414347 Excess Kurtosis: -0.993702 Ordinate at the mode: 0.7392 Parameters of VGGFR: 5 5 Std. Deviation: 0.399275 Excess Kurtosis: -0.983414 Ordinate at the mode: 0.768768 Parameters of VGGFR: 5 6 Std. Deviation: 0.386986 Excess Kurtosis: -0.975979 Ordinate at the mode: 0.794394 Parameters of VGGFR: 5 7 Std. Deviation: 0.376665 Excess Kurtosis: -0.970357 Ordinate at the mode: 0.817091 Parameters of VGGFR: 5 8 Std. Deviation: 0.367801 Excess Kurtosis: -0.965957 Ordinate at the mode: 0.837518 Parameters of VGGFR: 5 9 Std. Deviation: 0.360057 Excess Kurtosis: -0.962419 Ordinate at the mode: 0.856129 Parameters of VGGFR: 5 10 Std. Deviation: 0.353198 Excess Kurtosis: -0.959513 Ordinate at the mode: 0.873252 Parameters of VGGFR: 5 11 Std. Deviation: 0.347055 Excess Kurtosis: -0.957084 Ordinate at the mode: 0.889129 Parameters of VGGFR: 5 12 Std. Deviation: 0.341502 Excess Kurtosis: -0.955023 Ordinate at the mode: 0.903948 Parameters of VGGFR: 5.5 1 Std. Deviation: 0.504878 Excess Kurtosis: -1.094505 Ordinate at the mode: 0.590909 Parameters of VGGFR: 5.5 2 Std. Deviation: 0.467426 Excess Kurtosis: -1.054808 Ordinate at the mode: 0.644628 Parameters of VGGFR: 5.5 3 Std. Deviation: 0.442807 Excess Kurtosis: -1.034128 Ordinate at the mode: 0.683696 Parameters of VGGFR: 5.5 4 Std. Deviation: 0.424726 Excess Kurtosis: -1.021465 Ordinate at the mode: 0.714774 Parameters of VGGFR: 5.5 5 Std. Deviation: 0.410564 Excess Kurtosis: -1.012921 Ordinate at the mode: 0.740765 Parameters of VGGFR: 5.5 6 Std. Deviation: 0.398997 Excess Kurtosis: -1.006769 Ordinate at the mode: 0.763213 Parameters of VGGFR: 5.5 7 Std. Deviation: 0.389264 Excess Kurtosis: -1.002129 Ordinate at the mode: 0.783036 Parameters of VGGFR: 5.5 8 Std. Deviation: 0.380892 Excess Kurtosis: -0.998505 Ordinate at the mode: 0.800833 Parameters of VGGFR: 5.5 9 Std. Deviation: 0.373566 Excess Kurtosis: -0.995596 Ordinate at the mode: 0.817011 Parameters of VGGFR: 5.5 10 Std. Deviation: 0.367069 Excess Kurtosis: -0.99321 Ordinate at the mode: 0.831866 Parameters of VGGFR: 5.5 11 Std. Deviation: 0.361242 Excess Kurtosis: -0.991217 Ordinate at the mode: 0.845616 Parameters of VGGFR: 5.5 12 Std. Deviation: 0.355968 Excess Kurtosis: -0.989528 Ordinate at the mode: 0.858428 Parameters of VGGFR: 6 1 Std. Deviation: 0.509175 Excess Kurtosis: -1.106494 Ordinate at the mode: 0.583333 Parameters of VGGFR: 6 2 Std. Deviation: 0.474016 Excess Kurtosis: -1.072222 Ordinate at the mode: 0.631944 Parameters of VGGFR: 6 3 Std. Deviation: 0.450879 Excess Kurtosis: -1.054576 Ordinate at the mode: 0.667052 Parameters of VGGFR: 6 4 Std. Deviation: 0.433859 Excess Kurtosis: -1.043843 Ordinate at the mode: 0.694846 Parameters of VGGFR: 6 5 Std. Deviation: 0.420506 Excess Kurtosis: -1.036631 Ordinate at the mode: 0.718008 Parameters of VGGFR: 6 6 Std. Deviation: 0.409582 Excess Kurtosis: -1.031454 Ordinate at the mode: 0.737953 Parameters of VGGFR: 6 7 Std. Deviation: 0.400377 Excess Kurtosis: -1.027558 Ordinate at the mode: 0.755523 Parameters of VGGFR: 6 8 Std. Deviation: 0.392448 Excess Kurtosis: -1.02452 Ordinate at the mode: 0.771263 Parameters of VGGFR: 6 9 Std. Deviation: 0.385501 Excess Kurtosis: -1.022085 Ordinate at the mode: 0.785546 Parameters of VGGFR: 6 10 Std. Deviation: 0.379333 Excess Kurtosis: -1.02009 Ordinate at the mode: 0.798638 Parameters of VGGFR: 6 11 Std. Deviation: 0.373795 Excess Kurtosis: -1.018425 Ordinate at the mode: 0.810739 Parameters of VGGFR: 6 12 Std. Deviation: 0.368777 Excess Kurtosis: -1.017015 Ordinate at the mode: 0.821999 > > > > cleanEx() > nameEx("Gabbs") > ### * Gabbs > > flush(stderr()); flush(stdout()) > > ### Name: Gabbs > ### Title: Dimensions of fetal metacurpals and phalanges > ### Aliases: Gabbs > ### Keywords: datasets > > ### ** Examples > > data(Gabbs) > head(Gabbs) V15-29mm V30-44mm V45-59mm V60-74mm V75-89mm V90-104mm AR Metacarpal1 0.44 0.86 1.34 2.06 2.21 3.34 15 Metacarpal2 0.76 1.31 2.11 3.31 3.36 5.50 19 Metacarpal3 0.77 1.29 1.99 3.04 3.15 5.18 1 Metacarpal4 0.71 1.17 1.64 2.70 2.89 4.75 17 Metacarpal5 0.67 1.12 1.58 2.56 2.63 4.23 16 Proximal1 0.35 0.64 0.97 1.32 1.46 2.18 10 > > > > cleanEx() > nameEx("Insuhyper") > ### * Insuhyper > > flush(stderr()); flush(stdout()) > > ### Name: Insuhyper > ### Title: Role of insulin in endogenous hypertriglyceridemia > ### Aliases: Insuhyper > ### Keywords: datasets > > ### ** Examples > > data(Insuhyper) > head(Insuhyper) Age Sex Height Weight PI TG PG8am PG9am PG10am PGAvg PI8am PI9am PI10am G_L 45 1 170.0 68.0 12.62 243 88 164 112 121 10 177 131 A_T 51 1 154.9 44.7 14.32 316 84 198 154 145 14 172 136 E_A 61 1 154.9 57.9 12.13 367 121 265 309 232 14 112 160 H_G 76 1 158.0 70.1 11.62 420 331 525 585 480 36 62 82 E_D 58 1 152.4 68.2 11.30 495 179 471 465 372 26 117 153 P_M 47 1 162.0 111.0 10.18 509 236 488 533 419 41 112 91 PIAvg PF8am PF9am PF10am PFAvg G_L 106 375 204 190 256 A_T 107 343 188 147 226 E_A 95 417 217 211 282 H_G 60 1082 614 346 681 E_D 99 767 411 252 477 P_M 81 1020 689 330 680 > > > > cleanEx() > nameEx("Lambh") > ### * Lambh > > flush(stderr()); flush(stdout()) > > ### Name: Lambh > ### Title: Volcanic dust in the atmosphere > ### Aliases: Lambh > ### Keywords: datasets > > ### ** Examples > > data(Lambh) > head(Lambh) DVI Temp 1 113.0 9.043 2 24.5 9.059 3 30.5 9.230 4 142.0 8.861 5 52.0 9.108 6 18.0 9.159 > > > > cleanEx() > nameEx("Laudaher") > ### * Laudaher > > flush(stderr()); flush(stdout()) > > ### Name: Laudaher > ### Title: Data on the significance of lymphocytic infiltration > ### Aliases: Laudaher > ### Keywords: datasets > > ### ** Examples > > data(Laudaher) > head(Laudaher) Duration Infiltration 1 94 5.0 2 26 5.0 3 19 3.0 4 18 3.3 5 11 2.3 6 11 2.3 > > > > cleanEx() > nameEx("Locre") > ### * Locre > > flush(stderr()); flush(stdout()) > > ### Name: Locre > ### Title: Adult sex differences in crying behavior. > ### Aliases: Locre > ### Keywords: datasets > > ### ** Examples > > data(Locre) > head(Locre) Males Females Death_of_someone_close 78.9 86.0 Empathy_with_movie.TV.book 33.7 79.9 Break_up.of.romance 35.1 70.1 Someone.close.hurt 29.1 59.7 Death.of.pet 23.2 58.4 Fight.with.loved.one 25.6 57.5 > > > > cleanEx() > nameEx("Marozzi") > ### * Marozzi > > flush(stderr()); flush(stdout()) > > ### Name: Marozzi > ### Title: Trust in public institutions > ### Aliases: Marozzi > ### Keywords: datasets > > ### ** Examples > > data(Marozzi) > head(Marozzi) Ind_1 Ind_2 Ind_3 Ind_4 Ind_5 Ind_6 Ind_7 Albania 0.197 0.268 0.420 0.122 0.147 0.701 0.728 Belgium 0.442 0.437 0.665 0.313 0.296 0.466 0.529 Bulgaria 0.074 0.100 0.227 0.042 0.044 0.322 0.389 Switzerland 0.632 0.705 0.819 0.453 0.414 0.350 0.514 Cyprus 0.184 0.404 0.436 0.085 0.087 0.321 0.272 Czech Republic 0.186 0.297 0.425 0.140 0.138 0.276 0.382 > > > > cleanEx() > nameEx("Radiation") > ### * Radiation > > flush(stderr()); flush(stdout()) > > ### Name: Radiation > ### Title: Global radiation data > ### Aliases: Radiation > ### Keywords: datasets > > ### ** Examples > > data(Radiation) > head(Radiation) Mt Gr Msd January 25.32 17.61 7.86 February 26.80 21.11 9.15 March 28.81 23.07 8.77 April 30.92 23.55 9.11 May 32.91 21.94 8.24 June 32.35 20.60 6.77 > > > > cleanEx() > nameEx("Security") > ### * Security > > flush(stderr()); flush(stdout()) > > ### Name: Security > ### Title: Accessibility, security, and accuracy in statistical databases > ### Aliases: Security > ### Keywords: datasets > > ### ** Examples > > data(Security) > head(Security) OV AP MP 1 -5871.22 46265.37 -11993.33 2 -3851.81 -28798.95 -1929.97 3 -3851.81 -9120.43 -2865.00 4 -1206.41 -13267.11 -915.41 5 -248.34 -15618.57 -172.00 6 -49.19 -815242.00 -41.22 > > > > cleanEx() > nameEx("Sharpe") > ### * Sharpe > > flush(stderr()); flush(stdout()) > > ### Name: Sharpe > ### Title: Performance of mutual funds, 1954-1963 > ### Aliases: Sharpe > ### Keywords: datasets > > ### ** Examples > > data(Sharpe) > head(Sharpe) AVR VAR Affiliated 14.6 15.3 ABS 10.0 9.2 AxeHouA 10.5 13.5 AxeHouB 12.0 16.3 AxeHouSF 11.9 15.6 Boston 12.4 12.1 > > > > cleanEx() > nameEx("Starshi") > ### * Starshi > > flush(stderr()); flush(stdout()) > > ### Name: Starshi > ### Title: Policy-related determinants of health > ### Aliases: Starshi > ### Keywords: datasets > > ### ** Examples > > data(Starshi) > head(Starshi) Sm15F Sm15M Japan 13.3 60.4 Sweden 26.6 25.2 Canada 26.0 26.0 France 20.0 38.0 Australia 23.8 28.2 Spain 21.0 44.0 > > > > cleanEx() > nameEx("Viscoh") > ### * Viscoh > > flush(stderr()); flush(stdout()) > > ### Name: Viscoh > ### Title: Data on death measures for various conditions > ### Aliases: Viscoh > ### Keywords: datasets > > ### ** Examples > > data(Viscoh) > head(Viscoh) Deaths LLE EYLL All_causes 97299.5 13.53 13.16 All_disease 92483.4 12.58 11.64 Heart_disease 34956.3 10.56 3.69 All_cancer 21839.7 14.39 3.14 All_other_diseases 8188.6 11.97 0.98 Stroke 7238.9 9.54 0.69 > > > > cleanEx() > nameEx("Zoutus") > ### * Zoutus > > flush(stderr()); flush(stdout()) > > ### Name: Zoutus > ### Title: Data on time and dose of CT fluoroscopy-guided procedures > ### Aliases: Zoutus > ### Keywords: datasets > > ### ** Examples > > data(Zoutus) > head(Zoutus) LogTime LogDose 1 3.61 1.48 2 3.87 1.24 3 3.95 2.08 4 4.04 1.70 5 4.06 2.08 6 4.11 2.94 > > > > cleanEx() > nameEx("comprank") > ### * comprank > > flush(stderr()); flush(stdout()) > > ### Name: comprank > ### Title: Computes various rank correlation coefficients > ### Aliases: comprank > > ### ** Examples > > > data(Zoutus);attach(Zoutus);print(Zoutus) LogTime LogDose 1 3.61 1.48 2 3.87 1.24 3 3.95 2.08 4 4.04 1.70 5 4.06 2.08 6 4.11 2.94 7 4.19 2.24 8 4.20 1.85 9 4.32 2.84 10 4.32 3.93 11 4.42 3.03 12 4.42 3.23 13 4.45 3.87 14 4.50 3.55 15 4.52 2.81 16 4.57 4.07 17 4.58 4.44 18 4.61 3.16 19 4.74 4.19 > a<-comprank(LogDose,LogTime,"spearman","woodbury") ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110C: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC1791: s_wood_ (packages/tests-vg/pvrank/src/DealwT.f90:1559) ==19863== by 0x1ADC5AE3: woodbury_ (packages/tests-vg/pvrank/src/DealwT.f90:1604) ==19863== by 0x1ADC84AE: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1663) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC1791: s_wood_ (packages/tests-vg/pvrank/src/DealwT.f90:1559) ==19863== by 0x1ADC5AE3: woodbury_ (packages/tests-vg/pvrank/src/DealwT.f90:1604) ==19863== by 0x1ADC84AE: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1663) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC1791: s_wood_ (packages/tests-vg/pvrank/src/DealwT.f90:1559) ==19863== by 0x1ADC5AE3: woodbury_ (packages/tests-vg/pvrank/src/DealwT.f90:1604) ==19863== by 0x1ADC84AE: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1663) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC1791: s_wood_ (packages/tests-vg/pvrank/src/DealwT.f90:1559) ==19863== by 0x1ADC5AE3: woodbury_ (packages/tests-vg/pvrank/src/DealwT.f90:1604) ==19863== by 0x1ADC84AE: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1663) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110C: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC17A8: s_wood_ (packages/tests-vg/pvrank/src/DealwT.f90:1560) ==19863== by 0x1ADC5AE3: woodbury_ (packages/tests-vg/pvrank/src/DealwT.f90:1604) ==19863== by 0x1ADC84AE: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1663) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC17A8: s_wood_ (packages/tests-vg/pvrank/src/DealwT.f90:1560) ==19863== by 0x1ADC5AE3: woodbury_ (packages/tests-vg/pvrank/src/DealwT.f90:1604) ==19863== by 0x1ADC84AE: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1663) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC17A8: s_wood_ (packages/tests-vg/pvrank/src/DealwT.f90:1560) ==19863== by 0x1ADC5AE3: woodbury_ (packages/tests-vg/pvrank/src/DealwT.f90:1604) ==19863== by 0x1ADC84AE: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1663) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC17A8: s_wood_ (packages/tests-vg/pvrank/src/DealwT.f90:1560) ==19863== by 0x1ADC5AE3: woodbury_ (packages/tests-vg/pvrank/src/DealwT.f90:1604) ==19863== by 0x1ADC84AE: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1663) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== > cat(a$r,a$ities,"\n") 0.839473684 1 > a<-comprank(LogDose,LogTime,"kendall","woodbury") > cat(a$r,a$ities,"\n") 0.666666672 1 > detach(Zoutus) > ##### > # > # Yotopoulos, P. A. Nugent, J. B. (1973). A balanced-growth version of the > # linkage hypothesis: a test. The Quarterly Journal of Economics, 87, 157-171. > x<-1:18 > y<-c(7,1,4,8,3,9,2,5,10,6,17,13,14,12,11,16,15,18) > a<-comprank(x,y,"gini");cat(a$r,a$ities,"\n") 0.703703704 6 > ##### > # > data(Franzen);attach(Franzen) > op<-par(mfrow=c(1,1)) > plot(MECISSP,PPP, main="Environmental Attitudes in International Comparison", + xlab="Mean Environmental concern ISSP", ylab="Purchasing power parity", + pch=19, cex=0.8,col="salmon3") > abline(h=mean(PPP),col="darkred",lty=2,lwd=1) > abline(v=mean(MECISSP),col="darkred",lty=2,lwd=1) > par(op) > a<-comprank(MECISSP,PPP,"kendall","gh");cat(a$r,a$ities,"\n") ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110C: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC2149: giuls_ (packages/tests-vg/pvrank/src/DealwT.f90:690) ==19863== by 0x1ADC3111: ggh_ (packages/tests-vg/pvrank/src/DealwT.f90:716) ==19863== by 0x1ADC8456: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1665) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC2149: giuls_ (packages/tests-vg/pvrank/src/DealwT.f90:690) ==19863== by 0x1ADC3111: ggh_ (packages/tests-vg/pvrank/src/DealwT.f90:716) ==19863== by 0x1ADC8456: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1665) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC2149: giuls_ (packages/tests-vg/pvrank/src/DealwT.f90:690) ==19863== by 0x1ADC3111: ggh_ (packages/tests-vg/pvrank/src/DealwT.f90:716) ==19863== by 0x1ADC8456: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1665) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC2149: giuls_ (packages/tests-vg/pvrank/src/DealwT.f90:690) ==19863== by 0x1ADC3111: ggh_ (packages/tests-vg/pvrank/src/DealwT.f90:716) ==19863== by 0x1ADC8456: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1665) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== 0.570000014 2 > a<-comprank(MECISSP,PPP,"gini","wgh");cat(a$r,a$ities,"\n") ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110C: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC2149: giuls_ (packages/tests-vg/pvrank/src/DealwT.f90:690) ==19863== by 0x1ADC3111: ggh_ (packages/tests-vg/pvrank/src/DealwT.f90:716) ==19863== by 0x1ADC83F7: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1668) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC2149: giuls_ (packages/tests-vg/pvrank/src/DealwT.f90:690) ==19863== by 0x1ADC3111: ggh_ (packages/tests-vg/pvrank/src/DealwT.f90:716) ==19863== by 0x1ADC83F7: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1668) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC2149: giuls_ (packages/tests-vg/pvrank/src/DealwT.f90:690) ==19863== by 0x1ADC3111: ggh_ (packages/tests-vg/pvrank/src/DealwT.f90:716) ==19863== by 0x1ADC83F7: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1668) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC2149: giuls_ (packages/tests-vg/pvrank/src/DealwT.f90:690) ==19863== by 0x1ADC3111: ggh_ (packages/tests-vg/pvrank/src/DealwT.f90:716) ==19863== by 0x1ADC83F7: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1668) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== 0.334935897 3 > detach(Franzen) > ##### > # > data(Viscoh);attach(Viscoh) > Viscoh<-as.matrix(Viscoh) > a<-comprank(Viscoh,"spearman","gh",print=FALSE) Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 2 3 > print(a$r);cat(" method:", a$ities,"\n") [,1] [,2] [,3] [1,] 1.000000000 -0.676354730 0.975246305 [2,] -0.676354730 1.000000000 -0.548891672 [3,] 0.975246305 -0.548891672 1.000000000 method: 6 > b<-comprank(Viscoh,"r4","gh",print=FALSE) Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 2 3 > print(b$r);cat(" method:", b$ities,"\n") [,1] [,2] [,3] [1,] 1.000000000 -0.716649493 0.966682489 [2,] -0.716649493 1.000000000 -0.589805793 [3,] 0.966682489 -0.589805793 1.000000000 method: 6 > c<-comprank(Viscoh,"fy1","wgh",print=FALSE) Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 2 3 > print(c$r);cat("method:", c$ities,"\n") [,1] [,2] [,3] [1,] 1.000000000 -0.686428889 0.987974765 [2,] -0.686428889 1.000000000 -0.574442461 [3,] 0.987974765 -0.574442461 1.000000000 method: 6 > d<-comprank(Viscoh,"fy2","wgh",print=FALSE) Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 2 3 > print(d$r);cat(" method:", d$ities,"\n") [,1] [,2] [,3] [1,] 1.000000000 -0.522693865 0.777962940 [2,] -0.522693865 1.000000000 -0.374194971 [3,] 0.777962940 -0.374194971 1.000000000 method: 6 > d<-comprank(Viscoh,"sbz","wgh",print=FALSE) Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 2 3 > print(d$r);cat(" method:", d$ities,"\n") [,1] [,2] [,3] [1,] 1.000000000 -0.520689655 0.894334975 [2,] -0.520689655 1.000000000 -0.416748768 [3,] 0.894334975 -0.416748768 1.000000000 method: 6 > detach(Viscoh) > ##### > # > data(Laudaher);attach(Laudaher) > a1<-comprank(Duration,Infiltration,"gini","midrank") ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110C: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC152C: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:764) ==19863== by 0x1ADC83AC: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1672) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== by 0x506827: R_ReplConsole (svn/R-devel/src/main/main.c:308) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC152C: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:764) ==19863== by 0x1ADC83AC: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1672) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC152C: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:764) ==19863== by 0x1ADC83AC: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1672) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== by 0x506827: R_ReplConsole (svn/R-devel/src/main/main.c:308) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC152C: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:764) ==19863== by 0x1ADC83AC: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1672) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110C: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC153E: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:765) ==19863== by 0x1ADC83AC: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1672) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== by 0x506827: R_ReplConsole (svn/R-devel/src/main/main.c:308) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC153E: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:765) ==19863== by 0x1ADC83AC: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1672) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC153E: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:765) ==19863== by 0x1ADC83AC: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1672) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== by 0x506827: R_ReplConsole (svn/R-devel/src/main/main.c:308) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC153E: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:765) ==19863== by 0x1ADC83AC: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1672) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== > a2<-comprank(Duration,Infiltration,"gini","dubois") ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110C: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC152C: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:764) ==19863== by 0x1ADC835C: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1691) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== by 0x506827: R_ReplConsole (svn/R-devel/src/main/main.c:308) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC152C: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:764) ==19863== by 0x1ADC835C: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1691) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC152C: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:764) ==19863== by 0x1ADC835C: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1691) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== by 0x506827: R_ReplConsole (svn/R-devel/src/main/main.c:308) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC152C: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:764) ==19863== by 0x1ADC835C: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1691) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110C: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC153E: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:765) ==19863== by 0x1ADC835C: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1691) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== by 0x506827: R_ReplConsole (svn/R-devel/src/main/main.c:308) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC153E: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:765) ==19863== by 0x1ADC835C: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1691) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== ==19863== Conditional jump or move depends on uninitialised value(s) ==19863== at 0x1ADC110E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1401) ==19863== by 0x1ADC153E: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:765) ==19863== by 0x1ADC835C: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1691) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== by 0x506827: R_ReplConsole (svn/R-devel/src/main/main.c:308) ==19863== Uninitialised value was created by a heap allocation ==19863== at 0x4C2DB6B: malloc (/builddir/build/BUILD/valgrind-3.13.0/coregrind/m_replacemalloc/vg_replace_malloc.c:299) ==19863== by 0x1ADC109E: meanrank_ (packages/tests-vg/pvrank/src/DealwT.f90:1393) ==19863== by 0x1ADC153E: midrank_ (packages/tests-vg/pvrank/src/DealwT.f90:765) ==19863== by 0x1ADC835C: dealwt_ (packages/tests-vg/pvrank/src/DealwT.f90:1691) ==19863== by 0x49B287: do_dotCode (svn/R-devel/src/main/dotcode.c:1781) ==19863== by 0x4CC2C9: bcEval (svn/R-devel/src/main/eval.c:6772) ==19863== by 0x4DA82F: Rf_eval (svn/R-devel/src/main/eval.c:624) ==19863== by 0x4DC13E: R_execClosure (svn/R-devel/src/main/eval.c:1764) ==19863== by 0x4DA9BC: Rf_eval (svn/R-devel/src/main/eval.c:747) ==19863== by 0x4DE9BF: do_set (svn/R-devel/src/main/eval.c:2778) ==19863== by 0x4DABDB: Rf_eval (svn/R-devel/src/main/eval.c:699) ==19863== by 0x50644C: Rf_ReplIteration (svn/R-devel/src/main/main.c:258) ==19863== > a3<-comprank(Duration,Infiltration,"spearman","midrank") > a4<-comprank(Duration,Infiltration,"spearman","dubois") > cat("Coefficient","method","\n",a1$r,a1$ities,"\n",a2$r,a2$ities,"\n",a3$r,a3$ities, + "\n",a4$r,a4$ities,"\n") Coefficient method 0.6640625 4 0.661011698 5 0.738235289 4 0.738734924 5 > detach(Laudaher) > ##### > # Asymptotic confidence intervals. > r.cofint <- function(r, n, index, Level=.95) { + asd<-rep(0,4) + asd[1]<-1/sqrt(n-1) # Spearman + asd[2]<-sqrt((4*n+10)/(9*n*(n-1))) # Kendall + asd[3]<-1/sqrt(1.5*n) # Gini + asd[4]<-1/sqrt(1.00762*(n-1)) # r4 + # Fisher-Yates-means, Fisher-Yates-medians, symmetrical BZ + if (index<=4) {zse<-r*asd[index]} else {zse<-atan(r*(1-0.6/(n+8)))/sqrt(n-3)} + rlow <- r - zse * qnorm((1-Level)/2,lower.tail=FALSE);rlow<-max(-1,rlow) + rupp <- r + zse * qnorm((1-Level)/2,lower.tail=FALSE);rupp<-min(1,rupp) + out<-list(Lower.r=rlow, Upper.r=rupp) + return(out) + } > # > # Rajalingam S. and Zeya O. "Summative assessments by Spearman`s correlation > # coefficient: a case study. in enhancing learning: teaching & learning international > # conference. Nov 24-26 2011. Miri, Sarawak: Curtin University > > FA=c(39.09,39.77,35.62,34.69,34.42,35.57,35.94,38.68,38.41,36.4,37.95,37.03,38.03, + 35.5,38.55,30.68,26.3,36.28) > SA=c(27.5,34,24,24,17,17.5,16.5,26,25.5,28.5,26.5,13,12.5,9.5,28.5,23,24.5,22) > n<-length(FA) > op<-par(mfrow=c(1,1)) > plot(FA,SA, main="Academic progresses of students", + xlab="Formative assessment", ylab="summative assessment", + pch=19,col="steelblue4") > text(FA,SA,label=1:n,cex=0.8,pos=2) > abline(h=mean(SA),col="darkblue",lty=2,lwd=1) > abline(v=mean(FA),col="darkblue",lty=2,lwd=1) > par(op) > rct<-c("spearman","kendall","gini","r4", "fy1","fy2","sbz") > for (index in 1:7){ + r<-comprank(FA,SA,rct[index])$r + cir<-r.cofint(r,n,index,Level=0.99) + cat(rct[index],"Low:",cir$Lower,"Value:",r," Upp: ",cir$Upper.r,"\n") + } spearman Low: 0.195573956 Value: 0.521155831 Upp: 0.846737706 kendall Low: 0.18517681 Value: 0.333333337 Upp: 0.481489864 gini Low: 0.118288234 Value: 0.234567901 Upp: 0.350847568 r4 Low: 0.201190346 Value: 0.532762149 Upp: 0.864333953 fy1 Low: 0.20038111 Value: 0.505403038 Upp: 0.810424965 fy2 Low: 0.139745044 Value: 0.371001032 Upp: 0.60225702 sbz Low: 0 Value: 0 Upp: 0 > ##### > # > # Daniel, C. Wood, F. S. Fitting Equations to Data. New York: John Wiley, > # 1971, p. 45. Pilot-plant data > # The response variable (y) corresponds to the acid content determined by > # titration and the explanatory variable (x) is the organic acid content determined > # by extraction and weighting > y<-c(76, 70, 55, 71, 55, 48, 50, 66, 41, 43, 82, 68, 88, 58, 64, 88, 89, 88, + 84, 88) > x<-c(123, 109, 62, 104, 57, 37, 44, 100, 16, 28, 138, 105, 159, 75, 88, 164, + 169, 167, 149, 167) > out1<-comprank(x,y,"sp","woodbury") > out2<-comprank(x,y,"ke","woodbury") > out3<-comprank(x,y,"gi","woodbury") > out4<-comprank(x,y,"r4","woodbury") > out5<-comprank(x,y,"fy1","woodbury") > out6<-comprank(x,y,"fy2","woodbury") > out7<-comprank(x,y,"sbz","woodbury") > ind<-c(out1$r,out2$r,out3$r,out4$r,out5$r,out6$r,out7$r) > cat(out1$ities,"\n") 1 > print(round(ind,5)) [1] 0.98722 0.94211 0.47000 0.97821 0.98290 0.93177 0.00000 > ##### > # > data(DietFish);attach(DietFish) > op<-par(mfrow=c(1,1)) > plot(occ_D,occ_H, main="Comparison of the diets of banded killfish",xlab="occ F_diaphanus", + ylab="occ F_heteroclitus", pch=19, cex=0.8,col="darkcyan") > text(occ_D,occ_H, labels = rownames(DietFish), cex=0.5, pos=3) > abline(h=mean(occ_D),col="darkred",lty=2,lwd=1) > abline(v=mean(occ_H),col="darkred",lty=2,lwd=1) > par(op) > a<-comprank(occ_D,occ_H,"sp","woodbury");cat(a$r,a$ities,"\n") 0.270060484 1 > a<-comprank(occ_D,occ_H,"sp","gh");cat(a$r,a$ities,"\n") 0.271370979 2 > a<-comprank(occ_D,occ_H,"sp","wgh");cat(a$r,a$ities,"\n") 0.270491948 3 > a<-comprank(occ_D,occ_H,"sp","midrank");cat(a$r,a$ities,"\n") 0.275201624 4 > a<-comprank(occ_D,occ_H,"sp","dubois");cat(a$r,a$ities,"\n") 0.279985426 5 > detach(DietFish) > ##### > # > data(Radiation);attach(Radiation);Radiation<-as.matrix(Radiation) > r1<-comprank(Radiation,"spearman","midrank",print=FALSE);eigen(r1$r) eigen() decomposition $values [1] 2.041646997 0.905703491 0.052649512 $vectors [,1] [,2] [,3] [1,] 0.451507580 0.7970056605 0.401151945 [2,] 0.689919381 -0.0267422061 -0.723392081 [3,] 0.565819895 -0.6033795087 0.561943960 > r2<-comprank(Radiation,"kendall","midrank",print=FALSE);eigen(r2$r) eigen() decomposition $values [1] 1.793628844 0.944108421 0.262262734 $vectors [,1] [,2] [,3] [1,] -0.415773449 0.8329018558 0.365249144 [2,] -0.693778705 -0.0307925519 -0.719529657 [3,] -0.588050633 -0.5525634057 0.590652296 > r3<-comprank(Radiation,"gini","midrank",print=FALSE);eigen(r3$r) eigen() decomposition $values [1] 1.802508754 0.923890734 0.273600512 $vectors [,1] [,2] [,3] [1,] -0.419004951 0.8392048359 0.346655585 [2,] -0.688659922 -0.0448992825 -0.723693005 [3,] -0.591762083 -0.5419587603 0.596739759 > r4<-comprank(Radiation,"r4","gh",print=FALSE);eigen(r4$r) eigen() decomposition $values [1] 1.920456120 0.959041495 0.120502385 $vectors [,1] [,2] [,3] [1,] -0.445122772 0.7896473007 0.422282913 [2,] -0.698970313 -0.0116147043 -0.715056361 [3,] -0.559737634 -0.6134510895 0.557109991 > r5<-comprank(Radiation,"fy1","gh",print=FALSE);eigen(r5$r) eigen() decomposition $values [1] 1.960698182 0.935185891 0.104115927 $vectors [,1] [,2] [,3] [1,] -0.450370899 0.7904837506 0.415092151 [2,] -0.694549528 -0.0180432549 -0.719218599 [3,] -0.561041002 -0.6122171845 0.557156273 > r6<-comprank(Radiation,"fy2","gh",print=FALSE);eigen(r6$r) eigen() decomposition $values [1] 1.376258659 1.072696238 0.551045103 $vectors [,1] [,2] [,3] [1,] -0.203071406 0.891855693 0.404172522 [2,] 0.653206744 0.430891292 -0.622618378 [3,] 0.729440165 -0.137572227 0.670067854 > r7<-comprank(Radiation,"sbz","gh",print=FALSE);eigen(r7$r) eigen() decomposition $values [1] 1.811215612 0.933955298 0.254829090 $vectors [,1] [,2] [,3] [1,] -0.434317944 0.8154460334 0.382643032 [2,] -0.691661578 -0.0297787021 -0.721607573 [3,] -0.577037420 -0.5780666012 0.576946115 > detach(Radiation) > ##### > # > # Correlation matrix > data(Marozzi);attach(Marozzi) > Marozzi<-as.matrix(Marozzi) > cor1<-comprank(Marozzi,"spearman","midrank",print=FALSE) Tied scores are present in one or both rankings 1 2 Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 1 4 Tied scores are present in one or both rankings 1 5 Tied scores are present in one or both rankings 1 6 Tied scores are present in one or both rankings 1 7 Tied scores are present in one or both rankings 2 5 Tied scores are present in one or both rankings 2 6 Tied scores are present in one or both rankings 2 7 Tied scores are present in one or both rankings 3 5 Tied scores are present in one or both rankings 3 6 Tied scores are present in one or both rankings 3 7 Tied scores are present in one or both rankings 4 5 Tied scores are present in one or both rankings 4 6 Tied scores are present in one or both rankings 4 7 Tied scores are present in one or both rankings 5 6 Tied scores are present in one or both rankings 5 7 Tied scores are present in one or both rankings 6 7 > rownames(cor1$r)<-colnames(Marozzi);rownames(cor1$r)<-colnames(Marozzi) > print(round(cor1$r,3)) [,1] [,2] [,3] [,4] [,5] [,6] [,7] Ind_1 1.000 0.926 0.789 0.941 0.934 0.437 0.625 Ind_2 0.926 1.000 0.854 0.895 0.852 0.417 0.625 Ind_3 0.789 0.854 1.000 0.734 0.652 0.392 0.717 Ind_4 0.941 0.895 0.734 1.000 0.979 0.445 0.603 Ind_5 0.934 0.852 0.652 0.979 1.000 0.456 0.591 Ind_6 0.437 0.417 0.392 0.445 0.456 1.000 0.761 Ind_7 0.625 0.625 0.717 0.603 0.591 0.761 1.000 > cor1 <-comprank(Marozzi,"kendall","midrank",print=FALSE) Tied scores are present in one or both rankings 1 2 Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 1 4 Tied scores are present in one or both rankings 1 5 Tied scores are present in one or both rankings 1 6 Tied scores are present in one or both rankings 1 7 Tied scores are present in one or both rankings 2 5 Tied scores are present in one or both rankings 2 6 Tied scores are present in one or both rankings 2 7 Tied scores are present in one or both rankings 3 5 Tied scores are present in one or both rankings 3 6 Tied scores are present in one or both rankings 3 7 Tied scores are present in one or both rankings 4 5 Tied scores are present in one or both rankings 4 6 Tied scores are present in one or both rankings 4 7 Tied scores are present in one or both rankings 5 6 Tied scores are present in one or both rankings 5 7 Tied scores are present in one or both rankings 6 7 > print(round(cor1$r,3)) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.000 0.783 0.596 0.808 0.793 0.296 0.458 [2,] 0.783 1.000 0.700 0.729 0.680 0.310 0.468 [3,] 0.596 0.700 1.000 0.517 0.458 0.246 0.557 [4,] 0.808 0.729 0.517 1.000 0.901 0.315 0.419 [5,] 0.793 0.680 0.458 0.901 1.000 0.310 0.414 [6,] 0.296 0.310 0.246 0.315 0.310 1.000 0.562 [7,] 0.458 0.468 0.557 0.419 0.414 0.562 1.000 > cor1<-comprank(Marozzi,"gini","midrank",print=FALSE) Tied scores are present in one or both rankings 1 2 Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 1 4 Tied scores are present in one or both rankings 1 5 Tied scores are present in one or both rankings 1 6 Tied scores are present in one or both rankings 1 7 Tied scores are present in one or both rankings 2 5 Tied scores are present in one or both rankings 2 6 Tied scores are present in one or both rankings 2 7 Tied scores are present in one or both rankings 3 5 Tied scores are present in one or both rankings 3 6 Tied scores are present in one or both rankings 3 7 Tied scores are present in one or both rankings 4 5 Tied scores are present in one or both rankings 4 6 Tied scores are present in one or both rankings 4 7 Tied scores are present in one or both rankings 5 6 Tied scores are present in one or both rankings 5 7 Tied scores are present in one or both rankings 6 7 > print(round(cor1$r,3)) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.000 0.408 0.326 0.418 0.404 0.168 0.251 [2,] 0.408 1.000 0.724 0.771 0.698 0.360 0.543 [3,] 0.326 0.724 1.000 0.605 0.514 0.383 0.614 [4,] 0.418 0.771 0.605 1.000 0.900 0.357 0.471 [5,] 0.404 0.698 0.514 0.900 1.000 0.343 0.452 [6,] 0.168 0.360 0.383 0.357 0.343 1.000 0.676 [7,] 0.251 0.543 0.614 0.471 0.452 0.676 1.000 > cor1<-comprank(Marozzi,"r4","wgh",print=FALSE) Tied scores are present in one or both rankings 1 2 Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 1 4 Tied scores are present in one or both rankings 1 5 Tied scores are present in one or both rankings 1 6 Tied scores are present in one or both rankings 1 7 Tied scores are present in one or both rankings 2 5 Tied scores are present in one or both rankings 2 6 Tied scores are present in one or both rankings 2 7 Tied scores are present in one or both rankings 3 5 Tied scores are present in one or both rankings 3 6 Tied scores are present in one or both rankings 3 7 Tied scores are present in one or both rankings 4 5 Tied scores are present in one or both rankings 4 6 Tied scores are present in one or both rankings 4 7 Tied scores are present in one or both rankings 5 6 Tied scores are present in one or both rankings 5 7 Tied scores are present in one or both rankings 6 7 > print(round(cor1$r,3)) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.000 0.934 0.802 0.897 0.907 0.376 0.640 [2,] 0.934 1.000 0.865 0.861 0.864 0.352 0.654 [3,] 0.802 0.865 1.000 0.680 0.642 0.285 0.748 [4,] 0.897 0.861 0.680 1.000 0.974 0.463 0.551 [5,] 0.907 0.864 0.642 0.974 1.000 0.495 0.555 [6,] 0.376 0.352 0.285 0.463 0.495 1.000 0.685 [7,] 0.640 0.654 0.748 0.551 0.555 0.685 1.000 > cor1<-comprank(Marozzi,"fy1","wgh",print=FALSE) Tied scores are present in one or both rankings 1 2 Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 1 4 Tied scores are present in one or both rankings 1 5 Tied scores are present in one or both rankings 1 6 Tied scores are present in one or both rankings 1 7 Tied scores are present in one or both rankings 2 5 Tied scores are present in one or both rankings 2 6 Tied scores are present in one or both rankings 2 7 Tied scores are present in one or both rankings 3 5 Tied scores are present in one or both rankings 3 6 Tied scores are present in one or both rankings 3 7 Tied scores are present in one or both rankings 4 5 Tied scores are present in one or both rankings 4 6 Tied scores are present in one or both rankings 4 7 Tied scores are present in one or both rankings 5 6 Tied scores are present in one or both rankings 5 7 Tied scores are present in one or both rankings 6 7 > print(round(cor1$r,3)) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.000 0.927 0.783 0.905 0.912 0.399 0.618 [2,] 0.927 1.000 0.850 0.862 0.855 0.370 0.624 [3,] 0.783 0.850 1.000 0.688 0.636 0.319 0.715 [4,] 0.905 0.862 0.688 1.000 0.980 0.451 0.561 [5,] 0.912 0.855 0.636 0.980 1.000 0.470 0.561 [6,] 0.399 0.370 0.319 0.451 0.470 1.000 0.696 [7,] 0.618 0.624 0.715 0.561 0.561 0.696 1.000 > cor1<-comprank(Marozzi,"fy2","wgh",print=FALSE) Tied scores are present in one or both rankings 1 2 Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 1 4 Tied scores are present in one or both rankings 1 5 Tied scores are present in one or both rankings 1 6 Tied scores are present in one or both rankings 1 7 Tied scores are present in one or both rankings 2 5 Tied scores are present in one or both rankings 2 6 Tied scores are present in one or both rankings 2 7 Tied scores are present in one or both rankings 3 5 Tied scores are present in one or both rankings 3 6 Tied scores are present in one or both rankings 3 7 Tied scores are present in one or both rankings 4 5 Tied scores are present in one or both rankings 4 6 Tied scores are present in one or both rankings 4 7 Tied scores are present in one or both rankings 5 6 Tied scores are present in one or both rankings 5 7 Tied scores are present in one or both rankings 6 7 > print(round(cor1$r,3)) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.000 0.792 0.624 0.821 0.808 0.327 0.479 [2,] 0.792 1.000 0.848 0.855 0.850 0.360 0.617 [3,] 0.624 0.848 1.000 0.681 0.629 0.312 0.714 [4,] 0.821 0.855 0.681 1.000 0.976 0.447 0.548 [5,] 0.808 0.850 0.629 0.976 1.000 0.466 0.550 [6,] 0.327 0.360 0.312 0.447 0.466 1.000 0.686 [7,] 0.479 0.617 0.714 0.548 0.550 0.686 1.000 > cor1<-comprank(Marozzi,"sbz","wgh",print=FALSE) Tied scores are present in one or both rankings 1 2 Tied scores are present in one or both rankings 1 3 Tied scores are present in one or both rankings 1 4 Tied scores are present in one or both rankings 1 5 Tied scores are present in one or both rankings 1 6 Tied scores are present in one or both rankings 1 7 Tied scores are present in one or both rankings 2 5 Tied scores are present in one or both rankings 2 6 Tied scores are present in one or both rankings 2 7 Tied scores are present in one or both rankings 3 5 Tied scores are present in one or both rankings 3 6 Tied scores are present in one or both rankings 3 7 Tied scores are present in one or both rankings 4 5 Tied scores are present in one or both rankings 4 6 Tied scores are present in one or both rankings 4 7 Tied scores are present in one or both rankings 5 6 Tied scores are present in one or both rankings 5 7 Tied scores are present in one or both rankings 6 7 > print(round(cor1$r,3)) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1 0.000 0.000 0.000 0.000 0.000 0.000 [2,] 0 1.000 0.697 0.745 0.693 0.328 0.501 [3,] 0 0.697 1.000 0.566 0.490 0.293 0.584 [4,] 0 0.745 0.566 1.000 0.897 0.341 0.459 [5,] 0 0.693 0.490 0.897 1.000 0.341 0.446 [6,] 0 0.328 0.293 0.341 0.341 1.000 0.601 [7,] 0 0.501 0.584 0.459 0.446 0.601 1.000 > detach(Marozzi) > > > > graphics::par(get("par.postscript", pos = 'CheckExEnv')) > cleanEx() > nameEx("qrank") > ### * qrank > > flush(stderr()); flush(stdout()) > > ### Name: qrank > ### Title: Quantiles of exact and approximated null distribution of rank > ### correlations > ### Aliases: qrank > > ### ** Examples > > ##### > # > a<-qrank(0.10,61,"Ke","St") > a<-qrank(0.01,25,"Sp","Ga",FALSE,FALSE);a$Cq [1] -0.474863771 > a<-qrank(0.03,11,"fy1","Ga",FALSE,FALSE);a$Cq [1] -0.594759161 > a<-qrank(0.03,11,"fy2","Ga",FALSE,FALSE);a$Cq [1] -0.594759161 > a<-qrank(0.03,11,"sbz","Ga",FALSE,FALSE);a$Cq [1] -0.421921475 > a<-qrank(0.001,15,"r4","Ex",FALSE,FALSE);cat(a$Cq,a$Lq,"\n") -0.07444 -0.07443 > a<-qrank(0.01,14,"fy2","Ex",FALSE,FALSE);cat(a$Cq,a$Lq,"\n") -0.06144 -0.06143 > ##### > # > > > > cleanEx() > nameEx("rankor") > ### * rankor > > flush(stderr()); flush(stdout()) > > ### Name: rankor > ### Title: Critical values and p-values of some rank correlations > ### Aliases: rankor > > ### ** Examples > > > data(Lambh);attach(Lambh) > rankor(DVI,Temp,"spearman","ga","gh",FALSE,"less",TRUE) Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: Spearman's rho -- Approximation: Gaussian Alternative : less Contin. Adjust. : FALSE n= 22 --- Observed value: -0.67701 Est. p-value= 0.000960 $n [1] 22 $Statistic [1] "Spearman's rho" $Value [1] -0.67701 $Approx [1] "Gaussian" $Tails [1] "less" $Cpv [1] 0.000959503 $Lpv [1] 0.000959503 > rankor(DVI,Temp,"r4","ga","gh",FALSE,"less",TRUE) Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: r4 -- Approximation: Gaussian Alternative : less Contin. Adjust. : FALSE n= 22 --- Observed value: -0.77304 Est. p-value= 0.000188 $n [1] 22 $Statistic [1] "r4" $Value [1] -0.77304 $Approx [1] "Gaussian" $Tails [1] "less" $Cpv [1] 0.00018826 $Lpv [1] 0.00018826 > rankor(DVI,Temp,"fy1","ga","gh",FALSE,"less",TRUE) Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: Fisher-Yates - means -- Approximation: Gaussian Alternative : less Contin. Adjust. : FALSE n= 22 --- Observed value: -0.73850 Est. p-value= 0.000357 $n [1] 22 $Statistic [1] "Fisher-Yates - means" $Value [1] -0.7385 $Approx [1] "Gaussian" $Tails [1] "less" $Cpv [1] 0.000356855 $Lpv [1] 0.000356855 > detach(Lambh) > ##### > # > data(Security);attach(Security) > rankor(AP,OV,"sp","ga","gh",FALSE,"greater",TRUE) Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: Spearman's rho -- Approximation: Gaussian Alternative : greater Contin. Adjust. : FALSE n= 40 --- Observed value: 0.70834 Est. p-value= 0.000005 $n [1] 40 $Statistic [1] "Spearman's rho" $Value [1] 0.70834 $Approx [1] "Gaussian" $Tails [1] "greater" $Cpv [1] 4.852e-06 $Lpv [1] 4.852e-06 > rankor(MP,AP,"fy2","ga","gh",FALSE,"greater",TRUE) Method for breaking ties: untied rankings Statistic: Fisher-Yates - medians -- Approximation: Gaussian Alternative : greater Contin. Adjust. : FALSE n= 40 --- Observed value: 0.77278 Est. p-value= 0.000001 $n [1] 40 $Statistic [1] "Fisher-Yates - medians" $Value [1] 0.77278 $Approx [1] "Gaussian" $Tails [1] "greater" $Cpv [1] 6.96e-07 $Lpv [1] 6.96e-07 > rankor(OV,MP,"sbz","ga","gh",FALSE,"greater",TRUE) Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: Symmetric BZ rank correlation -- Approximation: Gaussian Alternative : greater Contin. Adjust. : FALSE n= 40 --- Observed value: 0.77828 Est. p-value= 0.000000 $n [1] 40 $Statistic [1] "Symmetric BZ rank correlation" $Value [1] 0.77828 $Approx [1] "Gaussian" $Tails [1] "greater" $Cpv [1] 0 $Lpv [1] 0 > detach(Security) > ##### > # > data(Dizytwin);attach(Dizytwin) > op<-par(mfrow=c(1,1), mgp=c(1.8,.5,0), mar=c(2.8,2.7,2,1),oma=c(0,0,0,0)) > plot(Latitude,DZT_Rate,main="Latitude and dizygotic twinning rates",xlab="Latitude", + ylab="DZT_Rate", pch=19,cex=0.9, col= "tan4") > text(Latitude,DZT_Rate,labels=rownames(Dizytwin),cex=0.6,pos=c(rep(3,10),1,3,1, + rep(3,4),1.3)) > abline(v=mean(Latitude),col="black",lty=2,lwd=1) > abline(h=mean(DZT_Rate),col="darkblue",lty=2,lwd=1) > par(op) > rankor(Latitude,DZT_Rate,"r4","vg","gh",FALSE,"two",TRUE) Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: r4 -- Approximation: Vianelli-GGFR Alternative : two-sided Contin. Adjust. : FALSE n= 19 Observed value: 0.69125 Est. p-value= 0.002527 $n [1] 19 $Statistic [1] "r4" $Value [1] 0.69125 $Approx [1] "Vianelli-GGFR" $Tails [1] "two-sided" $Cpv [1] 0.002527138 $Lpv [1] 0.002527138 $Lambda [1] 1.50000000 5.09346445 > rankor(Latitude,DZT_Rate,"sp","ex","gh",FALSE,"two",TRUE) Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: Spearman's rho value: 0.67368 Approximation: Exact_level Alternative : two-sided Corr. Cont. FALSE Conservative p-value: 0.002066 Liberal p-value: 0.002140 $n [1] 19 $Statistic [1] "Spearman's rho" $Value [1] 0.67368 $approx [1] "Exact_level" $tails [1] "two-sided" $Cpv [1] 0.002066199 $Lpv [1] 0.00214008 > rankor(Latitude,DZT_Rate,"ke","ex","gh",FALSE,"two",TRUE) Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: Kendall's tau value: 0.50292 Approximation: Exact_level Alternative : two-sided Corr. Cont. FALSE Conservative p-value: 0.001862 Liberal p-value: 0.002433 $n [1] 19 $Statistic [1] "Kendall's tau" $Value [1] 0.50292 $approx [1] "Exact_level" $tails [1] "two-sided" $Cpv [1] 0.001861923 $Lpv [1] 0.002433118 > detach(Dizytwin) > ##### > # > # Bland, J.M. and Mutoka, C. and Hutt, M.S.R. ``Kaposi's sarcoma in Tanzania". East African > # Journal of Medical Research, 4, 47--53 (1977). > # Data from a study of the geographical distribution of a tumor, Kaposi's sarcoma, in > # mainland Tanzania. > Region<-c("Tabora","Iringa","Coast","Mara","Ruvuma","Mbeya","Shinyanga","Kigoma", + "Singida","Dodoma","Morogoro","Westlake","Arusha","Mtwara","Mwanza","Tanga", + "Kilimanjara") > Popw10kHF<-c(1.8, 2.6, 4.0, 4.4, 6.6, 6.7, 9.0, 9.2, 10.8, 11.1, 11.7, 12.5, 13.7, 14.8, + 20.7, 23.0, 57.3) > CasePMY<-c(2.37, 8.46, 1.28, 4.29, 7.21, 2.06, 1.66, 4.22, 6.17, 2.6, 6.33, 6.6, 2.46, + 6.4, 8.54, 4.54, 6.65) > op<-par(mfrow=c(1,1), mgp=c(1.8,.5,0), mar=c(2.8,2.7,2,1),oma=c(0,0,0,0)) > plot(Popw10kHF,CasePMY,main="",pch=19,cex=0.9,cex.lab=0.9,cex.axis=0.8,col="navy") > text(Popw10kHF,CasePMY,labels=Region,pos=3,cex=0.7) > abline(v=mean(Popw10kHF),col="black",lty=2,lwd=1) > abline(h=mean(CasePMY),col="darkblue",lty=2,lwd=1) > par(op) > rankor(CasePMY,Popw10kHF,"sp","ga","dubois",FALSE,"greater",TRUE) Method for breaking ties: untied rankings Statistic: Spearman's rho -- Approximation: Gaussian Alternative : greater Contin. Adjust. : FALSE n= 17 --- Observed value: 0.37745 Est. p-value= 0.065547 $n [1] 17 $Statistic [1] "Spearman's rho" $Value [1] 0.37745 $Approx [1] "Gaussian" $Tails [1] "greater" $Cpv [1] 0.065546743 $Lpv [1] 0.065546743 > ##### > # > ##### > # > ##### > # > data(Locre);attach(Locre) > op<-par(mfrow=c(1,1)) > plot(Males,Females,main="Fer cryin' out loud - there is a sex difference", + xlab="Females",ylab="Males",pch=19,cex=0.8,col="steelblue") > text(Males,Females,labels=1:length(Females),cex=0.7,pos=2) > abline(h=mean(Females),col="black",lty=2,lwd=1) > abline(v=mean(Males),col="black",lty=2,lwd=1) > par(op) > out<-rankor(Males,Females,"gi","vg","gh",FALSE,"greater") Method for breaking ties: gh [1] "p-values of rank correlation tests are only indicative in presence of ties." [1] "The continuity correction should not be considered with ties present." Statistic: Gini's cograduation coefficient -- Approximation: Vianelli-GGFR Alternative : greater Contin. Adjust. : FALSE n= 20 Observed value: 0.40000 Est. p-value= 0.015473 > cat(out$Value,out$Cpv,"\n") 0.4 0.015472773 > detach(Locre) > ##### > # > ##### > # > ##### > # > > > > graphics::par(get("par.postscript", pos = 'CheckExEnv')) > cleanEx() > nameEx("ranktes") > ### * ranktes > > flush(stderr()); flush(stdout()) > > ### Name: ranktes > ### Title: Tests of independence between two rankings > ### Aliases: ranktes > > ### ** Examples > > > # G. P. Watkins (1933). An Ordinal Index of Correlation, Journal of the > # American Statistical Association, 28:182, 139-151. > # 20-item series for area and density have been made up to cover the > # original 13 states and the four others earliest admitted to the Union. > State<-c("Georgia","North_Carolina","New_York","Luisiana","Pennsylvania", + "Virginia","Tennessee","Ohio","Kentucky","Maine","South_Carolina", + "West_Virginia","Maryland","Vermont","New_Hampshire","Massachusetts", + "New_Jersey","Connecticut","Delaware","Rhode_Island") > Area<-c(1,2, 3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20) > Density<-c(17,13,5,18,6,14,12,7,11,20,15,10,8,19,16,2,3,4,9,1) > op<-par(mfrow=c(1,1)) > plot(Area,Density,main="",xlab="Area",ylab="Density",pch=19,cex=0.9, + col="darkgreen" ) > abline(h=mean(Area),col="black",lty=2,lwd=1) > abline(v=mean(Density),col="darkblue",lty=2,lwd=1) > par(op) > r<-comprank(Area,Density,"fy2","wgh")$r > ranktes(r, length(Area), "fy2", "ga",FALSE, "two", TRUE) Statistic: Fisher-Yates - medians -- Approximation: Gaussian Alternative : two-sided Contin. Adjust. : FALSE n= 20 --- Observed value: -0.33341 Est. p-value= 0.146137 $n [1] 20 $Statistic [1] "Fisher-Yates - medians" $Value [1] -0.33341 $Approx [1] "Gaussian" $Tails [1] "two-sided" $Cpv [1] 0.146136596 $Lpv [1] 0.146136596 > ##### > # > ##### > data(Sharpe);attach(Sharpe) > op<-par(mfrow=c(1,1)) > plot(AVR,VAR, type = "p",pch=19,cex=1.1,col="tomato",main="Mutual fund + performance") > text(AVR,VAR, labels = rownames(Sharpe), cex=0.5, pos=3) > abline(h=mean(AVR),col="black",lty=2,lwd=1) > abline(v=mean(VAR),col="black",lty=2,lwd=1) > par(op) > r<-comprank(AVR,VAR,"sbz","wgh")$r > ranktes(r, length(AVR), "sbz", "st",FALSE, "greater", TRUE) Student-t approx. not yet available. An approximate p-value is computed Statistic: Symmetric BZ rank correlation -- Approximation: Gaussian Alternative : greater Contin. Adjust. : FALSE n= 34 --- Observed value: 0.68897 Est. p-value= 0.000000 $n [1] 34 $Statistic [1] "Symmetric BZ rank correlation" $Value [1] 0.68897 $Approx [1] "Gaussian" $Tails [1] "greater" $Cpv [1] 3.3e-08 $Lpv [1] 3.3e-08 > detach(Sharpe) > ##### > # > ##### > # > ##### > # > ##### > # > > > > graphics::par(get("par.postscript", pos = 'CheckExEnv')) > ### *