0

我是 R 新手,我正在尝试为以下数字绘制直方图

输入.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) 坐标附近获得一个直方图条。任何帮助表示赞赏。

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 回答