2 Y 轴直方图（正常频率与相对频率）

Question

我需要你的帮助。

我有这两个地块，分别。一个是正​​常频率，另一个是具有完全相同数据的相对频率。

你能告诉我如何将它们加入一个带有 2 个 y 轴（频率和相对频率）的图中吗？

x<- AAA$starch 
h<-hist(x, breaks=40, col="lightblue", xlab="Starch ~ Corn", 
    main="Histogram with Normal Curve", xlim=c(58,70),ylim = c(0,2500),axes=TRUE)
xfit<-seq(min(x),max(x),length=40) 
yfit<-dnorm(xfit,mean=mean(x),sd=sd(x)) 
yfit <- yfit*diff(h$mids[1:2])*length(x)
lines(xfit, yfit, col="blue", lwd=3)

library(HistogramTools)
x<- AAA$starch 
c <- hist(x,breaks=10, ylab="Relative Frequency", main="Histogram with Normal Curve",ylim=c(0,2500), xlim=c(58,70), axes=TRUE)
PlotRelativeFrequency((c))

谢谢！！

编辑：

这只是我想要的示例图像......

Answer 1

我使用 latticeExtra 包中的 doubleYScale。

这是一个例子（我不确定相对频率计算）：

library(latticeExtra)

set.seed(42)
firstSet <- rnorm(500,4)

breaks = 0:10

#Cut data into sections
firstSet.cut = cut(firstSet, breaks, right=FALSE)
firstSet.freq = table(firstSet.cut) 

#Calculate relative frequency
firstSet.relfreq = firstSet.freq / length(firstSet) 

#Parse to a list to use xyplot later and assigning x values
firstSet.list <- list(x = 1:10, y = as.vector(firstSet.relfreq))

#Build histogram and relative frequency curve
hist1 <- histogram(firstSet, breaks = 10, freq = TRUE, col='skyblue', xlab="Starch ~ Corn", ylab="Frequency", main="Histogram with Normal Curve", ylim=c(0,40), xlim=c(0,10), plot=FALSE)
relFreqCurve <- xyplot(y ~ x, firstSet.list, type="l", ylab = "Relative frequency", ylim=c(0,1))

#Build double objects plot
doubleYScale(hist1, relFreqCurve, add.ylab2 = TRUE)

这是两个不同比例的 y 轴的结果：