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

[Sarah Goslee]

Hi,

The x value you want is lag from my.mgram$mgram

You can use lines() to add them to the plot.

Sarah

On Sun, Oct 16, 2011 at 7:10 PM, Nevil Amos <nevil.amos at gmail.com> wrote:
> 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

-- 
Sarah Goslee
http://www.functionaldiversity.org