(PR#7707) Re: [Rd] vcov on result of rlm() yields "-- please report!"

Prof Brian Ripley ripley at stats.ox.ac.uk
Sat Feb 26 08:56:04 CET 2005


This results from calling vcov() on an object for which it does not have a 
method.  That is allowed, but that you get an warning is not a bug: it 
should point up to you that the results are not reliable.  However, there 
is an error in the check!  (So the bug is not what you were asked to 
report: it seems strange to me that it asks for a report on something 
that is perfectly reasonable.)

I do not remember why vcov.lm calls summary.lm explicitly (I suspect it is 
to cope with aov objects), but that stops inheritance working correctly. 
Adding

vcov.rlm <- function (object, ...)
{
     so <- summary(object, corr = FALSE)
     so$stddev^2 * so$cov.unscaled
}

should circumvent this and also give the results you should expect.

I believe the results you get under S-PLUS are equally unreliable, for the 
same reason.  Nowhere that I can see is it stated that vcov works for rlm 
fits (do check its help page), and a leap of faith underlies this.


Finally, it is courtesy to acknowledge when you are using contributed 
software (I hope you _do_ give credit in your reports), and both the 
posting guide and FAQ do ask you to report problems with contributed 
packages to the maintainer.  (You seem unaware that rlm is part of the VR 
bundle.)  E.g.:

    Bug reports on contributed packages should be sent first to the package
    maintainer, and only submitted to the R-bugs repository by package
    maintainers, mentioning the package in the subject line.



On Fri, 25 Feb 2005 frank at trdlnk.com wrote:

> This is a multi-part message in MIME format.
> --------------020205020409020506080205
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
> Content-Transfer-Encoding: 7bit
>
> Dear r-bugs,
> I looked over the FAQ. Hope I'm reporting this correctly.
> I ran this on both solaris and windows. I've provided terminal snapshots
> which include how R was called from the command line, and the
> result of version at the R prompt.
>
> I have attached the .r file, and the data file and the output snapshots.
> Below also find everything except only a few lines of the data file.
>
> Note that it seems to work for lm() but not for rlm().
>
> Not that it matters for our purposes, but the same program(except for
> modifying read.csv()) runs without error under Splus Version 6.2.1on
> Linux 2.4.18.
>
> Thanks for looking at this.
> Of course let me know if I can provide any other information.
>
> Frank Hansen
> 312-264-2043
>
> -------- here is a cut and paste from the terminal window on solaris  ---
> brahe:/home/titan/frank $ R
>
> R : Copyright 2004, The R Foundation for Statistical Computing
> Version 1.9.1  (2004-06-21), ISBN 3-900051-00-3
>
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type 'license()' or 'licence()' for distribution details.
>
> R is a collaborative project with many contributors.
> Type 'contributors()' for more information and
> 'citation()' on how to cite R in publications.
>
> Type 'demo()' for some demos, 'help()' for on-line help, or
> 'help.start()' for a HTML browser interface to help.
> Type 'q()' to quit R.
>
> > version
>         _
> platform i386-pc-solaris2.8
> arch     i386
> os       solaris2.8
> system   i386, solaris2.8
> status
> major    1
> minor    9.1
> year     2004
> month    06
> day      21
> language R
> > source("test.r")
> [1] "before vcov( lm.out)"
> [1] "before vcov(rlm.out)"
> Error in if (rdf != z$df.residual) warning("inconsistent residual
> degrees of fre
> edom. -- please report!") :
>        missing value where TRUE/FALSE needed
> >
>
>
> ---here is a cut and paste of the terminal window on dos -----
>
> H:\Notes>R
>
> R : Copyright 2004, The R Foundation for Statistical Computing
> Version 2.0.1  (2004-11-15), ISBN 3-900051-07-0
>
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type 'license()' or 'licence()' for distribution details.
>
> R is a collaborative project with many contributors.
> Type 'contributors()' for more information and
> 'citation()' on how to cite R or R packages in publications.
>
> Type 'demo()' for some demos, 'help()' for on-line help, or
> 'help.start()' for a HTML browser interface to help.
> Type 'q()' to quit R.
>
> > version
>         _
> platform i386-pc-mingw32
> arch     i386
> os       mingw32
> system   i386, mingw32
> status
> major    2
> minor    0.1
> year     2004
> month    11
> day      15
> language R
> > source("test.r")
> [1] "before vcov( lm.out)"
> [1] "before vcov(rlm.out)"
> Error in if (rdf != z$df.residual) warning("inconsistent residual
> degrees of fre
> edom. -- please report!") :
>        missing value where TRUE/FALSE needed
> >
>
>
>
> here is the program
> ----------------- test.r --------
> library(MASS)
>
> dat <- read.csv( "data.for.r.help.no.I.csv", header=TRUE)
>
> attach( dat)
>
> log1pX <- log(1 + X)
> logZ <- log(Z)
>
> rlm.out <- rlm( Y ~ -1+ Factor + log1pX + logZ + W, method="M")
> lm.out <-   lm( Y ~ -1+ Factor + log1pX + logZ + W)
>
> print( "before vcov( lm.out)" )
> vcov( lm.out)
>
> print( "before vcov(rlm.out)" )
> vcov( rlm.out)
>
> print( "after vcov(rlm.out)" )
>
> ---------------
>
>
> ---------here is head on the data file --------
> "Y","Factor","X","Z","W"
> 8.73469387755,"A",-0.050335552279,3.49,1.30413295097
> 23.5315946724,"B",0.399664852202,4.59,0.914640733428
> 854.100293301,"C",6.65305970571,2.12,-0.127812880306
> 27.1833605167,"D",0.101059456576,3.52,-0.0648127731788
> 17.7277970012,"D",0.0446397114047,4.59,-0.137983116953
> 15.9722222222,"C",0.000378382282472,4.59,-0.495730019986
> 26.2627737226,"E",0.138801178471,3.55,-0.220289949486
> 21.9134615385,"A",0.387109096009,3.55,1.22120113878
> 16.5000165,"C",0.133098708566,3.52,-0.243375187344
>
>
>
>
> --------------020205020409020506080205
> Content-Type: text/x-comma-separated-values;
> name="data.for.r.help.no.I.csv"
> Content-Transfer-Encoding: 7bit
> Content-Disposition: inline;
> filename="data.for.r.help.no.I.csv"
>
> "Y","Factor","X","Z","W"
> 8.73469387755,"A",-0.050335552279,3.49,1.30413295097
> 23.5315946724,"B",0.399664852202,4.59,0.914640733428
> 854.100293301,"C",6.65305970571,2.12,-0.127812880306
> 27.1833605167,"D",0.101059456576,3.52,-0.0648127731788
> 17.7277970012,"D",0.0446397114047,4.59,-0.137983116953
> 15.9722222222,"C",0.000378382282472,4.59,-0.495730019986
> 26.2627737226,"E",0.138801178471,3.55,-0.220289949486
> 21.9134615385,"A",0.387109096009,3.55,1.22120113878
> 16.5000165,"C",0.133098708566,3.52,-0.243375187344
> 28.5052631579,"B",0.623439422305,3.52,-0.422811577042
> 25.8745622811,"D",0.229205420669,2.12,-0.200177307421
> 16.4680232558,"E",0.179418817508,3.51,-0.335809927225
> 11.8625,"F",-0.188194758578,3.63,-0.430203019688
> 28.75,"E",0.599979969827,3.55,-0.3972666749
> 11.0707456979,"A",0.0222998422651,3.55,-0.024764276122
> 19.5514705882,"C",0.114905693717,3.51,-0.429671828481
> 18.6742934051,"B",0.0297655497652,4.59,-0.474833455898
> 54.347826087,"E",0.0657718166744,3.63,-0.442615879212
> 8.61157024793,"A",-0.0679442165455,3.49,0.249656986243
> 55.9459459459,"H",0.402878814284,4.59,-0.523454051056
> 12.3357605624,"F",-0.174806348851,3.52,-0.484777106295
> 20.0869766087,"B",0.0997440131645,3.55,0.174734541472
> 12.3903345725,"B",0.0505165785386,3.49,-0.112688455668
> 26.0075,"A",0.157834271388,3.63,-0.132912730886
> 11.9518377693,"A",0.061384765068,4.59,-0.28334400572
> 21.7592592593,"C",0.119345555323,3.69,-0.539329198267
> 14.913409802,"E",0.102525055053,3.69,-0.419968044113
> 16.8187539333,"D",0.0850209681194,4.59,-0.152149322452
> 18.0620985011,"A",0.156630637908,3.51,0.264569349391
> 18.0671936759,"D",0.0804427049757,3.51,-0.26504156369
> 21.8181818182,"D",0.130936716684,3.52,-0.101664163159
> 20.0691823899,"A",0.375908602408,3.51,1.36043216402
> 74.2954545455,"A",1.41196236591,3.51,-0.427400156085
> 15.652173913,"E",0.0873392184994,3.5,-0.560279296469
> 15.3086419753,"C",0.0513013687336,4.59,-0.533314287835
> 20.1111111111,"C",0.0405328362591,3.52,-0.38391302595
> 28.7378947312,"B",0.586839294657,4.59,-0.40989782081
> 19.1723549488,"C",0.0639887969656,4.59,-0.527845508362
> 35.2649006623,"G",0.274574479051,4.59,-0.50384770061
> 0.719016779511,"D",-0.833129846672,3.55,-0.317205785189
> 30.085106383,"D",0.179936117396,3.55,-0.0870003808961
> 28.3265306122,"G",0.213416966543,3.49,-0.256214659806
> 31.4049586777,"C",0.144292683957,3.69,-0.37281444967
> 16.9541646131,"C",0.0880464420824,3.69,-0.197866625671
> 56.7961165049,"D",0.51303086006,2.12,-0.524345871293
> 33.2228915663,"G",0.341625506615,3.55,-0.41734819248
> 20.704787234,"D",0.135937099322,4.59,-0.398338602156
> 55.9024390244,"D",0.124499281523,3.63,0.0932087260108
> 12.2368421053,"A",0.102581097217,4.59,0.553328624412
> 11.9350282486,"F",-0.238121374793,4.59,-0.518984990035
> 18.3582915783,"D",0.146883469069,4.59,2.06751779427
> 638.75,"C",2.28305714236,4.59,-0.00459229003918
> 18.2692307692,"C",0.100688410774,4.59,-0.0855122155304
> 21.83,"A",0.283098380224,3.49,-0.085660368078
> 33.0893682589,"A",0.359909323826,2.12,-0.463869752385
> 35.3939393939,"A",0.914794869938,2.12,-0.502945505544
> 12.1089108911,"A",0.0861491202801,4.59,2.4615712013
> 34.8148148148,"B",0.607335568458,3.51,1.0114071716
> 21.9168900804,"D",0.124819185847,4.59,1.06303273202
> 16.8321916166,"B",-0.0680084359348,3.51,0.42731554537
> 23.6776859504,"E",0.375488475818,4.59,-0.375904488519
> 11.8466898955,"A",0.0381514439413,4.59,-0.442417097502
> 34.3509350935,"C",0.416855937657,4.59,-0.50460464808
> 16.0843373494,"A",0.15179876087,3.52,1.81905601218
> 21.3805253043,"C",0.165915546195,4.59,-0.503758477087
> 6.15492957746,"E",0.0127846528516,4.59,-0.490423917821
> 24.8085714286,"D",0.0867753888059,3.51,-0.21975087341
> 12.2097378277,"H",-0.0149331797033,3.51,0.0141799242148
> 29.476555477,"A",0.195826082472,3.63,-0.388398271703
> 16.5618448637,"J",0.0360542479099,4.59,0.451338667703
> 26.9806763285,"D",0.180334022977,4.59,-0.454131523587
> 23.6900893603,"A",0.369162188942,3.63,0.75107409603
> 6.36433218523,"E",-0.0650124215212,4.59,-0.535384688562
> 18.6930232558,"A",0.150747750833,3.52,-0.448016184817
> 20.6554350335,"A",-0.16903719864,4.59,-0.253885303365
> 23.1839359697,"B",0.330816121952,4.59,0.700759913906
> 14.2833333333,"A",0.0504728036037,3.63,-0.196384685202
> 11.401384083,"E",0.0588094263767,3.51,0.561548393346
> 21.1988304094,"D",0.126786476293,3.55,0.426740364892
> 21.902608808,"A",0.214354002576,3.63,-0.267421549293
> 9.85989492119,"A",0.0205611433351,3.55,2.05925569603
> 17.5729927007,"D",0.042278194822,4.59,-0.449556639318
> 28.1575919943,"G",0.0250474453866,3.55,-0.506641019371
> 16.0108765195,"B",0.0710275727451,4.59,-0.225880321946
> 14.232596571,"E",0.0563706340798,4.59,-0.209093849821
> 15.5168143974,"J",-0.0212597572544,4.59,3
> 19.3536673929,"C",0.17411115781,4.59,-0.40394017943
> 22.0322580645,"A",0.0746342213411,3.63,-0.463886767104
> 23.8358208955,"A",0.273664878954,3.63,-0.320777216069
> 14.0601851852,"A",0.0715585741236,3.63,-0.311936617412
> 21.8696883853,"E",0.155169033656,4.59,-0.466873887656
> 12.8799392097,"A",0.183681470544,4.59,-0.549320987867
> 15.7511581734,"A",0.0706974909567,4.59,-0.540393655646
> 35.9670530811,"H",0.302259030559,4.59,0.266881691114
> 15.0070126227,"H",0.143923838482,4.59,0.801914489443
> 15.5,"A",0.118286147591,3.52,-0.172178135905
> 12.6231884058,"E",-0.0304640285167,3.49,-0.551206301659
> 30.4935483871,"E",0.191880582349,3.63,-0.463242697765
> 20.7250755287,"B",0.116436155824,4.59,-0.563784328448
> 9.2807112069,"E",0.0947023744031,4.59,-0.551459447468
> 6.80956434264,"E",0.0291322117129,4.59,-0.537963455874
> 22.2639149468,"A",-0.014394485535,4.59,-0.539485235684
> 11.015435936,"C",-0.0393762420727,3.63,-0.485946556937
> 17.7212389381,"D",0.0681736812309,3.55,0.895551049429
> 65.5365853659,"G",0.392164002464,3.55,-0.423765231256
> 268.714285714,"A",3.09725517878,3.63,-0.171256851156
> 20.3508367484,"D",-0.133467466718,3.59,-0.462396941765
> 18.3097031458,"A",0.176576193731,4.59,-0.522722003174
> 81.619047619,"A",1.53521809901,3.51,-0.0691390765633
> 20.045045045,"D",0.0432678564067,4.59,0.213082396681
> 25.8842763909,"E",0.184626552532,4.59,-0.197630909573
> 32.005988024,"D",0.342763050869,3.55,-0.104045808734
> 29.1165334017,"E",0.367233759914,2.12,0.138132977357
> 47.0476190476,"H",0.664760815081,3.55,0.0819756919569
> 3.8962894337,"A",-0.382206694948,2.12,-0.557675629569
> 14.593639576,"B",0.117845114337,2.12,-0.412896561171
> 22.2195121951,"D",-0.137502209264,3.52,-0.372611518029
> 19.0950833256,"B",0.0821422178835,2.12,-0.451380119133
> 19.3523489933,"D",0.17937606215,3.51,-0.349684807547
> 16.8920471438,"F",0.186010553946,4.59,0.113705651715
> 34.7398477157,"J",-0.130566228444,4.59,-0.475698716574
> 12.8723404255,"A",0.0699942555247,3.55,-0.206506782646
> 14.0308370044,"D",0.0646825983561,3.55,-0.249744086917
> 16.0314579552,"C",-0.0236918978392,4.59,-0.498661946452
> 25.2911813644,"D",0.239699657685,3.59,-0.511894832402
> 23.4358047016,"D",0.111911505969,3.51,-0.367863166631
> 23.5654008439,"E",0.423832387139,3.51,-0.198764255812
> 16.9365671642,"A",0.0410579309873,3.49,-0.453447199915
> 19.8701298701,"H",0.157976291169,3.6,-0.342367233678
> 32,"E",0.153114522421,4.59,-0.169491885372
> 40.7386363636,"D",0.279737643567,4.59,-0.428702819521
> 43.5897435897,"E",0.356357430579,4.59,-0.497797515762
> 39.4897959184,"A",0.750261554798,3.55,1.37262092725
> 16.4838709677,"A",0.216883595037,2.12,1.52865792929
> 12.7197802198,"D",-0.130667668094,3.49,-0.513609169023
> 19.6818181818,"A",-0.172902116379,4.59,-0.543592007702
> 16.1052009456,"A",0.154350973583,4.59,-0.508886962192
> 10.4838709677,"H",-0.151054048736,4.59,-0.373613726439
> 10.8053482587,"A",0.0069533478648,3.59,0.225276969818
> 22.4045454545,"A",0.22021834298,3.51,-0.315822613085
> 13.1760154739,"A",0.176112074259,3.51,1.36775388782
> 20.0054436581,"D",0.284217740794,3.59,-0.48311132387
> 42.0470588235,"E",0.923441006603,3.51,1.3185203478
> 105.789473684,"H",0.384392057338,3.55,-0.374764917383
> 20,"D",0.12309069455,3.51,-0.331828483132
> 17.8890600924,"C",0.146768565037,3.55,-0.0926127479917
> 15.5873806559,"D",0.0318324392305,4.59,1.13956037022
> 36.9065100343,"H",0.0879124954014,4.59,-0.373212013089
> 14.7310513447,"A",0.153426149626,3.63,-0.104378218232
> 1.63905325444,"A",-0.644129253569,3.51,-0.298260103806
> 15.2627118644,"C",0.180879261613,3.51,0.46613151283
> 28.3611111111,"B",0.538548014995,3.49,-0.352794766066
> 27.1990018715,"C",0.304486236746,3.59,-0.548801001475
> 29.9590163934,"G",0.0977072342529,3.55,-0.391535619769
> 52.6451612903,"H",0.5902584797,3.51,3
> 19.1981747066,"D",0.0860864826947,4.59,-0.544536117074
> 15.3813559322,"A",0.152268293558,3.51,0.485330755015
> 13.7746710526,"E",0.0526889539257,4.59,-0.363267117712
> 126.315789474,"C",1.70314955678,3.49,0.41352573864
> 23.7654320988,"A",0.304469179209,3.49,-0.114049218143
> 20.5960264901,"C",0.0500650465503,4.59,-0.50853795297
> 27.7777777778,"D",0.126745552644,3.55,-0.310931089057
> 13.7538461538,"F",0.0759126251342,3.52,0.405647509052
> 22.0334779913,"F",0.0866180261048,3.51,0.42178078199
> 21.6451612903,"H",0.251392186018,3.5,-0.512703654013
> 9.77286312014,"E",-0.184666836952,3.69,-0.346289748747
> 23.0271839535,"H",0.0634428857159,4.59,-0.436288063959
> 11.1342592593,"A",-0.0725260731261,4.59,-0.266702366198
> 15.6127770535,"E",-0.117894837561,4.59,-0.468494850824
> 23.6014492754,"C",0.26191892066,3.55,-0.0351536291418
> 17.3876146789,"F",0.00820885323584,3.63,1.89576251242
> 19.5454545455,"F",0.137618132394,3.52,-0.451332394923
> 13.8489208633,"J",0.08691887475,3.63,3
> 21.3662790698,"H",-0.048573584644,3.69,-0.543103145795
> 11.8708609272,"F",-0.0337379748099,3.51,2.02929070206
> 13.3333333333,"C",-0.0774938771658,3.51,-0.579547842507
> 17.1666666667,"H",0.0855089129068,3.69,1.33901063362
> 30.80625,"H",0.121870647624,3.63,-0.50053481045
> 17.686533212,"E",0.215877973092,4.59,-0.401745695756
> 19.0769230769,"D",0.191872244057,3.52,-0.514102180862
> 25.6727828746,"C",0.118740372244,4.59,-0.419133077934
> 53.7398365677,"G",0.0955785834204,3.63,-0.591176365014
> 16.5805084746,"C",0.106557155003,3.52,-0.345995933611
> 22.7962085308,"C",0.12497828731,3.52,-0.257976720638
> 17.0238095238,"C",0.0908202826532,4.59,-0.501063096705
> 19.7086368366,"C",0.151933632034,3.55,-0.463464304097
> 26.8279569892,"C",0.214050870487,4.59,-0.557364799714
> 23.5935916546,"D",0.134426281626,3.51,-0.304631078344
> 17.6186046512,"C",-0.0716245747368,3.63,-0.276489979178
> 12.4597701149,"A",0.123721493188,3.51,0.886558563287
> 14.9221453287,"J",0.100871424252,3.55,3
> 16.9387755102,"A",0.166836722397,4.59,-0.337765374855
> 15.9454545455,"E",-0.0382307969164,3.55,-0.463398320189
> 29.5519713262,"D",0.310161625095,3.51,-0.43253777104
> 8.22366625872,"C",0.118889036586,3.63,-0.466267997686
> 12.4523809524,"A",0.062622541953,3.63,-0.434202308486
> 10.5663824604,"A",-0.0238781378493,3.69,-0.0915408207358
> 18.4360902256,"H",0.235178486996,3.55,2.60433630804
> 10.3611111111,"F",-0.100010061785,3.5,0.0748021207047
> 40.6984126984,"A",0.382261058889,3.63,1.15922398972
> 12.9638683856,"J",0.0378819860456,3.52,-0.474538395782
> 65.3594771242,"D",0.546131584251,3.69,-0.407604983765
> 16.4647887324,"F",0.0891782695302,3.52,-0.0535332598944
> 26.0534288075,"C",0.300760322538,2.12,-0.477424258011
> 21.6335540839,"D",0.224294183291,3.69,-0.498386806007
> 31.3909774436,"A",0.0919122427488,3.55,-0.337112590661
> 22.6646032639,"B",0.247232549994,2.12,-0.393807707158
> 13.1592827004,"E",-0.171981611631,4.59,-0.47575515564
> 21.9898989899,"E",0.10237130081,3.55,-0.460727424402
> 10.7835325365,"A",0.0432001236618,3.55,1.31669147308
> 19.2448512586,"D",0.0872286027114,4.59,-0.181966993865
> 50.2016129032,"H",0.451557633894,3.59,-0.514688566156
> 18.600028397,"D",0.0834025942902,4.59,-0.540474579307
> 16.9733333333,"D",0.0714769596527,4.59,3
> 2.3734939759,"G",-0.651163724667,3.49,-0.564940914303
> 24.4656776465,"E",0.15746444021,4.59,0.161204520442
> 35.386491969,"D",0.379737713576,4.59,-0.44124225195
> 13.5796766744,"A",0.149655068262,4.59,2.08756403743
> 14.5733173077,"D",-0.00421561648039,3.49,-0.300108068216
> 37.3758865248,"C",0.700855020694,3.51,-0.33301079856
> 16.2257281553,"D",-0.0231356758654,3.49,0.0869601744462
> 35.8024691358,"E",0.433001239815,4.59,-0.220729427211
> 26.5642775882,"E",0.145181111012,3.69,0.543715723539
> 35.9306496612,"A",0.0634360675793,4.59,-0.418075260445
> 13.431372549,"E",0.0761633310962,4.59,-0.529578934668
> 9.81481481481,"A",0.059674961812,3.51,-0.411268958113
> 16.4187116564,"C",0.0398966253484,4.59,-0.306076499417
> 14.7100840336,"C",-0.122146118665,3.55,-0.300106823237
> 16.091160221,"B",0.153101666693,3.69,-0.496541331556
> 23.1329552953,"C",0.193240449897,2.12,0.010088091949
> 91.5,"C",0.867904811444,3.51,-0.524764184369
> 80.25,"A",1.41771705334,3.51,-0.382457645041
> 18.5911283997,"A",0.142749463786,4.59,3
> 12.8681318681,"C",0.133500611982,3.5,-0.549932272749
> 79.2570718257,"A",1.46880889489,3.51,0.0501112744248
> 22.7941176471,"A",0.0558297967089,4.59,-0.477267805601
> 15.0298507463,"F",0.107463910668,3.52,0.35024551115
> 17.4358974359,"E",0.100630355518,3.52,-0.478049652658
> 14.9498151083,"B",0.123275916212,4.59,-0.355038218929
> 16.8370244179,"A",0.0979537283289,4.59,-0.54973100108
> 28.5833333333,"G",0.161631270384,3.52,-0.496276980932
> 70,"G",0.520211199442,3.51,-0.301673422304
> 43.0695652174,"J",0.61417819985,3.49,-0.524551707886
> 33.0543892603,"E",0.273685344545,4.59,-0.278226725429
> 40.8379888268,"A",0.149271959734,4.59,-0.26516357167
> 15.7894736842,"C",-0.0108066812806,4.59,-0.498171839565
> 18.334384858,"D",0.117215443414,4.59,0.213653842221
> 25.2991452991,"D",0.112479931635,3.51,0.29115712927
> 15.6873111782,"E",0.0924132422499,4.59,-0.485642781965
> 11.3830845771,"A",0.107334523199,3.49,2.09786665689
> 78.5123966942,"A",-0.27929098271,3.69,-0.189358436511
> 31.6652286454,"G",0.219034045239,4.59,-0.549154160629
> 19.3807753385,"F",0.092455401516,4.59,-0.408854943073
> 53.6098310292,"D",0.170117783271,3.59,-0.502660820257
> 10.5555555556,"C",-0.505707857414,4.59,-0.525631520011
> 10.6587301587,"C",-0.0825594349081,3.63,-0.427488134629
> 24.6998284734,"D",0.153680068972,4.59,-0.52737324618
> 31,"H",0.222957480864,4.59,-0.218224943669
> 31.8016759777,"E",0.268262926981,3.52,0.137670675009
> 16.2231030578,"E",-0.0171103304584,4.59,-0.571786640982
> 10.1031870428,"A",-0.099702851794,4.59,-0.554039044767
> 19.3030468893,"B",0.0915343879212,4.59,-0.434171598994
> 17.9284130504,"H",0.104955996928,4.59,-0.487734762336
> 17.4860335196,"A",0.219229215483,3.51,0.421632629443
> 14.7502356268,"F",0.0400977472099,4.59,-0.428364185127
> 15.8425348056,"E",0.11221498218,4.59,-0.478172490625
> 18.7735849057,"D",0.186047467153,3.5,-0.512582061026
> 16.9382022472,"E",0.129731544625,4.59,-0.0635748486708
> 23.2394366197,"D",0.126282194519,4.59,-0.506352184152
> 11.786389414,"C",-0.0366433260285,3.5,-0.401090836597
> 28.6923076923,"H",0.400782707028,3.49,0.325890809306
> 64.3571428571,"H",0.622211893857,2.12,-0.534120204494
> 23.7710219922,"A",0.0295035287388,4.59,-0.0814689374609
> 19.6341463415,"C",0.201182022518,3.55,0.124655245461
> 8.92045454545,"A",-0.121317221201,4.59,-0.0184368758627
> 38.3314329152,"A",0.109295679861,4.59,-0.361858631027
> 25.8883248731,"C",0.297105571241,3.51,-0.287330844724
> 8.70926243568,"A",-0.133096498571,4.59,1.64397663536
> 21.9054921075,"B",0.0896567343111,4.59,-0.502752948732
> 31.5789473684,"G",0.341957632314,4.59,-0.507104981691
> 25.7446808511,"B",0.342236933561,2.12,-0.487953048722
> 35.8054908337,"B",1.44209818133,3.55,0.016621328801
> 22.5916806714,"A",-0.00718018111553,4.59,-0.501414595886
> 91.3636363636,"E",0.831438414653,3.63,-0.471447526947
> 20.4584040747,"D",0.045593334453,3.59,-0.412785758005
> 28.2716666667,"D",0.097598385871,3.63,-0.201572099332
> 15.4305680794,"E",0.0902650689852,4.59,-0.47666316061
> 15.6596794081,"C",0.175068864892,3.59,0.952807651632
> 25.128,"C",0.470087515762,3.51,-0.367289231132
> 14.8765432099,"A",0.0827491302293,3.52,-0.464082643861
> 15.6976744186,"G",-0.00861267044142,2.12,-0.552899058638
> 18.5323943662,"A",-0.094098635783,3.63,-0.0474610804102
> 34.8705882353,"A",0.143232594316,3.63,-0.38516589021
> 26.9802631579,"D",0.36456013252,3.51,0.106685627921
> 27.68,"C",0.420455581847,3.5,-0.514044081824
> 8.94136807818,"E",-0.235223424932,4.59,-0.547269261831
> 28.0271280725,"C",0.115198398785,4.59,-0.518919006127
> 33.4384615385,"A",0.844951531571,3.63,-0.507836614579
> 14.9959742351,"E",0.0808356232679,4.59,-0.137900118327
> 12.4289099526,"E",0.0900984965715,3.55,-0.209497223143
> 26.14,"A",0.128160972557,3.51,-0.536928463007
> 21.5750915751,"H",0.094441500402,3.55,-0.17126141608
> 22.19,"H",0.0922393922579,3.55,-0.457189607965
> 28,"H",0.159120601844,3.63,-0.491988026929
> 16.9077306733,"H",0.0922129301505,3.51,-0.384886184841
> 22.8973214286,"D",0.0469911246895,3.51,-0.0872087074475
> 17.1564733916,"D",0.0145131704391,4.59,-0.169464080832
> 33.628125,"H",0.269777527266,3.63,0.0654162210743
> 73.5483870968,"G",0.417830009375,4.59,-0.514574028052
> 12.0839224728,"A",0.0231861611341,3.52,-0.467306725491
> 14.57244208,"A",0.137738308636,3.69,0.579362388461
> 22.9744525547,"G",0.0920845165283,3.55,-0.507153535887
> 19.7663551402,"D",0.0772217027542,2.12,3
> 33.4966088922,"C",0.136322284649,4.59,-0.518106034584
> 11.7216117216,"D",0.00657694006158,4.59,-0.553016501694
> 17.1363350126,"F",0.0612161665768,4.59,0.679595679249
> 25.8571428571,"E",0.0733892583774,3.52,-0.545014604153
> 13.8746854292,"J",0.0699062791292,3.63,0.228468681984
> 18.3177570093,"C",0.0939298957531,3.55,-0.394781281042
> 44.1489007728,"D",0.225753636962,2.12,-0.393755833017
> 18.7121212121,"H",0.0932385167682,3.5,2.48662392159
> 19.3791946309,"D",0.10255210046,2.12,3
> 21.3230355944,"D",0.0734596367224,3.59,0.195974719873
> 11.1613164284,"A",0.0803323258047,4.59,-0.344061235637
> 26.9696969697,"B",0.555298549404,3.63,-0.501548638668
> 13.0769230769,"D",-0.265104510599,3.49,-0.34097949665
> 23.8693467337,"C",0.00875563924,4.59,-0.550601656668
> 14.871897518,"F",0.138524039971,4.59,-0.539174405829
> 15.0225988701,"A",-0.0737659018802,3.49,-0.442706347714
> 17.0848267622,"E",0.168086216393,4.59,-0.288521045022
> 24.192139738,"D",0.0931991284747,3.59,-0.469108625664
> 11.6213592233,"G",-0.0719645449525,3.49,-0.558977463019
> 25.3666666667,"D",0.201318719852,3.51,-0.556741480032
> 15.2528089888,"A",0.168740799399,4.59,-0.20740482778
> 16.897260274,"J",0.256123547169,3.49,-0.229976304143
> 33.9958158996,"A",0.537416868522,3.69,-0.558643393549
> 43.4210526316,"H",0.290762456573,4.59,0.225154961837
> 26.3574660633,"D",0.0487592108549,3.55,1.50205396468
> 14.9861878453,"D",0.0747365731459,3.63,-0.322915260677
> 6.66666666667,"B",-0.0860887564069,3.5,-0.476260617273
> 23.8630806846,"H",0.104494164499,3.55,-0.320613708776
> 12.8429602888,"A",0.0381342040564,4.59,-0.472848958748
> 16.447944007,"C",0.011626949023,4.59,-0.563100834763
> 34.79,"H",-0.0123067277743,3.63,-0.459699486418
> 40.0943396226,"D",0.275472032489,2.12,-0.520140745903
> 9.54054054054,"C",-0.13038305523,3.49,0.850166570702
> 18.1481481481,"C",0.139393525489,4.59,-0.489943355776
> 61.9214285714,"A",1.09591439302,3.63,-0.475678381911
> 18.4525971667,"F",0.0596162323854,4.59,-0.516080038122
> 70.4751131222,"C",0.821192488851,4.59,-0.434791183738
> 8.82246376812,"A",-0.136068977525,3.5,-0.535683068623
> 14.9715909091,"A",0.131064525227,3.52,0.0157245286462
> 13.7242268041,"E",0.0962084205518,3.51,-0.115406660672
> 53.6470588235,"G",1.24686810786,3.5,-0.588524973904
> 16.2103174603,"E",0.0955623076034,3.55,-0.0785278811457
> 59.1111111111,"H",0.486090708063,3.52,-0.432458092359
> 22.7863474656,"A",0.347169217352,4.59,0.262687355545
> 31.3260847092,"H",0.319406357414,3.51,-0.447079960315
> 9.00493916396,"E",0.15402992877,4.59,-0.444010671123
> 119.636363636,"H",0.0964743009881,4.59,-0.213320139861
> 26.1304347826,"H",0.112114276072,3.55,-0.359394401819
> 14.4218942189,"E",0.139093332584,4.59,-0.476618341351
> 11.3636363636,"E",0.0641250839405,3.51,-0.455103022505
> 33.9364902607,"D",0.157791139452,3.49,-0.543053761612
> 20.4209302326,"A",0.152418003944,3.63,0.0227267077355
> 20.0439185068,"D",0.0805387261723,4.59,0.279703318871
> 14.142746315,"A",0.117575380713,4.59,3
> 17.9130434783,"A",0.0705695504351,3.49,-0.321861178126
> 75.6666666667,"H",0.581138830084,3.63,-0.432194571721
> 10.4063205418,"J",-0.144034861525,3.49,3
> 19.6428571429,"F",0.129494803089,3.52,-0.478568394071
> 37.3507462687,"H",0.0729417946823,4.59,-0.124944032795
> 38.4642857143,"H",0.197706901303,3.49,-0.188467031267
> 12.1939393939,"J",0.0540030590491,3.52,0.697326260745
> 22.9092752154,"C",0.181361536782,4.59,-0.420065567498
> 19.5220588235,"D",0.153094756099,3.55,-0.20950344804
> 20.1102364408,"B",0.284699481617,4.59,0.793791413909
> 35.8193277311,"C",0.45587341627,4.59,-0.437218063564
> 27.9421768707,"E",0.0540930653771,3.55,-0.29129278414
> 14.8177083333,"D",0.0259373534115,4.59,-0.480322570033
> 5.87959409909,"E",0.0224311204646,3.55,0.428481676067
> 23.9341085271,"D",0.134614201242,2.12,3
> 22.2038111019,"C",0.235470398358,4.59,-0.257220603154
> 55.0847457627,"D",0.795635914929,3.49,-0.411890202829
> 9.76941747573,"D",-0.119741766499,4.59,-0.387644229185
> 93.3066361556,"H",0.901021997809,4.59,-0.152179201958
> 22.9487179487,"B",0.936105441873,3.5,-0.549038377546
> 15.1459854015,"A",0.244449238392,2.12,0.307377550765
> 28.2307692308,"F",0.500251406939,3.51,0.801055868656
> 37.5519134495,"D",0.430848008848,4.59,0.0796085711411
> 17.8504672897,"A",0.194646676739,3.52,-0.301284158747
> 30.6488549618,"C",0.313646809052,3.55,-0.466225253394
> 16.015625,"A",0.274141868082,3.5,-0.305012872024
> 19.7260273973,"D",0.180448743615,3.55,3
> 17.6727272727,"A",0.16813176405,3.52,2.66873286705
> 16.1581920904,"C",0.0656754278457,3.69,-0.166892783396
> 29.6872037915,"G",0.138019931692,3.51,1.50696001347
> 23.0198019802,"D",0.062772285437,3.51,-0.549589488423
> 14.6839098815,"D",-0.219541601281,4.59,-0.227480535456
> 17.4679487179,"B",0.229805373126,3.69,-0.235817332453
> 23.4473684211,"D",0.111562244234,3.51,-0.0601470054137
> 16.2186379928,"H",0.0731589890462,2.12,0.418459591967
> 14.7203140334,"E",0.117872011306,3.69,-0.264549381837
> 26.9281045752,"B",0.322208491168,3.51,-0.520952887459
> 19.1935483871,"A",-0.0866827477673,3.69,-0.399762028593
> 50.1564455569,"A",0.447063402417,4.59,-0.432115308033
> 13.4441087613,"A",0.10398422871,3.69,-0.0637777803116
> 27.8606965174,"D",0.187170232819,3.69,-0.378796990638
> 18.4579444372,"D",0.101730414797,2.12,-0.229357964378
> 16.659242805,"C",0.0477698822753,4.59,-0.517904762916
> 21.3852376138,"G",0.0994125658227,4.59,-0.405228318107
> 16.2876254181,"C",0.11089118011,3.55,0.273653549016
> 23.8088642659,"D",0.172106562994,2.12,-0.37879989559
> 14.7008547009,"A",0.115845312932,3.69,0.0863505495377
> 12.5544740807,"J",-0.0458513772319,4.59,2.93518750629
> 19.1916241557,"D",0.169123435715,4.59,-0.397332658808
> 16.0261780105,"C",0.0852413054706,3.51,2.32167037715
> 14.4472361809,"E",0.0434478241911,4.59,-0.45884335559
> 12.6136363636,"C",-0.0769783785651,3.69,-0.406901570409
> 12.8126202386,"C",0.0380458001018,3.69,-0.390882835575
> 18.5580371959,"A",0.342975558177,2.12,-0.408075585975
> 23.9397321429,"A",0.0900455142857,4.59,-0.440129240374
> 21.6226415094,"C",0.134201200968,4.59,-0.237581053257
> 22.1382916909,"E",0.369944529726,4.59,-0.561494811348
> 25.1936218679,"D",0.0858645046926,4.59,-0.206964935062
> 21.2941847206,"D",0.204941222183,4.59,-0.535253550733
> 18.3368200837,"B",0.120862494786,3.51,-0.218412520564
> 16.2594202899,"A",0.0692642059047,3.63,0.50617378499
> 21.9377777778,"J",0.057552058083,3.63,-0.380831286963
> 15.0370245709,"F",0.0442051584781,4.59,0.0186950494736
> 18.32,"F",0.048801430205,3.63,-0.131998086026
> 27.8156996587,"D",0.0852590742528,4.59,-0.524143769639
> 4.94366197183,"B",-0.420434030053,3.69,-0.512083239283
> 18.0566182096,"F",0.107770684847,4.59,0.0130337131886
> 32.7198364008,"D",0.454490262013,4.59,-0.567479427282
> 22.0997395591,"A",0.177201824057,4.59,0.41815374203
> 21.9101123596,"A",0.176842921622,3.63,-0.510971472687
> 25.7954545455,"B",0.400971975906,3.5,-0.108790840187
> 9.30434782609,"J",-0.104398069793,3.63,0.86953720505
> 23.7337982707,"G",0.0603361078113,2.12,0.041630889805
> 51.3942307692,"D",0.238030417018,2.12,-0.468190245866
> 22.1722846442,"C",0.0600632462935,3.55,0.126006878086
> 18.7451133698,"F",0.0953967895979,4.59,-0.355437442321
> 17.9111111111,"E",0.0458283231398,3.55,-0.408307567134
> 14.7058823529,"D",-0.00737609375386,3.69,-0.449616813321
> 21.6666666667,"C",0.0961096453501,4.59,-0.569362666107
> 29.9010438968,"B",0.383530392662,2.12,-0.114098187332
> 38.8978494624,"D",0.141998639436,2.12,-0.0510507709883
> 9.97868420407,"A",-0.00102441760482,3.52,-0.4252426068
> 36.2333333333,"F",0.433981990821,3.55,0.568210693061
> 31.4814814815,"H",0.333422509364,3.63,-0.55155904582
> 67.5,"H",0.70433620754,3.63,-0.580792821898
> 16.1146496815,"D",0.0701493919164,4.59,-0.42722585897
> 22.3441615452,"H",0.103820121796,3.59,-0.213981638911
> 33.0909090909,"H",0.363031447646,3.52,3
> 22.5680933852,"H",0.222323797905,3.63,0.0525472841003
> 16.1538461538,"H",-0.0684763749641,3.69,-0.388038057666
> 27.3243243243,"H",0.178522586329,3.52,1.7496861755
> 27.8602150538,"H",0.152695895697,3.55,-0.315798958477
> 40.3846153846,"H",0.207080473386,4.59,-0.515952635231
> 23.3805309735,"G",0.139863258373,3.5,-0.497512830474
> 22.5583333333,"H",0.0926215416755,3.63,1.88550222226
> 21.8928571429,"H",0.0532990010981,3.63,1.93449672623
> 56.4074074074,"H",0.542324160096,3.49,-0.208033542372
> 9.55367231638,"B",-0.0126366360709,3.51,-0.146204130866
> 46.7714285714,"H",0.635336653851,3.52,-0.395713770606
> 15.1624548736,"H",0.0951298094925,3.69,0.544464371146
> 58.3047945205,"J",0.774955021618,4.59,-0.531539777209
> 20.5050505051,"E",0.0971440703532,3.55,-0.325810252755
> 16.8205003643,"H",0.124079046355,4.59,-0.455732152091
> 68.7679083095,"H",0.58834202541,4.59,-0.531985064838
> 34.2857142857,"C",0.290215786413,3.63,-0.57640509953
> 10.2560914968,"A",-0.0610637969297,3.59,-0.546943907216
> 28.6842105263,"C",0.525118632142,3.49,-0.0523795789918
> 30.0416666667,"H",0.156954082926,3.52,-0.472482519814
> 20.1554404145,"C",0.0943522550034,4.59,-0.432462657283
> 12.7670396745,"C",0.100871137529,3.69,-0.544531967143
> 20.9363636364,"F",0.0329810347639,3.63,-0.366829418743
> 19.9712643678,"D",0.143507877945,4.59,1.27487137035
> 8.79711425206,"E",0.169545007823,4.59,-0.479695515413
> 13.5026041667,"A",0.0560834440921,2.12,3
> 16.1768953069,"A",0.146732872812,3.63,0.887818482431
> 15.8369098712,"D",-0.039775301243,3.51,-0.323236050367
> 8.27361563518,"E",0.0406370104999,3.55,0.0423073286075
> 18.7425860024,"E",0.100508655177,2.12,-0.291833520189
> 21.326754386,"A",-0.0666644294356,3.55,-0.377558236144
> 16.5308118683,"D",0.00488001447699,4.59,0.550971878424
> 19.7704081633,"B",0.295825553381,3.51,-0.507839519531
> 16.7927631579,"D",0.0412655840118,3.49,0.96097762139
> 58.32,"G",0.36793901326,2.12,-0.543313547312
> 18.2540983607,"F",0.130183791089,3.52,-0.254570872016
> 21.5294117647,"A",0.216310718816,3.63,-0.44837847382
> 32.6153846154,"B",0.616852170827,3.5,-0.137872313787
> 12.100456621,"F",-0.0716737993601,4.59,-0.217552654798
> 96.2121212121,"C",1.43708718338,3.49,-0.562412776152
> 44.3225806452,"F",0.455679605168,3.55,-0.0195847468614
> 35.1515151515,"H",0.0634403289952,3.49,-0.301121066447
> 18.9349112426,"B",0.419806161243,3.49,-0.481343868127
> 13.6689038031,"E",0.129475482986,3.69,0.173172092336
> 13.4612676056,"E",0.038580753911,3.51,0.209952518491
> 17.2413793103,"F",0.0818016458158,4.59,-0.532234060716
> 32.7247290378,"C",0.223591135664,3.59,-0.522165497386
> 32.8036322361,"D",0.159738245449,3.59,-0.469923672172
> 14.0371621622,"A",0.00786263593081,3.49,-0.522361789137
> 18.6205016358,"C",0.0697282036396,4.59,-0.559036392043
> 22.5321637427,"C",0.253176165958,3.49,-0.462649257589
> 12.8644939966,"A",-0.257191950463,3.69,-0.535866080593
> 15.3495440729,"E",0.0686445282036,4.59,-0.559465909933
> 20.5068356119,"E",0.195081088435,4.59,-0.425675859628
> 6.37312126076,"E",0.0533894440023,4.59,-0.480184792314
> 16.2802641233,"E",0.105766795295,4.59,-0.107738832602
> 31.4438502674,"H",0.296937975502,4.59,-0.0697960106888
> 16.4912917271,"E",0.161621584567,4.59,-0.405132454694
> 50.1695168467,"B",1.43722620782,2.12,-0.0963983153273
> 53.9838337182,"H",0.377434883805,4.59,-0.361241121249
> 19.0052356021,"C",0.114790302635,3.55,-0.509974244194
> 10.2944111776,"A",0.117306642214,2.12,0.390494449889
>
> --------------020205020409020506080205
> Content-Type: text/plain;
> name="test.r"
> Content-Transfer-Encoding: 7bit
> Content-Disposition: inline;
> filename="test.r"
>
> library(MASS)
>
> dat <- read.csv( "data.for.r.help.no.I.csv", header=TRUE)
>
> attach( dat)
>
> log1pX <- log(1 + X)
> logZ <- log(Z)
>
> rlm.out <- rlm( Y ~ -1+ Factor + log1pX + logZ + W, method="M")
> lm.out <-   lm( Y ~ -1+ Factor + log1pX + logZ + W)
>
> print( "before vcov( lm.out)" )
> vcov( lm.out)
>
> print( "before vcov(rlm.out)" )
> vcov( rlm.out)
>
> print( "after vcov(rlm.out)" )
>
>
>
> --------------020205020409020506080205
> Content-Type: text/plain;
> name="test.unix.out"
> Content-Transfer-Encoding: 7bit
> Content-Disposition: inline;
> filename="test.unix.out"
>
> brahe:/home/titan/frank $ R
>
> R : Copyright 2004, The R Foundation for Statistical Computing
> Version 1.9.1  (2004-06-21), ISBN 3-900051-00-3
>
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type 'license()' or 'licence()' for distribution details.
>
> R is a collaborative project with many contributors.
> Type 'contributors()' for more information and
> 'citation()' on how to cite R in publications.
>
> Type 'demo()' for some demos, 'help()' for on-line help, or
> 'help.start()' for a HTML browser interface to help.
> Type 'q()' to quit R.
>
>> version
>         _
> platform i386-pc-solaris2.8
> arch     i386
> os       solaris2.8
> system   i386, solaris2.8
> status
> major    1
> minor    9.1
> year     2004
> month    06
> day      21
> language R
>> source("test.r")
> [1] "before vcov( lm.out)"
> [1] "before vcov(rlm.out)"
> Error in if (rdf != z$df.residual) warning("inconsistent residual degrees of freedom. -- please report!") :
>        missing value where TRUE/FALSE needed
>>
>
> --------------020205020409020506080205
> Content-Type: text/plain;
> name="test.dos.out"
> Content-Transfer-Encoding: 7bit
> Content-Disposition: inline;
> filename="test.dos.out"
>
> H:\Notes>R
>
> R : Copyright 2004, The R Foundation for Statistical Computing
> Version 2.0.1  (2004-11-15), ISBN 3-900051-07-0
>
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type 'license()' or 'licence()' for distribution details.
>
> R is a collaborative project with many contributors.
> Type 'contributors()' for more information and
> 'citation()' on how to cite R or R packages in publications.
>
> Type 'demo()' for some demos, 'help()' for on-line help, or
> 'help.start()' for a HTML browser interface to help.
> Type 'q()' to quit R.
>
>> version
>         _
> platform i386-pc-mingw32
> arch     i386
> os       mingw32
> system   i386, mingw32
> status
> major    2
> minor    0.1
> year     2004
> month    11
> day      15
> language R
>> source("test.r")
> [1] "before vcov( lm.out)"
> [1] "before vcov(rlm.out)"
> Error in if (rdf != z$df.residual) warning("inconsistent residual degrees of fre
> edom. -- please report!") :
>        missing value where TRUE/FALSE needed
>>
>
>
>
>
>
>
> --------------020205020409020506080205--
>
> ______________________________________________
> R-devel at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>
>

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
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



More information about the R-devel mailing list