meas.est: Measurement error variance estimation from repeated measures
In smatr: (Standardised) Major Axis Estimation and Testing Routines

Description Usage Arguments Details Value Author(s) References See Also Examples

Estimates the average variance matrix of measurement error for a set of subjects with repeated measures.

1	meas.est(datameas, id, data=NULL )

`datameas`	A data matrix containing the repeated measures of observations, with each variable in a different column and all observations on all subjects in different rows of the same column.
`id`	An id vector identifying the subject being measured for each observation in the data matrix.
`data`	Deprecated. Use with() instead (see Examples).

This function allows the estimation of measurement error variance, given a set of repeated measures on different subjects. Measurement error variance is estimated separately for each subject, then averaged across subjects. This provides terms that can be used to correct for measurement error in allometric analyses.

Any number of variables can be specified in the data matrix for measurement error calculation. If more than one variable is specified, the covariance of measurement error is estimated from the repeated measures as well as the variance.

As well as the estimated measurement error variance, a data matrix is returned which contains the averages of the repeated measures for each subject.

`V`	A matrix containing the average variances and average convariances of the repeated measures of subjects.
`dat.mean`	A matrix containing the values for each subject, averages across repeated measures. Subjects are in rows, variables in columns.

Warton, D. David.Warton@unsw.edu.au, translated to R by Ormerod, J. 2005-12-08

Warton, D.I., I.J. Wright, D.S. Falster and M. Westoby. 2006. Bivariate line-fitting methods for allometry. Biological Reviews 81, 259–291.

line.cis, slope.test

#load the individual level leaf example dataset
data(leafmeas)

#Estimate measurement error variance matrix, store in "meas.vr"
meas.vr <- meas.est(leafmeas[,3:4], leafmeas$spp)