8 lines
136 B
Text
8 lines
136 B
Text
dbenford <- function(x){
|
|
log10(1 + 1/x)
|
|
}
|
|
|
|
pbenford <- function(q){
|
|
cumprobs <- cumsum(dbenford(1:9))
|
|
return(cumprobs[q])
|
|
}
|