# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块数量1:R。隐藏# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #库(splicegear) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块2号:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #数据(spsite)打印(spsite)情节(spsite) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块3号:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #库文件名(XML) < -系统。文件(“数据”、“例子。xml”,包= " xml < splicegear”)——xmlTreeParse(文件名,asTree = TRUE) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块数量4:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # spsite <——buildSpliceSites (xml、verbose = FALSE)长度(spsite) (spsite[1:2]) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块5号:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #建SpliceSites图书馆(xml)文件名< -系统。文件(“数据”、“例子。xml”,包= " xml < - xmlTreeParse splicegear”)(文件名,asTree = TRUE) spsite < - buildSpliceSites (xml、verbose = FALSE) # #第二子集对象在我的列表中。spsite < - spsite[[2]] # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块6号:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #地块(my.spsites) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块7号:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #数据(spliceset) dataf <——as.data.frame (spliceset) colnames (dataf) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块8号:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # lm。小组< -函数(x, y,…){点p (x, y,…)。lm < - lm (y ~ x);abline (p.lm)} # #受制于情节探针强度值的位置探测器在mRNA: # # # #(注释掉,以避免一个警告)# # coplot(日志(exprs) ~材料|开始,数据= dataf面板= lm.panel) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块9号:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # 10 bp seq的窗口。长度< - as.integer(10) # #位置的外显子spsiteIpos < -矩阵(c(1、3.5、5、9、3、4、8、10),数控= 2)# #已知变量变异< -列表(= c (1、2、3、4), b = c (1、2、3), c = c (1、3、4)) # # n。外显子< - nrow (spsiteIpos) spvar < -新(“SpliceSitesGenomic spsiteIpos = spsiteIpos变异=的变异,seq.length = seq.length) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #块10号:# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # par (mfrow = c (3,1), 3 = c(3.1, 2.1, 2.1, 1.1))情节(spvar分裂= TRUE, col.exon =彩虹(n.exons))