### R代码来自Vignette Source'ivas.rnw'######################################################################################################################################################################## 图书馆(ivas)############################################# ###代码块数字3:加载表达式#################################################################################### ###代码块编号4:加载SNP数据#######################################################################################位置数据#############################数据(SamplesNPlocus)######################################################## ###代码块数字6:加载TXDB数据################################-system.file(“ extdata”,“ sampledb”,package =“ ivas”)sample.txdb <-loaddb(sampledb)######################################################################,calgene = null,ncor = 1,out.dir = null)asdb head(插槽(asdb,“ plicingmodel”)$“ ass”)##################################################################sample.Txdb,ASdb=ASdb,Total.expdata=sampleexp, CalIndex="ASS7",Ncor=1,out.dir=NULL) ASdb head(slot(ASdb,"Ratio")$"ASS") ################################################### ### code chunk number 9: sQTLsFinder function ################################################### ASdb <- sQTLsFinder(ASdb=ASdb,Total.snpdata=samplesnp, Total.snplocus=samplesnplocus,method="lm",Ncor=1) ASdb head(slot(ASdb,"sQTLs")$"ASS") ################################################### ### code chunk number 10: Multi-threds function (eval = FALSE) ################################################### ## ASdb <- Splicingfinder(GTFdb=sample.Txdb,calGene=NULL,Ncor=4) ## ASdb <- RatioFromFPKM(GTFdb=sample.Txdb,ASdb=ASdb,Total.expdata=sampleexp,Ncor=4) ## ASdb <- sQTLsFinder(ASdb=ASdb,Total.snpdata=samplesnp, ## Total.snplocus=samplesnplocus,method="lm",Ncor = 4) ## ASdb ################################################### ### code chunk number 11: saveBplot function ################################################### saveBplot(ASdb=ASdb,Total.snpdata=samplesnp,Total.snplocus=samplesnplocus, CalIndex="ASS7",out.dir="./result") ################################################### ### code chunk number 12: sessionInfo ################################################### sessionInfo()