similar to: qt with df<1

Sorry about the html-formatted message. Here it is again in plain text. Hello, The function qt returns NaN for degrees of freedom <1. For example: > qt(0.5,0.5) [1] NaN Warning message: In qt(p, df, lower.tail, log.p) : NaNs produced But qt(0.5,0.5) should be 0, since the distribution is symmetric. > pt(0,0.5) [1] 0.5 It actually fails with any value, as long as df<1. Is this a

Full_Name: Gerard Torrent Version: R version 2.8.0 (2008-10-20) OS: Linux 2.6.27.5-41.fc9.x86_64 #1 SMP Submission from: (NULL) (85.52.227.233) In some cases qt complains about NaNs and don't gives the correct result: > qt(0.1,, 0.1) [1] NaN Warning message: In qt(p, df, lower.tail, log.p) : NaNs produced But the result can be found: > pt(-1.60443e+06, 0.1) [1] 0.09999997 If I

Full_Name: Jens Lund Version: Version 1.7.0 (2003-04-16) OS: Win NT 4.0 SP 6 Submission from: (NULL) (193.3.225.210) As the T distribution is symmetrical around 0 qt(0.5,df) should return 0 for any df. However for df close to 1 it seems to have problems as seen by: qt(0.5,seq(1,1.1,by=0.0001)) For example: > qt(0.5,1.01) [1] -0.2300470 Higher df as df=1.1 (and up to approx df=1.7) seem

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

I would presume this is another manifestation of what I reported (reproduced below) on 2003-12-01. cajo.terbraak at wur.nl wrote: >Full_Name: cajo ter Braak >Version: 2.1.1 >OS: Windows >Submission from: (NULL) (137.224.10.105) > > ># I would like to attach the matrix C in the Rdata file; it is 50x50 and comes >from a geostatistical problem (spherical covariogram) >

Full_Name: cajo ter Braak Version: 2.1.1 OS: Windows Submission from: (NULL) (137.224.10.105) # I would like to attach the matrix C in the Rdata file; it is 50x50 and comes from a geostatistical problem (spherical covariogram) > rm(list=ls(all=TRUE)) > load(file= "test.eigen.Rdata") > ls() [1] "C" "eW" > > sym.check = max(abs(C - t(C))) # should

I was experimenting yesterday with a binomial make.link option for estimating student t binary response models, tentatively called gossit, and I noticed eventually that the R qt function doesn't like df < 1. Vaguely recalling that Splus didn't seem to mind such weirdness, I checked on our soon to be defunct Splus6.2 and sure enough, it produced plausible answers instead of R's

Dear R gurus, I am trying to work out the problem given in Nested design - Montgomery - Design of Experiments p.561 I have attached a pdf of the data as well the anova table. It is a mixed model with Supplier as fixed effect and batches within the supplier as random effects. I am able to work out the error stratums as below using aov. Which agrees perfectly with the book example

I define the following function to convert a t-value with degrees of freedom DF to another t-value with different degrees of freedom fullDF: tConvert <- function(tval, DF, fullDF) ifelse(DF>=1, qt(pt(tval, DF), fullDF), 0) It works as expected with the following case: > tConvert(c(2,3), c(10,12), 12) [1] 1.961905 3.000000 However, it gives me warning for the example below although the

help(qt) states that: "ncp non-centrality parameter delta; currently except for rt(), only for abs(ncp) <= 37.62" so I would expect that calling qt with non-centrality parameter exceeding 37.62 should fail, instead e.g. calling > mapply(function(x) qt(p = 0.9, df = 55, ncp = x),35:45) gives: [1] 40.21448 41.35293 42.49164 43.68862 44.82945 45.97048 47.11170 48.25310 [9]

Hi all, Recently our statistics students noticed that their Gibbs samplers were crashing due to some NaNs in some parameters. The NaNs came from mvrnorm (Ripley & Venables' MASS package multivariate normal sampling function) and with some more investigation it turned out that they were generated by function eigen, the eigenvalue computing function. The problem did not seem to happen

Hello All I am looking for high precision values for the normal distribution in the tail,(1e-10 and 1 - 1e-10) as the R package that I am using sets any number which is out of this range to these values and then calls the qnorm and qt function. What I have noticed is that the qnorm implementation in R is not symmetric when looking at the tails. This is quite surprising to me, as it is well known

Hi, Actually I am facing a similar problem. I would like to fit both an ordinary (symmetric) and a non-central t distribution to my (one-dimensional) data (quite some values.. > 1 mio.). For the symmetric one, fitdistr or funInfoFun (using fitdistr) from the qAnalyst package should do the job, and for the non-central one.. am I right to use gamlss(x ~ 1, family=GT()) ? Anyway, I am a little

I can't get confint.glm to work from within a function. Consider the following (using R 1.9.1, Windows 2000): # FIRST: SOMETHING THAT WORKS FROM A COMMAND PROMPT DF <- data.frame(y=.1, N=100) (fit <- glm(y~1, family=binomial, data=DF, weights=DF[,"N"])) Call: glm(formula = y ~ 1, family = binomial, data = DF, weights = DF[, "N"]) Coefficients:

Hello, I am trying to calculate and plot mean and confidence intervall for a set of data. This is the code that I am currently using: means <- sapply(data, mean, na.rm=TRUE) n <- sapply(data,length) stdev <- sqrt(sapply(data, var, na.rm=TRUE)) ciw <- qt(0.98, n) * stdev / sqrt(n) par(mgp=c(2,0.6,0), las=2, fin=c(7,3), mai=c(1,0.5,0.2,0.2), cex=0.8) plotCI(x=means, uiw=ciw,

On Thu, Mar 1, 2018 at 2:08 AM, Nuno Lopes <nunoplopes at sapo.pt> wrote: > We can do "add %x, undef" => "undef" because for any value of %x, we can > always find a value that when added to %x produces any value in the domain > of integers. > > This is not the case with floats since with some inputs, e.g., NaNs, we > are not able to produce some

CentOS Errata and Security Advisory 2012:0880 Moderate Upstream details at : https://rhn.redhat.com/errata/RHSA-2012-0880.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 50f94fe74c884b57a1a2101510283eeffdfd8d0af4af5b1246d61f18154919b5 phonon-backend-gstreamer-4.6.2-24.el6.i686.rpm

CentOS Errata and Bugfix Advisory 2012:1246 Upstream details at : https://rhn.redhat.com/errata/RHBA-2012-1246.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 317cf57e88f211adfb49da53292ccf6b97abc89424a89db7cb04680ce8d22dfe phonon-backend-gstreamer-4.6.2-25.el6.i686.rpm

CentOS Errata and Bugfix Advisory 2014:0201 Upstream details at : https://rhn.redhat.com/errata/RHBA-2014-0201.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 727dc3b50f032ae08c8d130621aceb876cced3a08e182011d6f9c2dae1ea7c73 phonon-backend-gstreamer-4.6.2-28.el6_5.i686.rpm

CentOS Errata and Bugfix Advisory 2016:1047 Upstream details at : https://rhn.redhat.com/errata/RHBA-2016-1047.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) x86_64: f63ab6aa5fdd59b19af69858bebd6828e8ecdf530873308fd9576493a633fead qt-4.8.5-12.el7_2.i686.rpm f26bb3e7c48329d902463e8ef2803617e858ade87dfc4ce9aa949b8169fc3161