frank@trdlnk.com
2005-Feb-25 23:04 UTC
[Rd] vcov on result of rlm() yields "-- please report!" (PR#7707)
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
Prof Brian Ripley
2005-Feb-26 08:56 UTC
(PR#7707) Re: [Rd] vcov on result of rlm() yields "-- please report!"
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@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@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>
>
--
Brian D. Ripley, ripley@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