Finding среза между двумя участками перекрывающихся плотности

Question

У меня есть эти данные выборки экспрессии двух различных вариантов гена:

value<-cbind(c(rnorm(100,500,90),rnorm(100,800,120))) 
genotype<-cbind(c(rep("A",100),rep("B",100))) 
df<-cbind(value,genotype) 
df<-as.data.frame(df) 
colnames(df)<-c("value","genotype") 
df$value<-as.numeric(as.character(df$value)) 

Я нанесенными эти два варианта генотипа их выражения и пытаюсь определить оптимальное значение отсечки анализ, который проводит различие между ними:

d <- density(value) 
plot(d, main="Genotypes A and B", ,type="n",xlim=c(200,1100),ylim=c(0,0.005),xlab="Units of expression",ylab="") 
d1 <- density(subset(value,genotype=="A")) 
polygon(d1, col = adjustcolor('gray', alpha.f = .40), border="black") 
d2 <- density(subset(value,genotype=="B")) 
polygon(d2, col = adjustcolor('gray', alpha.f = .40), border="black") 

Очевидно, что я могу использовать «abline» функцию, чтобы найти наилучшую отсечку между двумя плотностями, но есть аккуратнее способ определить обрезание?

Answer 1

Вот ответ, основанный на ссылке, предоставленной @Julius:

A.density <- density(subset(df, genotype == "A")$value, from = min(df$value), to = max(df$value), n = 2^10) 
B.density <- density(subset(df, genotype == "B")$value, from = min(df$value), to = max(df$value), n = 2^10) 
intersection.point <- A.density$x[which(diff((A.density$y - B.density$y) > 0) != 0) + 1]