Hist 与 R 中的线条

Question

我生成了4个大数据部分：cluster1（10000点），cluster2（15000点），cluster3（15000点）和throws（500点）。这是代码：

library('MASS')
library('fpc')
#library("dbscan")
library("factoextra")
library("clustertend")
library("boot")
library("stream")
set.seed(123)
mu1<-c(-5,-7)
mu1
sigma1<-matrix(c(4,-2,-2,2), nrow=2, ncol=2, byrow = TRUE)
sigma1
n<-10000
cluster1<-mvrnorm(n,mu1,sigma1)
cluster1
#cluster1<-as.data.frame(cluster1)
#cluster1
#c<-runif(10000,1,1000)
#c

phi <- runif(15000, max = 2*pi)
rho <- sqrt(runif(15000))
x <- sqrt(5)*rho*cos(phi) + 6
y <- sqrt(10/3)*rho*sin(phi) + 4
range(2*(x - 6)^2 + 3*(y - 4)^2)
#[1] 0.001536582 9.999425234
plot(x, y)
cluster2<-cbind(x,y)
cluster2

u <- runif(15000, max = 3)
v <- runif(15000, max = 2)
x <- u + v - 10
y <- v - u + 8
range(x + y)
#[1] -1.999774  1.999826
range(x - y + 15)
#[1] -2.999646  2.999692
plot(x, y)
cluster3<-cbind(x,y)
cluster3
#cluster3<-as.data.frame(cluster1)
#cluster3

x <- runif(500, -20, 20)
y <- runif(500, -20, 20)
#u <- runif(500, max = 20)
#v <- runif(500, max = 20)
#x <- u + v - 20
#y <- v - u
range(x)
range(y)
plot(x,y)
throws<-cbind(x,y)
throws

data<-rbind(cluster1,cluster2,cluster3,throws)
data<-as.data.frame(data)
data
plot(data)

然后我尝试用bootstrap的方法，构造一些H统计量的分布固定m，即来自生成点总数的7%（m=2835）。这是我执行此操作的代码：

B<-10#number of iterations
H<-NULL#value of Hopkins statistic
for(i in 1:B){
  N<-dim(data)[1]
  s<-sample(N,0.8*N)
  stat<-hopkins(data[s,], n=2835, byrow = TRUE)$H 
  H[i]<-stat
  #print(c(i, stat))
}

生成需要很长时间。然后我应该将这个结果与 beta 分布 - B(m,m) 进行比较。这是代码：

hist(H)
#(density(H), col="red")
#hist(distB)
X<-seq(min(H), max(H), 0.001)
X
lines(X, dbeta(X,2835,2835), type="l", col="red")

问题是 lined 没有根据 hist 绘制。谁能说出问题是什么？这是图像，我看到红线，但它并不完全正确。

Answer 1

dbeta() 绘制的 y 轴值太低，无法在提供的 y 轴上注册 (<0.0000001)。您需要叠加第二个图：

# sample data
H <- sample(seq(0.455,0.475,0.001), 1000, replace = TRUE)

#plot histogram
hist(H)

# prepare graphics to add second plot
par(new = TRUE)

# sample data for second plot
X <- seq(0.455,0.475, 0.001)
Y <- dbeta(X,2835,2835)

# plot second plot, remove axes
plot(X, dbeta(X,2835,2835), type="l", col="red", axes = FALSE)
axis(4, Y) # add axis on right side

Hist 与 R 中的线条

Hist with lines in R

r

plot

histogram

cluster-analysis

statistics-bootstrap