计算2正态分布的积分-解释

问题描述

我想问一下如何解释作为2个正态分布之间的整数给出的数字。

我已经在R中准备了一个示例:



library(tidyverse)
set.seed(1)


get_overlap_coef <- function(mu1,mu2,sd1,sd2){
  
  A = rnorm(500,mu1,sd1)
  B = rnorm(500,sd2)
  
  xs  <- seq(min(mu1 - 4*sd1,mu2 - 4*sd2),max(mu1 + 4*sd1,mu2 + 4*sd2),length.out = 5000)
  f1  <- dnorm(xs,mean=mu1,sd=sd1)
  f2  <- dnorm(xs,mean=mu2,sd=sd2)
  int <- xs[which.max(pmin(f1,f2))]
  l   <- pnorm(int,lower.tail = mu1>mu2)
  r   <- pnorm(int,sd2,lower.tail = mu1<mu2)
  res <- l+r
  
  
  df <- 
    A %>% as.data.frame() %>% mutate(CAT = "BEFORE CHANGE") %>% 
    rbind(.,B %>% as.data.frame() %>% mutate(CAT = "AFTER CHANGE"))
  df %>% 
    ggplot(aes(.)) + 
    geom_density(aes(fill = CAT),alpha = 0.2) + 
    labs(title = paste("Percentage of Overlapping:",round(res,2)))
}

get_overlap_coef(20,30,5,5)

t.test(rnorm(50,20,5),rnorm(50,5))

我的问题是32%的人具有什么解释? 这不是双面成对t检验的p值。

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)