Displaying 1 result from an estimated 1 matches for "logz".
Did you mean:
log
2005 Feb 25
1
vcov on result of rlm() yields "-- please report!" (PR#7707)
...sistent 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)"...