Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- c("00_11242T1.bmp", "00_7T.bmp", "01_677T.bmp", "106T.bmp", "106TV.bmp",
- "108T.bmp", "108TV.bmp", "124T.bmp", "124TV.bmp", "40T.bmp",
- "4497T.bmp", "44T.bmp", "44TV.bmp", "511T.bmp", "511TV.bmp",
- "514T.bmp", "514TV.bmp", "56T.bmp", "92_11145T.bmp", "94_10917T1.bmp",
- "95_549T.bmp", "97_12714T1.bmp", "97_7999T.bmp", "97_8073T2.bmp",
- "99_2221T.bmp", "99_6669T.bmp", "99_7417T1.bmp", "99_7417T2.bmp",
- "R01_80418T2.bmp", "R01_81197T.bmp", "R02_80456T2.bmp", "R03_80356T.bmp",
- "R03_80586T.bmp", "R04_80227T.bmp", "R04_80577T.bmp", "R04_80584T.bmp",
- "R04_81371T.bmp", "R04_81372T1.bmp", "R04_81449T.bmp", "R05_80479T.bmp",
- "R05_80481T.bmp", "R05_80611T.bmp")
- for( i in (1:ncol(data.lair))) {
- bmp(filename = bitlist[i], width = 1200, height = 1200, units = "px",bg = "white")
- par(mfrow=c(1,2))
- plot(density(assayData(data.lair)$lair.predicted[,i],na.rm=TRUE),main = bitlist[i])
- plot(density(assayData(data.lair)$predicted[,i],na.rm=TRUE),main = bitlist[i])
- cat(i,"n")
- dev.off()
- }
- i <- 2
- dat <- assayData(data.lair)$lair.predicted[ , i]
- any(!is.na(dat))
- # [1] FALSE
- for( i in (1:ncol(data.lair))) {
- bmp(filename = bitlist[i], width = 1200, height = 1200, units = "px",bg = "white")
- if (sum(!(is.na(assayData(data.lair)$lair.predicted[,i]))) > 1) {
- par(mfrow=c(1,2))
- plot(density(assayData(data.lair)$lair.predicted[,i],na.rm=TRUE))
- }
- plot(density(assayData(data.lair)$predicted[,i],na.rm=TRUE))
- cat(i,"n")
- dev.off()
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement