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