Displaying 2 results from an estimated 2 matches for "water2".
Did you mean:
water
2009 Apr 01
1
Noobie ANOVA intercept question
...dian 3Q Max
-1.23992 -0.18979 -0.01840 0.17471 0.74686
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.7765 0.2203 17.139 2.23e-16 ***
BED2 0.3185 0.1908 1.669 0.106242
BED3 -0.5044 0.1908 -2.643 0.013293 *
WATER2 0.7842 0.1908 4.109 0.000313 ***
WATER3 0.4489 0.1908 2.353 0.025905 *
SHADE2 0.1969 0.2203 0.894 0.379172
SHADE3 -0.2157 0.2203 -0.979 0.336068
SHADE4 -0.3673 0.2203 -1.667 0.106641
---
Signif. codes: 0 ?***? 0.001 ?**? 0....
2012 Mar 03
1
Problem running stepAIC within a function.
...on
if (opt==1){
fit <-
lm(stepAIC(lm(new~-1+SIC,s1,weights=w),scope=list(upper=~1+SIC+turn+exturn+new+exnew+hand+excise+merch+exmerch+tob+rub+water+turn1+exturn1+new1+exnew1+hand1+excise1+
merch1+exmerch1+tob1+rub1+water1+turn2+exturn2+new2+exnew2+hand2+excise2+merch2+exmerch2+tob2+rub2+water2,lower=~-1+SIC),
direction="both"))
sigdata <-
cbind(new=s1$new,SIC=s1$SIC,s1[,coeffs[which(summary(fit)$coeff[-c(1:128),4]
<= 0.05)]])
finalfit <- lm(new~.-1,data=sigdata,weights=w)
}
else {
fit <-
lm(stepAIC(lm(new~-1+SIC3+emp,data=s1),scope=list(upper=~1+SIC...