Re: [R] Orthogonal polynomial contrasts

From: Frank E Harrell Jr <>
Date: Thu, 10 Apr 2008 08:21:49 -0500

Peter Dalgaard wrote:
> Chris Bennett wrote:

>> How do you remove one of the terms from an ordered polynomial contrast in
>> your linear model.  For example, I have significant terms for linear and
>> cubic but not quadratic, how would i remove the quadratic term from
>> lm(response~treatment)
>> Cheers,
>> Chris

> Are you _sure_ you want to do that? The interpretation can get tricky,
> but if you insist, just set the contrast matrix for the term
> accordingly, e.g.,
> C(myterm, contr.poly(4)[,-2])

Related to Peter's note of caution, using statistical significance as a basis for this decision is asking for trouble. See

   author = {Grambsch, P. M. and {O'Brien}, P. C.},    year = 1991,
   title = {The effects of transformations and preliminary tests for

           non-linearity in regression},    journal = Stat in Med,
   volume = 10,
   pages = {697-709}


Frank E Harrell Jr   Professor and Chair           School of Medicine
                      Department of Biostatistics   Vanderbilt University

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Thu 10 Apr 2008 - 13:24:38 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 10 Apr 2008 - 13:30:27 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive