jasonelaw / nrsa

R package for calculating habitat metrics from the National Rivers and Streams Assessment protocol
3 stars 1 forks source link

metsChannelMorphology: not using transect K? #4

Open jasonelaw opened 11 years ago

jasonelaw commented 11 years ago

For the wetted width metrics 'n_w', 'xwidth', and 'sdwidth' in metsChannelMorphology.R the wetted width measurement for the K transect is not used in these calculations. I can see where the K transect gets subsetted out; around line 145 in this code bit:

wade <- subset (wade, !(wade$PARAMETER=='WETWID' & wade$UNITS=='NONE'))

You can verify by checking that n_w is 20 (stations 0 and 5/7 for transects A - J = 20) for the test data metrics in your test scripts. I'm getting n_w = 21 for some sites (with corresponding small differences in xwidth and sdwidth).

Why wouldn't all the available data be used?