# # shepard_diagram_1.r -- Shepard Diagram Program # # Reads in output from KYST program -- DATA DIST DHAT # and makes the classic Shepard Diagram graph # showing Distances and DHATS against the # Observed Data # T <- matrix(scan("D:/R_Files/color_shepard.txt",0),ncol=3,byrow=TRUE) TT <- T nrow <- length(T[,1]) ncol <- length(T[1,]) # plot(T[,1],T[,2], xlab="Observed Similarity",ylab="Distances/D-Hats From KYST") mtext(side=3,line=1.5,"Shepard Diagram \nColor Data",font=2) points(T[,1],T[,3],pch=16) lines(T[1:nrow,1],T[1:nrow,3],lty=1,lwd=3) #lines(lowess(T[,2],T[,1],f=.2),lwd=3) text(60,1.5,"Circles are Distances \nDots are D-Hats \nLine Shows Monotone Regression")