Warning message:In pt(q, df, lower.tail, log.p) : NaNs produced
Obviously your df=0 in pt(STATISTIC, df). Regards, Yihui -- Yihui Xie <xieyihui at gmail.com> Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.name School of Statistics, Room 1037, Mingde Main Building, Renmin University of China, Beijing, 100872, China On Tue, Jan 6, 2009 at 7:26 PM, Moumita Das
<das.moumita.online at gmail.com> wrote:
Hi friends, Any idea why do i get this warning?And also why all computed p-values are NaN. Have shown below what i did in Windows r-console.:--
df
c1 c2 1 1 50 2 NA NA 3 4 NA 4 7 6 5 NA 7 6 10 10
r<-cor(x=df,y=NULL,use="complete.obs",method=c("pearson"))
r
c1 c2 c1 1.0000000 -0.9148074 c2 -0.9148074 1.0000000
cor.p.values<- function(r, n)
+ {
+ df <- n - 2
+ ESTIMATE <- c(cor = r)
+ PARAMETER <- c(df = df)
+ STATISTIC <- c(sqrt(df) * r / sqrt(1 - r^2))
+ p <- pt(STATISTIC, df)
+ return(2 * pmin(p, 1 - p))
+ }
cor.p.values(r,2)
[1] NaN NaN NaN NaN
*Warning message:
In pt(q, df, lower.tail, log.p) : NaNs produced
*
Any help will be appreciated.. :)
--
Thanks
Moumita
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.