# Re: [R] test regression against given slope for reduced major axis regression (RMA)

From: david bird <bird.david_at_uqam.ca>
Date: Wed 26 Jul 2006 - 08:43:44 EST

Patrick Drechsler wrote on 11 Jul 2006 02:10:21 MET:

[...]
> I am now confronted with the problem that I have data which
> requires a modelII regression (also called reduced major axes
> regression (RMA) or geometric mean regression). For this I use
> the function "modelII" (see below).
>
> What would be a good way of adapting
> "test_regression_against_slope" for use with RMA regression?
>
> The question I am trying to answer is: "Does the slope acquired
> from experimental data differ significantly from theoretical
> predictions?"

JFTR: David Warton's "smatr" package solves the problem.

DFB: Note that there are conflicting solutions to this question. The confidence limits based on Pitman, EJG (1939) "A note on normal correlation". Biometrika 31: 9-12, as implemented in the smatr package, are liberal (i.e. narrow), compared to the alternative limits found in Tan, CY and B. Iglewicz (1999) Measurement-methods comparisons and linear statistical relationship, Technometrics 41: 192-201.

Note also that "equation error", if present, invalidates any confidence limits calculations. That is, in using the RMA, you are looking for the perfect underlying relationship between two variables. If the underlying relationship is not perfect, because missing variables are also important, the RMA is the wrong approach.

```--
View this message in context: http://www.nabble.com/test-regression-against-given-slope-for-reduced-major-axis-regression-%28RMA%29-tf1921881.html#a5494654
Sent from the R help forum at Nabble.com.

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
and provide commented, minimal, self-contained, reproducible code.
```
Received on Wed Jul 26 08:48:34 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 26 Jul 2006 - 11:19:46 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.