[BioC] Time course....makeConstrasts()
Khan, Sohail
khan at cshl.edu
Tue Mar 14 22:30:50 CET 2006
Thanks Sean.
My code and the sessionInfo:
Thanks Sean.
Here is my code and sessionInfo.
targets
FileName Target
1 DF_02.CEL 0hr
2 DF_04.CEL 2hr
3 DF_05.CEL 2hr
4 DF_07.CEL 4hr
5 DF_08.CEL 4hr
> lev <- c("0hr","2hr","4hr")
> lev
[1] "0hr" "2hr" "4hr"
> f <- factor(targets$Target, levels=lev)
> f
[1] 0hr 2hr 2hr 4hr 4hr
Levels: 0hr 2hr 4hr
> design <- model.matrix(~0+f)
> colnames(design) <- lev
> design
0hr 2hr 4hr
1 1 0 0
2 0 1 0
3 0 1 0
4 0 0 1
5 0 0 1
attr(,"assign")
[1] 1 1 1
attr(,"contrasts")
attr(,"contrasts")$f
[1] "contr.treatment"
> fit <- lmFit(eset, design)
> cont.wt <- makeContrasts(
+ "2hr-0hr",
+ "4hr-2hr",
+ levels=design)
Error in parse(file, n, text, prompt) : parse error
> sessionInfo()
R version 2.1.1, 2005-06-20, i386-pc-mingw32
attached base packages:
[1] "splines" "tools" "methods" "stats" "graphics" "grDevices" "utils" "datasets" "base"
other attached packages:
hgu133a2probe hgu133a2cdf limma affydata affyPLM gcrma matchprobes affy reposTools Biobase
"1.10.0" "1.10.0" "2.0.8" "1.4.1" "1.3.3" "1.1.4" "1.0.22" "1.6.7" "1.5.19" "1.5.12"
>
-Sohail
More information about the Bioconductor
mailing list