Script started on Thu 06 Sep 2007 08:22:36 PM EDT ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ fgrep Approval data | sed -e 's/^.*,//' | fgrep NA > approved ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ fgrep Approval data | sed -e 's/^.*,//' | fgrep NA > approvedpproveddpprovedipprovedspprovedapproved[1@D[1@i[1@s[1@a ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ wc -l *approved 1 approved 0 disapproved 1 total ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ cat less data [?1049h[?1h="","Vote","Status","CPI" "Finland","Abstention","P Member",9.6 "New Zealand","Disapproval","P Member",9.6 "Denmark","Disapproval","P Member",9.5 "Singapore","Approval with comments","P Member",9.4 "Switzerland","Approval with comments","P Member",9.1 "Norway","Disapproval","P Member",8.8 "Australia","Abstention","P Member",8.7 "Netherlands","Abstention","P Member",8.7 "Austria","Approval with comments","O Member",8.6 "Luxembourg","Abstention","O Member",8.6 "United Kingdom","Disapproval","P Member",8.6 "Canada","Disapproval","P Member",8.5 "Germany","Approval with comments","P Member",8 "Japan","Disapproval","P Member",7.6 "France","Disapproval","P Member",7.4 "Ireland","Disapproval","P Member",7.4 "Belgium","Abstention","P Member",7.3 "Chile","Abstention","O Member",7.3 "USA","Approval with comments","Secretariat",7.3 "Spain","Abstention","P Member",6.8 "Barbados","Approval",NA,6.7 "Portugal","Approval with comments","O Member",6.6 "Malta","Approval with comments","P Member",6.4 "Slovenia","Abstention","P Member",6.4 "Uruguay","Approval with comments","P Member",6.4 "United Arab Emirates","Approval",NA,6.2 "Qatar","Approval",NA,6 "Israel","Abstention","O Member",5.9 "Cyprus","Approval","P Member",5.6 "Jordan","Approval with comments",NA,5.3 data  [?1l>[?1049l]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ less data wc -l *approvedfgrep Disapproval data | sed -e 's/^.*,//' | fgrep NA > disapproved[1@A[1@p[1@p[1@r[1@o[1@v[1@a[1@lapprovedapprovedapproved[1@-[1@v[1@ ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ fgrep Approval data | sed -e 's/^.*,//' | fgrep -v NA > approved [1@D[1@i[1@s[1@adapprovediapprovedsapproved ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ wc =-l *approved 50 approved 18 disapproved 68 total ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ R R version 2.5.1 (2007-06-27) Copyright (C) 2007 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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. > v <- read.table("approved") [1@a > a <- read.table("approved") [1@ddapproved")iapproved")sapproved") > a V1 1 9.4 2 9.1 3 8.6 4 8.0 5 7.3 6 6.7 7 6.6 8 6.4 9 6.4 10 6.2 11 6.0 12 5.6 13 5.3 14 4.8 15 4.6 16 4.4 17 4.1 18 4.0 19 3.9 20 3.8 21 3.7 22 3.7 23 3.6 24 3.5 25 3.4 26 3.3 27 3.3 28 3.3 29 3.2 30 3.1 31 3.1 32 3.1 33 3.0 34 2.9 35 2.9 36 2.9 37 2.9 38 2.8 39 2.6 40 2.5 41 2.4 42 2.3 43 2.2 44 2.2 45 2.2 46 2.2 47 2.1 48 2.1 49 2.1 50 2.0 > d V1 1 9.6 2 9.5 3 8.8 4 8.6 5 8.5 6 7.6 7 7.4 8 7.4 9 5.1 10 4.8 11 4.6 12 3.6 13 3.3 14 3.3 15 3.3 16 2.7 17 2.5 18 2.3 > m1 = avg(mean(d) > m1 V1 5.716667 > m1 = mean(d)$)V)1) > m1 [1] 5.716667 > m1 = mean(d$V1) [1@0$V1)a$V1) > m0 [1] 4.116 > m1 [1] 5.716667 > m0 [1] 4.116 > n1 = count(V$V`1a$V1) Error: could not find function "count" > length(a$V1) [1] 50 > ma = mean(av$V1) > mbxbb = d = mean(d$V1) > na = lengtha(a$V1) > nd = length(a $vVd$V1) > n = ana + nd > y = c(a,d$V1,d$VC1) > y [1] 9.4 9.1 8.6 8.0 7.3 6.7 6.6 6.4 6.4 6.2 6.0 5.6 5.3 4.8 4.6 4.4 4.1 4.0 3.9 [20] 3.8 3.7 3.7 3.6 3.5 3.4 3.3 3.3 3.3 3.2 3.1 3.1 3.1 3.0 2.9 2.9 2.9 2.9 2.8 [39] 2.6 2.5 2.4 2.3 2.2 2.2 2.2 2.2 2.1 2.1 2.1 2.0 9.6 9.5 8.8 8.6 8.5 7.6 7.4 [58] 7.4 5.1 4.8 4.6 3.6 3.3 3.3 3.3 2.7 2.5 2.3 > n [1] 68 > std(y) Error: could not find function "std" > (ssdev(y) Error: could not find function "sdev" > var(s) Error in inherits(x, "data.frame") : object "s" not found > var(s))u)y)))y) [1] 5.168997 > var(y) svar(y) qvar(y) rvar(y)tvar(y)(var(y)) [1] 2.273543 > sqrt(var(y)) sqrt(var(y)) xsqrt(var(y))  sqrt(var(y)) sqrt(var(y)) sqrt(var(y)) ysqrt(var(y))  sqrt(var(y)) =sqrt(var(y))  sqrt(var(y))  > md -(m d - ma)/sy) [1@( * sqrt( (n1*n2)aa*nd) / (n (** (Nn - 1)) ) [1] 0.3129156 > ((md - ma)/sy) * sqrt( (na*nd) / (n * (n - 1)) ) [1@a[1@d [1] -0.3129156 > m,a Error: syntax error, unexpected ',', expecting '\n' or ';' in "m," > ma [1] 4.116 > md [1] 5.716667 > corr(x,y) Error: could not find function "corr" > cor(x.Y,y) Error in inherits(x, "data.frame") : object "x" not found > cor(x,y)$,y)V,y)1,y)$)V)1),)m)e)t)h)o)d)=)")P)e)a)r)s)o)n)") Error in match.arg(method) : 'arg' should be one of pearson, kendall, spearman > cor(x$V1,y$V1,method="Pearson")earson")pearson") Error in inherits(x, "data.frame") : object "x" not found In addition: Warning message: $ operator is deprecated for atomic vectors, returning NULL in: y$V1 > x Error: object "x" not found > axcor(x$V1,y$V1,method="pearson"),method="pearson") ,method="pearson") ,method="pearson") ,method="pearson") method="pearson") ,method="pearson") ,method="pearson") ,method="pearson") ,method="pearson") a,method="pearson") ,method="pearson") d,method="pearson")  Error in cor(x, y, na.method, method == "kendall") : incompatible dimensions > cor(a,d,method="pearson") quit function (save = "default", status = 0, runLast = TRUE) .Internal(quit(save, status, runLast)) > quit() Save workspace image? [y/n/c]: y ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ ls approved data disapproved typescript ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ cat aapawk '{prin t $1, " "\t" 1 }' approved 9.4 1 9.1 1 8.6 1 8 1 7.3 1 6.7 1 6.6 1 6.4 1 6.4 1 6.2 1 6 1 5.6 1 5.3 1 4.8 1 4.6 1 4.4 1 4.1 1 4 1 3.9 1 3.8 1 3.7 1 3.7 1 3.6 1 3.5 1 3.4 1 3.3 1 3.3 1 3.3 1 3.2 1 3.1 1 3.1 1 3.1 1 3 1 2.9 1 2.9 1 2.9 1 2.9 1 2.8 1 2.6 1 2.5 1 2.4 1 2.3 1 2.2 1 2.2 1 2.2 1 2.2 1 2.1 1 2.1 1 2.1 1 2 1 ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ awk '{print $1 "\t" 1 }' approved > approvedeverything ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ awk '{print $1 "\t" 1 }' approved > everything[1@0[1@d[1@i[1@s> everything ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ R R version 2.5.1 (2007-06-27) Copyright (C) 2007 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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. [Previously saved workspace restored] > e <- vtread.table("everything") > e V1 V2 1 9.4 1 2 9.1 1 3 8.6 1 4 8.0 1 5 7.3 1 6 6.7 1 7 6.6 1 8 6.4 1 9 6.4 1 10 6.2 1 11 6.0 1 12 5.6 1 13 5.3 1 14 4.8 1 15 4.6 1 16 4.4 1 17 4.1 1 18 4.0 1 19 3.9 1 20 3.8 1 21 3.7 1 22 3.7 1 23 3.6 1 24 3.5 1 25 3.4 1 26 3.3 1 27 3.3 1 28 3.3 1 29 3.2 1 30 3.1 1 31 3.1 1 32 3.1 1 33 3.0 1 34 2.9 1 35 2.9 1 36 2.9 1 37 2.9 1 38 2.8 1 39 2.6 1 40 2.5 1 41 2.4 1 42 2.3 1 43 2.2 1 44 2.2 1 45 2.2 1 46 2.2 1 47 2.1 1 48 2.1 1 49 2.1 1 50 2.0 1 51 9.6 0 52 9.5 0 53 8.8 0 54 8.6 0 55 8.5 0 56 7.6 0 57 7.4 0 58 7.4 0 59 5.1 0 60 4.8 0 61 4.6 0 62 3.6 0 63 3.3 0 64 3.3 0 65 3.3 0 66 2.7 0 67 2.5 0 68 2.3 0 > cor(e,$V1,e$V2) [1] -0.3129156 > cor(e$V1,e$V2) xcor(e$V1,e$V2)  cor(e$V1,e$V2) =cor(e$V1,e$V2)  cor(e$V1,e$V2) cor(e$V1,e$V2) cor(e$V1,e$V2)  x [1] -0.3129156 > x <- cor(e$V1,e$V2)  cor(e$V1,e$V2)  [1] -0.3129156 > cor(e$V1,e$V2)$ + p[pvalue] Error: syntax error, unexpected '[', expecting STR_CONST or SYMBOL in " cor(e$V1,e$V2)$" > [pvalue]  cor(e$V1,e$V2)$ [1] -0.3129156 > cor(e$V1,e$V2),)m)e)t)h)o)d)-)")))=)")s)e))p)e)a)r)m)a)n)") [1] -0.2863558 > cor(e$V1,e$V2,method="spearman")")")")")")")")")k")e")n")d")a")l")l") [1] -0.2375597 > cor(e$V1,e$V2,method="kendall")")")")")")")")p")e")a")r")s")o")n") [1] -0.3129156 > corhelpco cor Error: syntax error, unexpected SYMBOL, expecting '\n' or ';' in "help cor" > cor function (x, y = NULL, use = "all.obs", method = c("pearson", "kendall", "spearman")) { na.method <- pmatch(use, c("all.obs", "complete.obs", "pairwise.complete.obs")) method <- match.arg(method) if (is.data.frame(y)) y <- as.matrix(y) else stopifnot(is.atomic(y)) if (is.data.frame(x)) x <- as.matrix(x) else { stopifnot(is.atomic(x)) if (!is.matrix(x)) { if (is.null(y)) stop("supply both 'x' and 'y' or a matrix-like 'x'") x <- as.vector(x) } } if (method == "pearson") .Internal(cor(x, y, na.method, method == "kendall")) else if (na.method != 3) { Rank <- function(u) if (is.matrix(u)) apply(u, 2, rank, na.last = "keep") else rank(u, na.last = "keep") if (na.method == 2) { ok <- complete.cases(x, y) x <- if (is.matrix(x)) x[ok, ] else x[ok] if (!is.null(y)) y <- if (is.matrix(y)) y[ok, ] else y[ok] } x <- Rank(x) if (!is.null(y)) y <- Rank(y) .Internal(cor(x, y, na.method, method == "kendall")) } else { if (is.null(y)) { ncy <- ncx <- ncol(x) r <- matrix(0, nrow = ncx, ncol = ncy) for (i in seq(2, length = ncx - 1)) { for (j in seq_len(i - 1)) { x2 <- x[, i] y2 <- x[, j] ok <- complete.cases(x2, y2) x2 <- rank(x2[ok]) y2 <- rank(y2[ok]) r[i, j] <- .Internal(cor(x2, y2, na.method, method == "kendall")) } } r <- r + t(r) diag(r) <- 1 rownames(r) <- colnames(x) colnames(r) <- colnames(x) r } else { if (!is.matrix(x)) x <- matrix(x, ncol = 1) if (!is.matrix(y)) y <- matrix(y, ncol = 1) ncx <- ncol(x) ncy <- ncol(y) r <- matrix(0, nrow = ncx, ncol = ncy) for (i in 1:ncx) { for (j in 1:ncy) { x2 <- x[, i] y2 <- y[, j] ok <- complete.cases(x2, y2) x2 <- rank(x2[ok]) y2 <- rank(y2[ok]) r[i, j] <- .Internal(cor(x2, y2, na.method, method == "kendall")) } } rownames(r) <- colnames(x) colnames(r) <- colnames(y) r } } } > > corhelp cor  cor(e$V1,e$V2,method="pearson") r cor(e$V1,e$V2,method="pearson") e cor(e$V1,e$V2,method="pearson") p cor(e$V1,e$V2,method="pearson") o cor(e$V1,e$V2,method="pearson") r cor(e$V1,e$V2,method="pearson") t cor(e$V1,e$V2,method="pearson") 9 cor(e$V1,e$V2,method="pearson") ) cor(e$V1,e$V2,method="pearson")) cor(e$V1,e$V2,method="pearson")) (cor(e$V1,e$V2,method="pearson"))  Error: could not find function "report" > report(cor(e$V1,e$V2,method="pearson")) eport(cor(e$V1,e$V2,method="pearson")) port(cor(e$V1,e$V2,method="pearson")) ort(cor(e$V1,e$V2,method="pearson")) rt(cor(e$V1,e$V2,method="pearson")) t(cor(e$V1,e$V2,method="pearson")) (cor(e$V1,e$V2,method="pearson")) p(cor(e$V1,e$V2,method="pearson")) r(cor(e$V1,e$V2,method="pearson")) i(cor(e$V1,e$V2,method="pearson")) n(cor(e$V1,e$V2,method="pearson")) t(cor(e$V1,e$V2,method="pearson"))  [1] -0.3129156 > print(cor(e$V1,e$V2,method="pearson"))  [1] -0.3129156 > cor.prob <- function(X, dfr = nrow(X) - 2) { + R <- cor(X) + above <- row(R) < col(R) + r2 <- R[above]^2 + Fstat <- r2 * dfr / (1 - r2) + R[above] <- 1 - pf(Fstat, 1, dfr) + R + } > (c(cor.prob(e$#V1,$e$V2) Error in cor(X) : supply both 'x' and 'y' or a matrix-like 'x' > cor.prob(e) V1 V2 V1 1.0000000 0.00937335 V2 -0.3129156 1.00000000 > cor.prob(e)  x = c(1..2) Error: syntax error, unexpected NUM_CONST, expecting ',' in " x = c(1..2" > x = c(1..2)  c(1..2) < c(1..2) - c(1..2) Error: syntax error, unexpected NUM_CONST, expecting ',' in " x <- c(1..2" > x <- c(1..2)2)2):2) > x [1] 1 2 > x x <- c(1:2))1)0)0) > x <- c(1:100) [1@y[1@ :100)1:100)0:100)0:100))))2)0)0):200)1:200) > y <- c(101:200) x <- c(1:100) x x <- c(1:2)..2) = cor.prob(e))x),)y) Error in cor(X) : supply both 'x' and 'y' or a matrix-like 'x' > x [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 [19] 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 [37] 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 [55] 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 [73] 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 [91] 91 92 93 94 95 96 97 98 99 100 > y [1] 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 [19] 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 [37] 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 [55] 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 [73] 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 [91] 191 192 193 194 195 196 197 198 199 200 > matrix function (data = NA, nrow = 1, ncol = 1, byrow = FALSE, dimnames = NULL) { data <- as.vector(data) if (missing(nrow)) nrow <- ceiling(length(data)/ncol) else if (missing(ncol)) ncol <- ceiling(length(data)/nrow) x <- .Internal(matrix(data, nrow, ncol, byrow)) dimnames(x) <- dimnames x } > matrix (m <- matrix(c(c(x,y),(length(x)),2) > m [,1] [,2] [1,] 1 101 [2,] 2 102 [3,] 3 103 [4,] 4 104 [5,] 5 105 [6,] 6 106 [7,] 7 107 [8,] 8 108 [9,] 9 109 [10,] 10 110 [11,] 11 111 [12,] 12 112 [13,] 13 113 [14,] 14 114 [15,] 15 115 [16,] 16 116 [17,] 17 117 [18,] 18 118 [19,] 19 119 [20,] 20 120 [21,] 21 121 [22,] 22 122 [23,] 23 123 [24,] 24 124 [25,] 25 125 [26,] 26 126 [27,] 27 127 [28,] 28 128 [29,] 29 129 [30,] 30 130 [31,] 31 131 [32,] 32 132 [33,] 33 133 [34,] 34 134 [35,] 35 135 [36,] 36 136 [37,] 37 137 [38,] 38 138 [39,] 39 139 [40,] 40 140 [41,] 41 141 [42,] 42 142 [43,] 43 143 [44,] 44 144 [45,] 45 145 [46,] 46 146 [47,] 47 147 [48,] 48 148 [49,] 49 149 [50,] 50 150 [51,] 51 151 [52,] 52 152 [53,] 53 153 [54,] 54 154 [55,] 55 155 [56,] 56 156 [57,] 57 157 [58,] 58 158 [59,] 59 159 [60,] 60 160 [61,] 61 161 [62,] 62 162 [63,] 63 163 [64,] 64 164 [65,] 65 165 [66,] 66 166 [67,] 67 167 [68,] 68 168 [69,] 69 169 [70,] 70 170 [71,] 71 171 [72,] 72 172 [73,] 73 173 [74,] 74 174 [75,] 75 175 [76,] 76 176 [77,] 77 177 [78,] 78 178 [79,] 79 179 [80,] 80 180 [81,] 81 181 [82,] 82 182 [83,] 83 183 [84,] 84 184 [85,] 85 185 [86,] 86 186 [87,] 87 187 [88,] 88 188 [89,] 89 189 [90,] 90 190 [91,] 91 191 [92,] 92 192 [93,] 93 193 [94,] 94 194 [95,] 95 195 [96,] 96 196 [97,] 97 197 [98,] 98 198 [99,] 99 199 [100,] 100 200 > cor.prob(m) [,1] [,2] [1,] 1 0 [2,] 1 1 > cor.prob(m) m <- matrix(c(x,y),(length(x)),2) atrix yxcor.prob(x,y) xymatrix  <- matrix(c(x,y),(length(x)),2) cor.prob(m)  Save workspace image? [y/n/c]: lsn ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ ls approved data disapproved everything typescript ]0;abez@pinion: ~/projects/ooxmlabez@pinion:~/projects/ooxml$ vex data cat typescript Script done on Thu 06 Sep 2007 08:54:40 PM EDT