plot(c(-.5,1),c(0,5.5),pch=" ",xlab="log s", ylab="log n")
#title("Fraction Coil-state = for sigma=.001")
contour(log10(aaa.001$s),log10(aaa.001$n),aaa.001$th,add=T,draw=T,nlevels=7,levels=c(.1,.5,.9))
title(main=expression(paste("Fraction Coil-state = ",theta,", for ",sigma,"=.001")))
#
readline("? ")
#
persp(log10(aaa.001$s),log10(aaa.001$n),aaa.001$th,xlab="log s",ylab="log n", zlab="fraction coil",theta=-45,phi=45,scale=T,expand=.5)
title(main=expression(paste("Fraction Coil-state = ",theta," over n-s plane, for ",sigma,"=.001")))
#title("Fraction Coil-state over n-s plane, for sigma=.001")
#
readline("? ")
#
contour(log10(aaa.001$s),log10(aaa.001$n),log10(aaa.001$nu),nlevels=20, xlab="log s",ylab="log n")
contour(log10(aaa.001$s),log10(aaa.001$n),aaa.001$th,add=T,levels=c(.1,.9), draw=F, lty=3)
title(expression(paste("Log Coil-regions/chain = ",nu," over n-s plane"," for ",sigma,"=.001, with ",theta," overlay")))
#
readline("? ")
#
persp(log10(aaa.001$s),log10(aaa.001$n),aaa.001$dth,xlab="log s",ylab="log n", zlab="sharpness",theta=-20,expand=.5)
title(main = expression(paste("Sharpness = d",theta,"/dlns over n-s plane, for ",sigma,"=.001")))
#
readline("? ")
#
par(mfrow=c(3,1))
persp(log10(aaa.001$s),log10(aaa.001$n),aaa.001$th,xlab="log s",ylab="log n", zlab="fn coil",theta=-20,phi=20,scale=T,expand=.5)
title(main=expression(paste("Fraction Coil-state = ",theta," over n-s plane, for ",sigma,"=.001")))
persp(log10(aaa.01$s),log10(aaa.01$n),aaa.01$th,xlab="log s",ylab="log n", zlab="fn coil",theta=-20,phi=20,scale=T,expand=.5)
title(expression(paste(sigma,"=.01")))
persp(log10(aaa.999$s),log10(aaa.999$n),aaa.999$th,xlab="log s",ylab="log n", zlab="fn coil",theta=-20,phi=20,scale=T,expand=.5)
title(expression(paste(sigma,"=1")))
par(mfrow=c(1,1))
#
readline("? ")
#
par(mfrow=c(3,1))
persp(log10(aaa.001$s),log10(aaa.001$n),aaa.001$dth,xlab="log s",ylab="log n", zlab="sharpness",theta=-20,phi=20,scale=T,expand=.5)
title(main=expression(paste("Sharpness over n-s plane, for ",sigma,"=.001")))
persp(log10(aaa.01$s),log10(aaa.01$n),aaa.01$dth,xlab="log s",ylab="log n", zlab="sharpness",theta=-20,phi=20,scale=T,expand=.5)
title(expression(paste(sigma,"=.01")))
persp(log10(aaa.999$s),log10(aaa.999$n),aaa.999$dth,xlab="log s",ylab="log n", zlab="sharpness",theta=-20,phi=20,scale=T,expand=.5)
title(expression(paste(sigma,"=1")))
par(mfrow=c(1,1))
#
readline("? ")
#
contour(log10(aaa.001$s),log10(aaa.001$n),log10(matrix((aaa.001$data[,"theta"]*(aaa.001$data[,"n"]-3))/aaa.001$data[,"nu"],nrow=41,ncol=51)),nlevels=15,xlab="log s",ylab="log n")
contour(log10(aaa.001$s),log10(aaa.001$n),aaa.001$th,add=T,levels=c(.1,.9),draw=0,lty=3)
title(expression(paste("Log correlation length = <l> over n-s plane, for",sigma,"=.001")))