Вопросы Теги

Как построить гистограмму с чрезвычайно распространенными значениями / пробелами в графике или ggplot

Это вызвано тем, что есть файл класса, от которого зависит ваш код, и он присутствует во время компиляции, но не найден во время выполнения. Ищите различия в времени сборки и классах времени выполнения.

0
задан Davide Passaretti 5 March 2019 в 07:39
поделиться

1 ответ

Как насчет простого разбиения множества следующим образом: 

obs= c(rnorm(100, 0, 1), rnorm(100, sample(c(-1e6, 1e6), 1), 1))

histN <- function(x, n){
  clustRes <- kmeans(x,centers = n)
  par(mfrow=c(1,n))
  for(i in 1:n){
    hist(x[clustRes$cluster==i], xlab="", main = sprintf("subset = %d",i))
  }
}
dev.new()
histN(obs,2)

или как более симпатичная, но менее гибкая альтернатива: 

hist2 <- function(x){
  clustRes <- kmeans(x,centers = 2)
  parDefault <- par()
  layout(matrix(c(1,1,1,2,3,3,3),nrow=1))
  idx1 <- which.min(clustRes$centers)
  idx2 <- which.max(clustRes$centers)
  h1 <- hist(x[clustRes$cluster==idx1],plot=FALSE)
  h2 <- hist(x[clustRes$cluster==idx2],plot=FALSE)
  yRange <- c(min(c(h1$counts,h2$counts)),max(c(h1$counts,h2$counts)))
  par(mai = c(parDefault$mai[1:3],0))
  par(cex = parDefault$cex)
  hist(x[clustRes$cluster==idx1],ylim= yRange, xlab="", ylab ="", main = sprintf("subset = %d",1),axes=FALSE)
  axis(1)
  axis(2)
  par(mai = c(parDefault$mai[1],0,parDefault$mai[3],0))
  par(cex = 3)
  plot(1:10,type="n", axes=FALSE, xlab="...", ylab ="")
  text(x= 5, y= 5, "...")
  axis(1, tick=FALSE,labels = FALSE)
  par(mai = c(parDefault$mai[1],0,parDefault$mai[3:4]))
  par(cex = parDefault$cex)
  hist(x[clustRes$cluster==idx2],ylim= yRange, xlab="", ylab ="", main = sprintf("subset = %d",2),axes=FALSE)
  axis(1)
  axis(4)
}
dev.new()
hist2(obs)
0
ответ дан HolgerBarlt 5 March 2019 в 07:39
поделиться
Другие вопросы по тегам:

Похожие вопросы: