"VSS.plot" <-
function(x,title="Very Simple Structure",line=FALSE)
{ op <- par(no.readonly = TRUE) # the whole list of settable par's.
n=dim(x)
symb=c(49,50,51,52) #plotting sym
plot(x$cfit.1,ylim=c(0,1),type="b",ylab="Very Simple Structure Fit",xlab="Number of Factors",pch=49)
if (line) lines(x$fit)
title(main=title)
x$cfit.2[1]<-NA
x$cfit.3[1]<-NA
x$cfit.3[2]<-NA
x$cfit.4[1]<-NA
x$cfit.4[2]<-NA
x$cfit.4[3]<-NA
lines(x$cfit.2)
points(x$cfit.2,pch=50)
lines(x$cfit.3)
points(x$cfit.3,pch=symb[3])
lines(x$cfit.4)
points(x$cfit.4,pch=symb[4])
par(op)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.