Displaying 3 results from an estimated 3 matches for "upper2".
Did you mean:
upper
2010 Jun 13
1
using latticeExtra plotting confidence intervals
...f the formula. This is the example code:
library(lattice)
library(latticeExtra)
sample1<-rnorm(100,10,2)
sample2<-rnorm(100,50,3)
sample3<-rnorm(100,20,2)
sample4<-rnorm(100,40,1)
mu1<-mean(sample1)
ci.upper1<-mu1+2*2
ci.lower1<-mu1-2*2
mu2<-mean(sample2)
ci.upper2<-mu2+2*3
ci.lower2<-mu2-2*3
mu3<-mean(sample3)
ci.upper3<-mu3+2*2
ci.lower3<-mu3-2*2
mu4<-mean(sample4)
ci.upper4<-mu4+2*1
ci.lower4<-mu4-2*1
categories<-c("A","B")
mu<-cbind(mu1,mu2,mu3,mu4)
ci.upper<-cbind(ci.upper1,ci.upper2,ci.up...
2008 May 13
1
R help: problems with step function
...AICs (* and **) are different.
I noticed that the Df for rs13482096:rs13483699 is 4, while I think
Df should be 6, 2 from rs13483699 and 4 from interactions. When I ran
add1 directly, I got Df=6 and AIC 848.75.
> step2.bic.out <- step(step.bic.out, scope=list(lower=scope.lower2,
upper=scope.upper2),
+ direction="both", k=log(length(step.bic.out$y)),
trace=1)
Start: AIC=841.18
pheno.dat ~ rs13479085 + rs13480057 + rs13482096 + rs8254221
Df Deviance AIC
+ rs13482096:rs13483699 4 216.63 840.63 (*)
<none>...
2008 May 14
0
Problems with step function
...AICs (* and **) are different.
I noticed that the Df for rs13482096:rs13483699 is 4, while I think
Df should be 6, 2 from rs13483699 and 4 from interactions. When I ran
add1 directly, I got Df=6 and AIC 848.75.
> step2.bic.out <- step(step.bic.out, scope=list(lower=scope.lower2,
upper=scope.upper2),
+ direction="both", k=log(length(step.bic.out$y)),
trace=1)
Start: AIC=841.18
pheno.dat ~ rs13479085 + rs13480057 + rs13482096 + rs8254221
Df Deviance AIC
+ rs13482096:rs13483699 4 216.63 840.63 (*)
<none>...