绘制置信区间

时间:2015-08-01 19:18:37

标签: r graphics plot bioconductor confidence-interval

我正在策划Gviz手册示例:

library(Gviz)
data(twoGroups)

dTrack <- DataTrack(twoGroups, name = "uniform")

plotTracks(dTrack, groups = rep(c("control", "treated"),
each = 3), type = c("a", "p", "confint"))

以上代码应该绘制手册(http://www.bioconductor.org/packages/release/bioc/vignettes/Gviz/inst/doc/Gviz.pdf)第21页中的数字:

enter image description here

但是,我得到的是一个缺少阴影区域的图,它给出了每个组的置信区间('control'和'case'): enter image description here

“confint”类型不能再用于DataTrack功能了吗?

非常感谢!

更新@SabDeM评论:R和/或Gviz包的版本可能是原因。请检查您的版本:

[R

> R.Version()
$platform
[1] "x86_64-w64-mingw32"

$arch
[1] "x86_64"

$os
[1] "mingw32"

$system
[1] "x86_64, mingw32"

$status
[1] ""

$major
[1] "3"

$minor
[1] "1.1"

$year
[1] "2014"

$month
[1] "07"

$day
[1] "10"

$`svn rev`
[1] "66115"

$language
[1] "R"

$version.string
[1] "R version 3.1.1 (2014-07-10)"

$nickname
[1] "Sock it to Me"

会话信息:

> sessionInfo()
R version 3.1.1 (2014-07-10)
Platform: x86_64-w64-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252   
[3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C                          
[5] LC_TIME=English_United States.1252    

attached base packages:
 [1] stats4    parallel  grid      stats     graphics  grDevices utils     datasets  methods  
[10] base     

other attached packages:
[1] installr_0.16.0      BiocInstaller_1.16.5 data.table_1.9.4     Gviz_1.10.11        
[5] GenomicRanges_1.18.4 GenomeInfoDb_1.2.5   IRanges_2.0.1        S4Vectors_0.4.0     
[9] BiocGenerics_0.12.1 

0 个答案:

没有答案