Re: [R] p-value in changepoint package

2023-02-15 Thread Ivan Krylov
В Wed, 15 Feb 2023 14:20:16 + Nick Wray пишет: > how to get the p values for changepoints in the "changepoint" > package? The documentation of this package is not always precise in describing what the code actually does. It's possible to find out by reading the source code that with

[R] p-value in changepoint package

2023-02-15 Thread Nick Wray
Hello Does anyone know how to get the p values for changepoints in the "changepoint" package? There is guidance in the CRAN vignette https://cran.r-project.org/web/packages/changepoint/changepoint.pdf but I can't make it work to get the p value out... Thanks Nick Wray [[alternative HTML

Re: [R] P value of VAR model from Portmanteau Test

2017-05-23 Thread Rui Barradas
Hello, It seems to be right. You have Chi-squared = 23.724, df = 8, p-value = 0.002549. So try the R function ?pchisq: pchisq(23.724, df = 8, lower = FALSE) [1] 0.002549054 Hope this helps, Rui Barradas Em 23-05-2017 13:08, Dhivya Narayanasamy escreveu: Hi, I am working with bivariate

[R] P value of VAR model from Portmanteau Test

2017-05-23 Thread Dhivya Narayanasamy
Hi, I am working with bivariate time series data. I used VAR model to fit and forecast. But the "*p*" value from seria.test (Portmanteau Test) gives values *p<< 0.05*. Is that okay? > var1 = VAR(datax.ts, p= 8) > serial.test(var1, lags.pt=10, type = "PT.asymptotic") Portmanteau Test

[R] p-value=0 running log-rank test

2017-05-18 Thread Anne Karin da Mota Borges
Dear all, I have a question concerning the p-value. When running log-rank test I get a p-value = 0. What is it mean? Can this be true? Why aren´t there decimal points? Is there a way to find out the exact p-value? Here is the output: > survdiff(Surv(tempo2,status)~tphist, data=base,rho=0) Call:

[R] p-value or the t-test value

2016-04-08 Thread Alaa Sindi
Hello, How can I print out the p-value or the t-test value for coefficients estimated from a likelihood function. Thanks __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] P value from jaccard's index matrix

2015-07-24 Thread sreenath
My table having 40 raw and 4 columns, in that 4 columns first column belongs to one group and the remaining constitute the other group. using following commands for calculating jaccard's index x - read.csv(file name,header=T, sep= ) jac - vegdist(x,method=jaccard) from this out file(jac) how can

Re: [R] P-value from Matching

2015-02-16 Thread Rob Wood
Hi Peter, That was my first port of call before I posted this thread. Unfortunately, it does not seem to explicitly state which test is used or how the p-value is calculated. Thanks, Rob. -- View this message in context:

Re: [R] P-value from Matching

2015-02-16 Thread peter dalgaard
On 16 Feb 2015, at 00:31 , Rob Wood robert.w...@adelphigroup.com wrote: Hi all, When using the match command from the matching package, the output reports the treatment effect, standard error, t-statistic and a p-value. Which test is used to generate this p-value, or how us it generated?

[R] P-value from Matching

2015-02-15 Thread Rob Wood
Hi all, When using the match command from the matching package, the output reports the treatment effect, standard error, t-statistic and a p-value. Which test is used to generate this p-value, or how us it generated? Thanks, Rob. -- View this message in context:

Re: [R] p value of trends for odds ratios (or hazard ratios)

2014-01-06 Thread zhu yao
Thanks Michael. I have read the original article and I guess whether the following formula could resolve my question? *for adjusted HR of each BMI category,* fit-coxph(Surv~factor(BMI-category)+..+covariates) *for trend, * fit-coxph(Surv~as.numeric(BMI-category)+..+covariates) *# p trend may be

Re: [R] p value of trends for odds ratios (or hazard ratios)

2014-01-05 Thread zhu yao
Thanks peter. *Yao Zhu* *Department of UrologyFudan University Shanghai Cancer CenterShanghai, China* 2014/1/5 peter dalgaard pda...@gmail.com On 04 Jan 2014, at 13:56 , zhu yao mailzhu...@gmail.com wrote: Thanks for the suggestion. The results is presented in following table. The

Re: [R] p value of trends for odds ratios (or hazard ratios)

2014-01-05 Thread Michael Friendly
On 1/4/2014 12:39 AM, zhu yao wrote: Dear Sir Many papers calculated the p value of trends for odds ratios of ordered category variables. I have found the tabodds command in Stata. But how to do it in R? Thanks *Yao Zhu* If what you are looking at is a 2 x 2 x k table, where you want the odds

Re: [R] p value of trends for odds ratios (or hazard ratios)

2014-01-04 Thread Uwe Ligges
On 04.01.2014 06:39, zhu yao wrote: Dear Sir Many papers calculated the p value of trends for odds ratios of ordered category variables. I have found the tabodds command in Stata. But how to do it in R? Depends on the method you want to use ... and most ladies and gents on this list won't

Re: [R] p value of trends for odds ratios (or hazard ratios)

2014-01-04 Thread peter dalgaard
On 04 Jan 2014, at 12:53 , Uwe Ligges lig...@statistik.tu-dortmund.de wrote: On 04.01.2014 06:39, zhu yao wrote: Dear Sir Many papers calculated the p value of trends for odds ratios of ordered category variables. I have found the tabodds command in Stata. But how to do it in R?

Re: [R] p value of trends for odds ratios (or hazard ratios)

2014-01-04 Thread peter dalgaard
On 04 Jan 2014, at 13:56 , zhu yao mailzhu...@gmail.com wrote: Thanks for the suggestion. The results is presented in following table. The authors calculated p value for linearity and trend and they stated in the methods: Linear and nonlinear trends of BMI associated with each mortality

[R] p value of trends for odds ratios (or hazard ratios)

2014-01-03 Thread zhu yao
Dear Sir Many papers calculated the p value of trends for odds ratios of ordered category variables. I have found the tabodds command in Stata. But how to do it in R? Thanks *Yao Zhu* *Department of UrologyFudan University Shanghai Cancer CenterShanghai, China* [[alternative HTML

Re: [R] p value for mu: anova()

2013-12-06 Thread Ben Bolker
Rosario Garcia Gil M.Rosario.Garcia at slu.se writes: Hello I have run an anova analysis for the fallowing model: H_obs=mu+REGION+MANAGEMENT + e When I run it in ASRelm I get the p-value for mu, and, of course also for the two dependent variables (REGION and MANAGEMENT) When I run

[R] p value for mu: anova()

2013-12-05 Thread Rosario Garcia Gil
Hello I have run an anova analysis for the fallowing model: H_obs=mu+REGION+MANAGEMENT + e When I run it in ASRelm I get the p-value for mu, and, of course also for the two dependent variables (REGION and MANAGEMENT) When I run it in R, I do not get the pvalue for mu. Can some one help me

[R] P value discrepancy plots

2012-10-04 Thread bilelsan
Hello R-users, I am trying to plot P value discrepancy plots as describe by Russell Davidson and James G. MacKinnon in Graphical Methods for Investigating the Size and Power of Hypothesis Tests (1998). I know ecdf () to obtain the empirical cumulative distribution frequencies and knots() to

Re: [R] p value from lmekin()

2012-09-06 Thread Terry Therneau
On 09/06/2012 05:00 AM, r-help-requ...@r-project.org wrote: Hi, R experts I am currently using lmekin() function in coxme package to fit a mixed effect model for family based genetic data. How can I extract the p value from a lmekin object? When I print the object in R console, I can

[R] p-value for the fitted parameters in linear models

2012-06-24 Thread Li SUN
Hi All, when using the linear model lm() to fit data of the form y = k * x + b, is it possible to know the p-value for the parameters k and b? i.e. can we find the result of the form (k, Δk; b, Δb)? Thanks in advance! Li Sun __ R-help@r-project.org

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread Uwe Ligges
On 24.06.2012 17:47, Li SUN wrote: Hi All, when using the linear model lm() to fit data of the form y = k * x + b, is it possible to know the p-value for the parameters k and b? i.e. can we find the result of the form (k, Δk; b, Δb)? If you explain what Δk means in a linear model, we may be

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread Li SUN
Sorry for the confusion. Let me state the question again. I missed something in my original statement. When using the linear model lm() to fit data of the form y = k * x + b, where k, b are the coefficients to be found, and x is the variable and has an error bar (uncertainty) Δx of the same

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread David Winsemius
On Jun 24, 2012, at 1:21 PM, Li SUN wrote: Sorry for the confusion. Let me state the question again. I missed something in my original statement. When using the linear model lm() to fit data of the form y = k * x + b, where k, b are the coefficients to be found, and x is the variable and

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread Prof Brian Ripley
On 24/06/2012 18:39, David Winsemius wrote: On Jun 24, 2012, at 1:21 PM, Li SUN wrote: Sorry for the confusion. Let me state the question again. I missed something in my original statement. When using the linear model lm() to fit data of the form y = k * x + b, where k, b are the

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread Li SUN
Thanks David and Brian. But what if x is exact while y has some uncertainty Δy, in the relation y = k * x + b? Now I need to fit some data like x = 1, 2, 3, 4, 5 y±Δy = 1.1±0.1, 2.0±0.2, 3.1±0.2, 4.1±0.1, 5.0±0.2 Is there any mechanism to pass x, y and

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread Uwe Ligges
On 24.06.2012 20:35, Li SUN wrote: Thanks David and Brian. But what if x is exact while y has some uncertainty Δy, in the relation y = k * x + b? Now I need to fit some data like x = 1, 2, 3, 4, 5 y±Δy = 1.1±0.1, 2.0±0.2, 3.1±0.2, 4.1±0.1, 5.0±0.2

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread Li SUN
2012/6/24 Uwe Ligges lig...@statistik.tu-dortmund.de: On 24.06.2012 20:35, Li SUN wrote: Thanks David and Brian. But what if x is exact while y has some uncertainty Δy, in the relation y = k * x + b? Now I need to fit some data like x       = 1,          2,          3,          4,      

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread Bert Gunter
Dear Li Sun: You appear to have a good deal of statistical confusion. R-help is not a statistical consulting service. You should consult a local statistician or, if that's not possible,post on a statistical help list like stats.stackexchange.com -- Bert On Sun, Jun 24, 2012 at 12:04 PM, Li SUN

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread S Ellison
But what if x is exact while y has some uncertainty Δy, in the relation y = k * x + b? Again, no: this is not a linear model. Assumption in a linear model is that the errors are identically distributed. Surely not; errors in a linear model do not need to be homoscedastic. lm handles

Re: [R] p-value for the fitted parameters in linear models

2012-06-24 Thread Li SUN
Dear All, Thanks for all your explanations and sorry for the confusion. I will need to consult some statistician for help. Sincerely, Li Sun 2012/6/24 S Ellison s.elli...@lgcgroup.com: But what if x is exact while y has some uncertainty Δy, in the relation y = k * x + b? Again, no: this

[R] p-value using survival ROC

2012-03-28 Thread Anamaria Crisan
It is possible to calculate the c-index for time dependent outcomes (such as disease) using the survivalROChttp://cran.r-project.org/web/packages/survivalROC/survivalROC.pdf package in R. My question is : is it possible to produce a p-value for the c-index that is calculated (at a specific point

[R] p-value and mlogit

2012-03-18 Thread Ville Iiskola
Hi How does mlogit count the p-values of the variables? Ville [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] p-value of the pooled Z score

2012-03-13 Thread cheba meier
Hello, I have to compute the pooled z-value and I would like to know which way is more appropriate b - c( -0.205,1.040,0.087) s - c(0.449,0.167,0.241) n - c(310, 342, 348) z - b/s Z - sum(z)/sqrt(length(n)) P - 2*(1-pnorm(abs(Z))) P w - sqrt(n) Zw - sum(w * z)/sqrt(sum(w^2)) Pw - 1 -

Re: [R] p-value of the pooled Z score

2012-03-13 Thread Thomas Lumley
On Wed, Mar 14, 2012 at 3:41 AM, cheba meier cheba.me...@googlemail.com wrote: Hello, I have to compute the pooled z-value and I would like to know which way is more appropriate b - c( -0.205,1.040,0.087) s - c(0.449,0.167,0.241) n - c(310, 342, 348) z - b/s Z - sum(z)/sqrt(length(n))

Re: [R] p-value from GLM

2012-03-05 Thread Dunia Scheid
Dear Peter, Many thanks for the reply. Sure, this is just an example and it makes zero sense! Regards, Dunia 2012/3/4 peter dalgaard pda...@gmail.com On Mar 4, 2012, at 12:21 , Dunia Scheid wrote: Dear all, I am fitting a GLM similar to library(MASS) anorex.1 -

[R] p-value from GLM

2012-03-04 Thread Dunia Scheid
Dear all, I am fitting a GLM similar to library(MASS) anorex.1 - glm(Treat~Postwt+Prewt,family = binomial, data = anorexia) I have found two ways of computing the p-value of the fitted model: pval1 - 1-pchisq(anorex.1$deviance,anorex.1$df.residual) pval2 - 1-pchisq(anorex.1$null.deviance -

Re: [R] p-value from GLM

2012-03-04 Thread peter dalgaard
On Mar 4, 2012, at 12:21 , Dunia Scheid wrote: Dear all, I am fitting a GLM similar to library(MASS) anorex.1 - glm(Treat~Postwt+Prewt,family = binomial, data = anorexia) I hope that's just for an example. The actual analysis makes zero sense to me... I have found two ways of

Re: [R] p-value for hazard ratio in Cox proportional hazards regression?

2011-12-16 Thread Joshua Wiley
did) seems reasonable to me at least to ask. Cheers, Josh Best, Thierry - Ursprüngliche Mail - Von: Joshua Wiley jwiley.ps...@gmail.com An: Thierry Julian Panje tpa...@stanford.edu CC: r-help@r-project.org Gesendet: Samstag, 10. Dezember 2011 13:53:00 Betreff: Re: [R] p

[R] p-value for hazard ratio in Cox proportional hazards regression?

2011-12-10 Thread Thierry Julian Panje
Hi, I'm new to R and using it for Cox survival analysis. Thanks to this great forum I learned how to compute the HR with its confidence interval. My question would be: Is there any way to get the p-value for a hazard ratio in addition to the confidence interval? Thanks, Thierry --

Re: [R] p-value for hazard ratio in Cox proportional hazards regression?

2011-12-10 Thread Joshua Wiley
Hi Thierry, Could you give us an example of what exactly you are doing (preferablly reproducible R code)? I may be misunderstanding you, but if you are fitting cox proportional hazard models using the coxph() function from the surival package, summary(yourmodel) should give the SE, p-value based

Re: [R] p-value for hazard ratio in Cox proportional hazards regression?

2011-12-10 Thread Joshua Wiley
) seems reasonable to me at least to ask. Cheers, Josh Best, Thierry - Ursprüngliche Mail - Von: Joshua Wiley jwiley.ps...@gmail.com An: Thierry Julian Panje tpa...@stanford.edu CC: r-help@r-project.org Gesendet: Samstag, 10. Dezember 2011 13:53:00 Betreff: Re: [R] p-value

Re: [R] p value in R - beginners question

2011-10-19 Thread Ken Hutchison
The p-value as I understand it, is the probability of seeing a phenomenon at least as extreme as the one which you observed assuming the null hypothesis is true. In other words, if you assume there is no difference (which is the NULL for you here I believe), and there is an 81% chance of seeing

Re: [R] p value in R - beginners question

2011-10-19 Thread niki
Thanks for your help. I was not familiar with this form of notating the p value. -- View this message in context: http://r.789695.n4.nabble.com/p-value-in-R-beginners-question-tp3915873p3918140.html Sent from the R help mailing list archive at Nabble.com.

[R] p value in R - beginners question

2011-10-18 Thread niki
Dear all, i have done some regression analyses but i do not understand the p value. These are the results t-value p value geno.1 -0.229 0.978 -0.234 8.15e-01 geno.50.647 1.146 0.565 5.73e-01 stress:geno.5-1.337 1.022 -1.307

Re: [R] p value in R - beginners question

2011-10-18 Thread B77S
This is just scientific notation, so 8.15e-01 is the same as: 8.15*10^-1 [1] 0.815 niki wrote: Dear all, i have done some regression analyses but i do not understand the p value. These are the results t-value p value geno.1

[R] p-value for non linear model

2011-09-15 Thread Tatiana Donnay
Hello, I want to understand how to tell if a model is significant. For example I have vectX1 and vectY1. I seek first what model is best suited for my vectors and then I want to know if my result is significant. I'am doing like this: model1 - lm(vectY1 ~ vectX1, data= d), model2 - nls(vectY1

Re: [R] p-value for non linear model

2011-09-15 Thread Bert Gunter
Tatiana: It sounds like you are in way over your head statistically. This is not a statistics tutorial site (though sometimes good folks do help witjh this). I suggest you try http://stats.stackexchange.com/ . Better yet, work with your local statistician. Cheers, Bert On Thu, Sep 15, 2011 at

Re: [R] p-value for y in non-parametric distribution, Y

2011-02-08 Thread Sascha Vieweg
On 11-02-07 19:21, Sascha Vieweg wrote: Hello, knowing that some index y, with y(341)=2, SE=3, is t-distributed, I (think I) can find an appropriate (left/one-sided) p-value with R: y - 2 R: yse - 3 R: (p - 1-pt(y/yse, df=341)) Now, some simulation resulted in the non-parametric

[R] p-value for y in non-parametric distribution, Y

2011-02-07 Thread Sascha Vieweg
Hello, knowing that some index y, with y(341)=2, SE=3, is t-distributed, I (think I) can find an appropriate (left/one-sided) p-value with R: y - 2 R: yse - 3 R: (p - 1-pt(y/yse, df=341)) Now, some simulation resulted in the non-parametric distribution, Y, of my index, y: R: Y -

Re: [R] p value for joint probability

2011-02-01 Thread Petr Savicky
On Mon, Jan 31, 2011 at 09:42:27PM +0100, moleps wrote: Dear all, Given rr-data.frame(r1-rnorm(1000,10,5),r2-rnorm(1000,220,5)) Hello. There is already an answer to your question. However, i think that the above command works in a different way than you expect. The embedded assignments

Re: [R] p value for joint probability

2011-02-01 Thread moleps
My terminology is probably way off. I´ll try again in plain english. I´d like to generate a scatter plot of r1 r2 and color code each pair according to the probability of observing the pair given that the two samples (r1 r2) are drawn from two independent normal distributions.

Re: [R] p value for joint probability

2011-02-01 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of moleps Sent: Tuesday, February 01, 2011 11:32 AM To: Peter Ehlers Cc: r-help@r-project.org Subject: Re: [R] p value for joint probability My terminology is probably way off

Re: [R] p value for joint probability

2011-02-01 Thread Mike Marchywka
Date: Tue, 1 Feb 2011 12:10:31 -0800 From: nord...@dshs.wa.gov To: r-help@r-project.org Subject: Re: [R] p value for joint probability -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of moleps Sent: Tuesday

Re: [R] p value for joint probability

2011-02-01 Thread David Winsemius
On Feb 1, 2011, at 2:31 PM, moleps wrote: My terminology is probably way off. I´ll try again in plain english. I´d like to generate a scatter plot of r1 r2 and color code each pair according to the probability of observing the pair given that the two samples (r1 r2) are drawn from two

Re: [R] p value for joint probability

2011-02-01 Thread Petr Savicky
On Tue, Feb 01, 2011 at 08:31:57PM +0100, moleps wrote: My terminology is probably way off. I´ll try again in plain english. I´d like to generate a scatter plot of r1 r2 and color code each pair according to the probability of observing the pair given that the two samples (r1 r2) are

Re: [R] p value for joint probability

2011-02-01 Thread moleps
Allright.. Appreciate the input on non-zero terminology (:-). What I wanted was: rr-data.frame(r1=rnorm(1000,10,5),r2=rnorm(1000,220,5)) with(rr,plot(r1,r2)) r3-kde2d(r1,r2,lims=c(2,18,200,240)) filled.contour(r3) //M On 1. feb. 2011, at 21.26, David Winsemius wrote: On Feb 1, 2011, at

[R] p value for joint probability

2011-01-31 Thread moleps
Dear all, Given rr-data.frame(r1-rnorm(1000,10,5),r2-rnorm(1000,220,5)) How can I add a column (rr$p) for the joint probability of each r1 r2 pair? I know how to add the column.. I just dont know how to compute the p value for joint probabilities given the two samples. //M

Re: [R] p value for joint probability

2011-01-31 Thread Peter Ehlers
On 2011-01-31 12:42, moleps wrote: Dear all, Given rr-data.frame(r1-rnorm(1000,10,5),r2-rnorm(1000,220,5)) How can I add a column (rr$p) for the joint probability of each r1 r2 pair? If you take the values in each pair to be observations from two independent Normal distributions, it's

[R] p-value from regsubsets

2010-11-10 Thread James Stegen
Hi, does anyone know if there is a way to easily extract p-values from the regsubsets() function? Thanks, James Stegen -- James C. Stegen NSF Postdoctoral Fellow in Bioinformatics University of North Carolina Chapel Hill, NC 919-962-8795 ste...@email.unc.edu

Re: [R] p-value from regsubsets

2010-11-10 Thread Greg Snow
10, 2010 8:31 AM To: r-help@r-project.org Subject: [R] p-value from regsubsets Hi, does anyone know if there is a way to easily extract p-values from the regsubsets() function? Thanks, James Stegen -- James C. Stegen NSF Postdoctoral Fellow in Bioinformatics University of North Carolina

[R] p-value from fisher.test odds ratio 'estimate' (without 2x2 table)

2010-08-17 Thread dschruth
Hello, I'm trying to generate a p.value directly from the output of a fisher exact test odds ratio statistic without resorting to performing another test on the 2x2 table itself. The motivation for this stems from a need to modify the test statistic before calculating it's new p- value. I've

[R] p-VALUE calculation

2010-07-22 Thread jd6688
Here is my dataframe with 1000 rows: employee_id weigth p-value 100 150 101 200 102 300 103 180 . My question: how can I calculate the p-value in R for each employee? the

Re: [R] p-VALUE calculation

2010-07-22 Thread Hadley Wickham
What is your null hypothesis? What is your alternate hypothesis? What is the test statistic? Why do you want a p-value? Hadley On Thu, Jul 22, 2010 at 5:40 PM, jd6688 jdsignat...@gmail.com wrote: Here is my dataframe with 1000 rows: employee_id         weigth       p-value 100            

Re: [R] p-VALUE calculation

2010-07-22 Thread Bill.Venables
103180 0.6722391 (You might want to think about the spelling of 'weigth' a bit more.) -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of jd6688 Sent: Friday, 23 July 2010 8:41 AM To: r-help@r-project.org Subject: [R] p-VALUE

Re: [R] p value

2010-05-18 Thread Greg Snow
-project.org Subject: Re: [R] p value runif(1) Bert Gunter Genentech Nonclinical Statistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On Behalf Of Soham Sent: Saturday, May 15, 2010 9:05 AM To: r-help@r-project.org Subject: [R

[R] p value

2010-05-16 Thread Soham
How to compute the p-value of a statistic generally? -- View this message in context: http://r.789695.n4.nabble.com/p-value-tp2217867p2217867.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] p value

2010-05-16 Thread Tal Galili
Hi Soham, I don't feel your question is well defined. But an equally ill defined answer would be: Through a permutation test. Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com

Re: [R] p value

2010-05-16 Thread Bert Gunter
runif(1) Bert Gunter Genentech Nonclinical Statistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Soham Sent: Saturday, May 15, 2010 9:05 AM To: r-help@r-project.org Subject: [R] p value How to compute the p-value

Re: [R] p-value calculation on a joint distribution

2009-11-16 Thread Greg Snow
Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On Behalf Of Leon Yee Sent: Tuesday, November 10, 2009 11:23 PM To: R help Subject: [R] p-value calculation on a joint distribution

[R] p-value 1

2009-11-13 Thread amor Gandhi
Dear all, I am trying to use SAMr-library(samr), it gives me p-value = 1.001, any idea why? Many thanks, Amor __ z gegen Massenmails. [[alternative HTML version deleted]] __

Re: [R] p-value 1

2009-11-13 Thread Peter Dalgaard
amor Gandhi wrote: Dear all, I am trying to use SAMr-library(samr), it gives me p-value = 1.001, any idea why? Well, the package maintainer might be a better place to ask. (Bring an example if you do.) However, as a generic matter, the result is obviously nonsensical, and could be due to

Re: [R] p-value 1

2009-11-13 Thread amor Gandhi
the p-value of 0.004332? Sorry for these questions. Best wishes, Amor --- Peter Dalgaard p.dalga...@biostat.ku.dk schrieb am Fr, 13.11.2009: Von: Peter Dalgaard p.dalga...@biostat.ku.dk Betreff: Re: [R] p-value 1 An: amor Gandhi amorigan...@yahoo.de CC: r-h...@stat.math.ethz.ch Datum

Re: [R] p-value 1

2009-11-13 Thread Ista Zahn
Dalgaard p.dalga...@biostat.ku.dk schrieb am Fr, 13.11.2009: Von: Peter Dalgaard p.dalga...@biostat.ku.dk Betreff: Re: [R] p-value 1 An: amor Gandhi amorigan...@yahoo.de CC: r-h...@stat.math.ethz.ch Datum: Freitag, 13. November 2009, 15:24 amor Gandhi wrote: Dear all, I am trying to use SAMr

[R] p-value calculation on a joint distribution

2009-11-10 Thread Leon Yee
Dear R users, For a uni-variable distribution represented in a numerical vector, we can obtain a distribution function using 'ecdf', and then calculate corresponding p-values. But if I have a 2-column dataframe representing a bi-variable joint distribution, given a pair of values, how can

[R] P-value and R-squared variable selection criteria

2009-09-24 Thread Lucas Sevilla García
Hi R community I have a question. I'll explain my situation. I have to build a climate model to obtain monthly and annual temperature from 2004 to 2008 from a specif area in Almeria (Spain). To build this climate model, I will use Multiple regression. My dependant variable will be monthly and

Re: [R] P-value and R-squared variable selection criteria

2009-09-24 Thread JLucke
Lucas This problem is very old --- older than keypunches. There are several methods for selecting variables (forward, backwards, both, all subsets) using a variety of criteria (p-values, R^2, adjusted R^2, Cp, AIC, BIC, and more). Be sure you understand the methods, especially the tendency to

Re: [R] P-value and R-squared variable selection criteria

2009-09-24 Thread Daniel Malter
- -Ursprüngliche Nachricht- Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Im Auftrag von Lucas Sevilla García Gesendet: Thursday, September 24, 2009 8:49 AM An: r-help@r-project.org Betreff: [R] P-value and R-squared variable selection criteria

[R] p value from F value and dfs

2009-09-23 Thread Sascha Wolfer
Dear List, is there an easy and fast way to compute the p value from a given F value and given degrees of freedom for an effect and the dfs for the residuals? I think of a function like this: compute.p(F.value, numerator.dfs, denominator.dfs) which returns the p value. Thanks, Sascha

Re: [R] p value from F value and dfs

2009-09-23 Thread Jorge Ivan Velez
Hi Sascha, Take a look at ?pf HTH, Jorge On Wed, Sep 23, 2009 at 11:34 AM, Sascha Wolfer wrote: Dear List, is there an easy and fast way to compute the p value from a given F value and given degrees of freedom for an effect and the dfs for the residuals? I think of a function like this:

[R] p-value for the parameter in ARIMA model with R

2009-06-17 Thread FMH
Dear All, I decided to use an AR(1) model for the residual series and trying to find the p-value for each parameter by using arima command in R, but i never find it from the output. The output gives me the parameter and mean's value, its standard error, estaimated variance, AIC and

[R] p-value=0 running coxph

2008-10-30 Thread GSt
Dear all, I have a question concerning the p-value. When running coxph I get a p-value = 0. :confused: Can this be true? Why aren´t there decimal points? Is there a way to find out the exact p-value? Here is the output:

Re: [R] p-value=0 running coxph

2008-10-30 Thread C.H.
Because your test statistics (z=-9.16) is too small. (well, z 1.96 or Z -1.96 is p=0.05) If you really want the exact p, you can try pnorm(-9.16) if I were you I will quote p0.001 as well as the test statistics. (if you are going to submit it to medical journals) C On Thu, Oct 30, 2008 at

Re: [R] p-value=0 running coxph

2008-10-30 Thread Thomas Lumley
On Thu, 30 Oct 2008, GSt wrote: Dear all, I have a question concerning the p-value. When running coxph I get a p-value = 0. :confused: Can this be true? Why aren??t there decimal points? Is there a way to find out the exact p-value? The p-values are rounded, and of course they are not

Re: [R] p-value for Nonmetric Multidimentional Scaling?

2008-07-05 Thread Gavin Simpson
On Fri, 2008-07-04 at 21:14 -0700, Michael Denslow wrote: Dear R-helpers, I am running metaMDS in the vegan package, which uses isoMDS in MASS, to perform Nonmetric Multidimentional Scaling (NMDS). I have seen some authors report a p-value for the NMDS ordination based on randomization of

[R] p-value for Nonmetric Multidimentional Scaling?

2008-07-04 Thread Michael Denslow
Dear R-helpers, I am running metaMDS in the vegan package, which uses isoMDS in MASS, to perform Nonmetric Multidimentional Scaling (NMDS). I have seen some authors report a p-value for the NMDS ordination based on randomization of the dataset. As I understand it this is meant to compare the

Re: [R] p-value for Nonmetric Multidimentional Scaling?

2008-07-04 Thread Bill.Venables
July 2008 2:14 PM To: r-help@r-project.org Subject: [R] p-value for Nonmetric Multidimentional Scaling? Dear R-helpers, I am running metaMDS in the vegan package, which uses isoMDS in MASS, to perform Nonmetric Multidimentional Scaling (NMDS). I have seen some authors report a p-value for the NMDS

[R] p-value

2008-06-12 Thread mohamed nur anisah
Dear R User, say I have this sample of data ( attach with). What i'm going to do is to test whether this data is uniformly distributed or not by finding the p-value. I've tried using the punif command but it gave me the value of 1 of all the data. Any suggestion on R command to find the

Re: [R] p-value

2008-06-12 Thread John Gant
Here's a sample: unif_rand_1 - runif(1000); unif_rand_2 - runif(1000); ks.test(unif_rand_1,unif_rand_2); Two-sample Kolmogorov-Smirnov test data: unif_rand_1 and unif_rand_2 D = 0.021, p-value = 0.9802 alternative hypothesis: two-sided So in your case: ks.test( runif( length(

Re: [R] p-value

2008-06-12 Thread Richardson, Patrick
(midpoint) HTH, Patrick -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of mohamed nur anisah Sent: Thursday, June 12, 2008 9:56 AM To: r-help@r-project.org Subject: [R] p-value Dear R User, say I have this sample of data ( attach with). What i'm going

Re: [R] p-value

2008-06-12 Thread Felipe Carrillo
Not sure if this is what you are looking for but you can get the p-value with something like this: # Create a vector mydata- c(132968364, 135945080, 156539568, 157817896, 162399496, 168344072, 173146584, 176302744, 182878168, 183946152,

[R] p-value in Spearman rank order

2008-02-28 Thread Anne-Katrin Link
Dear R-helpers, I would like to do a Spearman rank order test, and used the cor() function with the method spearman. It gives me a number (correlation coefficient?) , but how can I get the p-value? Thank you for the help in advance! Regards, Anne-Katrin -- [[alternative HTML

Re: [R] p-value in Spearman rank order

2008-02-28 Thread Gavin Simpson
On Thu, 2008-02-28 at 16:17 +0100, Anne-Katrin Link wrote: Dear R-helpers, I would like to do a Spearman rank order test, and used the cor() function with the method spearman. It gives me a number (correlation coefficient?) , but how can I get the p-value? Thank you for the help in

Re: [R] p-value in Spearman rank order

2008-02-28 Thread Liviu Andronic
On 2/28/08, Anne-Katrin Link [EMAIL PROTECTED] wrote: I would like to do a Spearman rank order test, and used the cor() function with the method spearman. It gives me a number (correlation coefficient?) , but how can I get the p-value? You're probably looking for rcorr() from Hmisc. It

Re: [R] p-value in Spearman rank order

2008-02-28 Thread Daniel Malter
-- -Ursprüngliche Nachricht- Von: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Im Auftrag von Anne-Katrin Link Gesendet: Thursday, February 28, 2008 10:17 AM An: r-help@r-project.org Betreff: [R] p-value in Spearman rank order Dear R-helpers, I would like to do a Spearman rank order test, and used the cor

[R] p-value for fixed effect in generalized linear mixed model

2008-02-20 Thread Martijn Vandegehuchte
Dear R-users, I am currently trying to switch from SAS to R, and am not very familiar with R yet, so forgive me if this question is irrelevant. If I try to find the significance of the fixed factor spikes in a generalized linear mixed model, with site nested within zone as a random factor, I

Re: [R] p-value for fixed effect in generalized linear mixed model

2008-02-20 Thread Douglas Bates
The lme4 package is still in development. I can't guarantee the evaluation of the deviance for the quasipoisson family. There is a certain amount of mystery about the role of the dispersion parameter and the null deviance in families like quasipoisson and quasibinomial. For the time being don't

[R] P-value

2007-10-03 Thread Emre Unal
How can i print only the P-Value of the kolmogorov smirnov test? ks.test(VeriSeti1, VeriSeti2) Two-sample Kolmogorov-Smirnov test data: VeriSeti1 and VeriSeti2 D = 0.5, p-value = 0.4413 alternative hypothesis: two-sided This expression gives me the whole test results, i need only

Re: [R] P-value

2007-10-03 Thread Chung-hong Chan
ks.test(x,y)[2] On 10/3/07, Emre Unal [EMAIL PROTECTED] wrote: How can i print only the P-Value of the kolmogorov smirnov test? ks.test(VeriSeti1, VeriSeti2) Two-sample Kolmogorov-Smirnov test data: VeriSeti1 and VeriSeti2 D = 0.5, p-value = 0.4413 alternative hypothesis:

  1   2   >