similar to: p values from GLM

Bert Gunter wrote on 01.04.2016 23:46: > ... of course, whether one **should** get them is questionable... > > http://www.nature.com/news/statisticians-issue-warning-over-misuse-of-p-values-1.19503#/ref-link-1 > This paper repeats the common place statement that a small p-value does not necessarily indicate an important finding. Agreed, but maybe I overlooked examples of important

... of course, whether one **should** get them is questionable... http://www.nature.com/news/statisticians-issue-warning-over-misuse-of-p-values-1.19503#/ref-link-1 Cheers, Bert 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 Fri, Apr

On 01/04/2016 6:14 PM, John Sorkin wrote: > How can I get the p values from a glm ? I want to get the p values so I can add them to a custom report > > > fitwean<- glm(data[,"JWean"]~data[,"Group"],data=data,family=binomial(link ="logit")) > summary(fitwean) # This lists the coefficeints, SEs, z and p values, but I can't

On 01/04/2016 6:46 PM, Bert Gunter wrote: > ... of course, whether one **should** get them is questionable... They're just statistics. How could it hurt to look at them? Duncan Murdoch > > http://www.nature.com/news/statisticians-issue-warning-over-misuse-of-p-values-1.19503#/ref-link-1 > > > Cheers, > Bert > > > > Bert Gunter > > "The trouble

How can I get the p values from a glm ? I want to get the p values so I can add them to a custom report fitwean<- glm(data[,"JWean"]~data[,"Group"],data=data,family=binomial(link ="logit")) summary(fitwean) # This lists the coefficeints, SEs, z and p values, but I can't isolate the pvalues. names(summary(fitwean)) # I see the coefficients,

> On Apr 1, 2016, at 5:01 PM, Duncan Murdoch <murdoch.duncan at gmail.com> wrote: > > On 01/04/2016 6:46 PM, Bert Gunter wrote: >> ... of course, whether one **should** get them is questionable... > > They're just statistics. How could it hurt to look at them? Like Rolf, I thought that this utterance on April 1 deserved fortune enshrinement. It reminded me of one

Dear Heinz, ---------------------------------------------- > On Dec 17, 2018, at 10:19 AM, Heinz Tuechler <tuechler at gmx.at> wrote: > > Dear All, > > do you think that use of a data argument is best practice in the example below? No, but it is *normally* or *usually* the best option, in my opinion. Best, John > > regards, > > Heinz > > ###

[I tried to send this message privately, but the return address bounced.] I think this has been fixed in R-patched, but I doubt if the fix has been tested in Win98. Could you please download a copy from <http://cran.r-project.org/bin/windows/base/rpatched.html> and confirm that it has been fixed? Duncan Murdoch On Sat, 27 Nov 2004 23:31:23 +0100, Heinz Tuechler <tuechler at gmx.at>

Following p.206 of "Statistical Models in S", I wish to change the code for summary.glm() so that it estimates the dispersion for binomial & poisson models when the parameter dispersion is set to zero. The following changes [insertion of ||dispersion==0 at one point; and !is.null(dispersion) at another] will do the trick: "summary.glm" <- function(object, dispersion =

I've changed the heading because this really is another thread. I think it inevitable that there will, in the course of time, be other lists that are devoted, in some shape or form, to the concerns of practitioners (at all levels) who are using R. One development I'd not like to see is fracture along application area lines, allowing those who are comfortable in coteries whose

OK, I see now that I was supposed to twig that the reference was to putting the ?.Rnw' files back into the vignettes directory from the inst/doc directory where they?d been placed in the course of creating the tar.gz file. I am still trying to work out what I need to put into ?.Rinstignore? so that ?.install_extras? is not installed. John Maindonald email: john.maindonald at

Sorry. This, and the description in the ?Writing R Extensions? manual, leaves me completely mystified. Is it that I have to remove the PDFs that are created when I run ?R CMD build?, and somehow ensure that they are rebuilt when the package is installed? Do I need a Makefile? John Maindonald email: john.maindonald at anu.edu.au<mailto:john.maindonald at anu.edu.au> phone :

Confirmation from John Maindonald. -- Brian D. Ripley, ripley@stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272860 (secr) Oxford OX1 3TG, UK Fax: +44 1865 272595 ---------- Forwarded message ---------- Date: Fri, 3

Dear All, is there a simple way to print a data.frame without its row.names? example: datum <- as.Date(c("2004-01-01", "2004-01-06", "2004-04-12")) content <- c('Neujahr', 'Hl 3 K.', 'Ostern') df1 <- data.frame(datum, content) print(df1) datum content 1 2004-01-01 Neujahr 2 2004-01-06 Hl 3 K. 3 2004-04-12 Ostern Can I get

Hi, I am trying to estimate the sample size needed for the comparison of two groups on a certain measurement, given some previous data at hand. I find that the data collected does not follow a normal distribution, so I would like to use a non-parametric option for sample size calculation. I found the pwr package but I don't think it has this option and on the internet found that

From a quick look at the paper in the SAS proceedings, the simulations seem limited to nested designs. The major problems are with repeated measures designs where the error structure is not compound symmetric, which lme4 does not at present handle (unless I have missed something). Such imbalance as was investigated was not a serious issue, at least for the Kenward and Roger degree of freedom

You've omitted a comma. races2000 is a data frame, which for purposes of extracting rows behaves like a 2-dimenional object. The following works fine: hills2000 <- races2000[races2000$type == 'hill', ] Additionally, you might like to ponder > type <- races2000[names(races2000)=="type"] > type[1:4] Error in "[.data.frame"(type, 1:4) :

In the following titletxt <- substitute(paste(tx, tilde(y) == sqrt(y)), list(tx="Replace y by ")) I'd like to be able to replace e.g. sqrt(y) by a symbol or symbolic expression that has the effect of inserting sqrt(y) into the place that it occupies. This would allow me to put an arbitrary expression into that position. [One can use text() or mtext() etc. to place titletxt

I have posted, at http://wwwmaths.anu.edu.au/~john/r/newsets/ image (.rda) files, and first stabs at .Rd files for various data on deaths in London from 1629 to 1939. (There are of course gaps.) The sources (Guy 1882 & Stocks 1942) are documented in the .Rd files: (1) poxetc: measles, smallpox & total deaths: 1629-1881 [I have deliberately left several inconsistencies that were in

<soap-box> While I agree with John about abuse of methodology, I can't subscribe to the "failure to read the manual" proposition. I don't mean to lecture, but it would seem to me that people who, like me, will probably gain far more from this list than we will be able to contribute (at least in the near term) owe it to those who will likely contribute far more than they