Displaying 3 results from an estimated 3 matches for "mod6".
Did you mean:
mod
2010 Mar 06
1
Robust SE for lrm object
I'm trying to obtain the robust standard errors for a multinomial ordered logit model:
mod6 <- lrm(wdlshea ~ initdesch + concap + capasst + qualrat + terrain,data=full2)
The model is fine but when I try to get the RSE I get an error.
coeftest(mod6, vcov = vcovHAC(mod6))
Error in match.arg(type) :
'arg' should be one of “ordinary”, “score”, “score.binary”, “pearson”, “d...
2010 Nov 15
1
Executing Command on Multiple R Objects
...want to print a number of results from lme function objects out to a txt
file. How could I do this more efficiently than what you see here:
out2 <- capture.output(summary(mod2a))
out3 <- capture.output(summary(mod3))
out4 <- capture.output(summary(mod5))
out5 <- capture.output(summary(mod6))
out6 <- capture.output(summary(mod7))
cat(out2,file="out.txt",sep="\n",append=TRUE)
cat(out3,file="out.txt",sep="\n",append=TRUE)
cat(out4,file="out.txt",sep="\n",append=TRUE)
cat(out5,file="out.txt",sep="\n",appen...
2002 Dec 18
6
Can I build an array of regrssion model?
...:cuts[1]]);
mod2 =
lm(rawData[[2]][cuts[1]:cuts[2]]~rawData[[1]][cuts[1]:cuts[2]]);
mod3 =
lm(rawData[[2]][cuts[2]:cuts[3]]~rawData[[1]][cuts[2]:cuts[3]]);
mod4 =
lm(rawData[[2]][cuts[3]:cuts[4]]~rawData[[1]][cuts[3]:cuts[4]]);
mod5 =
lm(rawData[[2]][cuts[4]:cuts[5]]~rawData[[1]][cuts[4]:cuts[5]]);
mod6 =
lm(rawData[[2]][cuts[5]:cuts[6]]~rawData[[1]][cuts[5]:cuts[6]]);
mod7 =
lm(rawData[[2]][cuts[6]:cuts[7]]~rawData[[1]][cuts[6]:cuts[7]]);
plot(rawData[[1]],rawData[[2]],type='l', col="green", xlab="Da/z",
ylab="m/z");
abline(mod1, lty="1", col="...