- > plotData[1:40,1]
- [1] 0.020 1.000 1.000 0.132 NA NA NA NA 0.000 0.072 1.000 0.000 0.000 0.999 1.000 0.751 NA 0.015
- [19] 0.903 0.000 0.040 0.001 0.995 0.054 0.010 0.058 1.000 0.010 NA 0.967 1.000 0.004 0.040 0.946 1.000 0.000
- [37] 0.020 0.629 1.000 0.414
- Score.value <- matrix(plotData[1:40,1],byrow=TRUE,40/4,4)
- z.value <- matrix(0,10,4)
- for(i in 1:4){
- for(j in 1:nrow(Score.value)){
- if(!is.na(Score.value[j,i])){
- if(Score.value[j,i]<=0.1){
- z.value[1,i] <- z.value[1,i]+1
- }
- if(Score.value[j,i]<=0.2 & Score.value[j,i]>0.1){
- z.value[2,i] <- z.value[2,i]+1
- }
- if(Score.value[j,i]<=0.3 & Score.value[j,i]>0.2){
- z.value[3,i] <- z.value[3,i]+1
- }
- if(Score.value[j,i]<=0.4 & Score.value[j,i]>0.3){
- z.value[4,i] <- z.value[4,i]+1
- }
- if(Score.value[j,i]<=0.5 & Score.value[j,i]>0.4){
- z.value[5,i] <- z.value[5,i]+1
- }
- if(Score.value[j,i]<=0.6 & Score.value[j,i]>0.5){
- z.value[6,i] <- z.value[6,i]+1
- }
- if(Score.value[j,i]<=0.7 & Score.value[j,i]>0.6){
- z.value[7,i] <- z.value[7,i]+1
- }
- if(Score.value[j,i]<=0.8 & Score.value[j,i]>0.7){
- z.value[8,i] <- z.value[8,i]+1
- }
- if(Score.value[j,i]<=0.9 & Score.value[j,i]>0.8){
- z.value[9,i] <- z.value[9,i]+1
- }
- if(Score.value[j,i]<=1 & Score.value[j,i]>0.9){
- z.value[10,i] <- z.value[10,i]+1
- }
- }
- }
- }
- #-----------------------------------------------------------------
- # plot and output plot
- #-----------------------------------------------------------------
- hist3D(1:nrow(z.value)/nrow(z.value),1:ncol(z.value),z = z.value/nrow(Score.value)*100, theta=25, phi=15,scale = TRUE, expand = 0.5,
- lwd = 2,ticktype="detailed", xlab=paste0("\n\n\n", "Score Range"),bty = "b2",facets = NA,colkey = FALSE,
- ylab=paste0("\n\n", "Score Classification"),zlab=paste0("\n\n", "Frequency"),axes=TRUE,alpha=1)
想将图中的Score classification中的1、2、3、4替换成score_LOD、score_LOH、score_SV、score_CNV
各位有修改过的吗?


雷达卡







京公网安备 11010802022788号







