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