similar to: p values from GLM

... 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: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

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

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

How small does a p-value need to be to warrant attention, however? Witness Fisher?s comment that: ?. . . we may, if we prefer it, draw the line at one in fifty (the 2 per cent point), or one in a hundred (the 1 per cent point). Personally, the writer prefers to set a low standard of significance at the 5 per cent point, and ignore entirely all results which fail to reach this level. A scientific

> 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

R 2.7.2 Windows XP I am using apply to compute a series of Student's t-test from two matrices, sample1 and sample2. boo<-apply(sample1,1,t.test,sample2) I want to pick of the p-values from the tests, but can't seem to get it to work. I have tried several methods to get the values including: boo<-apply(sample1,1,t.test$t.test,sample2) boo<-apply(sample1,1,t.test,sample2)$t.test

Widows XP R 2.3.1 I have been trying to make a data structure that will contain both the coefficients from a linear regression along with column and row titles AND the call, i.e. myreg<-lm(y~x+y+z) whatIwant<-cbind(c(summary(myreg)$call,"",""),summary(myreg)$coefficients) Neither the statement above, nor any one of twenty variations I have tried work. I would appreciate

(1) I know that \n when used in cat, e.g. cat("\n") produces a line feed (i.e. skips to the next line). Is there any escape sequence that will go to the top of the next page? (2) I know that control L will clear the console. Is there an equivalent function or other means that can be used in R code to clear the console? Thanks, John John David Sorkin M.D., Ph.D. Chief, Biostatistics

I am looking for a package that will allow me to choose R (a real number) that minimizes the correlation of y and x^R, i.e. find R such that corr(y,x^R) is minimized. Any suggestions for packages I might look at would be helpful. Thanks, John John David Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics University of Maryland School of Medicine Division of Gerontology Baltimore VA

??I would like to apply a function, fract, to the columns of a dataframe. I tried the following apply(data5NonEventEpochs,2,fract) but, no surprise it did not work as apply works on matrices not data frames. How can I apply a fuction to the columns of a data frame? (I can't covert data5NonEventsEpochs to a matrix as it contains character data). Thank you, John John David Sorkin M.D., Ph.D.

An apparatus exists whereby a collection of balls is displaced to the top of a stack by suction. A top level (Level 1) each ball is shifted 1 unit to the left or 1 unit to the right at random with equal probability. The ball then drops down to level 2. At Level 2, each ball is again shifted 1 unit to the left or 1 unit to the right at random. The process continues for 15 levels and the balls are

R 2.8.1 Windows XP I am trying to analyze repeated measures data (the data are listed at the end of this Email message) and I need help to make sure that I have properly specified my model, and would like to know why lmer does not return a p value for Group, my fixed effect. My subjects are divided into two groups (variable GROUP), individual subjects are indicated by the variable SS, Value is

I'm looking for an icon to represent an R package. Perhaps something like http://cdn2.iconfinder.com/data/icons/DarkGlass_Reworked/128x128/apps/package.png but with the R logo rather than KDE. -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Street Web:

I am having trouble selecting rows of a dataframe that will be included in a regression. I am trying to select those rows for which the variable Meno equals PRE. I have used the code below: difffitPre<-lm(data[,"diff"]~data[,"Age"]+data[,"Race"],data=data[data[,"Meno"]=="PRE",]) summary(difffitPre) The output from the summary indicates that

I would like to select rows if a row contains any one of several values. I can do the selection as follows: result[,"Subject"]=="JEFF" | result[,"Subject"]=="BG" But this is very unwieldily if one wishes to select many, many rows as one has to continuously repeat the source: result[,"Subject"]=="JEFF" |

I am trying to learn more about how to write functions. I would like to pass a data frame (or matrix) and depending on the parameters passed to the function work with a given column of the data frame or matrix. My function, learnfn is given below as are two calls to the function. The first call is an attempt to print the x column from the data frame, the second call is an attempt to print the y

windows Vista R 2.10.1 Is it possible to get p values from gee? Summary(geemodel) does not appear to produce p values.: > fit4<- gee(y~time, id=Subject, data=data.frame(data)) Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27 running glm to get initial regression estimate (Intercept) time 1.1215614 0.8504413 > summary(fit4) GEE: GENERALIZED LINEAR MODELS FOR

I would like to plot 3-dimensional data on a two-dimensional scatter-plot. Is there a way I can automatically modify the plot symbol (e.g. changing size or color) to indicate the value of a third variable? E.g. How can I plot weight vs. age and indicate the value of muscle mass for each value weight-age pair by making the plot point proportional to the subject's muscle mass? Thanks, John

R 2.8.1 windows XP I am trying to learn how to use the boot function to perform a bootstrap of a regression. I have written a short trial program, shown below. Clearly I have done something wrong as the output of each of the 100 bootstrap values for the regression are exactly the same - there does not appear to be any bootstrap respampling!. What have I done wrong? # Define function to be run.