R-help doesn't foward attached data files like this, but Brian kindly forwarded it to me.
You need to restrict X so that it is full rank, it now has rank 19 and column dimension 29 (with intercept). See for example svd(cbind(1,x)).
I'll add some better checking for this, but it will basically amount to setting singular.ok = FALSE in lm() and forcings users to do the rank reduction themselves.
url: www.econ.uiuc.edu/~roger Roger Koenker email firstname.lastname@example.org Department of Economics vox: 217-333-4558 University of Illinois fax: 217-244-6678 Champaign, IL 61820
On Jun 7, 2006, at 3:05 PM, Mu Tian wrote:
> I attached the data file here. I restarted the PC but it still
> happens. It
> says a memory address could not be written. I am not sure it is a
> problem of
> R or quantreg but I plot without problems before I load quantreg.
> Thank you.
> On 6/7/06, Prof Brian Ripley <email@example.com> wrote:
>> Without y and x we cannot reproduce this.
>> PLEASE do read the posting guide!
>> On Wed, 7 Jun 2006, Mu Tian wrote:
>> > I forgot to mention my R version is 2.3.1 and quantreg is the most
>> > too.
>> It has a version number, which the posting guide tells you how to
>> > On 6/7/06, Mu Tian <firstname.lastname@example.org> wrote:
>> >> I was trying "quantreg" package,
>> >> lm1 <- lm(y~x)
>> >> rq1 <- rq(y~x)
>> >> plot(summary(rq1)) #then got a warning says singular value,
>> etc. but
>> >> line can be omited
>> >> plot(lm1) #crash here
>> >> It happened every time on my PC, Windows XP Pro Serv. Pack 1,
>> >> 3.00G.
>> > [[alternative HTML version deleted]]
>> Brian D. Ripley, email@example.com
>> Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
>> University of Oxford, Tel: +44 1865 272861 (self)
>> 1 South Parks Road, +44 1865 272866 (PA)
>> Oxford OX1 3TG, UK Fax: +44 1865 272595
> Rfirstname.lastname@example.org mailing list
> PLEASE do read the posting guide! http://www.R-project.org/posting-
Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Archive generated by hypermail 2.1.8, at Thu 08 Jun 2006 - 08:10:47 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.