example(nlme)
fm2
intervals(fm2)
confint(fm2)
?confint.nlme
profile(fm2)
do()
do()
do()
do()
do()
do()
do()
xl
xu
yl
yu
summary(trapmat0
summary(trapmat)

summary(trapmat)
do()
do()
?do_admb
do(list(sigma=200,lam0=0.1))
do(list(sigma=200,lam0=0.1))
dim(trapmat)
trapmat
do(list(sigma=200,lam0=0.1))
do()
do()
unlink("spatialmr_admb_gen.tpl")
head(IndBYtrap)
head(indivs)
head(IndBYtrap)
head(traplocs)
head(traplocs)
names(i2)
head(i2)
head(traplocs)
head(x2)
aggregate(x2,by="Trap")
aggregate(x2,by="Trap",fun=sum)
head(aggregate(x2,by="Trap",FUN=sum))
(aggregate(x2,by="Trap",FUN=sum)
)
aggregate(x2,by=list("Trap"),FUN=sum)
)
aggregate(x2,by=list("Trap"),FUN=sum)
head(x2)
head(i2)
aggregate(i2,Trap~value,FUN=sum)
aggregate(Trap~value,data=i2,FUN=sum)
aggregate(value~Trap,data=i2,FUN=sum)
head(traptot)
head(x2)
head(trap.n)
head(i2)
head(trap.n)
head(trap.n)
?trap.n
?aggregate
head(trap.n)
tail(trap.n)
head(i2)
tail(i2)
head(trap.n)
?merge
head(x2)
str(traplocs[,-1])
str(traplocs)
str(read.csv("traplocs.csv"))
q()
n
?do_admb
detach("package:R2admb")
library(R2admb)
do_admb
q()
n
debug(do_admb)
debug(read_admb)
undebug(do_admb)
debug(read_pars)
nsdpar
sd_dat
ncorpar
cor_dat
corat
cormat
parnames
Q
q()
n
m1
debug(read_pars)
par_dat
par_dat
sd_dat
sd_dat
sd_dat[,2]
sd_dat[-(1:npar),3]
Q
q()
n
Sweave("R2admb.Rnw")
traceback()
options(error=recover)
2
pnames
gsub("r_","",tpldat$info$sdnums$vname)
tpldat$info$dnums
tpldat$info$sdnums
debug(read_psv)
read_psv(fn,names=pnames)
dim(ans)
names
ls()
Q
?mcmc.control
0
2
mcmc.opts$mcmcpars
tpldat$info
summary(ans)
colMeans(ans)
fn
fn <- "toy1"
paste(fn,"std",sep=".")
fn <- "toy1_gen"
   sdinfo <- read.table(paste(fn,"std",sep="."))
head(sdinfo)
pnames
rep_pars
sdinfo[,2]
split(parnames,factor(parnames))
pnames <- sdinfo[,2]
split(pnames,factor(pnames))
rep_pars
pnames
pnames
?grep
length(pnames)
pnames
grep("^r_.*",pnames,invert=TRUE)
grep("^r_.*",pnames,value=TRUE,invert=TRUE)
length(grep("^r_.*",pnames,value=TRUE,invert=TRUE))
Q
q()
n
