npde
asks the user the name
and structure of the files containing the data, using
pdemenu
, while autonpde
takes these
variables and others as arguments.
autonpde(namobs, namsim, iid, ix, iy, imdv = 0, icens = 0, icov = 0, iipred = 0, boolsave = TRUE, namsav = "output", type.graph = "eps", verbose = FALSE, calc.npde=TRUE, calc.pd=TRUE, decorr.method = "cholesky", cens.method = "cdf", units = list(x="",y=""), detect=FALSE, ties=TRUE)
iid
,ix
and iy
below)npde.cens.method
) defaults to cdfnpde.decorr.method
) defaults to choleskyNpdeObject
npde
is an interactive function
whereas autonpde
takes all required input as
arguments.When the computation of npde fails because of numerical problems, error messages are printed out, then pd are computed instead and graphs of pd are plotted so that the user may evaluate why the computation failed.
The function also prints out the characteristics of the distribution of the npde (mean, variance, skewness and kurtosis) as well as the results of the statistical tests applied to npde. In addition, if boolsave is TRUE, two files are created:
npde.graphs
, gof.test
data(theopp)
data(simtheopp)
# Calling autonpde with dataframes
x<-autonpde(theopp,simtheopp,1,3,4,boolsave=FALSE)
x
# Calling autonpde with names of files to be read from disk
write.table(theopp,"theopp.tab",quote=FALSE,row.names=FALSE)
write.table(simtheopp,"simtheopp.tab",quote=FALSE,row.names=FALSE)
x<-autonpde(namobs="theopp.tab", namsim="simtheopp.tab", iid = 1,
ix = 3, iy = 4, imdv=0, boolsave = FALSE)
head(x["results"]["res"])
Run the code above in your browser using DataCamp Workspace