#\donttest{
mu.undPE2D(1.2)
mu.undPE2D(1.2,ugraph="r")
rseq<-seq(1.01,5,by=.05)
lrseq<-length(rseq)
muOR = muAND <- vector()
for (i in 1:lrseq)
{
muOR<-c(muOR,mu.undPE2D(rseq[i]))
muAND<-c(muAND,mu.undPE2D(rseq[i],ugraph="r"))
}
plot(rseq, muOR,type="l",xlab="r",ylab=expression(mu(r)),lty=1,
xlim=range(rseq),ylim=c(0,1))
lines(rseq,muAND,type="l",lty=2,col=2)
legend("bottomright", inset=.02,
legend=c(expression(mu[or](r)),expression(mu[and](r))),
lty=1:2,col=1:2)
#}
#\donttest{
asy.var.undPE2D(1.2)
asy.var.undPE2D(1.2,ugraph="r")
rseq<-seq(1.01,5,by=.05)
lrseq<-length(rseq)
avarOR<-avarAND<-vector()
for (i in 1:lrseq)
{
avarOR<-c(avarOR,asy.var.undPE2D(rseq[i]))
avarAND<-c(avarAND,asy.var.undPE2D(rseq[i],ugraph="r"))
}
oldpar <- par(mar=c(5,5,4,2))
plot(rseq, avarAND,type="l",lty=2,col=2,xlab="r",
ylab=expression(paste(sigma^2,"(r)")),xlim=range(rseq))
lines(rseq,avarOR,type="l")
legend(3.75,.02,
legend=c(expression(paste(sigma["underlying"]^"2","(r)")),
expression(paste(sigma["reflexivity"]^"2","(r)")) ),
lty=1:2,col=1:2)
par(oldpar)
#}
Run the code above in your browser using DataLab