source(coda_analyze.R)
gtwd()
getwd()
source("coda_analyze.R")
source("coda_analyze.R")
source("coda_analyze.R")
source("coda_analyze.R")
source("coda_analyze.R")
source("coda_analyze.R",echo="verbose")
source("coda_analyze.R",echo="true")
source("coda_analyze.R",echo=true)
source("coda_analyze.R",echo=True)
source("coda_analyze.R",echo=1)
source("coda_analyze.R",echo=TRUE)
source("coda_analyze.R",echo=TRUE,prompt=">")
source("coda_analyze.R",echo=TRUE,prompt=">>")
source("coda_analyze.R",echo=TRUE,prompt=">>")
fg
source("coda_analyze.R",echo=TRUE,prompt=">>")
FG
source("coda_analyze.R",echo=TRUE,prompt=">>")
source("coda_analyze.R",echo=TRUE,prompt=">>")
fg
 library(coda)
help(coda)
??coda
library(help=coda)
source("coda_analyze.R",echo=TRUE,prompt=">>")
data<-read.table(bayesallfile, header=T)
bfile<-scan(n=1)
bfile<-scan(n=1)
bayesallfile<-"bayesallfile"
data<-read.table(bayesallfile, header=T)
data
data.list<-split(data,data$Replicate)
data.list.1<-lapply(data.list,subset,select=c(9:34)
)
data.list.1<-lapply(data.list,subset,select=c(1:100))
data.list.1<-lapply(data.list,subset,select=c(1:7))
data.list.2<-lapply(data.list.1,mcmc)
data.list.3<-mcmc.list(data.list.2)
plot(data.list.3, ask=T)
q()
