{MeasurementError.cor}R Documentation

A function to calculate measurement error model estimates for correlation coefficient betweem two variables


Given the observed value of two variables and their respective standard error, the measurement error estimate for their correlation coefficient is returned

Usage, se1, exp2, se2)


exp1 observed value for vector 1
se1 estimated standard error for vector 1
exp2 observed value for vector 2
se2 estimated standard error for vector 2



estimate Vecotr containing the estimates from the measurement error model, i.e., cor.true, mu1, mu2, s1, s2 which are correlation for the measurement error distribution of the two variables, true correlation between the two variables, and true mean and standard deviation estimates for the two variables respectively
count numer of function and gradient evaluation

{ 0 if converged. See optim() for details}


Most applicable for microarray expression data where standard errors are readily estimated by most low level analysis softwares. Hence variables can be thought of as genes. One also need to differentiate between and cor.true: the first one being the correlation between the measurement error distributions of the two genes whereas the second one is the quantity of interest, i.e true correlation between the two gene expression profiles.\

The function involves using quasi-newton for linear optimization, "BFGS" is the only implemented method now.


Beiying Ding


Ding, B.Y. and Gentleman, R. (2003) Measurement Error Model for correlation coefficient estimation and its application in microarray analysis

See Also


 exp <- matrix(abs(rnorm(200,1000,20)),ncol=10)
 se <- matrix(abs(rnorm(200,50,5)),ncol=10)[1,],se[1,],exp[2,],se[2,])

[Package MeasurementError.cor version 1.0.2 Index]