(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