For Duplicates - Generate Mean Values

[Hirindu Kawshala]

To xtset, I want to remove duplicates. I have duplicates with the same GVKEY and datadate but with different tickers. I want to merge those duplicates with the same tickers together and for the variables, I want to have a mean value for those duplicates only. To do so, how can I complete or revise the following codes. Please help me with this. Thank you.

Error message:
. xtset GVKEY datadate
repeated time values within panel
r(451);


Codes:
sort GVKEY datadate
quietly by GVKEY datadate: gen dup = cond(_N == 1, 0, _n)


destring GVKEY, replace
xtset GVKEY datadate
gen returns = (prccd - L.prccd) / L.prccd



----------------------- copy starting from the next line -----------------------

* Example generated by -dataex-. To install: ssc install dataex
clear
input long(datadate GVKEY) str8 tic str10 cusip str70 conm
15342 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15343 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15344 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15347 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15348 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15349 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15350 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15351 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15354 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15355 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15356 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15357 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15358 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15362 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15363 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15364 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15365 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15368 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15369 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15370 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15371 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15372 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15375 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15376 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15377 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15378 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15379 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15382 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15383 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15384 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15385 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15386 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15390 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15391 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15392 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15393 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15396 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15397 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15398 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15399 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15400 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15403 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15404 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15405 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15406 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15407 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15410 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15411 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15412 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15413 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15414 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15417 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15418 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15419 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15420 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15421 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15424 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15425 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15426 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15427 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15431 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15432 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15433 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15434 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15435 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15438 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15439 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15440 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15441 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15442 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15445 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15446 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15447 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15448 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15449 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15452 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15453 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15454 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15455 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15456 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15459 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15460 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15461 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15462 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15463 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15466 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15467 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15468 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15469 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15470 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15473 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15474 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15475 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15476 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
15477 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15480 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15481 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15482 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15483 1076 "PRG"   "74319R101" "PROG HOLDINGS INC"
15484 1076 "AAN.2" "002535201" "PROG HOLDINGS INC"
end
format %td datadate

------------------ copy up to and including the previous line ------------------

[Clyde Schechter]

collapse (mean) prccd, by(GVKEY datadaate)
xtset GVKEY datadate
gen returns = (prccd-L.prccd)/L.prccd

[Hirindu Kawshala]

Thank you Clyde Schechter. It worked.

[Clyde Schechter]

egen tag = tag(GVKEY datadate)
by GVKEY datadate, sort: egen mean_ret = mean(prccd)
replace mean_ret = . if !tag

Note: In the future, please be sure your example data includes all of the variables necessary to implement the solution to your problem. Your data example in this thread lacks the prccd variable. The example also contains no instances of duplicate observations for any combination of GVKEY and datadate. To be truly helpful, the example data needs to exhibit the problem that needs to be solved and contain all relevant variables.

[Hirindu Kawshala]

Thank you Clyde Schechter.
These are my total codes based on daily CRSP data. My purpose is to calculate ILLIQUIDITY and TOTAL_VOLATILITY based on daily data and collapse them into means to convert quarterly and merge them with quarterly dataset using GVKEY and datadate. Please can you please correct my codes.

CODES:
sort GVKEY datadate
quietly by GVKEY datadate: gen dup = cond(_N == 1, 0, _n)

collapse (mean) prccd, by(GVKEY datadate fyrc)
destring GVKEY, replace
xtset GVKEY datadate
gen returns = (prccd-L.prccd)/L.prccd

gen fiscal_date = datadate
replace fiscal_date = fiscal_date + (12 - fyrc)*30 if month(datadate) > fyrc
replace fiscal_date = fiscal_date - fyrc*30 if month(datadate) <= fyrc

gen fyearq = year(fiscal_date)
gen fqtr = quarter(fiscal_date)

drop if fyearq < 2001 | fyearq > 2023
egen quarter_id = group(fyearq fqtr)

gen daily_return = (prccd - prcod) / prcod
gen abs_return = abs(daily_return)
gen amihud_numerator = abs_return / cshtrd

egen quarterly_amihud_numerator = sum(amihud_numerator), by(GVKEY quarter_id)
egen quarterly_volume = sum(cshtrd), by(GVKEY quarter_id)
gen ILLIQUIDITY = quarterly_amihud_numerator / quarterly_volume
collapse (mean) ILLIQUIDITY, by(GVKEY quarter_id)

collapse (sd) sd_returns_by_fqtr = returns (first) fyearq fqtr, by(GVKEY quarter_id)
gen log_sd_returns_by_fqtr = 2*log(sd_returns_by_fqtr)
gen log_sd = log(sd_returns_by_fqtr)
rename log_sd TOTAL_VOLATILITY

[Hirindu Kawshala]

For the above codes, I revised them as follows. Can you please confirm whether my below codes are correct (I will merge based on GVKEY and quarter_id):

gen fiscal_date = datadate
replace fiscal_date = fiscal_date + (12 - fyrc)*30 if month(datadate) > fyrc
replace fiscal_date = fiscal_date - fyrc*30 if month(datadate) <= fyrc

gen fyearq = year(fiscal_date)
gen fqtr = quarter(fiscal_date)

drop if fyearq < 2001 | fyearq > 2023
egen quarter_id = group(fyearq fqtr)

gen daily_return = (prccd - prcod) / prcod
gen abs_return = abs(daily_return)
gen amihud_numerator = abs_return / cshtrd

egen quarterly_amihud_numerator = sum(amihud_numerator), by(GVKEY quarter_id)
egen quarterly_volume = sum(cshtrd), by(GVKEY quarter_id)
gen ILLIQUIDITY = quarterly_amihud_numerator / quarterly_volume

egen sd_returns_by_fqtr = sd(daily_return), by(GVKEY quarter_id)
gen TOTAL_VOLATILITY = log(sd_returns_by_fqtr)

collapse (mean) ILLIQUIDITY TOTAL_VOLATILITY, by(GVKEY quarter_id)

tostring GVKEY, replace
replace GVKEY = "0" + GVKEY if length(GVKEY) == 5
replace GVKEY = "00" + GVKEY if length(GVKEY) == 4
replace GVKEY = "000" + GVKEY if length(GVKEY) == 3
replace GVKEY = "0000" + GVKEY if length(GVKEY) == 2
replace GVKEY = "00000" + GVKEY if length(GVKEY) == 1

[Clyde Schechter]

I'm sorry, but I cannot critique this code for several reasons. First, it contains several variables that are not in your example data and whose meaning I do not understand. You also calculate new variables amihud and illiquidity. As for amihud, I have no idea what that is. As for illiquidity, I have a general understanding of the term, but I do not know how it is supposed to be related to the variables you are dealing with. I don't understand what the manipulations you are doing with the fiscaldate variable are trying to accomplish. I don't work in finance.

So you would need to post a data example that includes all of the variables (except the ones calculated de novo in the code) mentioned, along with an explanation of what each one is (as the meanings of several of the variable names are obscure, at least to me). And it would be best if it were reviewed by somebody who is familiar with the formal definitions of illiquidity and amihud in finance.

[Hirindu Kawshala]

Thank you, Clyde Schechter, and sorry for the confusion. I need to calculate quarterly figures for ILLIQUIDITY and TOTAL_VOLATILITY based on daily data and then merge them with a quarterly dataset. This is the code I want to construct: I want to know the correctness of my codes (specially red colour codes)

STEP1
use "D: Daily_data.dta"

destring GVKEY, replace

gen fiscal_date = datadate
replace fiscal_date = fiscal_date + (12 - fyrc)*30 if month(datadate) > fyrc
replace fiscal_date = fiscal_date - fyrc*30 if month(datadate) <= fyrc

gen fyearq = year(fiscal_date)
gen fqtr = quarter(fiscal_date)

drop if fyearq < 2001 | fyearq > 2023
egen quarter_id = group(fyearq fqtr)

Calculate variables using daily data for illiquidity
sort GVKEY datadate
gen volume_dollars = cshtrd * prccd
gen daily_return = D.prccd / L.prccd - 1
gen daily_illiquidity = abs(daily_return) / volume_dollars

Calculate variables using daily data for total volatility
egen sd_returns_by_fqtr = sd(returns), by(GVKEY quarter_id)

collapse (mean) ILLIQUIDITY = daily_illiquidity (sd) sd_returns_by_fqtr, by(GVKEY quarter_id)
gen TOTAL_VOLATILITY = 2*log(sd_returns_by_fqtr)

Keep ILLIQUIDITY TOTAL_VOLATILITY GVKEY quarter_id

STEP2
use "D: Quarterly_data.dta"
merge 1:1 GVKEY quarter_id using "D:\Daily_data.dta"

[Hirindu Kawshala]

I need to calculate quarterly figures for ILLIQUIDITY and TOTAL_VOLATILITY based on daily data and then merge them with a quarterly dataset. Please correct my codes. I tried the above codes in another way:

STEP1
use "D: Daily_data.dta"

destring GVKEY, replace

gen fiscal_date = datadate
replace fiscal_date = fiscal_date + (12 - fyrc)*30 if month(datadate) > fyrc
replace fiscal_date = fiscal_date - fyrc*30 if month(datadate) <= fyrc

gen fyearq = year(fiscal_date)
gen fqtr = quarter(fiscal_date)

drop if fyearq < 2001 | fyearq > 2023
egen quarter_id = group(fyearq fqtr)

sort GVKEY datadate
gen volume_dollars = cshtrd * prccd


sort GVKEY datadate
quietly by GVKEY datadate:gen dup=cond(_N==1,0,_n)
drop if dup>1
drop dup

tsset GVKEY datadate
gen daily_return = D.prccd / L.prccd - 1
gen daily_illiquidity = abs(daily_return) / volume_dollars

egen sd_returns_by_fqtr = sd(daily_return), by(GVKEY quarter_id)
save "D:\Daily_Data1.dta", replace


use "D:\Daily_Data1.dta"
count if missing(daily_illiquidity)
drop if missing(daily_illiquidity)
collapse (mean) ILLIQUIDITY = daily_illiquidity, by(GVKEY quarter_id)
save "D:\Daily_Data2.dta", replace


use "D:\Daily_Data1.dta"
collapse (sd) sd_returns_by_fqtr = returns (first) fyearq fqtr, by(GVKEY quarter_id)
gen TOTAL_VOLATILITY = 2*log(sd_returns_by_fqtr)
save "D:\Daily_Data3.dta", replace


use "D:\Daily_Data3.dta"
merge 1:1 GVKEY quarter_id using "D:Daily_Data2.dta"
save "D:\Daily_Data4.dta", replace

STEP2
use "D: Quarterly_data.dta"
merge 1:1 GVKEY quarter_id using "D:\Daily_Data4.dta"

[Clyde Schechter]

I explained in #7 why I can't review the correctness of this code. Your posts in #8 and #9 do not respond to any of the obstacles I have pointed out, so I still can't help you.

[Hirindu Kawshala]

Clyde Schechter. Please disregard my posts 6 & 8 because I changed my variables. According to post 9, I want to compute two main variables from a daily dataset and derive QUARTERLY values to merge with a quarterly dataset. Please review my codes in post 9 only.

Variables:
1. ILLIQUIDITY




2. TOTAL_VOLATILITY
The Quarterly standard deviation of stock returns.

Data Example:
----------------------- copy starting from the next line -----------------------

* Example generated by -dataex-. To install: ssc install dataex
clear
input long(GVKEY datadate) double(cshtrd prccd prcod) float(quarter_id volume_dollars daily_return daily_illiquidity sd_returns_by_fqtr)
1004 15342  58000   9.2 . 3  533600          .             .  .02689223
1004 15343 122800  9.66 . 3 1186248       -.95  8.008443e-07  .02689223
1004 15344 110100  9.84 . 3 1083384  -.9813665  9.058344e-07  .02689223
1004 15347  53500   9.7 . 3  518950          .             .  .02689223
1004 15348  86100  9.85 . 3  848085  -.9845361 1.1608932e-06  .02689223
1004 15349  97200  9.85 . 3  957420         -1 1.0444737e-06  .02689223
1004 15350  42500   9.7 . 3  412250 -1.0152284  2.462652e-06  .02689223
1004 15351 115100  9.55 . 3 1099205  -1.015464  9.238167e-07  .02689223
1004 15354 155200  9.09 . 3 1410768          .             .  .02689223
1004 15355  83300  8.92 . 3  743036 -1.0187019 1.3709994e-06  .02689223
1004 15356 150800  8.28 . 3 1248624 -1.0717489   8.58344e-07  .02689223
1004 15357  56500  8.49 . 3  479685  -.9746377 2.0318287e-06  .02689223
1004 15358  44400  8.17 . 3  362748 -1.0376914   2.86064e-06  .02689223
1004 15362 105200  8.15 . 3  857380          .             .  .02689223
1004 15363  57000  8.28 . 3  471960  -.9840491 2.0850264e-06  .02689223
1004 15364  39600  8.25 . 3  326700 -1.0036232 3.0720025e-06  .02689223
1004 15365  46800  8.16 . 3  381888 -1.0109091  2.647135e-06  .02689223
1004 15368  27200  8.22 . 3  223584          .             .  .02689223
1004 15369 155700  8.23 . 3 1281411  -.9987835  7.794404e-07  .02689223
1004 15370  61000  8.25 . 3  503250  -.9975699  1.982255e-06  .02689223
1004 15371 129400   8.6 . 3 1112840  -.9575757  8.604792e-07  .02689223
1004 15372  97800  8.45 . 3  826410 -1.0174419 1.2311588e-06  .02689223
1004 15375  42200  8.22 . 3  346884          .             .  .02689223
1004 15376  64800  8.15 . 3  528120 -1.0085158  1.909634e-06  .02689223
1004 15377  75400  8.05 . 3  606970   -1.01227  1.667743e-06  .02689223
1004 15378 115700  7.97 . 3  922129 -1.0099379 1.0952241e-06  .02689223
1004 15379  79300  8.24 . 3  653432   -.966123 1.4785364e-06  .02689223
1004 15382  56600  8.34 . 3  472044          .             .  .02689223
1004 15383 139600  8.06 . 3 1125176 -1.0335732   9.18588e-07  .02689223
1004 15384 173100   7.9 . 3 1367490 -1.0198511  7.457832e-07  .02689223
1004 15385 200800   7.6 . 3 1526080 -1.0379747  6.801575e-07  .02689223
1004 15386 118900  7.41 . 3  881049     -1.025  1.163386e-06  .02689223
1004 15390 135900   7.3 . 3  992070          .             .  .02689223
1004 15391  63300  7.51 . 3  475383  -.9712329 2.0430534e-06  .02689223
1004 15392  94500  7.29 . 3  688905 -1.0292943  1.494102e-06  .02689223
1004 15393 109200  7.47 . 3  815724  -.9753087 1.1956356e-06  .02689223
1004 15396 103600  7.46 . 3  772856          .             .  .02689223
1004 15397  48600  7.35 . 3  357210 -1.0147454  2.840753e-06  .02689223
1004 15398 114900  7.45 . 3  856005  -.9863946 1.1523234e-06  .02689223
1004 15399 265700   7.4 . 4 1966180 -1.0067114 5.1201386e-07 .024723677
1004 15400  75300  7.44 . 4  560232  -.9945946 1.7753263e-06 .024723677
1004 15403 194900  7.55 . 4 1471495          .             . .024723677
1004 15404 322400  7.95 . 4 2563080  -.9470199  3.694851e-07 .024723677
1004 15405 290000   8.4 . 4 2436000  -.9433962  3.872727e-07 .024723677
1004 15406 187900  8.88 . 4 1668552  -.9428571  5.650751e-07 .024723677
1004 15407 148000   9.2 . 4 1361600   -.963964  7.079641e-07 .024723677
1004 15410 157500   9.2 . 4 1449000          .             . .024723677
1004 15411 124600  9.41 . 4 1172486  -.9771739  8.334205e-07 .024723677
1004 15412 124700  9.75 . 4 1215825  -.9638682  7.927688e-07 .024723677
1004 15413 192600  9.91 . 4 1908666  -.9835898  5.153284e-07 .024723677
1004 15414  93900   9.8 . 4  920220 -1.0110999  1.098759e-06 .024723677
1004 15417 136100  9.85 . 4 1340585          .             . .024723677
1004 15418 200900  10.1 . 4 2029090  -.9746193  4.803233e-07 .024723677
1004 15419 220600  10.2 . 4 2250120   -.990099  4.400205e-07 .024723677
1004 15420 112100 10.05 . 4 1126605 -1.0147059  9.006758e-07 .024723677
1004 15421 102000 10.11 . 4 1031220  -.9940299  9.639358e-07 .024723677
1004 15424 141800 10.26 . 4 1454868          .             . .024723677
1004 15425 183700 10.54 . 4 1936198  -.9727095  5.023812e-07 .024723677
1004 15426 252400  10.7 . 4 2700680  -.9848197  3.646562e-07 .024723677
1004 15427 353600 10.77 . 4 3808272   -.993458 2.6086843e-07 .024723677
1004 15431 144600    11 . 4 1590600          .             . .024723677
1004 15432 197900    11 . 4 2176900         -1  4.593688e-07 .024723677
1004 15433 113100 11.07 . 4 1252017  -.9936364  7.936285e-07 .024723677
1004 15434  69400 11.08 . 4  768952  -.9990966 1.2992965e-06 .024723677
1004 15435  94500  11.1 . 4 1048950  -.9981949  9.516135e-07 .024723677
1004 15438 124900  11.7 . 4 1461330          .             . .024723677
1004 15439 359000 12.39 . 4 4448010  -.9410256 2.1156104e-07 .024723677
1004 15440 571700 13.09 . 4 7483553  -.9435028 1.2607686e-07 .024723677
1004 15441 323100 13.29 . 4 4293999  -.9847212 2.2932497e-07 .024723677
1004 15442 376900  13.5 . 4 5088150  -.9841986 1.9342956e-07 .024723677
1004 15445  97300 13.11 . 4 1275603          .             . .024723677
1004 15446  96200 13.48 . 4 1296776  -.9717773  7.493794e-07 .024723677
1004 15447 234900 13.65 . 4 3206385  -.9873887  3.079445e-07 .024723677
1004 15448  96400  13.2 . 4 1272480 -1.0329671  8.117747e-07 .024723677
1004 15449  86700 13.15 . 4 1140105 -1.0037879  8.804346e-07 .024723677
1004 15452 155900    13 . 4 2026700          .             . .024723677
1004 15453  74300    13 . 4  965900         -1 1.0353039e-06 .024723677
1004 15454 143000 13.35 . 4 1909050  -.9730769  5.097179e-07 .024723677
1004 15455  79100 13.19 . 4 1043329 -1.0119851  9.699578e-07 .024723677
1004 15456  60400  13.1 . 4  791240 -1.0068233 1.2724626e-06 .024723677
1004 15459 441700 12.54 . 4 5538918          .             . .024723677
1004 15460 361000 12.74 . 4 4599140   -.984051 2.1396414e-07 .024723677
1004 15461 336800  12.6 . 4 4243680 -1.0109891 2.3823405e-07 .024723677
1004 15462 413100 12.33 . 4 5093523 -1.0214286  2.005348e-07 .024723677
1004 15463 232600 12.27 . 4 2854002 -1.0048661  3.520902e-07 .024723677
1004 15466  75300 12.16 . 4  915648          .             . .024723677
1004 15467  58400 11.95 . 4  697880 -1.0172697  1.457657e-06 .024723677
1004 15468 497000 12.28 . 4 6103160  -.9723849 1.5932483e-07 .024723677
1004 15469 217100  11.9 . 4 2583490 -1.0309446  3.990511e-07 .024723677
1004 15470  69900 11.85 . 4  828315 -1.0042017 1.2123427e-06 .024723677
1004 15473 226600  11.7 . 4 2651220          .             . .024723677
1004 15474 325600  12.5 . 4 4070000  -.9316239 2.2890023e-07 .024723677
1004 15475  72600 12.47 . 4  905322    -1.0024 1.1072304e-06 .024723677
1004 15476  80700 12.14 . 4  979698 -1.0264635 1.0477346e-06 .024723677
1004 15477  58300 11.98 . 4  698434 -1.0131795 1.4506446e-06 .024723677
1004 15480  55100 11.89 . 4  655139          .             . .024723677
1004 15481 187900 11.85 . 4 2226615 -1.0033642  4.506231e-07 .024723677
1004 15482  63800    12 . 4  765600  -.9873418 1.2896313e-06 .024723677
1004 15483  30800 11.98 . 4  368984 -1.0016667 2.7146616e-06 .024723677
1004 15484  32200 11.76 . 4  378672  -1.018364 2.6893035e-06 .024723677
end
format %td datadate

------------------ copy up to and including the previous line ------------------

THANKS A LOT!