A
on the left and the rotated version of B on the right. Below is a code fragment
showing how to do this:
#
par(mfrow=c(1,2))
#
# mar
# A numerical vector of the form c(bottom, left, top, right) which
# gives the number of lines of margin to be specified on the four
# sides of the plot. The default is c(5, 4, 4, 2) + 0.1.
# Note that you might want to make the right-hand side margin narrower so that the plots are closer together
# You may want to experiment with the parameters until the plot looks nice!
# This puts more white space
# on the Right-Hand-Side Margin
#
par(mar=c(4.1,5.1,4.1,5.1))
#
#
plot(dimension1,dimension2,type="n",asp=1,
main="",
xlab="",
ylab="",
xlim=c(-2.0,2.0),ylim=c(-2.0,2.0))
points(dimension1,dimension2,pch=16,col="red")
axis(1,font=2)
axis(2,font=2,cex=1.2)
# Main title
mtext("Morse Code Signals Lower Half Matrix",side=3,line=1.00,cex=1.2,font=2)
# x-axis title
mtext("Best Interpretation!!",side=1,line=2.75,cex=1.2)
# y-axis title
mtext("Best Interpretation!!",side=2,line=2.5,cex=1.2)
#
# pos -- a position specifier for the text. Values of 1, 2, 3 and 4,
# respectively indicate positions below, to the left of, above and
# to the right of the specified coordinates
#
namepos <- rep(4,nrow)
#
#
text(dimension1,dimension2,names,pos=namepos,offset=00.00,col="blue")
#
# Right hand Plot
#
plot(dimension11,dimension22,type="n",asp=1,
main="",
xlab="",
ylab="",
xlim=c(-2.0,2.0),ylim=c(-2.0,2.0))
points(dimension11,dimension22,pch=16,col="red")
axis(1,font=2)
axis(2,font=2,cex=1.2)
# Main title
mtext("Morse Code Signals Upper Half Matrix",side=3,line=1.00,cex=1.2,font=2)
# x-axis title
mtext("Best Interpretation!!",side=1,line=2.75,cex=1.2)
# y-axis title
mtext("Best Interpretation!!",side=2,line=2.5,cex=1.2)
#
# pos -- a position specifier for the text. Values of 1, 2, 3 and 4,
# respectively indicate positions below, to the left of, above and
# to the right of the specified coordinates
#
namepos <- rep(4,nrow)
#
#
text(dimension11,dimension22,names,pos=namepos,offset=00.00,col="blue")
Use the Orthogonal Procrustes Rotation function
from MCMCpack to do the rotation. You can
see how to do it by inspecting this R program:
smacof_nonmetric_nations_Procrustes_rotation.r
Turn in the resulting two panel plot and a listing of your R program.