From frank at trdlnk.com  Fri Feb 25 23:04:29 2005
From: frank at trdlnk.com (frank@trdlnk.com)
Date: Fri Feb 25 23:04:38 2005
Subject: [Rd] vcov on result of rlm() yields "-- please report!" (PR#7707)
Message-ID: <20050225220429.60023BAFE@slim.kubism.ku.dk>

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--


[Rd] vcov on result of rlm() yields "-- please report!" (PR#7707)
frank@trdlnk.com