Barras rosas -> dados butland (995 instâncias) Barras roxas -> todos os 5 milhões
NÚMERO DE BINS = 20
fprob = read.csv("resumeprob.txt", header = FALSE, sep = "\t")
butprob = read.csv("prob_but.txt", header = FALSE, sep = "\t")
p1 <- hist(fprob$V3, breaks = 20, plot = FALSE)
p2 <- hist(butprob$V3, breaks = 20, plot = FALSE)
p1$density = p1$counts/sum(p1$counts) * 100
p2$density = p2$counts/sum(p2$counts) * 100
plot(p1, freq = F, col = rgb(0, 0, 1, 1/4), xlab = "Probability of being AGG",
ylab = "Density", main = "Histogram", ylim = c(0, 50))
plot(p2, freq = F, col = rgb(1, 0, 0, 1/4), add = T)
AUMENTEI O NÚMERO DE BINS = 100
p1 <- hist(fprob$V3, breaks = 100, plot = FALSE)
p2 <- hist(butprob$V3, breaks = 100, plot = FALSE)
p1$density = p1$counts/sum(p1$counts) * 100
p2$density = p2$counts/sum(p2$counts) * 100
plot(p1, freq = F, col = rgb(0, 0, 1, 1/4), xlab = "Probability of being AGG",
ylab = "Density", main = "Histogram")
plot(p2, freq = F, col = rgb(1, 0, 0, 1/4), add = T)
