[R] How to plot CI's (llim ulim) on ecodist mgram

Nevil Amos nevil.amos at gmail.com
Mon Oct 17 01:10:47 CEST 2011


I would like to put confidence intervals on a mantel corellogram
they are already calculated in the pmgram object but I am unsure how I 
get the x value in order to plot them?

package(ecodist)
X<-1:100
Y<-rnorm(1:100)
Z<-rnorm(1:100)
XY<-dist(data.frame(X,Y))
YX<-dist(data.frame(Y,X))
my.mgram<-mgram(XY,XZ)
plot(my.mgram)
print(my.mgram)
 > print(my.mgram)
$mgram
             lag ngroup      mantelr  pval        llim         ulim
  [1,]  3.770055    672  0.500012737 0.001  0.49689923  0.504301550
  [2,] 11.310165    691  0.383960457 0.001  0.38000201  0.387324434
  [3,] 18.850274    584  0.232086251 0.001  0.22670074  0.237501735
  [4,] 26.390384    587  0.114097397 0.001  0.10243901  0.122973735
  [5,] 33.930494    463 -0.003113351 0.835 -0.01928101  0.008839295
  [6,] 41.470603    468 -0.106354446 0.001 -0.12682280 -0.089539628
  [7,] 49.010713    357 -0.181250278 0.001 -0.20154017 -0.164863572
  [8,] 56.550823    348 -0.266397615 0.001 -0.28498271 -0.251134864
  [9,] 64.090933    252 -0.298705798 0.001 -0.31421396 -0.284154643
[10,] 71.631042    228 -0.353134525 0.001 -0.36468910 -0.341422330
[11,] 79.171152    147 -0.337181781 0.001 -0.34854961 -0.322161075
[12,] 86.711262    108 -0.334465576 0.001 -0.35500933 -0.309763543
[13,] 94.251371     43 -0.238965642 0.001 -0.26437371 -0.196038662

$resids
[1] NA

attr(,"class")
[1] "mgram"


Thanks

Nevil Amos



More information about the R-help mailing list