R使用分布曲线创建直方图

Question

我想在同一窗口中使用分布曲线制作两个不同的直方图，其中包含来自txt文件的“data1”和“data2”。我可以创建一个直方图，其中包含“data1”和“data2”列的数据，但不是分开的。我如何分离我的数据？感谢。

data1  data2
155      130
195      10
21       26
15       210
190      15
2        205
182      50
115      55
170      1
17       56

Data = as.matrix( read.table( "c:\\Data.txt",header=TRUE ) )
attach( Data )

par( mfrow=c(1,2) ) #c(rows,columns)

hist(Data ,plot=TRUE, col=c("red","blue"),
     main = "Histogram of Data1",
     xlab="X-Axis", ylab="Y-Axis", cex.lab= 1, col.lab="blue" )

#Curve is not working even if data is combined
curve(dnorm(x, mean=mean(Data), sd=sd(Data)), add=TRUE, col="blue", lwd=2)

Answer 1

试试这个

par( mfrow=c(1,2) ) 
invisible(lapply(1:ncol(Data), function(i){
  x <- Data[,i]
  hist(x, col=c("red","blue"),
       main = paste0("Histogram of ", names(Data)[i]),
       xlab="X-Axis", ylab="Y-Axis", 
       cex.lab= 1, col.lab="blue",
       prob=TRUE)

  curve(dnorm(x, mean=mean(x), sd=sd(x)), 
        add=TRUE, col="blue", lwd=2)

}))

par( mfrow=c(1,1))

Answer 2

或试试这个：

Data <- structure(list(data1 = c(155, 195, 21, 15, 190, 2, 182, 115, 
                                 170, 17), data2 = c(130, 10, 26, 210, 15, 205, 50, 55, 1, 56)), .Names = c("data1", 
                                                                                                            "data2"), row.names = c(NA, -10L), class = "data.frame")
par( mfrow=c(1,2) ) 
invisible(lapply(1:ncol(Data), function(i){
  x <- Data[,i]
  h <- hist(x, col=c("red","blue"),
            main = paste0("Histogram of ", names(Data)[i]),
            xlab="X-Axis", ylab="Y-Axis", breaks=5,
            cex.lab= 1, col.lab="blue")
  xfit<-seq(min(x),ceiling(210/50)*50,length=40) 
  yfit<-dnorm(xfit,mean=mean(x),sd=sd(x)) 
  yfit <- yfit*diff(h$mids[1:2])*length(x) 
  lines(xfit, yfit, col="green", lwd=2)
}))

（source）