=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= DENDROCHRONOLOGY PROGRAM LIBRARY Update number 99 10 FEB 1998 =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= SOME NEW FEATURES OF THE PROGRAMS * [This README.TXT file is updated when a new version of one or more programs or routines becomes available. It describes new features in the programs, with the most recent changes first.] =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= DISCLAIMER > The following disclaimer applies to the Dendrochronology Program Library, Program ARSTAN, all programs in the Dendroecology Program Library, and any other software distributed by myself directly, via Anonymous FTP, or by any other means. > All software has received only limited testing and neither the author(s), the Laboratory of Tree-Ring Research, nor the University of Arizona can accept responsibility for its compatibility with other computer systems or user needs. The user must judge its efficacy for the intended use. > Although every effort is made to ensure that the software functions properly, no guarantee is made to this effect and we cannot be responsible for any problems with it, including adaptation of the software to any particular computer system. > The software is made available free of charge. Any payment accepted is for actual costs of diskettes, copying, printing and/or postage. =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= * If you find problems or have suggestions about any of the programs, please contact: Richard Holmes, c/o Laboratory of Tree-Ring Research, University of Arizona, Tucson, Arizona 85721 USA E-mail: RHolmes@LTRR.Arizona.EDU * If you have detected a problem, it is possible that it has already been found and the corrected version has been made available. Make sure you have the latest versions of the programs; see ANONYMOUS FTP below. =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= * In Programs ARSTAN and CRONOL ring indices will not be negative, and statistical tables will be printed correctly regardless of option selected. [10FEB98] * In Program SCRPLT (graphics on the screen) placement is more accurate of the lines at zero, the mean, and plus and minus one standard deviation. [07FEB98] * The negative exponential curve fit in Program ARSTAN and CRONOL may go negative, that is, the curve may have a negative asymptote. The user may override this provision and revert to linear regression fit. This is in Option 2, "Detrending Method", when curve fit type "1" is selected. When the series of measurements is divided by the curve, however, a constant is added to raise both series to all positive values, the division is performed, and the constant is subtracted again. [02FEB98] * Program ARSTAN tree summaries summaries are truncated if required to include only continuous time spans; they do not include voids between non-continuous radial series. [29JAN98] * Program DPL/YUX (produces spreadsheet files in columns) now runs faster [29JAN98] * Unnecessary lines in Program OUTBREAK are not printed. [29JAN98] * Some changes were made to the code, invisible to the user, to make it more efficient, and all programs were recompiled. [22JAN98] * DPL.FMT now allows you to fit a spline curve of rigidity proportional to the length of the series. For example, a rigidity of "-67" will fit a spline of 50% frequency response of 67% of the series length. You may also save both spline and indices in a single run. [01DEC97] * The negative exponential curve fit in Programs ARSTAN, CRONOL and DPL/CRN was modified so that the constant which adjusts the level of the exponential curve will not be negative, that is, the curve will never have a negative asymptote. If it is negative, the programs will default to a linear regression or line through the mean. In ARSTAN the user may override this provision and allow negative asymptotes. This is in Option 2, "Detrending Method". (But see later note on ARSTAN above.) [29JUN97] * The date computation in all programs was corrected so that there will be no problem when the year 2000 arrives. Version was set to 3.00 in most programs. [14JUL97] * The algorithm in Programs ARSTAN, CRONOL and DPL/CRN for computing the negative exponential curve fit was extensively modified and dramatically improved by Ed Cook of the Tree Ring Laboratory at Lamont-Doherty Earth Observatory. The new routine is now included in the public versions of ARSTAN, CRONOL and DPL/CRN. [03JUN97] * A minor problem in calculating the mean correlations in DPL/MAT has been corrected. [08JUL97] * Programs ARSTAN, CRONOL and DPL/CRN produce a new tab-delimited file with extension ".COE" for use in graphics programs. This file contains the coefficients of each deterministic curve fit to the measurement series. [22MAY97] * In Programs ARSTAN, CRONOL and DPL/CRN the routine which fits the modified negative exponential curve did not check for curves that have a negative asymptote. This has now been corrected. A linear regression line will be fit instead of a negative exponential when the asymptote is negative. Program ARSTAN checks for indices over 4.0 and lists them by series; this is one tool you may use to verify the curve fit. [08APR97] * Program ARSTAN prints the total number of rings in the data set and number and percent absent, on the first page of the output, as well as at the head of the chronology listing. [07MAR97] * In running COFECHA (DPL/COF) on a PC and selecting the parts to print, the program occasionally caused an error; this has been corrected. [27FEB97] * In COFECHA (DPL/COF), if you have a file of undated series, you may now list the adjustments of segments in order of the adjustment as well as by highest correlation. [07FEB97] * Program CRONOL (DPL/CRN), you may now detrend by subtracting the curve values from the measurements as well as by dividing the measurements by the curve values. [07FEB97] * Bar plotting routines DPL/BAR and DPL/BPG now allow the user to select the rigidity of the spline fit to the series before plotting. * If your data are in ITRDB Index format, you may select a new output format in Program DPL/FMT, with three columns: year, value and number of samples. [31OCT96] * You may now set your favorite default colors for text and screen. You may have one color combination during execution of programs and another when you are outside them. Use either DPL and type (~), or use program SETCLR. If the file ANSI.SYS is not activated in your CONFIG.SYS file, use SETCLR and respond with a slash (/). All programs were recompiled today. [24OCT96] * A new version of COFECHA, called COF10K, is now available. It will work with time spans of up to 10240 years in length. [08OCT96] * Program EVENT (Superposed Epoch Analysis) now will handle a maximum window size of 102 years -- those before plus those after plus key year. [11SEP96] * Program CLR may be used to change the colors and brightness of the text and the screen on the monitor. [03SEP96] * ITRDB ("Tucson") measurement format may now be successfully read by the programs even when the file has one or more end-of-file marks on the last line. These are ASCII 26 or Control-Z characters (^Z), and are inserted by some editing and file-transfer software. [03SEP96] * ITRDB ("Tucson") index or chronology format may now be successfully read by the programs even if it varies slightly from the current ITRDB standard in two respects. (1) If the date of the first decade line of the data gives the decade year rather than the actual first year of the chronology: WILSTD17609990 09990 09990 09990 09990 09990 09990 09990 01048 11023 1 Date=^^^^ First data=\_____/ The ITRDB standard specifies that the actual first year of the chronology should appear as the date of the dirst decade: WILSTD17689990 09990 09990 09990 09990 09990 09990 09990 01048 11023 1 Date=^^^^ First data=\_____/ (2) If the last year of the chronology is a year ending in "9", with the data in the last place on the line, a following line may be omitted, if it is also the last data line in the file: WILSTD1980 940 26 854 261422 261643 26 940 26 854 261422 261643 26 940 26 854 26 The ITRDB standard specifies that the actual last year of the chronology should be followed by filling out the line with null values ("9990 0") even if this means inserting an entire null line: WILSTD1980 940 26 854 261422 261643 26 940 26 854 261422 261643 26 940 26 854 26 WILSTD19909990 09990 09990 09990 09990 09990 09990 09990 09990 09990 0 Null data=\_____/\_____/\_____/\_____/\_____/\_____/\_____/\_____/\_____/\_____/ [16AUG96] * Programs CRONOL and DPL/CRN now finish correctly the file for printing. Data files have been correct, but the file for printing was not quite complete. [13MAY96] In the common interval analysis Program ARSTAN shows the signal strength of samples of different size, and compares it to a hypothetical infinite ("population") sample, and the theoretical number of trees needed to attain levels of signal strength equal to 75%, 80%, 85% and 90% of this infinite sample. [24APR96] * Mean sensitivity is now calculated even for series that contain negative values, according to the formula: MS = SUM(I=2,N)/(N-1) 2. * |(Y(I) - Y(I-1))| / ( |Y(I)| + |Y(I-1)| ) if ( |Y(I)| + |Y(I-1)| ) is non-zero Where MS is mean sensitivity of the series Y is the time series of values N is the number of values in the series With the caveat: If there are any negative values in a series, the decimal point for mean sensitivity is printed as an asterisk. In such cases mean sensitivity cannot strictly be calculated correctly. [24APR96] * ARSTAN now gives the mean correlation of radii vs mean as one of the statistics on the last page, "Summary of Chronology Statistics". [24APR96] * Program ARSTAN has yet another protection against dividing by zero. When there are two or more identical time series, their correlation is 1.0 and the Fisher Z-transform then would require a divide by zero. When two series are found to be identical a warning is printed and the correlation is set to 0.99999 for purposes of the statistical analysis. [07MAR96] * Program DPL/FMT now allows you to set the first or last date of all series to a user-specified year. [05MAR96] * Program SCRPLT (High-resolution plots on screen) will plot two sets of time series even when the pairs do not coincide in time span. SCRPLT also handles all data formats used in the DPL, and shows in the lower left corner of the screen the sequential series number as well as the identification and time span. * Program RESPO (correlation and response function) now permits you to save the table of correlation and response function coefficients in a file ready for importation into a spreadsheet program. The file extension is .COL * A new program, EDTRIM, will copy time series from one file to another, allowing you to trim data values from either or both ends, or to copy the series intact. The new file will have the same name with the last letter changed to `Z`, or if the original name ends in `Z` the last letter will be `Y`. Each series is plotted in high resolution on the screen. Touch `Y` to copy the series without changes, or to see the first and last twenty indices in the series and trim values from the start or end of the series. Touching will cause no change to be made. A log file is kept, noting changes made to the data. * In opening files two new convenient features are added. The program spells out the format it finds for the data file; you may correct it if it is incorrect. After seeing the first eight lines of the data file, you may now close the file with a backslash (\) and open another file as desired. Data series in Compact format may be interspersed among those in ITRDB Measurement or Index format, if the first series is not in Compact format. * In reading about algorithms I found that sophisticated and efficient sorting algorithms can on rare occasions cause problems by not sorting correctly. The recommendation is to use the simpler routines except in the case of extremely large data series, when a more efficient (but less safe) algorithm may be used. For this reason the more efficient routines were replaced on 5 February 1996 by the simpler and safer routines. * DPL/MAT, the correlation matrix routine, now gives you the option of marking the correlations over 99% and 95% significance (default) or marking those above a user-supplied level, as before. Mean correlations for each variable are now printed for matrices with two files as well as for those with one file as before. * EDRM has been added as a stand-alone version of DPL/EDT, the routine to edit ring measurement files. * DPL/YUX now eliminates spaces between columns (unless you have specified a space as the column delimiter). This yields a smaller file which will work better in Excel and other graphics programs. * Responding at any prompt with a reverse apostrophe (`) [ASCII 96] will allow you to give one DOS command while still running the program. You may thus do a Directory or any DOS other command, including running other programs. This is to save you trouble with file names, etc., in the midst of running a program. After the DOS command you are returned to the same prompt for a response. You may use this feature to access the ITRDB Program Library help (if installed on your PC) by giving the DOS command ITRHELP. On normal termination of a program the screen will be returned to default colors -- gray on black. If you quit with a slash the screen will also be cleared. * All routines give you a choice of combinations of screen colors for color monitors, most of which are also OK for monochrome monitors. Default for mainframe versions is bright green on black with yellow typing. Default for PC versions is gray on blue with yellow typing. On a monochrome monitor this becomes color on black with bold typing. Type a tilde (~) at any prompt to change screen colors. Use option X if ANSI.SYS is not enabled on your PC. * Program ARSTAN will process up to 95 ring measurement series covering a time span of up to 4096 years. If you have specified that AR modeling be done on each series to its own order, and the pooled autoregression order is zero, ARSTAN will now go ahead and make a residual chronology At the bottom of the tables of statistics for the original measurements, detrended indices, and AR-modeled series a line appears showing the means for each column. =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= ANONYMOUS FTP The Dendrochronology Program Library and programs in the Dendroecology Program Library are in the ANONYMOUS FTP directory /SYS/HOME/ANONYMOU/PUB/DPL They are not compressed, but alive and well and ready to work. Each time you unload programs, be sure to use the most updated version of the memory manager program DOSXMSF.EXE. Instructions below show how to obtain the programs. =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= HOW TO GET PROGRAMS VIA ANONYMOUS FTP Following are instructions to obtain the latest versions of computer programs available from the Tree-Ring Laboratory of the University of Arizona. First is a table of computer prompts, your responses and an explanation. After the table is a more complete record of a typical FTP session. TRANSFERRING EXECUTABLE FILES TO YOUR PC Be sure to set your FTP client to binary mode to transfer executable files to your PC. Put executable files in a subdirectory listed in the PATH statement in the AUTOEXEC.BAT file. You may then run any program simply by typing the name, e.g.: "DPL" or "EVENT" -------------------- ------------------------- --------------------------- Prompt_by_computer Your_input Explanation -------------------- ------------------------- --------------------------- $> [Your computer] ftp kati.ltrr.arizona.edu Connect to server KATI.LTRR.ARIZONA.EDU> login anonymous Anonymous login Password: rholmes@ltrr.arizona.edu Your E-mail address KATI.LTRR.ARIZONA.EDU> binary Set to BINARY transfer KATI.LTRR.ARIZONA.EDU> cd pub Change to PUBLIC directory KATI.LTRR.ARIZONA.EDU> cd dpl Change to DPL directory KATI.LTRR.ARIZONA.EDU> mget dosxmsf.exe Get Memory manager * KATI.LTRR.ARIZONA.EDU> mget dpl.exe Get Program Library KATI.LTRR.ARIZONA.EDU> mget arstan.exe Get Program ARSTAN KATI.LTRR.ARIZONA.EDU> mget [any other programs] Get other programs KATI.LTRR.ARIZONA.EDU> bye Exit connection to server FTP> exit [may vary by system] Exit FTP * This is a new version of the memory manager, required for any version of Programs DPL, ARSTAN and others dated after 24 May 1994. =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Following is a nearly complete record of a typical FTP session. A blank line comes before a line where you respond, and the remaining lines are printed on the screen by the computer. Your response follows the symbol ">" or in the case of the password, ":". -------------------------------------------------------------------------------- $> ftp kati.ltrr.arizona.edu KATI.LTRR.ARIZONA.EDU>login anonymous binary Type: Image, Structure: File, Mode: Stream KATI.LTRR.ARIZONA.EDU>cd pub cd dpl mget dosxmsf.exe mget dpl.exe mget arstan.exe bye exit =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=