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