# NOT RUN {
thStrt <- list(
list(beta=c(26.07,48808),sigsq=1.1573,lambda=0.33333333),
list(beta=c(23.48,32387),sigsq=1.8730,lambda=0.33333333),
list(beta=c(-0.0597,20760),sigsq=0.2478,lambda=0.33333333)
)
kfit <- mixreg(y ~ x,ncomp=3,data=kilnAoneOut,thetaStart=thStrt)
plot(kfit,pch=8,col=c("red","green","blue","black"))
plot(kfit,pch=8,col=c("green","blue","black"),cMeth="d")
plot(kfit,pch=8,col=c("green","blue","black"),cMeth="p")
afit <- mixreg(plntsInf ~ aphRel,data=aphids,ncomp=2)
plot(afit,cMeth="p",col=c("blue","red"),pch=20)
# Separates the points into two groups incredibly clearly!
npfit <- mixreg(plntsInf ~ 1, data=aphids,ncomp=2)
plot(npfit,cMeth="p",col=c("blue","red"),pch=20)
# }
Run the code above in your browser using DataLab