similar to: R-alpha: signif( small , d) gives NA

The Intel compiled version also fails the below test: > ###------------ Very big and very small > umach <- unlist(.Machine)[paste("double.x", c("min","max"), sep='')] > xmin <- umach[1] > xmax <- umach[2] > tx <- unique(outer(-1:1,c(.1,1e-3,1e-7)))# 7 values (out of 9) > tx <- unique(sort(c(outer(umach,1+tx))))# 11 values

Hello, I am all new here. Thanks for the job done, R really helped me in my thesis lately. However, I am kind of new in statistics, coming from mecanical engineering, and I mostly teached myself with "The R Book", so I may do silly things some time. PLease tell me if you think so. Anyway, I've just build up a piecewise linear model to fit some data, including some interaction and i

I use a 32 bit machine. For those using 64 bit machines, what is the .Machine$double.xmin for for machines? regards. Stephen. ----- Original Message ----- From: Achim Zeileis <Achim.Zeileis at wu-wien.ac.at> Date: Tuesday, July 12, 2005 10:51 am Subject: Re: [R] exact values for p-values - more information. > On Tue, 12 Jul 2005, S.O. Nyangoma wrote: > > > > If they

2017-06-22 19:49 GMT+02:00 Uwe Ligges <ligges at statistik.tu-dortmund.de>: > On 22.06.2017 17:11, Bernd Funovits wrote: >> >> Hello, >> >> I experienced some unexpected behaviour while determining the rank of matrices (sometimes 1x1 matrices): >> base::qr(matrix(1e-20))$rank returns 1 (incorrect) >> base::qr(diag(c(1, 1e-20)))$rank returns 2

printCoefmat(signif.legend =FALSE) does not work properly. The option "signif.legend = FALSE" is ignored as shown in the example below. cmat <- cbind(rnorm(3, 10), sqrt(rchisq(3, 12))) cmat <- cbind(cmat, cmat[,1]/cmat[,2]) cmat <- cbind(cmat, 2*pnorm(-cmat[,3])) colnames(cmat) <- c("Estimate", "Std.Err", "Z value", "Pr(>z)") #

Dear list, I'm trying to write a class for Gaussian error propagation of measured values and their (estimated) errors, > setClass("sec", representation(val="numeric", err="numeric")) I've already successfully implemented basic arithmetics using mostly the "Arith" group generics. But I'm running into trouble when trying to get signif() to

PROBLEM `signif' does change to scientic notation at different levels depending on the number of significant digits in the input. This can generate tables where figures change ``irregularly'' from normal to scientific notation. PROPOSAL The change to the scientific notation should be made only if the figure in scientific notation - with potentially as

Hi R fans, As a newbie following the five-hour rule (after hitting my head against the wall for five hours, post to this list), I am appealing for some help understanding geom_rect() in ggplot2. What I want to do is very simple. I want to generate a plot of rectangles. Each one represents a business cycle. The x-values will be pairs representing the start and end of each cycle. The y-values

Dear R-Devel, As I am sure many of you know, a special issue of The American Statistician just came out, and its theme is the [mis]use of P values and the many common ways in which they are abused. The lead editorial in that issue mentions the 2014 ASA guidelines on P values, and goes one step further, by now recommending that the words "statistically significant" and related simplistic

Dear R users, Consider the following example function: f = function(a,b) { g = function(x) a*x + b h = function(x) g(x)^2 + x^2 opt = optimize(h,lower = -1, upper = 1) x.min = opt$minimum h.xmin = opt$objective g.xmin = g(x.min) return(c(x.min, h.xmin, g.xmin)) } In my real problem the function that plays the role of "g" is costly to compute. Now, to

This is too bad (and an old bug at least since 0.49); it also breaks some code I have been using in S-plus: > tst <- function(...) ls(...) > tst()#--- should really give the same as ls() !! character(0) > ls() [1] "%in%" "As.integer" "Dir" [4] "Dir0" "Fn12"

Hi all, I get two different results when using sum() and the switch na.rm. The result is correct when na.rm=FALSE. Linux Redhat 7.3, R version 1.6.1. I've had no luck searching the mail archives, so I was hoping somebody could explain/check this one for me. I will need to apply the function to missing data, simple as it is. Code: x<-matrix(runif(20,0,5)%/%1,4,5) # random matrix

2017-06-22 20:31 GMT+02:00 Uwe Ligges <ligges at statistik.tu-dortmund.de>: > > > On 22.06.2017 20:09, I?aki ?car wrote: >> >> 2017-06-22 19:49 GMT+02:00 Uwe Ligges <ligges at statistik.tu-dortmund.de>: >>> >>> On 22.06.2017 17:11, Bernd Funovits wrote: >>>> >>>> >>>> Hello, >>>> >>>> I

I'm still trying to install R on my Irix machine. Now I have a new problem that crops up during the checks. I've found the root cause, and it's that R is returning zero for certain things for reasons I don't understand. 2.225073859e-308, entered directly into R, responds "2.225074e-308". 2.225073858e-308 responds "0". Their negative values respond

I'm still trying to install R on my Irix machine. Now I have a new problem that crops up during the checks. I've found the root cause, and it's that R is returning zero for certain things for reasons I don't understand. 2.225073859e-308, entered directly into R, responds "2.225074e-308". 2.225073858e-308 responds "0". Their negative values respond

Hi, I recently spent quite a bit of time trouble shooting a function that I had written only to discover that the problem I was having was with the comparison operator. I assumed that the following would return TRUE: > testMean <- 82.8 + 0.1 > testMean [1] 82.9 > testMean == 82.9 [1] FALSE Apparently this has to do with deciml places. Look: > newTest <- 82.0 > newTest [1]

[Environment: Win Xp, Miktex 2.7, R 2.9.2] In an Sweave document, I'm displaying the results of car:::Anova() tests, that look like this in the generated .tex file: \begin{Soutput} Type III MANOVA Tests: Pillai test statistic Df test stat approx F num Df den Df Pr(>F) (Intercept) 1 0.86 90.38 4 60 <2e-16 *** --- Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1

###Dear R users ###I have been using SensoMineR package from CRAN for most of my work in sensory data analysis and from my usage experience, I encountered some areas for improvement and considered ###modifying the function in SensoMineR package for my personal use. I felt that it could be useful to share this to the community for enabling adoption by other users where they might require a

Hi R help, I am trying to determine how nls() generates a function based on the self-starting SSlogis and what the formula for the function would be. I've scoured the help site, and other literature to try and figure this out but I still am unsure if I am correct in what I am coming up with. ************************************************************************** dat <-

I have a xml exported by Manifold GIS but I'm not being able to import it into R using XLM package. The file have this structure: <?xml version="1.0" encoding="UTF-8" ?> - <layout> <name>Layout 2</name> <pagesByX>1</pagesByX> <pagesByY>1</pagesByY> - <elements> <legend