GoFy: Goodness-of-Fit in the forward direction (y)

In david-borchers/LT2D: Line transect estimation in two dimensions

Description

Calculates the Goodness-of-Fit in the forward direction (y)

Plot f(y) and forward distance distribution resulting from a call of fityx.

Calculates the Goodness-of-Fit in the forward direction (y)

Plot f(y) and forward distance distribution resulting from a call of fityx.

Usage

GoFy(fit, plot = FALSE, dotitle = FALSE)

GoFy(fit, plot = FALSE, dotitle = FALSE)

Arguments

fit	object resulting from a call of fityx.
plot=FALSE	boolean TRUE - display Q-Q plot.
fit	object resulting from a call of fityx.
plot=FALSE	boolean TRUE - display Q-Q plot.

Details

Calculates goodness-of-fit for forward distances and calculates p-values using kolomogarov and Cramer-von Mises.

Calculates goodness-of-fit for forward distances and calculates p-values using kolomogarov and Cramer-von Mises.

Value

data frame with these elements $p.cvm = Cramer-von Mises p-value. $D.kolomogarov = x value of Kolmogarov-distributed random variable $p.kolomogarov = kolomogarov p-value.

data frame with these elements $p.cvm = Cramer-von Mises p-value. $D.kolomogarov = x value of Kolmogarov-distributed random variable $p.kolomogarov = kolomogarov p-value.

See Also

fityx,

fityx p.kolomogarov GoFx

fityx,

fityx p.kolomogarov GoFx

Examples

## Not run: 
ystart=0.05;w=0.03
logphi=c(0.0180552, -4.4215995)
b=c(-23.725809,  -3.136638  , 2.122910)
N=200 #true number of animals
#generate some observations
simDat=simXY(N=N,pi.x=pi.norm,logphi=logphi,
hr=ip1,b=b,w=w,ystart=ystart)

x=simDat$locs$x; y=simDat$locs$y 
est.yx=fityx(y,x,b,hr=ip1,ystart,pi.x=pi.norm,logphi,w) 
plotfit.y(y,x,est.yx,nclass=10)
GoFy(fit=est.yx,plot=TRUE)

## End(Not run)
## Not run: 
ystart=0.05;w=0.03
logphi=c(0.0180552, -4.4215995)
b=c(-23.725809,  -3.136638  , 2.122910)
N=200 #true number of animals
#generate some observations
simDat=simXY(N=N,pi.x=pi.norm,logphi=logphi,
hr=ip1,b=b,w=w,ystart=ystart)

x=simDat$locs$x; y=simDat$locs$y
est.yx=fityx(y,x,b,hr=ip1,ystart,pi.x=pi.norm,logphi,w)
plotfit.y(y,x,est.yx,nclass=10)
GoFy(fit=est.yx,plot=TRUE)

## End(Not run)

david-borchers/LT2D documentation built on Aug. 17, 2020, 1:37 a.m.