0

私はRを初めて使用し、次の数値のヒストグラムをプロットしようとしています。

Input.txt

a b c 1.111 
d e f 2.3433 
g h i 9.87667 
............ 
............

合計レコード-2444910

最大値-142701.38999976

だから、これまで私はこれを試しました

setwd ("/path/")
data <- read.table ("input.txt", sep="\t", header=F)
x <-data$V4
bins = seq(0,150000,by=10000)
h <- hist(x, breaks=bins, probability=TRUE)
s = sd(x)
m = mean(x)
curve(dnorm(x, mean=m, sd=s), add=TRUE, lwd=3, col="red")
lines(density(x), col="blue")
abline(v=mean(x),col="blue")
mtext(paste("mean ", round(mean(x),1), "; sd ", round(sd(x),1), "; N ",length(x),sep=""), side=1, cex=.75)
dev.off()

しかし、(0,0)座標の近くに1つのヒストグラムバーを取得できました。どんな助けでも大歓迎です。

4

1 に答える 1

0
dat <- data.frame(V4=c(runif(10000, 1,5000), 149999))
hist(dat$V4, breaks=c(seq(0,5000,100), 150000) )

外れ値を除外するには:

hist(dat$V4[dat$V4 <5000], breaks=c(seq(0,5000,100)) )
于 2012-06-12T15:54:10.740 に答える