找回密码
 注册
查看: 449|回复: 0

R语言 WaveCD包 PlotBreakPoints()函数中文帮助文档(中英文对照)

[复制链接]
发表于 2012-10-1 16:52:36 | 显示全部楼层 |阅读模式
PlotBreakPoints(WaveCD)
PlotBreakPoints()所属R语言包:WaveCD

                                        Plot the series and detected jump points
                                         绘制一系列检测到的跳跃点

                                         译者:生物统计家园网 机器人LoveR

描述----------Description----------

A lattice style plot is used to show break points for each chromosome
的格子风图是用来显示每个染色体的破发点,


用法----------Usage----------


PlotBreakPoints(x, Chromosome, layout = c(length(unique(Chromosome)),1), ...)



参数----------Arguments----------

参数:x
any series of observations
任何系列的意见


参数:Chromosome
chromosome corresponding to data in x
在x中的数据相应的染色体


参数:layout
the layout of the lattice plot
的晶格图的布局


参数:...
arguments that can be passed to plot command
可以通过绘制命令的参数


Details

详细信息----------Details----------

Uses MODWT method described in Islam (2008).
使用MODWT伊斯兰教(2008年)中描述的方法。


值----------Value----------

lattice plot
格子图


(作者)----------Author(s)----------


M. Shahidul Islam



参考文献----------References----------

M.S. Islam (2008). Periodicity, Change Detection and Prediction in Microarrays. Ph.D. Thesis, University of Western Ontario.

参见----------See Also----------

BreakPoints
BreakPoints


实例----------Examples----------


# Example-1: Data simulated from AR(1) process with phi=0.8[例1:数据AR(1)过程的模拟与披= 0.8]
phi<-0.8
SigE<- 0.12
ex.c1<- c(rep(0, 80), rep(-0.7,30), rep(0,40))+ SimulateError(phi, SigE, 150)
ex.c2<- c(rep(0,40), rep(0.7, 30), rep(0, 50))+ SimulateError(phi, SigE, 120)
ex.c<-c(ex.c1, ex.c2)
PlotBreakPoints(ex.c, c(rep(1, 150),rep(2, 120)))


# Example-2: seven jump points[例2:7个跳跃点]

fx<-function(x){ifelse(x<50,0,
        ifelse(x>=50 &amp;&amp; x<60,1,
          ifelse(x>=60 &amp;&amp; x<92,0,
            ifelse(x>=92 &amp;&amp; x<106,-1,
              ifelse(x>=106 &amp;&amp; x<145,0,
                ifelse(x>=145 &amp;&amp; x<169,1,
                  ifelse(x>=169 &amp;&amp; x<=180,2, 0)))))))+rnorm(1,0,0.1)}


set.seed(123)
y<-sapply(1:200,fx)
PlotBreakPoints(y, c(rep(1, length(y))), type="l", lwd=1.7)


#[]
# CGH Data[CGH数据]
#[]
data(arrayCGH)
# Plotting first 12 chromosomes [绘制第12条染色体]
PlotBreakPoints(arrayCGH[1:765,1], arrayCGH[1:765,2], type="l",layout=c(3,4))
# Plotting rest of the chromosomes[绘制其余染色体]
PlotBreakPoints(arrayCGH[766:1040,1], arrayCGH[766:1040,2], type="l",layout=c(3,4))


转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。


注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|小黑屋|生物统计家园 网站价格

GMT+8, 2024-11-25 21:48 , Processed in 0.035755 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表