RDX2 X  .grid.started  .Random.seed ^b> .required MASS mass .Traceback min(x) range(..., na.rm = na.rm) range(fd) biplot.correspondence source >function (obj, type = c("symmetric", "rows", "columns"), ...) {  if (length(obj$cor) < 2) 2 stop("biplot is only possible if nf >= 2")  type <- match.arg(type)  X <- obj$rscore[, 1:2]  if (type != "columns") % X <- X %*% diag(obj$cor[1:2])  colnames(X) <- rep("", 2)  Y <- obj$cscore[, 1:2]  if (type != "rows") % Y <- Y %*% diag(obj$cor[1:2]) D switch(type, symmetric = biplot.default(X, Y, var.axes = FALSE, D ...), rows = biplot.bdr(X, Y, ...), columns = biplot.bdr(Y,  X, ...))  invisible() } obj type c symmetric rows columns ... { if < length $ cor@ stop "biplot is only possible if nf >= 2 <- match.arg X [ rscore :?@  != columns %*% diag cor?@ colnames rep @ Y cscore?@  rows cor?@ switch symmetric biplot.default var.axes   rows biplot.bdr  columns!  invisible: @function (x, y, var.axes = TRUE, col, cex = rep(par("cex"), 2), F xlabs = NULL, ylabs = NULL, expand = 1, xlim = NULL, ylim = NULL,  arrow.len = 0.1, ...) {  n <- nrow(x)  p <- nrow(y)  if (missing(xlabs)) { ! xlabs <- dimnames(x)[[1]]  if (is.null(xlabs))  xlabs <- 1:n  } xlabs <- as.character(xlabs) 0 dimnames(x) <- list(xlabs, dimnames(x)[[2]])  if (missing(ylabs)) { ! ylabs <- dimnames(y)[[1]]  if (is.null(ylabs)) & ylabs <- paste("Var", 1:p)  } ylabs <- as.character(ylabs) 0 dimnames(y) <- list(ylabs, dimnames(y)[[2]])  if (length(cex) == 1)  cex <- c(cex, cex)  if (missing(col)) {  col <- par("col")  if (!is.numeric(col)) ( col <- match(col, palette())  col <- c(col, col + 1)  }  else if (length(col) == 1)  col <- c(col, col) > unsigned.range <- function(x) c(-abs(min(x)), abs(max(x))) $ rangx1 <- unsigned.range(x[, 1]) $ rangx2 <- unsigned.range(x[, 2]) $ rangy1 <- unsigned.range(y[, 1]) $ rangy2 <- unsigned.range(y[, 2]) ( if (missing(xlim) && missing(ylim)) A xlim <- ylim <- rangx1 <- rangx2 <- range(rangx1, rangx2)  else if (missing(xlim))  xlim <- rangx1  else ylim <- rangx2 5 ratio <- max(rangy1/rangx1, rangy2/rangx2)/expand  on.exit(par(oldpar))  oldpar <- par(pty = "s") @ plot(x, type = "n", xlim = xlim, ylim = ylim, col = col[1], ...) 3 text(x, xlabs, cex = cex[1], col = col[1], ...)  par(new = TRUE) I plot(y, axes = FALSE, type = "n", xlim = xlim * ratio, ylim = ylim * 7 ratio, xlab = "", ylab = "", col = col[1], ...)  axis(3, col = col[2])  axis(4, col = col[2])  box(col = col[1]) < text(y, labels = ylabs, cex = cex[2], col = col[2], ...)  if (var.axes) ? arrows(0, 0, y[, 1] * 0.8, y[, 2] * 0.8, col = col[2],  length = arrow.len)  invisible() } x y  col cex par cex@ xlabs ylabs expand? xlim ylim arrow.len?   n nrow$ p0%  missing) ) [[ dimnames$?  is.null))?/) as.character)4$ list)34$@ 2* *34%? 5** paste Var?1*6*4%7*34%@  =='?' '' 2& &( col  ! is.numeric&& match& palette& & +&? 9&?& && unsigned.range function$  - abs min$B max$ (function(x) c(-abs(min(x)), abs(max(x))) rangx1?$? rangx2?$@ rangy1?%? rangy2?%@  &&2,2-,-EF rangeEF 2,,E-F ratio /DLGELHF+ on.exit( oldparN( pty s plot$ n,,--&&?  text$)''?&&? ( new P% axes  n, *,K-T-K xlab  ylab &&?  axis@&&@W@&&@ box&&?Q% labels*''@&&@   arrowsT%??陙T%@?陙&&@.# circles1.corresp?m)??i(@?Ж7@x2 ?0?%x&͠40Z0Zk0Zf0Zj&͠-0Zk0Zi0Zi0Zi?Jq70ZiÜ9?,KQֿ0Zm2Ӵ0ń0Zi!TX~!TX~!TX~?Cš }?Cš }*z?,KQ?ʱ2e?3r?%п0Zi?,KQ?,KQ?,KQ?,KQ?%?)W?,KQ?,KQֿMt?q8zٕ߼?%ӿ2?7Z!TX~ Zg>!TX~?Cš }?,KQֿ <ǿ9m?,KQ?,KQֿ!TX~?,KQֿEץ?,KQ?u d0Zi&͠-Lͱ9?,KQֿ0Zi?`T jO`::?%?,KQֿ)ÿ:?,KQֿnמ7Ri?%?,KQ?,KQֿ&͠-?,KQֿ&͠-&͠-6&￿J%}+?,KQֿ!TX~!TX~EM?,KQֿ!TX~!TX~?,KQֿ0Zi?,KQ?,KQ?,KQֿ&͠)?,KQֿ:!TX~!TX~?,KQ?,KQ?,KQ?,KQ?,KQ?,KQ?,KQ?,KQ?*;dhwX>HF?RSD@?6?OW^?OW?OW?OW??OW?OW?OW?OW?I$̎ npgnpg̎ ?OW̎ ̎ ̎ ?̎ ?Qпnpgnpg̎ ̎ ̎ ̎ ̎ ̎ ̎ ̎  dim n4n Ca1 Ca2 Ca3 Ca4 Cb5 Cb6 Cc7 Cc8 Cc9 Cc10 Cc11 Cc12 Cc13 Cc14 Cc15 Cc16 Cc17 Cc18 Cc19 Cc20 Cd21 Cd22 Cd23 Cd24 Cd25 Cd26 Cd27 Ce28 Ce29 Cf30 Cf31 Cf32 Cf33 Cf34 Cf35 Cf36 Cf37 Cf38 Cg39 Cg40 Ch41 Ch42 Ci43 Ci44 Cj45 Ck46 Ck47 Cm49 Cm50 Cm51 Cm52 Cm53 Cm54 Cm55 Cm56 Cm57 Cm58 Cn59 Cn60 Cn61 Cn62 Cn63 Cn64 Cn65 Cn66 Cn67 Cn68 Cn69 Cn70 Cn71 Co72 Co73 Co74 Co75 Co76 Co77 Co78 Co79 Cp80 Cq81 Cq82 Cq83 Cq84 Cq85 Cq86 Cq87 Cq88 Cq89 Cq90 Cq91 Cq92 Cq93 Cq94 Cq95 Cq96 Cq97 Cr98 Cr99 Cr100 Cr101 Cs102 Cs103 Cs104 Cs105 Cs106 Ct107 C108 C109 C110 C111QDA+2 ?읺&ۿuW=D?Q׋.c@?M^%)KٹTc?@ m; C?&4^?٥CX=Hc%ceA?F#CG?% /q4(g?o/}5?VX}?(ֺ1\  4 ENB Fiji Micro Manus Milne Morobe Sialum NC NI NWG Sogeri Tonga WNB????????????????@@?@??????????@?????????????@?????@????@???????????????????????????????????????????????????????????@???@???????@??@????@@?@?@??@?@????@???????????@@??@?@@??@???@??\ n 4n Ca1 Ca2 Ca3 Ca4 Cb5 Cb6 Cc7 Cc8 Cc9 Cc10 Cc11 Cc12 Cc13 Cc14 Cc15 Cc16 Cc17 Cc18 Cc19 Cc20 Cd21 Cd22 Cd23 Cd24 Cd25 Cd26 Cd27 Ce28 Ce29 Cf30 Cf31 Cf32 Cf33 Cf34 Cf35 Cf36 Cf37 Cf38 Cg39 Cg40 Ch41 Ch42 Ci43 Ci44 Cj45 Ck46 Ck47 Cm49 Cm50 Cm51 Cm52 Cm53 Cm54 Cm55 Cm56 Cm57 Cm58 Cn59 Cn60 Cn61 Cn62 Cn63 Cn64 Cn65 Cn66 Cn67 Cn68 Cn69 Cn70 Cn71 Co72 Co73 Co74 Co75 Co76 Co77 Co78 Co79 Cp80 Cq81 Cq82 Cq83 Cq84 Cq85 Cq86 Cq87 Cq88 Cq89 Cq90 Cq91 Cq92 Cq93 Cq94 Cq95 Cq96 Cq97 Cr98 Cr99 Cr100 Cr101 Cs102 Cs103 Cs104 Cs105 Cs106 Ct107 C108 C109 C110 C111 ENB Fiji Micro Manus Milne Morobe Sialum NC NI NWG Sogeri Tonga WNB names Row Column] cor rscore cscore Freq class correspondence circtab1????????????????@@?@??????????@?????????????@?????@????@???????????????????????????????????????????????????????????@???@???????@??@????@@?@?@??@?@????@???????????@@??@?@@??@???@??\ n 4n Ca1 Ca2 Ca3 Ca4 Cb5 Cb6 Cc7 Cc8 Cc9 Cc10 Cc11 Cc12 Cc13 Cc14 Cc15 Cc16 Cc17 Cc18 Cc19 Cc20 Cd21 Cd22 Cd23 Cd24 Cd25 Cd26 Cd27 Ce28 Ce29 Cf30 Cf31 Cf32 Cf33 Cf34 Cf35 Cf36 Cf37 Cf38 Cg39 Cg40 Ch41 Ch42 Ci43 Ci44 Cj45 Ck46 Ck47 Cm49 Cm50 Cm51 Cm52 Cm53 Cm54 Cm55 Cm56 Cm57 Cm58 Cn59 Cn60 Cn61 Cn62 Cn63 Cn64 Cn65 Cn66 Cn67 Cn68 Cn69 Cn70 Cn71 Co72 Co73 Co74 Co75 Co76 Co77 Co78 Co79 Cp80 Cq81 Cq82 Cq83 Cq84 Cq85 Cq86 Cq87 Cq88 Cq89 Cq90 Cq91 Cq92 Cq93 Cq94 Cq95 Cq96 Cq97 Cr98 Cr99 Cr100 Cr101 Cs102 Cs103 Cs104 Cs105 Cs106 Ct107 C108 C109 C110 C111 ENB Fiji Micro Manus Milne Morobe Sialum NC NI NWG Sogeri Tonga WNB cmdplot function(yleg=0.55){ Moldpar <- par(cex.axis=1.25, mgp=c(3,.5,0), mar=c(4.6,4.6,5.1,2.1), cex=1.25) on.exit(par(oldpar)) ,tricts<-pacific$District[-c(47,54,60,63,92)] tab <- -sort(-table(tricts)) lev<-names(tab) !tricts<-factor(tricts,levels=lev) ;eqscplot(pacific.cmdscale$points, xlab="MDS1", ylab="MDS2", H pch=(1:15)[as.numeric(tricts)], col=(1:15)[unclass(tricts)], cex=1.25) xran <- par()$usr[1:2] %xran <- xran+diff(xran)*c(-.035,.075) 'fourths <- xran[1] + (0:3)*diff(xran)/4 lev[12] <- "Micronesia" par(xpd=T) Qlegend(x=fourths[1], y=yleg, legend=lev[1:4], pch=1:4, bty="n", cex=0.6, col=1:4) Glegend(x=fourths[1], y=yleg, legend="", pch=1, bty="n", cex=0.6, col=1) Qlegend(x=fourths[2], y=yleg, legend=lev[5:8], pch=5:8, bty="n", cex=0.6, col=5:8) Tlegend(x=fourths[3], y=yleg, legend=lev[9:12], pch=9:12, bty="n", cex=0.6, col=9:12) Wlegend(x=fourths[4], y=yleg, legend=lev[13:15], pch=13:15, bty="n", cex=0.6, col=13:15) } yleg?ᙙ N( cex.axis? mgp @? mar @ffffff@ffffff@ffffff@'?M(N tricts pacific DistrictA @G@K@N@O@W tabA sortA tablee lev]he factore levelsk eqscplot pacific.cmdscale pointsU MDS1V MDS2 pch (?@. as.numerice&r?@. unclasse'? xran( usr?@u>uT diffu A?Q?333333 fourths>u?LTr@wu@k@( Micronesia( xpd T legend$x?%a{k?@q?@ bty n'?333333&?@{$x?%a{ q?| n'?333333&?{$x@%a{k@@ q@@ | n'?333333&@@ {$x@%a{k@"@(q@"@(| n'?333333&@"@({$x@%a{k@*@.q@*@.| n'?333333&@*@. compact1m]@@@@@"@&@(@*@,@0@1@2@3@4@5@6@:@<@=@>@?@@@@@A@B@B@C@D@E@E@F@G@H@I@I@J@P@@P@P@Q@R@S@S@@U@@U@U@V@V@@V@W@W@X@X@Y@Y@@Y@Y@Z@@[@@[@[@\@\@@\@\@]@^@^@@^@^@_@_@@_@_@`@`@@``@`@`@`@`@a@a@b@b@@b`@b@b@c @c@@c`@c@d] Sogeri Sogeri Sogeri Sogeri Sogeri Sogeri Sogeri Sogeri Sogeri Sogeri Sogeri Central Milne Bay Milne Bay Milne Bay Milne Bay Milne Bay Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe (Sialum) Morobe Manus Province Manus Province Manus Province Manus Province West New Britain West New Britain West New Britain West New Britain East New Britain East New Britain East New Britain New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland New Ireland Bougainville Province Bougainville Province Bougainville Province Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Northwest Guadalcanal Fiji Fiji Fiji Fiji Fiji Fiji Fiji Fiji Fiji Tonga Micronesia New Caledonia] E P P E E P P E P P P P E E E E P P P P P P P P P P P P P P P E P P P P E E E E E E E E E E E E E P E P P E E E E E E E E P P P E E E E E E E E E E E E E E E E E E E E E E E E E P E E E]????????]?????????]???????????????]????????????]??]???????????]?????]?]????]???????????????]???????????????????]?????????????]??????????????]????]??]?]?]?]?]????????????]????]?]??]??]???]???????????]??]?]?]?]??????????]?????]???]?????????]??]???????]????]?]??]?]??????????????????????]?????????]???????????????]???]???]????]????]????????????????]????????????????]???]?????????????]?]??????]??]?????]???]??]?????]?]?]?????????????????]?????????????????????????????????]???]?]??????]?????????????????]????????????????]???????????????]??????????????]??]????????????]??]?????]?]?]?]?]?]?]????]?]?]??]??????????????????]?]???????????]?]?]?]?]?]??????????]????]???]??????]????]???????????????]??]??????????????]????]???????]??]???]??]?]?]m Site.No. District Technique Ca4 Cb5 Cb6 Cc7 Cc8 Cc9 Cc10 Cc11 Cc12 Cc13 Cc14 Cc15 Cc16 Cc17 Cc18 Cc19 Cc20 Cd21 Cd22 Cd23 Cd24 Cd25 Cd26 Cd27 Ce28 Ce29 Cf30 Cf31 Cf32 Cf33 Cf34 Cf35 Cf36 Cf37 Cf38 Cg39 Cg40 Ch41 Ch42 Ci43 Ci44 Cj45 Ck46 Ck47 Cl48 Cm49 Cm50 Cm51 Cm52 Cm53 Cm54 Cm55 Cm56 Cm57 Cm58 Cn59 Cn60 Cn61 Cn62 Cn63 Cn64 Cn65 Cn66 Cn67 Cn68 Cn69 Cn70 Cn71 Co72 Co73 Co74 Co75 Co76 Co77 Co78 Co79 Cp80 Cq81 Cq82 Cq83 Cq84 Cq85 Cq86 Cq87 Cq88 Cq89 Cq90 Cq91 Cq92 Cq93 Cq94 Cq95 Cq96 Cq97 Cr98 Cr99 Cr100 Cr101 Cs102 Cs103 Cs104 Cs105 Cs106 Ct107 C108 C109^ data.frame row.names] 2 6 7 7.1 9 11 12 13 14 16 17 18 19 20 21 22 26 28 29 30 31 32 33 35 36 37 39 40 42 43 45 46 49 50 51 52 65 66 67 71 75 76 77 85 86 87 88 89 90 94 95 96 99 100 101 102 103 105 109 110 111 112 113 114 115 119 120 121 122 123 124 125 126 127 128 130 131 132 133 134 135 141 143 144 146 147 148 150 153 154 155 159 160 compact1.cmdscale?`Z?he??_?E+?p█?:y$.9+&D?{#?^얱?\P?'"C3!z^U3[,?/0Ղ+?}0?#꽘?O#?Mcy?]; -?Mv#[?6u1Ɏ?. a?z20,C?fh8V?e?.UD?/?f3:?ygl?g_))?p?\rhCݒ?U>Oe?7Jщ*.50H