nv
n
c
nv
fs
(fs-isize)/dsize
37000*8
Q
c
c
c
c
c
c
q()
n
getwd()
setwd("owls_mcmcdir0")
debug(read_psv)
debug(glmmADMB:::read_psv)
r <- read_psv("glmmadmb")
r <- glmmADMB:::read_psv("glmmadmb")
dim(ans)
glmmADMB:::read_psv
getwd()
debug(glmmadmb)
run
debug
list.files()
list.files()
list.files()
getwd()
setwd("..")
tmp
nrow(tmp)
rep_pars(tmp[,1])
rep_pars(tmp[,1])
tmp[,1]
tmp[,2]
rep_pars(tmp[,2])
rep_pars(tmp[,2])
Q
Q
q()
n
debug(glmmadmb)
library(glmmADMB)
q()
n
getwd()
names(test$mcmc)
colnames(test$mcmc)
getwd()
fit0
head(fit0$mcmc)
plotTrace(fit0)
library(scapeMCMC)
plotTrace(as.mcmc(fit0$mcmc))
mcmcControl(mcmc=50000)
getwd()
setwd("~/students/mccoy")
library(glmmADMB)
subset(X1,ttt=="LF")
head(s1)
library(ggplot2)
g_mussel5
ranef(g_mussel5)
?curve3d
curve3d
?colorRampPalette
?rgl
?rgl.material
pop3d()
pop3d()
pop3d()
pop3d()
pop3d()
pop3d()
pop3d()
pop3d()
rgl.off()
rgl.close()
?spheres3d
pop3d()
pop3d()
pop3d()
pop3d()
pop3d()
pop3d()
rgl.close()
q()
n
library(glmmADMB)
head(HPDinterval(m))
?geweke.diag
?geweke.diag
?coda
apropos(".diag$")
?heidel.diag
?autocorr.diag
?ls.diag
?apropose
?apropos
apropos(".diag$",where=TRUE)
search()
str(gg)
head(m)
t(apply(m,2,quantile,c(0.025,0.975)))
fit_zinbinom2_bs_mcmc <- glmmadmb(NCalls~(FoodTreatment+ArrivalTime)*SexParent+
                                     BroodSize+(1|Nest),
                                     data=Owls,
                                     zeroInflation=TRUE,
                                     family="nbinom1",
                                  mcmc=TRUE)
head(fit0)
head(fit_zinbinom2_bs_mcmc)
library(glmmADMB)
head(fit0$mcmc)
head(fit_zinbinom2_bs_mcmc$mcmc)
pnames
colnames(fit0$mcmc)
ls()
L <- load("owls_mcmcbatch.RData")
L
colnames(fit_zinbinom1_bs_mcmc)
colnames(fit_zinbinom1_bs_mcmc$mcmc)
tracePlot(fit_zinbinom1_bs_mcmc)
library(scapeMCMC)
plotTrace(fit_zinbinom1_bs_mcmc)
plotTrace(fit0)
plotTrace(fit0$mcmc)
plotTrace(fit_zinbinom1_bs_mcmc$mcmc)
dim(fit_zinbinom1_bs_mcmc$mcmc)
mcmctime
mcmctime[3]
mcmctime[3]/(60*60)
mcmctime[1]/(60*60)
mcmctime["elapsed"]/(60*60)
mcmctime["user"]/(60*60)
str(mcmctime)
mcmctime["user.self"]/(60*60)
mcmctime["user.child"]/(60*60)
?Sys.time
?system.time
?proc.time
sum(mcmctime[1:2])
str(sum(mcmctime[1:2]))
sum(mcmctime[1:2])
c(mcmctime)
mcmctime[4:5]
sum(mcmctime[4:5])
sum(mcmctime[4:5])/(60*60)
)
q()
n
