[R] p-Value

2007-08-27 Thread amna khan
Hi Sir When we use Kendall Package to obtain Kendall's Tau statistic. Then we also get two-sided p value. What does two-sided p-value mean? The word two-sided is confusing to understand. Kindly provide help in this regard. -- AMINA SHAHZADI Department of Statistics GC University Lahore,

Re: [R] p-Value

2007-08-27 Thread Daniel Lakeland
On Mon, Aug 27, 2007 at 11:49:19AM +0500, amna khan wrote: Hi Sir When we use Kendall Package to obtain Kendall's Tau statistic. Then we also get two-sided p value. What does two-sided p-value mean? The word two-sided is confusing to understand. Two-sided is sometimes also called

[R] p value statistic for lm object

2007-08-13 Thread Arjun Bhandari
Hi, I conduct a univariate regression with lm function. I would like to get the p value for the regression. Is there a method that would enable me to extract the p value into a variable. Thanks. Arjun Bhandari

Re: [R] p value statistic for lm object

2007-08-13 Thread Dieter Menne
Arjun Bhandari arb.em at adia.ae writes: I conduct a univariate regression with lm function. I would like to get the p value for the regression. Is there a method that would enable me to extract the p value into a variable. # From lm docs (Sorry, Annette, for misusing your data) ctl -

Re: [R] p-value from survreg

2007-07-12 Thread Terry Therneau
The question was how to get the p-value from the fit below, as an S object sr-survreg(s~groups, dist=gaussian) Coefficients: (Intercept) groups -0.02138485 0.03868351 Scale= 0.01789372 Loglik(model)= 31.1 Loglik(intercept only)= 25.4 Chisq= 11.39 on 1 degrees of freedom, p=

Re: [R] p-value from survreg

2007-07-12 Thread hadley wickham
On 7/12/07, Terry Therneau [EMAIL PROTECTED] wrote: The question was how to get the p-value from the fit below, as an S object sr-survreg(s~groups, dist=gaussian) Coefficients: (Intercept) groups -0.02138485 0.03868351 Scale= 0.01789372 Loglik(model)= 31.1 Loglik(intercept only)=

[R] p-value from survreg, library(survival)

2007-07-11 Thread vladosr
dear r experts: I would appreciate advice on how to get the p-value from the object 'sr' created with the function survreg() as given below. vlad sr-survreg(s~groups, dist=gaussian) Coefficients: (Intercept) groups -0.02138485 0.03868351 Scale= 0.01789372 Loglik(model)= 31.1

[R] p-value from survreg(), library(survival)

2007-07-11 Thread Vlado Sremac
dear r experts: It seems my message got spam filtered, another try: i would appreciate advice on how to get the p-value from the object 'sr' created with the function survreg() as given below. vlad sr-survreg(s~groups, dist=gaussian) Coefficients: (Intercept) groups -0.02138485 0.03868351

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread hadley wickham
str(survreg(s~groups, dist=gaussian)) is probably a good place to start. Hadley On 7/11/07, Vlado Sremac [EMAIL PROTECTED] wrote: dear r experts: It seems my message got spam filtered, another try: i would appreciate advice on how to get the p-value from the object 'sr' created with the

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread Henrique Dallazuanna
Try also: pchisq(summary(sr)$chi, degrees_freedom, lower=FALSE) *You need know your degrees of freedom -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 11/07/07, Vlado Sremac [EMAIL PROTECTED] wrote: dear r experts: It seems my message got spam filtered, another

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread Marc Schwartz
Actually, in this case, looking at the code for: survival:::print.survreg would be better, as the p value is calculate there, rather than being part of the survreg object. As with many R functions, the p value is calculated in the print method for the object. In this case, it is a pretty

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread hadley wickham
On 7/11/07, Marc Schwartz [EMAIL PROTECTED] wrote: Actually, in this case, looking at the code for: survival:::print.survreg would be better, as the p value is calculate there, rather than being part of the survreg object. As with many R functions, the p value is calculated in the print

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread Marc Schwartz
On Wed, 2007-07-11 at 16:41 +0200, hadley wickham wrote: On 7/11/07, Marc Schwartz [EMAIL PROTECTED] wrote: Actually, in this case, looking at the code for: survival:::print.survreg would be better, as the p value is calculate there, rather than being part of the survreg object. As

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-12 Thread Carmen Meier
Benilton Carvalho schrieb: Well, AFAIK, the definition of a p-value is the probability of observing something at least as extreme as the observed data. If you observed z, and Z follows a std-normal p-value = P( Z -abs(z) ) + P( Z abs(z) ) = 2*P ( Z abs(z) ) = 2*pnorm(z,

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Carmen Meier
I got an answer for the other question (thank you) But there is another question (I am afraid this is a basic question ...) In this tread there is a hint hwo to calculate the p-vlue of an GEE: _http://finzi.psych.upenn.edu/R/Rhelp02a/archive/74150.html_ Then, get the P values using a normal

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Benilton Carvalho
the recommendation was to use lower.tail=FALSE. b On Jun 11, 2007, at 11:21 AM, Carmen Meier wrote: I got an answer for the other question (thank you) But there is another question (I am afraid this is a basic question ...) In this tread there is a hint hwo to calculate the p-vlue of

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Carmen Meier
Benilton Carvalho schrieb: the recommendation was to use lower.tail=FALSE. b O but then the results are significant and this does not match the observation. The results are matching the observations if the formula is

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Robert A LaBudde
At 11:21 AM 6/11/2007, Carmen wrote: snip In this tread there is a hint hwo to calculate the p-vlue of an GEE: _http://finzi.psych.upenn.edu/R/Rhelp02a/archive/74150.html Then, get the P values using a normal approximation for the distribution of z: / 2 *

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Benilton Carvalho
Well, AFAIK, the definition of a p-value is the probability of observing something at least as extreme as the observed data. If you observed z, and Z follows a std-normal p-value = P( Z -abs(z) ) + P( Z abs(z) ) = 2*P ( Z abs(z) ) = 2*pnorm(z, lower.tail=FALSE) try z=0 (you should

[R] p-value from GEE

2007-06-06 Thread carmei3
Hi to all, I found in the R-help archive how to calculate the p-value for a gee result: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/74150.html but there are two questions (I am afraid they are basic questions ...) 1. why is the result multiplicated with 2 2. how could I decide between

[R] p-value gee

2007-06-04 Thread Web de
Hi to all, I found in the R-help archive how to calculate the p-value for a gee result: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/74150.html but there are two questions (I am afraid they are basic questions ...) 1. why is the result multiplicated with 2 2. how could I decide between

[R] p-value and significance code

2007-04-10 Thread Paul Suckling
Dear R list members. Is there anyone out there that might be able to explain the difference between p-values and significance codes for terms in a linear regression model (as given by summary.lm for example) for me please? Do they both give the probability that the observed variation in the

[R] p value for coefficients in multinomial model

2007-04-05 Thread Xingwang Ye
Dear all, 1)how can I easily get p value for the coefficients of factors in a multinomial model? 2)why the p values for type III test with Anova are not identical to that from SAS? for example: A,B and C are categorical variables,but the proportions of each level in each categorical

Re: [R] p value for coefficients in multinomial model

2007-04-05 Thread John Fox
://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Xingwang Ye Sent: Thursday, April 05, 2007 9:01 AM To: r-help Subject: [R] p value for coefficients in multinomial model Dear all, 1)how can I

[R] p-value of Logrank-Test

2006-01-03 Thread Verena Hoffmann
Hello! I want to compare two Kaplan-Meier-Curves by using the Logrank-Test: logrank(Surv(time[b], status[b]) ~ group[b]) This way I only get the value of the test-statistic, but not the p-value. Does anybody know how I can get the p-value? Thanks in advance! Verena Hoffmann

Re: [R] p-value of Logrank-Test

2006-01-03 Thread Thomas Lumley
On Tue, 3 Jan 2006, Verena Hoffmann wrote: Hello! I want to compare two Kaplan-Meier-Curves by using the Logrank-Test: logrank(Surv(time[b], status[b]) ~ group[b]) This way I only get the value of the test-statistic, but not the p-value. Does anybody know how I can get the p-value? You

[R] p-value calculation

2005-10-18 Thread richard mendes
hello everybody i'm very new at using R so probably this is a very stupid question. I have a problem calculating a p-value. When i do this with excel i can use the method CHIDIST for 1.2654 with 1 freedom degree i get the answer 0.261 i just want to do the same thing in R but i can't find a

Re: [R] p-value calculation

2005-10-18 Thread Marc Schwartz (via MN)
On Tue, 2005-10-18 at 17:13 +0200, richard mendes wrote: hello everybody i'm very new at using R so probably this is a very stupid question. I have a problem calculating a p-value. When i do this with excel i can use the method CHIDIST for 1.2654 with 1 freedom degree i get the answer 0.261

Re: [R] p-value calculation

2005-10-18 Thread Robert Kinley
pchisq(1.2654,1) [1] 0.7393686 looks like excel gave you the other tail bob kinley richard mendes [EMAIL PROTECTED] Sent by: [EMAIL PROTECTED] 18/10/2005 16:13 To R-help@stat.math.ethz.ch cc Subject [R] p-value calculation hello everybody i'm very new

Re: [R] p-value calculation

2005-10-18 Thread ronggui
pchisq(1.2654,df=1,low=F) [1] 0.2606314 === 2005-10-18 23:13:20 您在来信中写道:=== hello everybody i'm very new at using R so probably this is a very stupid question. I have a problem calculating a p-value. When i do this with excel i can use the method CHIDIST for 1.2654 with 1 freedom

Re: [R] p-value calculation

2005-10-18 Thread Peter Dalgaard
richard mendes [EMAIL PROTECTED] writes: hello everybody i'm very new at using R so probably this is a very stupid question. I have a problem calculating a p-value. When i do this with excel i can use the method CHIDIST for 1.2654 with 1 freedom degree i get the answer 0.261 i just want

Re: [R] p-value calculation

2005-10-18 Thread Nordlund, Dan
Look at ?pchisq Daniel Nordlund Bothell, WA -Original Message- From: [EMAIL PROTECTED] [mailto:r-help- [EMAIL PROTECTED] On Behalf Of richard mendes Sent: Tuesday, October 18, 2005 8:13 AM To: R-help@stat.math.ethz.ch Subject: [R] p-value calculation hello everybody i'm very

Re: [R] p-value calculation

2005-10-18 Thread Don MacQueen
1-pchisq(1.2654,1) [1] 0.2606314 -Don At 5:13 PM +0200 10/18/05, richard mendes wrote: hello everybody i'm very new at using R so probably this is a very stupid question. I have a problem calculating a p-value. When i do this with excel i can use the method CHIDIST for 1.2654 with 1 freedom

[R] p-value for non-linear variable in overdispersed glm()

2005-09-30 Thread tsnall
Dear all, I am fitting an nonlinear glm() using optim() by first minimising glm(resp~ var1 + var2, family=binomial, data=data)$deviance where var1= exp(-a1*dist1), and var2= exp(-a2*dist2), where a1 and a2 are parameters and dist1 and dist2 are independent variables. Next, I calculate the

[R] p-value in lrm

2005-08-30 Thread donghu
Hi, How can I get the p-value, S.E., and Z from the result of lrm fit? Thank you. Donglei __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] p-value in lrm

2005-08-30 Thread Frank E Harrell Jr
[EMAIL PROTECTED] wrote: Hi, How can I get the p-value, S.E., and Z from the result of lrm fit? Thank you. Donglei That's automatic. The print method for lrm in the Design package does that. Better though is the anova method which gives you the right pooled tests when you have

[R] p-value 1 in fisher.test()

2005-06-03 Thread Yi-Xiong Zhou
The following contingency table generates p-value 1 from fisher.test() ff = c(0,10,250,5000); dim(ff) = c(2,2); fhisher.test(ff)$p.value Sean [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

Re: [R] p-value 1 in fisher.test()

2005-06-03 Thread Peter Dalgaard
Yi-Xiong Zhou [EMAIL PROTECTED] writes: The following contingency table generates p-value 1 from fisher.test() ff = c(0,10,250,5000); dim(ff) = c(2,2); fhisher.test(ff)$p.value On this system (PIII, Fedora Core 3): ff = c(0,10,250,5000); dim(ff) = c(2,2); fisher.test(ff)$p.value -1

Re: [R] p-value 1 in fisher.test()

2005-06-03 Thread Don MacQueen
At 10:01 PM +0200 6/3/05, Peter Dalgaard wrote: Yi-Xiong Zhou [EMAIL PROTECTED] writes: The following contingency table generates p-value 1 from fisher.test() ff = c(0,10,250,5000); dim(ff) = c(2,2); fhisher.test(ff)$p.value On this system (PIII, Fedora Core 3): ff =

FW: [R] p-value 1 in fisher.test()

2005-06-03 Thread Yi-Xiong Zhou
The problem is from the most recent version. Here is the platform: R 2.1.0 2005-04-18, on winXP SP2, Sean -Original Message- From: Douglas Grove [mailto:[EMAIL PROTECTED] Sent: Friday, June 03, 2005 12:49 PM To: Yi-Xiong Zhou Subject: Re: [R] p-value 1 in fisher.test() You need

Re: [R] p-value 1 in fisher.test()

2005-06-03 Thread Ted Harding
On 03-Jun-05 Peter Dalgaard wrote: Yi-Xiong Zhou [EMAIL PROTECTED] writes: The following contingency table generates p-value 1 from fisher.test() ff = c(0,10,250,5000); dim(ff) = c(2,2); fhisher.test(ff)$p.value On this system (PIII, Fedora Core 3): ff = c(0,10,250,5000);

Re: [R] p-value 1 in fisher.test()

2005-06-03 Thread Ted Harding
On 03-Jun-05 Ted Harding wrote: And on mine (A: PII, Red Had 9, R-1.8.0): ff - c(0,10,250,5000); dim(ff) - c(2,2); 1-fisher.test(ff)$p.value [1] 1.268219e-11 (B: PIII, SuSE 7.2, R-2.1.0beta): ff - c(0,10,250,5000); dim(ff) - c(2,2); 1-fisher.test(ff)$p.value [1]

[R] p-value for nonparamentric multiple comparison testing

2004-11-02 Thread Antonio Olinto
Dear R users, I wrote a function to perform a nonparametric multiple comparison test. The function below solves the example 11.10 (pg. 228) from Zar’s Biostatistical Analysis (3rd ed.). Nevertheless I couldn’t find a function to get the p-value associated to Q, so I still have to consult Table

[R] P-value from the joint cumulative distribution of an n-dimensional order statistic

2004-07-30 Thread Jeroen Van Goey
Hello, I want to compute the P-value from the joint cumulative distribution of an n-dimensional order statistic in R, using the formula found on http://cmgm.stanford.edu/%7Ekimlab/multiplespecies/Supplement/methods_network.html My data consists of three different techniques (G2D, POCUS and

Re: [R] p-value from chisq.test working strangely on 1.8.1

2003-12-09 Thread Torsten Hothorn
Hello everybody, I'm seeing some strange behavior on R 1.8.1 on Intel/Linux compiled with gcc 3.2.2. The p-value calculated from the chisq.test function is incorrect for some input values: chisq.test(matrix(c(0, 1, 1, 12555), 2, 2), simulate.p.value=TRUE) Pearson's

Re: [R] p-value from chisq.test working strangely on 1.8.1

2003-12-09 Thread Marc Schwartz
On Tue, 2003-12-09 at 02:23, Torsten Hothorn wrote: Hello everybody, I'm seeing some strange behavior on R 1.8.1 on Intel/Linux compiled with gcc 3.2.2. The p-value calculated from the chisq.test function is incorrect for some input values: chisq.test(matrix(c(0, 1, 1, 12555),

Re: [R] p-value from chisq.test working strangely on 1.8.1

2003-12-09 Thread John Dougherty
It happens with suse 9.0 as well. ... chisq.test(matrix(c(0, 1, 1, 12556), 2, 2), simulate.p.value=TRUE) Pearson's Chi-squared test with simulated p-value (based on 2000

Re: [R] p-value from chisq.test working strangely on 1.8.1

2003-12-09 Thread Jeffrey Chang
On Dec 9, 2003, at 9:37 AM, Peter Dalgaard wrote: Marc Schwartz [EMAIL PROTECTED] writes: Confirmed on Fedora Core 1 with R Version 1.8.1 Patched (2003-12-07) compiled using gcc (GCC) 3.3.2 20031107 (Red Hat Linux 3.3.2-2). chisq.test(matrix(c(0, 1, 1, 12555), 2, 2), simulate.p.value=TRUE) ...

[R] p-value from chisq.test working strangely on 1.8.1

2003-12-08 Thread Jeffrey Chang
Hello everybody, I'm seeing some strange behavior on R 1.8.1 on Intel/Linux compiled with gcc 3.2.2. The p-value calculated from the chisq.test function is incorrect for some input values: chisq.test(matrix(c(0, 1, 1, 12555), 2, 2), simulate.p.value=TRUE) Pearson's Chi-squared

[R] p value in MANOVA

2003-11-20 Thread Tu Yu-Kang
Dear R users, Can anyone tell me how to get the p value out of the output from summary.manova? I tried all the methods I can think of, but failed. Many thanks Yu-Kang _ MSN Mobile MSN Hotmail http://msn.com.tw/msnmobile

Re: [R] p value in MANOVA

2003-11-20 Thread Prof Brian Ripley
Hint: print.summary.manova manages it. Try example(summary.manova) summary(fit, test=Wilks)$stats[,Pr(F)] On Thu, 20 Nov 2003, Tu Yu-Kang wrote: Can anyone tell me how to get the p value out of the output from summary.manova? I tried all the methods I can think of, but failed. --