Notes:
Changes: Fixes: May-June 05 1) Fixed error in binned data analysis; create.model.frame was changing binned to numeric 0/1 instead of T/F because it was thrown into a matrix. This caused indexing errors in accessing data and it screwed up fitting for binned data with covariates. 2) Ordering of partials and v-c matrix of exponent and scale parameters in hazard rate didn't agree which created incorrect variance estimation for hazard rate. 3) If distance was not included in mrmodel for method=io.fi or trial.fi, the results were incorrect because the analytical integral for the logistic was not correct. 4) Made changes to gof, qqplot and plot routines to handle left truncation properly. Changes: May-June 05 1) Prevent use of "observer" in model for trial data. 2) plot and gof functions now all contain a breaks= field and use it if provided. 3) predict functions now all contain compute and newdata argument and use them as specified. 4) if an object is not detected by either observer an error is issued and the analysis stops 5) if fields other than observer or detected are different between observers a warning is issued that it may cause analysis to fail. 6) bins used for an analysis are automatically used in plot and gof without specifying in the arguments. A binned analysis is easily created with a call to make.bins prior to the analysis. Change: Aug 05 1) modified ddf.ds to restrict exponent of haz rate to be >1 2) modified gstdint and tablecgf to use a standardized set of grid points (exponentially spaced with more points closest to 0) for computing integrals to create a spline for interpolation. calls to tablecgf in integratedetfct and flt.lnl were modfied. 3) ddf.ds and flt were also changed such that the option showit now has 3 levels: 0 - no messages; 1- basic messages about refitting; 2 - all iterations are shown 4) Also in ddf.ds the parscale command in optim is used to resolve problems with scaling continuous covariates - need to test this more but appeared to be working well. It uses the max value of each column in the design matrix as the scalar for the parameter. This puts each parameter on the same footing in the optimization. 5) added functions model.table, collect.models and collect.model.names; these are presently not documented and aren't complete for all model types but does work for method="ds". It collects sets of results and builds a table for model selection. 6) dht.se was modified to return the component v-c matrices (vc1 for detection prob and vc2 for er). dht was modified to pass vc,vc1 and vc2 as a list in its return arguments. The list elements under vc are named total, detection and er. Also it computes the average detection probability and passes it back as average.p Change: Aug 05 1) Above changes integrated. 2) New structure to overall program. Now instead of detfct being a kind of pointer to another function we call it with options and it works out what it is. getpar() has been re-written to be less messy to use (1 call rather than 2). 3) Groundwork for adjustment terms put in.
Copyright © 2010 Geeknet, Inc. All rights reserved. Terms of Use