.Rhistory

fit
fit
fit
fit
fit
fit
fit
U48i
fit
traceback()
traceback()
list(...)
UPb
weightedmean(UPb)
traceback()
q()
n
out$mean
out$mean
alpha
fit$sigma
Q
out$mean
nfact(alpha/2)*out$mean['s[x]']
nfact(alpha)
Q
dev.copy2pdf('/home/pvermees/Desktop/randomeffects.pdf')
dev.copy2pdf(file='/home/pvermees/Desktop/randomeffects.pdf')
dev.copy2pdf(file='/home/pvermees/Desktop/ludwig.pdf')
fract
i
j
valid
fract
i
j
valid
valid
valid
valid
fract
out$fract
valid
i
j
fract
i
j
Q
10:3
pat
plat
Q
i
j
valid
valid
i
j
i
j
j
j
i
valid
j
j
i
i
j
fract
i
j
valid
i
j
Y
sY
rep(TRUE,j-i+1)
i
j
X1,1
X[1,1]
Y[1]
YsY[1,1]
YsY[1,2]
random.effects
valid
Y
Q
fit
imin
fit
X
sX
Q
Q
valid
valid
i
i
j
Q
i
j
i
j
ns
Q
Y
sY
valid
rep(TRUE,j-i+1)
valid
Q
X
sX
fit
mu
prob
ns
Q
X
X
valid
prob
mu
prob
fit
prob
ns
minp
X
Q
minp
mu
sigma
Q
X
sX
valid
valid
X
sX
mu
sigma
minp
ns
Q
mu
valid
Q
minp
X
sigma
prob
X
mu
sigma
sX
X
sX
sigma
mu
prob
Q
X
sX
valid
prob
mu
sigma
Q
X
sX
valid
prob
mu
sigma
Q
sY
YsY
i
j
Q
X
sX
Q
tt
Q
tt
Q
traceback
traceback()
q()
n
traceback()
d
search.range
intersection.misfit.york(search.range[1],a=a,b=b,d=d)
intersection.misfit.york(search.range[2],a=a,b=b,d=d)
intersection.misfit.york(1,a=a,b=b,d=d)
midpoint
a
b
Q
intercept
slope
plot(xy[,'X'],xy[,'Y'],type='p')
lines(range(xy[,'X'],intercept+slope~*xy[,'X'])
lines(range(xy[,'X'],intercept+slope*xy[,'X'])
)
lines(range(xy[,'X'],intercept+slope*range(xy[,'X'])))
lines(range(xy[,'X']),intercept+slope*range(xy[,'X']))
ta0bo
ta0b0
q()
n
q()
n
q()
n
common.Pb
Q
common.Pb
Q
traceback()
traceback()
type
tt<cutoff.76
tt
tt
Q
tt
k
fit
Q
k
k
fit
Q
k
k
Q
k
i2i
x
xommon.Pb
common.Pb
common.Pb
Q
common.Pb
Q
traceback()
traceback()
Q
tt
tijd
xy
Q
?quote
?quote
fit
fit
fit
fit$w
w
wl
wu
ll
us
ul
Q
traceback()
traceback()
J11
J12
J21
J22
E11
J22
E12
J21*J21*E11
J21*J22*E12
J21*J22*E12
J22*J22*E22
J22
E22
J21
J22
Q
J11
J12
J21
x
Q
plot(y[,1],y[,2])
data2york(x,inverse=FALSE)
Y <- data2york(x,inverse=FALSE)
plot(Y[,1],Y[,2])
plot(Y[,1],Y[,3])
plot(y[,1],y[,3])
yi
fit$a
fit
1/fit$a[1]
yi
foo <- c(1,2)
subset(foo,select=1)
subset(foo,subset=1)
subset(foo,select=1)
?select
?subset
?subset
?I
foo <- c(1,2)
foo[,1]
subset
foo
subset(foo,subset=c(TRUE,FALSE))
initial
dat[,'Y']
Q
y
fit
yi
y
plot(y[,'X'],y[,'Y'])
Q
x
y
Q
initial
dat[,'sY']
dat
scatterplot(dat)
Q
scatterplot(dat)
dev.new()
scatterplot(dat)
Q
scatterplot(dat)
dev.new()
scatterplot(dat)
Q
x.ess.source('/home/pvermees/Dropbox/Programming/R/IsoplotRtest.R', visibly = FALSE, output = TRUE)
x
y
yi
fit$b
Q
.ess.source('/home/pvermees/Dropbox/Programming/R/IsoplotRtest.R', visibly = FALSE, output = TRUE)
Q
.ess.source('/home/pvermees/Dropbox/Programming/R/IsoplotRtest.R', visibly = FALSE, output = TRUE)
Q
.ess.source('/home/pvermees/Dropbox/Programming/R/IsoplotRtest.R', visibly = FALSE, output = TRUE)
dat
scatterplot(dat)
initial$y0
dat
Q
q()
n
scatterplot(y)
Q
k.ess.source('/home/pvermees/Dropbox/Programming/R/IsoplotRtest.R', visibly = FALSE, output = TRUE)
Q
.ess.source('/home/pvermees/Dropbox/Programming/R/IsoplotRtest.R', visibly = FALSE, output = TRUE)
k
k
k
Q
q()
n
k
k
fit
x
k
x
BIC
newBIC
newBIC
newBIC
newBIC
traceback()
hess
hess
k
hess
k
hess
solve(hess)
Q
hess
hess
blockinverse(Au,Bu,t(Bu),Cu,doall=TRUE)
hess
    blockinverse(Au,Bu,t(Bu),Cu,doall=TRUE)
hess
    blockinverse(Au,Bu,t(Bu),Cu,doall=TRUE)
Q
traceback()
k
BIC_fit(x,5,log=log,type=type,cutoff.76=cutoff.76,
 cutoff.disc=cutoff.disc,detritus=detritus,commonPb=common.Pb)
Q
pii
piu
blockinverse(Au,Bu,t(Bu),Cu,doall=TRUE)
Q
out
piu
k
k
k
piu
aiu
biu
dist(aiu)
Q
q()
n
out$mean
fit$mean
Q
tt
R
get.KCa.age(R[1],R[2],exterr=TRUE)
Q
DP
sDP
tt
st
exterr
DP
sDP
Q
R
tt
st
Q
tt
st
exterr
R
sR
bratio
Q
tt
st
L
L[2]/L[1]
sR/R
sR
exterr
exterr
st/tt
sR/R
sDP/DP
st/tt
exterr
exterr
L[2]/L[1]
sDP/DP
Q
q()
n
Q
errconvert(X,'other',format,ierr)
errconvert(X,'other',format,ierr)
format
ierr
X
errconvert(X,'other',format,ierr)
Q
i
gc
format
ierr
Q
x
i
Q
i
Q
format
Q
format
Q
i
q()
n
traceback()
X
out
X$x
Q
q()
n