acts=apply(matrix(rep(letters[1:4],7),nrow=2),2,paste,collapse=" ")[1:13]
reps=apply(matrix(rep(letters[5:8],7),nrow=2),2,paste,collapse=" ")[1:13]
grn=data.frame("Target"= LETTERS[1:26] ,"coact"=c(acts,reps),"corep"= c(reps,acts),"R2"=runif(26),stringsAsFactors=FALSE)
co=coregnet(grn)
samples= paste("S",1:100,sep="")
expression=matrix(rnorm(3400),ncol=100)
dimnames(expression) = list(c(grn$Target,names(regulators(co))),samples)
TFA = regulatorInfluence(co,expression,minTarg=4)
colnames(TFA) = samples
if(interactive()){
display(co,TFA=TFA,expressionData=expression)
}
CNA =matrix( sample(-2:2,800,replace=TRUE),ncol=100)
dimnames(CNA) = list(names(regulators(co)),samples)
if(interactive()){
display(co,TFA=TFA,expressionData=expression,alteration=CNA)
}
clinicGrp = factor(paste("grp",sample(1:3,100,replace=TRUE),sep=""))
names(clinicGrp) =samples
if(interactive()){
display(co,TFA=TFA,expressionData=expression,alteration=CNA,clinicalData=clinicGrp)
}
Run the code above in your browser using DataLab