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 [EMAIL PROTECTED] 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 Rsource("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 Rsource("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@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
-- Brian D. Ripley, [EMAIL PROTECTED] 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
______________________________________________ R-devel@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
