R: genera dati da una distribuzione di densità di probabilità

Question

Dire che ho un array semplice, con una distribuzione di probabilità corrispondente.

library(stats)  
data <- c(0,0.08,0.15,0.28,0.90) 
pdf_of_data <- density(data, from= 0, to=1, bw=0.1) 

C'è un modo per generare un altro insieme di dati utilizzando la stessa distribuzione. Poiché l'operazione è probabilistica, non è necessario che corrisponda esattamente alla distribuzione iniziale, ma verrà generata da essa.

Ho avuto successo nel trovare una soluzione semplice per conto mio. Grazie!

Answer 1

La cosa migliore è quello di generare la funzione di densità cumulativa empirica, approssimare l'inverso, e poi trasformare l'ingresso.

L'espressione composto si presenta come

random.points <- approx(
    cumsum(pdf_of_data$y)/sum(pdf_of_data$y), 
    pdf_of_data$x, 
    runif(10000) 
)$y 

rendimenti

hist(random.points, 100) 

Answer 2

Dagli esempi nella documentazione di ?density si ottiene (o quasi) la risposta.

Quindi, qualcosa di simile a questo dovrebbe farlo:

library("stats")  
data <- c(0,0.08,0.15,0.28,0.90) 
pdf_of_data <- density(data, from= 0, to=1, bw=0.1) 

# From the example. 
N <- 1e6 
x.new <- rnorm(N, sample(data, size = N, replace = TRUE), pdf_of_data$bw) 

# Histogram of the draws with the distribution superimposed. 
hist(x.new, freq = FALSE) 
lines(pdf_of_data) 

È appena possibile rifiutare l'pareggi fuori della vostra intervallo come in rejection sampling. In alternativa, è possibile utilizzare l'algoritmo descritto nel link.

Answer 3

Per disegnare dalla curva:

sample(pdf_of_data$x, 1e6, TRUE, pdf_of_data$y)