Displaying 1 result from an estimated 1 matches for "cultb".
Did you mean:
cult
2011 Aug 06
1
How set lm() to don't return NA in summary()?
...>
> # plot
> require(lattice)
> xyplot(y~fert|cult, da)
>
> # adjust
> m0 <- lm(y~cult*fert, da)
> summary(m0)
...
Coefficients: (1 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 7.55401 10.18956 0.741 0.469
cultB -8.04486 13.54672 -0.594 0.561
cultC -3.83644 6.74848 -0.568 0.578
fert 0.87486 0.08265 10.586 1.24e-08 ***
cultB:fert 0.13509 0.11688 1.156 0.265
cultC:fert NA NA NA NA
...
> require(gmodels)
> estimable(m0, cm=c(1,0,0,...