Description Usage Arguments Details Value Author(s) References See Also Examples
Computes the standard weight equation using the regressionlinepercentile method when given the log(a) and b values for a population of lengthweight regression equations.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29  rlp(log.a, b, min, max, w = 10, qtype = 8, probs = 0.75, digits = NULL)
## S3 method for class 'rlp'
plot(x, what = c("both", "raw", "log"), col.pop = c("rich",
"cm", "default", "grey", "gray", "heat", "jet", "rainbow", "topp", "terrain"),
lwd.pop = 1, lty.pop = 1, order.pop = TRUE, col.Ws = "black",
lwd.Ws = 3, lty.Ws = 1, ...)
## S3 method for class 'rlp'
anova(object, ...)
## S3 method for class 'rlp'
coef(object, ...)
## S3 method for class 'rlp'
predict(object, ...)
## S3 method for class 'rlp'
summary(object, ...)
## S3 method for class 'rlp'
fitPlot(object, pch = 16, col.pt = "black", col.Ws = "red",
lwd.Ws = 3, lty.Ws = 1, xlab = "log10(midpt Length)",
ylab = paste("log10(", 100 * object$prob,
" Percentile of Predicted Weight)", sep = ""), main = "RLP Equation Fit",
...)
## S3 method for class 'rlp'
residPlot(object, ...)

log.a 
A numeric vector that contains the log_{10}(a) values for the population of lengthweight regression equations. 
b 
A numeric vector that contains the b values for the population of lengthweight regression equations 
min 
A number that indicates the midpoint value of the smallest Xmm length category. 
max 
A number that indicates the midpoint value of the largest Xmm length category. 
w 
A number that indicates the widths for which to create length categories. 
qtype 
Type of quantile method to use. See details. 
probs 
A number that indicates the probability of the quantile. Must be between 0 and 1. 
digits 
Number of digits to round predicted weights. 
x 
An object saved from the 
object 
An object saved from 
what 
A string that indicates the type of plot to produce. See details. 
col.pop 
A string that indicates the type of color or palette to use for the population of lengthweight regression lines. See details. 
order.pop 
A logical that indicates whether the populations should be plotted from the smallest to largest weight in the initial length category. See details. 
lwd.pop 
A numeric that indicates the width of the line to use for the population of lengthweight regression lines. 
lty.pop 
A numeric that indicates the type of line to use for the population of lengthweight regression lines. 
col.Ws 
A string that indicates the type of color to use for the standard lengthweight regression line. 
lwd.Ws 
A numeric that indicates the width of the line to use for the standard lengthweight regression line. 
lty.Ws 
A numeric that indicates the type of line to use for the standard lengthweight regression line. 
pch 
A single numeric that indicates what plotting characther codes should be used for the points in the fitPlot. 
col.pt 
A string used to indicate the color of the plotted points. 
xlab 
A label for the xaxis of fitPlot. 
ylab 
A label for the yaxis of fitPlot. 
main 
A label for the main title of fitPlot. 
... 
Additional arguments for methods. 
The main function follows the steps of the regressionlinepercentile method detailed in Murphy et al. (1990). In summary, a predicted weight is constructed for each 1cm length class from each population from the given log_{10}(a) and b values, the predicted weight at the prob
th percentile (wq) is identified, and a linear regression equation is fit to the log_{10}(wq) and log_{10}(midpoint length) data.
Note that log_{10}(a) and b must be from the regression of log_{10}(W) on log_{10}(L) where W is measured in grams and L is the total length measured in mm.
It appears that Murphy et al. (1990) used qtype=6
in their SAS program. Types of quantile calculation methods are discussed in the details of of quantile
.
The plot
, coef
, and summary
methods are used to construct a plot (see below), extract the coefficients of the standard weight equation, and find summary results of the lm
object returned by the main function. The what
argument in the plot
method can be set to "both"
, "log"
, or "raw"
. The "raw"
plot plots lines on the lengthweight scale for each population represented in the log.a
and b
vectors with the resultant standard weight equation superimposed in red. The "log"
plot constructs a similar plot but on the log_{10}(weight)log_{10}(length) scale. The "both"
option produces both plots sidebyside.
If the col.pop
argument is set equal to one of these palettes – “rich”, “cm”, “default”, “grey”, “gray”, “heat”, “jet”, “rainbow”, “topo”, or “terrain” – and the order.pop=TRUE
then the populations plotted should form a general color gradient from smallest to largest weight in the initial length category. This will make it easier to identify populations that “cross over” other populations.
A list is returned with the following items:
log.a
is a numeric vector of the observed log_{10}(a) values sent in the log.a
argument.
b
is a numeric vector of the observed b values sent in the b
argument.
data.pred
is a matrix of the predicted weight at length for all populations.
data.reg
contains a data frame with the prob
th quartile of predicted weights and the midpoint lengths.
Ws
is an lm
object that contains the results of the regression of log_{10}(wq) on log_{10}(midpoint length).
Derek H. Ogle, [email protected]
Murphy, B.R., M.L. Brown, and T.A. Springer. 1990. Evaluation of the relative weight (Wr) index, with new applications to walleye. North American Journal of Fisheries Management, 10:8597.
emp
, FroeseWs
, and wsValidate
; and quantile
in stats
1 2 3 4 5 6 7 8 9 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.