Displaying 9 results from an estimated 9 matches for "treat3".
Did you mean:
treat
2011 Aug 03
1
Coefficient names when using lm() with contrasts
...er=T)
> dt
ID VarA VarB VarC
1 1 cat 2 23
2 2 dog 3 56
3 3 rat 5 35
4 4 cat 2 43
5 5 cat 7 51
6 6 dog 3 31
7 7 dog 4 65
8 8 rat 1 18
9 9 rat 6 49
10 10 dog 3 28
> dt$VarA=relevel(dt$VarA,ref="dog")
> treat3=matrix(-1/3,ncol=2,nrow=3); for (i in 1:2) {treat3[i+1,i]=2/3}
> contrasts(dt$VarA)=treat3
> levels(dt$VarA)
[1] "dog" "cat" "rat"
> results=lm(formula=VarC~VarA+VarB, data=dt)
> names(results$coef)
[1] "(Intercept)" "VarA1" "...
2003 May 22
1
Experimental Design
...uld be greatly appreciated.
Isaac
A representation of the experiments follows:
Experiment 1
Control1 Treat1 Treat2
Blk1 s1 s2 s3
Blk2 s4 s5 s6
Blk3 s7 s8 s9
Experiment 2
Control2 Treat3 Treat4
Blk1 s1a s2a s3a
Blk2 s4a s5a s6a
Blk3 s7a s8a s9a
Control1 and Control2 I are the same control cell line. I would like to
compare Treat1 to Treat3 and Treat 4 and also I would like to compare
Treat2 to...
2007 Nov 12
0
Resid() and estimable() functions with lmer
...e (Intercept) 71.585 8.4608
Residual 49.856 7.0608
number of obs: 323, groups: initial.size, 292
Fixed effects:
Estimate Std. Error t value
(Intercept) 12.8846 1.3028 9.890
infl.treat1 -0.4738 1.1819 -0.401
def.treat2 -3.5522 1.6022 -2.217
def.treat3 -2.1757 1.6461 -1.322
def.treat4 -2.1613 1.7003 -1.271
Correlation of Fixed Effects:
(Intr) infl.1 df.tr2 df.tr3
infl.treat1 -0.413
def.treat2 -0.616 0.013
def.treat3 -0.641 0.002 0.493
def.treat4 -0.638 0.028 0.469 0.524
#Coefficients#
mean.infl.treat0=c(1,0,...
2004 Jan 27
2
Probability for ANOVA
Hi all!
I have 4 treatments on 5 animals
Treat1 Treat2 Treat3 Treat4
Animal1 36 37 35 39
Animal2 33 34 36 37
Animal3 37 35 33 38
Animal4 34 36 34 35
Animal5 35 36 33 36
I use an Anova and i try to verify calcul
So i retrieve:...
2007 Nov 01
2
F distribution from lme()?
...al" as random factor
> incub.lme2<-lme(egg.temp~kjday+treat,random=~1|ind,data=incub.df)
Fixed effects: egg.temp ~ kjday + treat
Value Std.Error DF t-value p-value
(Intercept) 24.937897 6.662475 11 3.743038 0.0032
kjday 0.108143 0.152540 7 0.708945 0.5013
treat3 -1.506605 0.485336 7 -3.104254 0.0172
#generating an anova table to get the F-distribution
> anova(incub.lme2)
numDF denDF F-value p-value
(Intercept) 1 11 1176.6686 <.0001
kjday 1 7 5.7060 0.0483
treat 1 7 9.6364 0.0172
>...
2007 Oct 31
0
Problems with generating F-distr from lme()
...l" as random factor
> incub.lme2<-lme(egg.temp~kjday+treat,random=~1|ind,data=incub.df)
Fixed effects: egg.temp ~ kjday + treat
Value Std.Error DF t-value p-value
(Intercept) 24.937897 6.662475 11 3.743038 0.0032
kjday 0.108143 0.152540 7 0.708945 0.5013
treat3 -1.506605 0.485336 7 -3.104254 0.0172
#generating an anova table to get the F-distribution
> anova(incub.lme2)
numDF denDF F-value p-value
(Intercept) 1 11 1176.6686 <.0001
kjday 1 7 5.7060 0.0483
treat 1 7 9.6364 0.0172
>...
2012 Nov 19
0
glht function in multcomp gives unexpected p=1 for all comparisons
...I tried the glht function in the multcomp package, but I get surprising results, with p=1 for all comparisons.?
# My data:
surv.data <- data.frame(
Site=c(rep("Site1", 9), rep("Site2", 9)),
Treatment=rep(c(rep("Treat1", 3), rep("Treat2", 3), rep("Treat3", 3)), 2),
survival=c(0.9, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.8, 0.4, 0.6, 0.5, 1.0, 0.7, 0.2, 0.2, 0.4)
)
# Binomial GLM:
glm.out <- glm(survival~Site*Treatment, data=surv.data, family="binomial", weights=rep(10, nrow(surv.data)))
anova(glm.out, test="Chisq")???...
2004 Jun 07
0
dfs in lme
...create response variable
data.df$zv <- rnorm (80, 2)
# create potential explanatory variables
data.df$explan <- rnorm (80, 2)
data.df$treat <- as.factor (sample (1:3, 80, T))
data.df$treat1 <- as.factor (sample (1:4, 80, T))
data.df$treat2 <- as.factor (sample (1:5, 80, T))
data.df$treat3 <- as.factor (sample (1:6, 80, T))
# with each of the explanatory variables
withoutInt <- lme (zv ~ explan, data= data.df, random= ~1 | gruppe)
withInt <- lme (zv ~ explan, data= data.df, random= ~ explan | gruppe)
anova (withoutInt)
anova (withInt)
anova (withoutInt, withInt)
There are...
2007 Oct 29
3
Strange results with anova.glm()
...poisson)
Deviance Residuals:
Min 1Q Median 3Q Max
-6.1541 -2.1872 -0.7045 0.9926 7.5629
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.8788185 0.1093306 17.185 < 2e-16 ***
infl.treat1 -0.2456909 0.0924677 -2.657 0.00827 **
def.treat2 -0.1478382 0.1277139 -1.158 0.24788
def.treat3 -0.0780282 0.1207796 -0.646 0.51871
def.treat4 -0.2581576 0.1221538 -2.113 0.03532 *
initial.size 0.0013834 0.0000993 13.931 < 2e-16 ***
---
Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
(Dispersion parameter for quasipoisson family taken to be 6.845172)
Null deviance: 3261.5 o...