Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Fox, John
Dear Eric, Have you looked at your data? -- for example: plot(log(Moons) ~ Volume, data = moon_data) text(log(Moons) ~ Volume, data = moon_data, labels=Name, adj=1, subset = Volume > 400) The negative-binomial model doesn't look reasonable, does it? After you eliminate Jupiter

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Jim Lemon
Hi Eric, When I run your code (using the MASS library) I find that rstudent(fit2) also returns NaN in the seventh position. Perhaps the problem is occurring there and not in the "influence" function. Jim On Wed, Apr 3, 2019 at 9:12 AM Eric Bridgeford wrote: > > I agree the influence documentatio

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Eric Bridgeford
I agree the influence documentation suggests NaNs may result; however, as these can be manually computed and are, indeed, finite/existing (ie, computing the held-out influence by manually training n models for n points to obtain n leave one out influence measures), I don't possibly see how the func

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Eric Bridgeford
How can I add attachments? The following two files were attached in the initial message On Tue, Apr 2, 2019 at 3:34 PM Bert Gunter wrote: > Nothing was attached. The r-help server strips most attachments. Include > your code inline. > > Also note that > > > 0/0 > [1] NaN > > so maybe something l

[R] Free financial data - equities, equity options and ETFs - for quantmod package (or other packages)

2019-04-02 Thread H
I am relatively new to analyzing financial data but have some experience with R. I understand that the data available from Yahoo Finance via its API is often questionable in quality and Google Finance is no longer available. Although Googling pointed me to some other sources such as Quandl etc.,

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Bert Gunter
Also, I suggest you read ?influence which may explain the source of your NaN's . Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Apr 2, 2019 at 1:29 PM Ber

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Bert Gunter
I told you already: **Include code inline ** See ?dput for how to include a text version of objects, such as data frames, inline. Otherwise, I believe .txt text files are not stripped if you insist on *attaching* data or code. Others may have better advice. Bert Gunter "The trouble with having

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Bert Gunter
Nothing was attached. The r-help server strips most attachments. Include your code inline. Also note that > 0/0 [1] NaN so maybe something like that occurs in the course of your calculations. But that's just a guess, so feel free to disregard. Bert Gunter "The trouble with having an open mind

[R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Eric Bridgeford
Hi R core team, I experienced the following issue with the attached data/code snippet, where the studentized residual for a single observation appears to be NaN given finite predictors/responses, which appears to be driven by the glm.influence method in the stats package. I am curious to whether t

Re: [R] Pairwise testing with pairwise.prop.test

2019-04-02 Thread Jenny Liu
Thank you Stephen and Peter for your code and your answers! Stephen, your NaN explanation makes sense - thank you for putting it so clearly. Cheers, Jenny On Tue, Apr 2, 2019 11:10 AM, peter dalgaard pda...@gmail.com wrote: >> My questions:1) Why is there a "-" instead of a numerical result

Re: [R] Pairwise testing with pairwise.prop.test

2019-04-02 Thread peter dalgaard
>> My questions:1) Why is there a "-" instead of a numerical result for pairs >> 1-2, >> 1-16, and 2-16? > > When the difference between the pair is zero, the p.value is NaN (not a > number). > Not quite: When both groups have 0 successes (or both 0 failures), the test stat has a divide-by-ze