Displaying 7 results from an estimated 7 matches for "edadysexo".
2012 Jan 11
2
problems with glht for ancova
I've run an ancova, edadysexo is a factor with 3 levels,and log(lcc) is the
covariate (continous variable)
I get this results
> ancova<-aov(log(peso)~edadysexo*log(lcc))
> summary(ancova)
Df Sum Sq Mean Sq F value Pr(>F)
edadysexo 2 31.859 15.9294 803.9843 <2e-16 ***
log(lcc)...
2011 Jul 24
1
GLM different results with the same factors
...ow can i avoid this
problem.
Why the order of the factors give different results? I suppose it's because
the order of the factors, i've just changed "lcc" from the first position to
the last in the model, and the significance change completely
> modo<-glm(prevalencia~*lcc*+edadysexo/lcc+edadysexo/mes,binomial)
> anova(modo,test="Chisq")
Df Deviance Resid. Df Resid. Dev P(>|Chi|)
NULL 524 206.97
lcc 1 10.5715 523 196.40 0.001148 **
edadysexo 2 1.0725 521 195.32...
2012 Jan 10
1
different results from fligner.test
I've made fligner test with the same data, changing the orders of the
variables, and this what i get
> fligner.test(rojos~edadysexo*zona*ano*estacion)
Fligner-Killeen test of homogeneity of variances
data: rojos by edadysexo by zona by ano by estacion
Fligner-Killeen:med chi-squared = 15.7651, df = 2, p-value = 0.0003773
> fligner.test(rojos~ano*edadysexo*zona*estacion)
Fligner-Killeen test of homogenei...
2011 Jan 28
1
Problems making subsets with [] or "subset"
...cted individual
0.- Non infected
So what i'm trying to do is a subset in a model like this,
model<-aov(intensity~ageandsex*month*zone*year,subset=(prevalence=="1"))
Is a correct way?
I prefer to make the subsets with [], but doesn't work, why
> model<-aov(hemogregarinas~edadysexo*mes*zona*ano,[prevalencia=="1"])
Error: inesperado '[' en
"model<-aov(hemogregarinas~edadysexo*mes*zona*ano,["
And finally, look at this, i don't know why I get different results
> summary(hemogregarinas[prevalencia=="1"])
Min. 1st Qu. Median...
2012 Jan 12
0
glht (multicomparisons) with an interaction factor
Hi,
i was working with this model
> mq<-glm(rojos~edadysexo*zona*estacion,quasipoisson)
and i get this minimal adequate model
> anova(mq5,test="F")
Df Deviance Resid. Df Resid. Dev F Pr(>F)
NULL 518 64799
edadysexo 2 1556.5 516 63243 8.94...
2011 Jan 13
2
send commands from script to console
Hi,
I'm working with R in windows and I wonder if there is any command (of the
kind CTRL+ ) to transfer the commands I've worked with (like: >
model<-glm(prevalencia~edadysexo*mes*zona,binomial)) to a script
automatically, without the results I received from R after execute them.
Another question, what's the combination of CTRL+ for transfer the orders in
the other way, fom the script to the console?
Thanks very much
Sorry for the explanation, but I don't know t...
2012 Feb 13
0
pairwise comparisons with multcomp package
Hi,
I've got this model and following Hothorn et al advices, I used glht for a
post hoc comparison
> modezqM<-glm(rojos~estacion*zona3,quasipoisson,subset=(edadysexo=="M"))
> anova(modezqM,test="F")
Df Deviance Resid. Df Resid.
Dev F Pr(>F)
NULL 293 41148
estacion 1 3647.0 292...