Displaying 2 results from an estimated 2 matches for "treeatment".
Did you mean:
treatment
2006 Jul 28
0
tests performed by anova
...the columns coding for factor b
a:b : the removal of the columns coding for the interaction of factors a and b
BUT :
print(anova(mod))
does not provide the expected results. Moreover, the results of
"anova" do not depend on the contrasts, since setting
options(contrasts=c("contr.treeatment", "contr.poly"))
leads to the same result.
I manage to obtain the desired results by creating the submodels
without the corresponding columns, and by comparing them with "anova"
, for example :
xsub = mod$x[,-c(2)]
submod = lm(Y ~ xsub-1)
print(anova(submod,mod))
So m...
2011 Mar 25
2
A question on glmnet analysis
Hi,
I am trying to do logistic regression for data of 104 patients, which
have one outcome (yes or no) and 15 variables (9 categorical factors
[yes or no] and 6 continuous variables). Number of yes outcome is 25.
Twenty-five events and 15 variables mean events per variable is much
less than 10. Therefore, I tried to analyze the data with penalized
regression method. I would like please some of the