### This analysis of the language test data.
### $Revision: 1.1.1.1 $
### $Date: 2005/01/06 14:58:12 $

###Parameters:

asp.user <- "ralmond"
asp.dbname <- "Language"

### Set things up
library("RMySQL")
asp.driver <- MySQL()
source("access.q")

###Connect to database
asp.dbConnect()

###List available tables
dbListTables(asp.con)

r1.model.name <- "lm_R1_1_SkillRead1"

###List available parameters
asp.listParameters(r1.model.name)

###Check number of chains
asp.getChains(paste(r1.model.name,"_Parameters",sep=""))

###Load parameters
paramNames <- asp.listParameters(r1.model.name)
r1.smp <- asp.getParameter(r1.model.name,paramNames)

postscript("reading1a.ps")
summary(r1.smp)
plot(r1.smp,ask=FALSE)
plot(r1.smp[[1]],ask=FALSE)
plot(r1.smp[[2]],ask=FALSE)
plot(r1.smp[[3]],ask=FALSE)
plot(r1.smp[[4]],ask=FALSE)
plot(r1.smp[[5]],ask=FALSE)
autocorr(r1.smp)
autocorr.plot(r1.smp,ask=FALSE)
gelman.diag(r1.smp)
gelman.plot(r1.smp,ask=FALSE)
dev.off()



###When done, disconnect
##asp.dbDisconnect()