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 && x<60,1,
ifelse(x>=60 && x<92,0,
ifelse(x>=92 && x<106,-1,
ifelse(x>=106 && x<145,0,
ifelse(x>=145 && x<169,1,
ifelse(x>=169 && 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:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|