r 可视化信令数据的条件概率分布

Posted 2021-05-24

data(tcells, package="bninfo")
library(bnlearn)
data(tcells, package="bninfo")
.tcells <- tcells$raw_data$cd3cd28icam2
.dtcells <- discretize(.tcells, method = "hartemink", breaks = 3, ibreaks = 50)
cdplot(.dtcells$Erk ~ .tcells$Mek, xlab = "Mek", ylab = "Erk", main = "original scale")
.dtcells <- discretize(log(.tcells), method = "hartemink", breaks = 3, ibreaks = 50)
cdplot(.dtcells$Erk ~ log(.tcells$Mek), xlab = "Mek", ylab = "Erk", main = "log scale")