Chapter 4 issues

[bonnfire]

When introducing rbinom() and runif(), should I specify the arguments (i.e. x=, prob=, etc.)? In the code chunk or in the text after the code chunk?

[nicholasjhorton]

Please name these arguments. After usage, you can call x x and size n for rbinom(). You'll also want to add in a simple example of xpbinom(). Nick

On Jul 17, 2018, at 3:57 PM, Bonnie Lin ***@***.***> wrote: When introducing rbinom() and runif(), should I specify the arguments (i.e. x=, prob=, etc.)? In the code chunk or in the text after the code chunk? — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or mute the thread.

Nicholas Horton Professor of Statistics Department of Mathematics and Statistics, Amherst College PO Box 5000, AC #2239 Amherst, MA 01002-5000

[bonnfire]

Thanks for the response! I will specify the arguments in the function call.

[bonnfire]

Another question that I had was regarding the Benford's law histogram. Would you like me to use the Benford's base R function to plot the probability histogram?

[nicholasjhorton]

That sounds like a good plan. Nick

On Jul 17, 2018, at 4:25 PM, Bonnie Lin ***@***.***> wrote: Another question that I had was regarding the Benford's law histogram. Would you like me to use the Benford's base R function to plot the probability histogram? — You are receiving this because you commented. Reply to this email directly, view it on GitHub, or mute the thread.

Nicholas Horton Professor of Statistics Department of Mathematics and Statistics, Amherst College PO Box 5000, AC #2239 Amherst, MA 01002-5000

[bonnfire]

Goal is to recreate these two plots:

For the first one, I've tried runif and rbinom, specifying 2000 trials, selecting 1 each trial, each with a prob of success = 1/9 or 0.1111.. that doesn't seem to work?

For the second one, I've written this function
benlaw <- function(d){
log10(1 + 1 / d)
}
digits <- 1:9
data <- data.frame(Outcome = benlaw(digits))

However, neither the bar or histogram plots match the one in the textbook

[bonnfire]

Also, would you mind clarifying how you would like me to include xpbinom here? I don't quite understand how to apply this function to this particular question with the cumulative number of heads in a series of tosses?
I think I'm mainly confused about the first argument in the function.

xpbinom(c(1, 10),size=1,prob=0.5)

[nicholasjhorton]

We may need to defer these questions until my return unless Margaret and Shukry can assist (I presume that you've spoken with them).

But see:

xpbinom(c(1,10), size=15, prob=0.5)

for a better example

[nicholasjhorton]

For Benford, I'd suggest just displaying the probabilities from https://en.wikipedia.org/wiki/Benford%27s_law Nick

On Jul 18, 2018, at 11:23 AM, Bonnie Lin ***@***.***> wrote: Goal is to recreate these two plots: For the first one, I've tried runif and rbinom, specifying 2000 trials, selecting 1 each trial, each with a prob of success = 1/9 or 0.1111.. that doesn't seem to work? For the second one, I've written this function benlaw <- function(d){ log10(1 + 1 / d) } digits <- 1:9 data <- data.frame(Outcome = benlaw(digits)) However, neither the bar or histogram plots match the one in the textbook — You are receiving this because you commented. Reply to this email directly, view it on GitHub, or mute the thread.

Nicholas Horton Professor of Statistics Department of Mathematics and Statistics, Amherst College PO Box 5000, AC #2239 Amherst, MA 01002-5000

[bonnfire]

So I'll just include an image from the WIki page then?

[nicholasjhorton]

No. Instead why don't you use the probabilities from the table. Nick

On Jul 18, 2018, at 12:33 PM, Bonnie Lin ***@***.***> wrote: So I'll just include an image from the WIki page then? — You are receiving this because you commented. Reply to this email directly, view it on GitHub, or mute the thread.

Nicholas Horton Professor of Statistics Department of Mathematics and Statistics, Amherst College PO Box 5000, AC #2239 Amherst, MA 01002-5000

[nicholasjhorton]

Any updates on this front?